[R-br] Simulação de dados de um modelo linear

[Greice Laureano]

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)))
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110922/7ceea08a/attachment.html>