[R-br] Blocos com algumas repetições do material

[Eder David Borges da Silva]

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