Exigiu apenas uma modificação do código porque você tá desdobrando um fator dentro da combinação dos níveis de dois fatores,
# blocos devem vir antes nas fórmulas
av <- aov(resp~anos*ambientes/(blocos+genotipos))
anova(av) # nem todo F tá com denominador correto
names(coef(av))
pattern <- c(outer(levels(anos), levels(ambientes),
function(x,y) paste("anos",x,":ambientes",y,":",sep="")))
des.tab <- sapply(pattern, simplify=FALSE,
grep, x=names(coef(av)[av$assign==5]))
summary(av, split=list("anos:ambientes:genotipos"=des.tab))
À disposição.
Walmes.
==========================================================================Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmeslinux user number: 531218
==========================================================================