Uma maneira alternativa

grupo<-ex11
attach(grupo)
grupo.av<-by(grupo,local,function(x) anova(aov(resposta~factor(trat)+factor(bloco),data=x)));grupo.av

grupo.av[1] # uma ANOVA

grupo.av[[1]]$Mean[3] # um QMRes

# Todos os QMRes
qm<-numeric()
for(i in 1:4) qm[i]<-grupo.av[[i]]$Mean[3];qm