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 . <[email protected]> > 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 > [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. >
_______________________________________________ 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.
