<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Fernando,<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Acho que é assim:<br><span style="font-family:courier new,monospace"><br>
## fatores cruzados<br>xtabs(~local+trat, dados)<br><br>## fatores aninhados<br>xtabs(~trat+eras, dados)<br><br>##  local<br>## +bloco dentro de local<br>## +tratamento<br>## +interação local tratamento<br>mf1 <- lm(terms(resp~local/rep+trat+local:trat, keep=TRUE),<br>
          data=dados)<br>anova(mf1)<br><br>##  local<br>## +bloco dentro de local<br>## +eras<br>## +tratamento dentro de eras<br>## +interação local eras<br>## +interação local tratamento dentro de eras<br>mf2 <- lm(terms(resp~local/rep+eras/trat+local:eras/trat, keep=TRUE),<br>
          data=dados)<br>anova(mf2)<br><br>require(lme4)<br><br>mm1 <- lmer(resp~(1|local/rep)+(1|trat)+(1|local:trat), data=dados)<br>summary(mm1)<br><br>mm2 <- lmer(resp~(1|local/rep)+(1|eras/trat)+<br>            (1|local:eras)+(1|local:eras:trat), data=dados)<br>
summary(mm2)</span><br><br>À disposição.<br>Walmes.</div><br></div>