augustofilho — Jun 3, 2014, 11:24 AM
banco1=function(){
#Semente fixa igual a 2014
set.seed(2014)
#Tamanho da amostra
cat("Digite o tamanho da amostra:\n")
n=scan( ,nlines=1)
#Simulando o alpha
alpha=seq(0.2,4,length=n)
#Simulando o delta
delta=seq(-2.4,-1.6,length=n)
#Restricoes para os tau's:
#artigo do J.S. Roberts, p. 178
#1) tau0=0 (por definicao)
#2) tau(z) = tau(-z) (simetrico em torno do zero)
#3) soma(tau_{k=0 ate M)=0
tau=rep(0,3) #cria um vetor com zeros de tamanho 3 para
#receber os valores tau a esquerda do zero
tau[1]=runif(1,min=-3,max=-2)
tau[2]=runif(1,min=-2,max=-1.33)
tau[3]=runif(1,min=-1,max=-0.67)
#Vetor final de tau's
tau=c(tau,0,-tau)
#Agrupando a saida de dados
cat("Parametros dos itens utilizados na primeira simulacao\n")
cat("\n")
saida=cbind(alpha,delta,tau[1],tau[2],tau[3],tau[4],tau[5],tau[6],tau[7])
colnames(saida)=c("alpha","delta","tau1","tau2","tau3","tau4","tau5","tau6","tau7")
return(saida)
}