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.