Thiago,

Confesso que nunca vi alguém apresentar em um artigo uma comparação múltipla multivariada. Acredito que você pode usar um algoritmo k-means classificar os centróides em grupos homogeneos e tals, da mesma forma como SkottKnott faz. De qualquer forma, a comparação de modelos que eu me referi está representada no código abaixo

da <- data.frame(tr=gl(3,8))
da$y1 <- rnorm(da$tr)
da$y2 <- rnorm(da$tr)
da$tr12vs3 <- as.factor(ifelse(as.numeric(da$tr)<=1,1,2))

m0 <- aov(cbind(y1,y2)~1, da)
m1 <- aov(cbind(y1,y2)~tr12vs3, da)
m2 <- aov(cbind(y1,y2)~tr, da)

anova(m0, m1)
anova(m1, m2)

À 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
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================