Simulação de dados de um modelo linear

Bom dia, 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) 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. 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 ao gerar os mesmos a cada vez que montar um banco de dados!!! Obrigado! Saudações Greice. Abaixo esta o programa no R # Tamanho de amostra n.af=2 n.suj=29 rep=3 repl=5 ntot=174 var.af=0.915544 var.sj=51.987741 var.sjaf=1.215483 var.erro=0.8563 theta = -1.24795 ### Dados sujeito = rep(1:n.suj,each=(n.af*rep)) aferidor = rep(rep(1:n.af,each=rep),n.suj) interacao = rep(1:(n.af*n.suj),each=rep) ### Efeitos aleatórios set.seed(123) ef.aferidor = rnorm(n.af,0,sqrt(var.af)) ef.sujeito = rnorm(n.suj,0,sqrt(var.sj)) ef.inter = rnorm((n.suj*n.af),0,sqrt(var.sjaf)) ef.suj.rep = rep(ef.sujeito,each=(n.af*rep)) ef.af.rep = rep(rep(ef.aferidor,each=rep),n.suj) ef.inter.rep = rep(ef.inter,each=rep) ### modelo linear media = theta + ef.suj.rep + ef.af.rep + ef.inter.rep ### Geração dos dados #### y = replicate(repl,rnorm(ntot,media,sqrt(var.erro)))

Greice, Depende dos seus objetivos. Acho mais apropriado gerar novos efeitos aleatórios porque você está replicando em computador um experimento de campo. Aqui eu considero que cada vez que fizer o experimento de campo, você terá outra amostra de sujeitos e aferidores sorteadas respectivas populações para compor o experimento. Então, os efeitos aleatórios de experimento para experimento vão mudar. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Ok! Muito Obrigado! 2011/9/22 Walmes Zeviani <walmeszeviani@gmail.com>
Greice,
Depende dos seus objetivos. Acho mais apropriado gerar novos efeitos aleatórios porque você está replicando em computador um experimento de campo. Aqui eu considero que cada vez que fizer o experimento de campo, você terá outra amostra de sujeitos e aferidores sorteadas respectivas populações para compor o experimento. Então, os efeitos aleatórios de experimento para experimento vão mudar.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (2)
-
Greice Laureano
-
Walmes Zeviani