> Boa tarde Pessoal, > > Tenho experimentos em blocos(Locais) e dentro deles gostaria de por > algumas de repetições de alguns tratamentos para estimar a variância dentro > do bloco, que desconfio ser grande. Por exemplo, tenho 5 Cultivares > gostaria de repetir 3 delas dentro do bloco, tendo ao todo 8 parcelas, o > ideal seria ter uma repetição de tudo, porem na pratica isso não > é possível. > Minhas duvidas são: > 1) Estou definindo de maneira correta essas repetições dentro do modelo? > 2) Como o método de calculo dos dois modelos abaixo é igual para os > efeitos fixos, por que o intercepto esta com valores diferentes? > 3) Neste dois casos(m0 e m1, abaixo) o BLUP seria a média mais o efeito da > cultivar? > Obrigado pela ajuda > > require(lme4) > ### Simulando dados > dados <- > expand.grid(Variedade=factor(c(letters[1:3],letters[1:5])),Local=factor(1:5)) > dados$Rep <- factor(rep(c(2,2,2,1,1,1,1,1),t=5)) > dados$Resp <- sort(rnorm(nrow(dados))) > dados > > ### Modelos > # Sem efeito dentro do local > m0 <- lmer(Resp~1+(1|Local)+(1|Variedade),dados) > summary(m0) > ranef(m0) > # Com efeito dentro do local > m1 <- lmer(Resp~1+(0+Rep|Local)+(1|Variedade),dados) > summary(m1) > ranef(m1) >
_______________________________________________ R-br mailing list [email protected] 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.
