<div>Bom dia,<br><br>Tenho uma duvida a respeito de uma simulação de um modelo linear Two-way com os efeitos fatorias aleatorios (sujeito, aferidor e interação)</div>

<div><br></div><div>A questão é: Quero gerar vários bancos de 
dados (y)  para realizar a simulação de diversos métodos de estimação para o coeficiente de correlação intraclasse.<br>  <br></div><div>Eu gerei os efeitos fatoriais uma vez e a partir dessa geração geramos os y do modelo linear... Minha dúvida é se preciso gerar os efeitos 
fatoriais novamente a cada vez que gero um y ou se posso fazer como disse anteriormente, minha maior preocupação é perder informação dos efeitos fatorias<br>ao gerar os mesmos a cada vez que montar um banco de dados!!!<br>
</div><div><br></div><div>Obrigado!</div><div>Saudações</div><div>Greice.</div>

<div><br></div><div>Abaixo esta o programa no R</div><div><br></div><div><font color="#993399"># Tamanho de amostra</font></div><div><div><font color="#993399"><a href="http://n.af/" target="_blank">n.af</a>=2</font></div>


<div><font color="#993399">n.suj=29</font></div><div><font color="#993399">rep=3</font></div><div><font color="#993399">repl=5</font></div><div><font color="#993399">ntot=174</font></div>
<div><font color="#993399"><br></font></div><div><font color="#993399"><a href="http://var.af/" target="_blank">var.af</a>=0.915544</font></div><div><font color="#993399"><a href="http://var.sj/" target="_blank">var.sj</a>=51.987741</font></div>


<div><font color="#993399">var.sjaf=1.215483</font></div><div><font color="#993399">var.erro=0.8563</font></div><div><font color="#993399">theta = -1.24795</font></div>
<div><font color="#993399"><br></font></div><div><font color="#993399">### Dados</font></div><div><font color="#993399"><br></font></div><div><font color="#993399">sujeito = rep(1:n.suj,each=(<a href="http://n.af/" target="_blank">n.af</a>*rep))  </font></div>


<div><font color="#993399">aferidor = rep(rep(1:<a href="http://n.af/" target="_blank">n.af</a>,each=rep),n.suj)</font></div><div><font color="#993399">interacao = rep(1:(<a href="http://n.af/" target="_blank">n.af</a>*n.suj),each=rep)</font></div>


<div><font color="#993399"><br></font></div><div><font color="#993399">### Efeitos aleatórios</font></div><div><font color="#993399"><br></font></div>
<div><font color="#993399">set.seed(123)</font></div><div><font color="#993399"><br></font></div><div><font color="#993399">ef.aferidor = rnorm(<a href="http://n.af/" target="_blank">n.af</a>,0,sqrt(<a href="http://var.af/" target="_blank">var.af</a>))</font></div>


<div><font color="#993399">ef.sujeito = rnorm(n.suj,0,sqrt(<a href="http://var.sj/" target="_blank">var.sj</a>))</font></div><div><font color="#993399">ef.inter = rnorm((n.suj*<a href="http://n.af/" target="_blank">n.af</a>),0,sqrt(var.sjaf))</font></div>


<div><font color="#993399"><br></font></div><div><font color="#993399">ef.suj.rep = rep(ef.sujeito,each=(<a href="http://n.af/" target="_blank">n.af</a>*rep))</font></div><div><font color="#993399">ef.af.rep = rep(rep(ef.aferidor,each=rep),n.suj)</font></div>


<div><font color="#993399">ef.inter.rep = rep(ef.inter,each=rep)</font></div><div><font color="#993399"><br></font></div><div><font color="#993399">### modelo linear</font></div>
<div><font color="#993399"><br></font></div><div><font color="#993399">media = theta + ef.suj.rep + ef.af.rep + ef.inter.rep</font></div><div><font color="#993399"><br>
</font></div><div><font color="#993399">### Geração dos dados ####</font></div><div><font color="#993399"><br></font></div><div><font color="#993399">y = replicate(repl,rnorm(ntot,media,sqrt(var.erro)))</font></div>
<div><font color="#993399"><br></font></div><div><font color="#993399"><br></font></div></div>