<div dir="ltr"><div>Sem entrar no mérito da análise, outra ideia que parece atender sua necessidade...</div><div><br></div><div><font face="monospace, monospace">### <code r></font></div><div><font face="monospace, monospace">### dados fictícios!</font></div><div><font face="monospace, monospace">df <- data.frame(UF=rep(c("AC","AL","AM","MS","PR"), each=4), ano=rep(2003:2006, 5), </font></div><div><font face="monospace, monospace"> V1=rnorm(20), V2=rnorm(20), v3=rnorm(20), v4=rnorm(20))</font></div><div><font face="monospace, monospace">df</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">df2 <- aggregate(.~UF, data=df, mean)</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">### criando rownames únicos a partir dos dados...</font></div><div><font face="monospace, monospace">rownames(df2) <- df2$UF</font></div><div><font face="monospace, monospace">df2</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">PCA <- prcomp(df[3:6],scale=T) ### informar as variáveis númericas</font></div><div><font face="monospace, monospace">hc <- hclust(dist(df[3:6]), method='ward.D')</font></div><div><font face="monospace, monospace">### </code></font></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div></div>
<br><div class="gmail_quote">Em 15 de abril de 2015 00:24, Andre Oliveira <span dir="ltr"><<a href="mailto:andreolsouza@yahoo.com.br" target="_blank">andreolsouza@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:bookman old style,new york,times,serif;font-size:16px"><div dir="ltr"><span>Exatamente Éder,</span></div><div dir="ltr"><span>você chegou exatamente onde eu gostaria, mas continuando seu código ao plotar plot(as.dendrogram(hc),main = "")</span> cria os grupos e separa os <span style="font-family:monospace,monospace;font-size:13px">ACs, os Als e os AMs no dendograma e eu gostaria de apenas um Estado no agrupamento! Pensei o seguinte! </span></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px"><br></span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">V1=tapply(df$V1,df$UF, sum) # Ou trocar pelo mean</span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">V2=tapply(df$V2,df$UF, sum) # Ou trocar pelo mean</span></font></div><div dir="ltr"><span style="font-family:monospace,monospace;font-size:13px">dadostotal=cbind(V1,V2)</span><br></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px"><br></span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">hc<-hclust(dist(dadostotal),method='ward.D')</span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">plot(as.dendrogram(hc),main = "Cluster Dendrograma")</span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">rect.hclust(hc, k = 2, border = "red")</span></font></div><div dir="ltr"><br></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">require(bpca) </span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">summary(bpca(dadostotal, scale=TRUE))</span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">plot(bpca(dadostotal, scale = T,d=1:2))</span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px"><br></span></font></div><div dir="ltr"><span style="font-size:13px;font-family:monospace,monospace">resultado=prcomp(dadostotal, scale=T)</span><br></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">summary(resultado)</span></font></div><div dir="ltr"><font face="monospace, monospace"><span style="font-size:13px">biplot(resultado)</span></font></div><div dir="ltr"><br></div><div dir="ltr">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. </div><div dir="ltr"><br></div><div dir="ltr">obrigado a todos! </div><span class=""><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr"><span style="font-family:monospace,monospace;font-size:13px"><br></span></div><div dir="ltr"><span style="font-family:monospace,monospace;font-size:13px"><br></span></div><div dir="ltr"><span style="font-family:monospace,monospace;font-size:13px"><br></span></div><div dir="ltr"><br></div><div></div><div> </div><div><div>
<br>
<br>André Oliveira Souza.
<br>Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES
<br>
<br></div></div> <br><div><br><br></div></span><div style="display:block"> <div style="font-family:bookman old style,new york,times,serif;font-size:16px"> <div style="font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"><div><div class="h5"> <div dir="ltr"> <font size="2" face="Arial"> Em Terça-feira, 14 de Abril de 2015 23:10, Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>> escreveu:<br> </font> </div> <br><br> </div></div><div><div><div class="h5"><div><div><div dir="ltr"><div>Andre, boa noite!</div><div><br clear="none"></div>Uma ideia que talvez possa ajudar...<div><br clear="none"><div><div><font face="monospace, monospace">### <code r></font></div><div><font face="monospace, monospace">### dados fictícios!</font></div><div><font face="monospace, monospace">df <- data.frame(UF=rep(c("AC","AL","AM"), each=4), ano=rep(2003:2006, 3), V1=rnorm(12), V2=rnorm(12))</font></div><div><font face="monospace, monospace">df</font></div><div><font face="monospace, monospace"><br clear="none"></font></div><div><font face="monospace, monospace">### criando rownames únicos a partir dos dados...</font></div><div><font face="monospace, monospace">rownames(df) <- paste0(df$UF, df$ano)</font></div><div><font face="monospace, monospace">df</font></div><div><font face="monospace, monospace"><br clear="none"></font></div><div><font face="monospace, monospace">PCA <- prcomp(df[3:4],scale=T) ### informar as variáveis númericas</font></div><div><font face="monospace, monospace">hc <- hclust(dist(df[3:4]), method='ward.D')</font></div></div><div><font face="monospace, monospace">### </code></font></div><div><br clear="all"><div><div><div dir="ltr">Éder Comunello <<a rel="nofollow" shape="rect" href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a rel="nofollow" shape="rect" href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br clear="none">Dourados, MS - [22 16.5'S, 54 49'W]<div><br clear="none"></div></div></div></div><div>
<br clear="none"><div><br clear="none"></div></div></div></div></div></div></div><br></div></div><span class=""><div>_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div><br><br></span></div> </div> </div> </div></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>