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 <comunello.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 <comunello.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.