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.