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: [email protected] twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
