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 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================