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>