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