[R-br] Cluster + ANOVA
Guilherme Moraes Ferraudo
guiferraudo em gmail.com
Sexta Julho 29 14:59:44 BRT 2011
Vinicius segue o código. Lembrando que no meu código o objeto "dados"
é um data.frame
E MANOVA significa "Multivariate Analysis of Variance" já que vc tem
mais que uma variável resposta.
#### Código R ####
paises_kmeans7 = kmeans(dados, centers = 7, nstart = 1111) # nstart:
semente aleatória
dados_kmeans7_final = data.frame(dados,cluster =
as.factor(paises_kmeans7$cluster))
##Ajustando o modelo multivariado considerando os grupos formados pelo
método K-means
mod_manova = manova(as.matrix(dados) ~ dados_kmeans7_final$cluster)
E = mod_manova$res #matriz de resíduos (erros)
#por default já é o Pillai
summary(mod_manova)
#ou
summary(mod_manova, test = "Pillai")
summary(mod_manova, test = "Wilks")
summary(mod_manova, test = "Hotelling")
summary(mod_manova, test = "Roy")
##
Se surgir mais dúvidas vai falando.
Abs
--
Guilherme Moraes Ferraudo
Jaboticabal/Campinas - SP
http://lattes.cnpq.br/2096118558794430
Mais detalhes sobre a lista de discussão R-br