[R-br] DESDOBRAMENTO DE UM FATORIAL TRIPLO

[Walmes Zeviani]

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 em ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130108/640e3c2e/attachment.html>