[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