[R-br] Parametrização com lme4

[FHRB Toledo]

Adriano e Walmes,

Obrigado pelas sugestões, e por compartilhar o livro!

O modelo ajustou corretamente, perfeito!

Obrigado mesmo.

abraço,
FH


2013/11/20 walmes . <walmeszeviani em gmail.com>

> Fernando,
>
> Acho que é assim:
>
> ## fatores cruzados
> xtabs(~local+trat, dados)
>
> ## fatores aninhados
> xtabs(~trat+eras, dados)
>
> ##  local
> ## +bloco dentro de local
> ## +tratamento
> ## +interação local tratamento
> mf1 <- lm(terms(resp~local/rep+trat+local:trat, keep=TRUE),
>           data=dados)
> anova(mf1)
>
> ##  local
> ## +bloco dentro de local
> ## +eras
> ## +tratamento dentro de eras
> ## +interação local eras
> ## +interação local tratamento dentro de eras
> mf2 <- lm(terms(resp~local/rep+eras/trat+local:eras/trat, keep=TRUE),
>           data=dados)
> anova(mf2)
>
> require(lme4)
>
> mm1 <- lmer(resp~(1|local/rep)+(1|trat)+(1|local:trat), data=dados)
> summary(mm1)
>
> mm2 <- lmer(resp~(1|local/rep)+(1|eras/trat)+
>             (1|local:eras)+(1|local:eras:trat), data=dados)
> summary(mm2)
>
> À disposição.
> Walmes.
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131120/438b1466/attachment.html>