
Sem entrar no mérito da análise, outra ideia que parece atender sua necessidade... ### <code r> ### dados fictícios! df <- data.frame(UF=rep(c("AC","AL","AM","MS","PR"), each=4), ano=rep(2003:2006, 5), V1=rnorm(20), V2=rnorm(20), v3=rnorm(20), v4=rnorm(20)) df df2 <- aggregate(.~UF, data=df, mean) ### criando rownames únicos a partir dos dados... rownames(df2) <- df2$UF df2 PCA <- prcomp(df[3:6],scale=T) ### informar as variáveis númericas hc <- hclust(dist(df[3:6]), method='ward.D') ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 15 de abril de 2015 00:24, Andre Oliveira <andreolsouza@yahoo.com.br> escreveu:
Exatamente Éder, você chegou exatamente onde eu gostaria, mas continuando seu código ao plotar plot(as.dendrogram(hc),main = "") cria os grupos e separa os ACs, os Als e os AMs no dendograma e eu gostaria de apenas um Estado no agrupamento! Pensei o seguinte!
V1=tapply(df$V1,df$UF, sum) # Ou trocar pelo mean V2=tapply(df$V2,df$UF, sum) # Ou trocar pelo mean dadostotal=cbind(V1,V2)
hc<-hclust(dist(dadostotal),method='ward.D') plot(as.dendrogram(hc),main = "Cluster Dendrograma") rect.hclust(hc, k = 2, border = "red")
require(bpca) summary(bpca(dadostotal, scale=TRUE)) plot(bpca(dadostotal, scale = T,d=1:2))
resultado=prcomp(dadostotal, scale=T) summary(resultado) biplot(resultado)
Ando estudando um pouco de multivariada sozinho aqui e me veio desde o primeiro poster a pergunta: Se eu teria e suporte técnico para argumentar o que fiz? Já que não consigo analisar com banco de dados repetidos (Com estas duas propostas e que não são as melhores!). Gostaria aqui de ouvir dos que conhecem do assunto, pois sempre tem dados da biologia com esta estrutura e muitas vezes o cara diz tal programa eu faço e eu no R não consigo andar.
obrigado a todos!
André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES
Em Terça-feira, 14 de Abril de 2015 23:10, Éder Comunello < comunello.eder@gmail.com> escreveu:
Andre, boa noite!
Uma ideia que talvez possa ajudar...
### <code r> ### dados fictícios! df <- data.frame(UF=rep(c("AC","AL","AM"), each=4), ano=rep(2003:2006, 3), V1=rnorm(12), V2=rnorm(12)) df
### criando rownames únicos a partir dos dados... rownames(df) <- paste0(df$UF, df$ano) df
PCA <- prcomp(df[3:4],scale=T) ### informar as variáveis númericas hc <- hclust(dist(df[3:4]), method='ward.D') ### </code>
Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.