[R-br] Análise de Fatores PCA e Dendograma

Andre Oliveira andreolsouza em yahoo.com.br
Quarta Abril 15 01:24:56 BRT 2015


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 meanV2=tapply(df$V2,df$UF, sum)  # Ou trocar pelo meandadostotal=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 em 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úmericashc  <- hclust(dist(df[3:4]), method='ward.D')### </code>
Éder Comunello <comunello.eder em gmail.com> 
Dourados, MS - [22 16.5'S, 54 49'W]



_______________________________________________
R-br mailing list
R-br em 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.

  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150415/b02305a1/attachment-0001.html>


Mais detalhes sobre a lista de discussão R-br