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

[Éder Comunello]

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 em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]

Em 15 de abril de 2015 00:24, Andre Oliveira <andreolsouza em 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 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úmericas
> hc  <- hclust(dist(df[3:4]), method='ward.D')
> ### </code>
>
> Éder Comunello <c <comunello.eder em gmail.com>omunello.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.
>
>
>
> _______________________________________________
> 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/1eced0a8/attachment.html>