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@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@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.