<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
André,<br>
<br>
Assim como o Cesar comentou, não consegui identificar os data.frame
por ano que você menciona.<br>
Agora se você tem esses data.frames por ano, será que uma análise de
agrupamento ou PCA é o ideal? Uma abordagem longitudinal não seria
melhor? E também não está claro o que deseja. É calcular médias
anuais?<br>
O pacote prcomp() extrai sim as PCA's <br>
<br>
No seu comando se colocar o argumento scale=T, ele vai criar as
PCA's de acordo com a matriz de correlação. Veja:<br>
PCA <- prcomp(USArrests,scale=T)<br>
summary(PCA)<br>
names(PCA)<br>
PCA$sdev #### Ríaz quadrada dos autovalores da matriz de
correlação<br>
PCA$rotation #### Matriz de autovetores normalizados<br>
<br>
Comparando com os autovalores e autovetores da matriz de correlação
terá valores como os extraídos ao acima.<br>
<br>
mcor <- cor(USArrests)<br>
mcor<br>
vv <- eigen(mcor)<br>
vv<br>
<br>
Logo desta forma terá suas PCA's. Por exemplo a PCA1:<br>
Y1
<-vv$vectors[1,1]*USArrests[,1]+vv$vectors[2,1]*USArrests[,2]+vv$vectors[3,1]*USArrests[,3]+vv$vectors[4,1]*USArrests[,4]<br>
Y1<br>
roY1x1 <- (vv$vectors[1,1]*sqrt(vv$values[1]))/sqrt(mvar[1,1]) ##
Correlação da PCA1 com Murder<br>
roY1x1<br>
roY1x2 <- (vv$vectors[2,1]*sqrt(vv$values[1]))/sqrt(mvar[2,2]) ##
Correlação da PCA1 com Assault<br>
roY1x2<br>
roY1x3 <- (vv$vectors[3,1]*sqrt(vv$values[1]))/sqrt(mvar[3,3]) ##
Correlação da PCA1 com UrbanPop<br>
roY1x3<br>
roY1x4 <- (vv$vectors[4,1]*sqrt(vv$values[1]))/sqrt(mvar[4,4]) ##
Correlação da PCA1 com Rape<br>
roY1x4<br>
<br>
<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">On 12-04-2015 22:00, Andre Oliveira
wrote:<br>
</div>
<blockquote
cite="mid:606918295.1445473.1428886857730.JavaMail.yahoo@mail.yahoo.com"
type="cite">
<div style="color:#000; background-color:#fff; font-family:bookman
old style, new york, times, serif;font-size:16px">
<div id="yui_3_16_0_1_1428885607258_5993"><span
style="font-family: 'Helvetica Neue', 'Segoe UI', Helvetica,
Arial, 'Lucida Grande', sans-serif; font-size: 13px;"
class="" id="yui_3_16_0_1_1428885607258_6318">Cesar Rabak,</span><span></span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
style="font-family: 'Helvetica Neue', 'Segoe UI', Helvetica,
Arial, 'Lucida Grande', sans-serif; font-size: 13px;"
class="" id="yui_3_16_0_1_1428885607258_6319">obrigado pela
resposta. Esta é a questão como organizar os três bancos de
dados levando em conta a variável ano, pois, pelo que tenho
olhado as funções </span>hclust() e prcomp() não extrai os
PCAs e nem o dendograma replicando ano. Todos os exemplos do R
trabalham com dados de um ano apenas. O objetivo é agrupar as
cidades e os bancos de dados <span style="font-family:
'Helvetica Neue', 'Segoe UI', Helvetica, Arial, 'Lucida
Grande', sans-serif; font-size: 13px;" class=""
id="yui_3_16_0_1_1428885607258_10868">USArrests2007,
USArrests2009 e USArrests2010 são fictícios e estão ai
apenas para ilustrar que tenho 3 bancos com mesma estrutura
do </span><span class=""
id="yui_3_16_0_1_1428885607258_10868" style="font-family:
'Helvetica Neue', 'Segoe UI', Helvetica, Arial, 'Lucida
Grande', sans-serif; font-size: 13px;">USArrests. </span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"><br>
</span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><font
id="yui_3_16_0_1_1428885607258_14979" face="Helvetica Neue,
Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span
style="font-size: 13px;"
id="yui_3_16_0_1_1428885607258_14978">Dado que tenho
apenas o </span></font><span class=""
id="yui_3_16_0_1_1428885607258_10868" style="font-family:
'Helvetica Neue', 'Segoe UI', Helvetica, Arial, 'Lucida
Grande', sans-serif; font-size: 13px;">USArrests, tudo pode
ser resolvido com o CMR com pequenas variações de escala e
padronização das variáveis se desejar.</span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"><br>
</span></div>
<div id="yiv6363743822yui_3_16_0_1_1428692662702_21288"
dir="ltr" class="" style=""><span class=""
id="yiv6363743822yui_3_16_0_1_1428692662702_21091" style="">hc=hclust(dist(</span><span
style="font-family: 'Helvetica Neue', 'Segoe UI', Helvetica,
Arial, 'Lucida Grande', sans-serif; font-size: 13px;"
class="" id="yui_3_16_0_1_1428885607258_13872">USArrests)</span>)</div>
<div id="yiv6363743822yui_3_16_0_1_1428692662702_21288"
dir="ltr" class="" style="">plot(hc)</div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"></span></div>
<div id="yiv6363743822yui_3_16_0_1_1428692662702_21303"
dir="ltr" class="" style=""><span class=""
id="yiv6363743822yui_3_16_0_1_1428692662702_21091" style="">summary(prcomp(</span><span
style="font-family: 'Helvetica Neue', 'Segoe UI', Helvetica,
Arial, 'Lucida Grande', sans-serif; font-size: 13px;"
class="" id="yui_3_16_0_1_1428885607258_13873">USArrests</span>))</div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"><br>
</span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;" id="yui_3_16_0_1_1428885607258_12964">Obrigado </span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"><br>
</span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"><br>
</span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><span
class="" style="font-family: 'Helvetica Neue', 'Segoe UI',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-size:
13px;"><br>
</span></div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><br>
</div>
<div id="yui_3_16_0_1_1428885607258_5993" dir="ltr"><br>
</div>
<div id="yui_3_16_0_1_1428885607258_5992"> </div>
<div id="yui_3_16_0_1_1428885607258_5983">
<div id="yui_3_16_0_1_1428885607258_5982">
<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 class="qtdSeparateBR"><br>
<br>
</div>
<div class="yahoo_quoted" 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 dir="ltr"> <font face="Arial" size="2"> Em Domingo,
12 de Abril de 2015 18:14, Cesar Rabak
<a class="moz-txt-link-rfc2396E" href="mailto:cesar.rabak@gmail.com"><cesar.rabak@gmail.com></a> escreveu:<br>
</font> </div>
<br>
<br>
<div class="y_msg_container">
<div id="yiv2538766343">
<div>
<div dir="ltr">Embora você cite os procedimentos de
PCA e plotar dendrograma e "análise dos fatores",
você não dá muita informação de porque você
precisa fazer essa análise e a que conclusões
deseja chegar ou quais hipóteses avaliar.
<div><br clear="none">
</div>
<div>Não sei onde estão os dataframes
USArrests2007, USArrests2009 e USArrests2010,
tendo nos meus defaults apenas o USArrests.
Considerando este último, vê-se que a estrutura
dele é 50 linhas com cada uma um estado
americano por quatro variáveis (três crimes e
uma com a população), se você for anexar bancos
com estatísticas de diferentes anos você terá
que considerar como colocar a variável ano em
jogo e como analisar os estados americanos para
manter apenas as cinquenta linhas...</div>
<div><br clear="none">
</div>
<div>A análise de componentes principais já faz
todos os cálculos necessários, assim calcular as
médias em separado só faria sentido para alguma
análise descritiva (incidentalmente você pode
[dependendo do pacote que use para fazer a PCA]
escolher se a análise será centrada na média de
cada variável ou não).<br clear="none">
<div><br clear="none">
</div>
<div>Quanto aos comandos R que você colocou,
eles não formam um CMR e por isso não dá para
comentar a respeito.</div>
</div>
<div><br clear="none">
</div>
<div>HTH</div>
<div>--</div>
<div>Cesar Rabak</div>
</div>
<div class="yiv2538766343gmail_extra"><br
clear="none">
<div class="yiv2538766343gmail_quote">2015-04-10
16:39 GMT-03:00 Andre Oliveira <span dir="ltr"><<a
moz-do-not-send="true" rel="nofollow"
shape="rect"
ymailto="mailto:andreolsouza@yahoo.com.br"
target="_blank"
href="mailto:andreolsouza@yahoo.com.br">andreolsouza@yahoo.com.br</a>></span>:<br
clear="none">
<blockquote class="yiv2538766343gmail_quote"
style="margin:0 0 0 .8ex;border-left:1px #ccc
solid;padding-left:1ex;">
<div class="yiv2538766343yqt5422653120"
id="yiv2538766343yqt35273">
<div>
<div
style="color:#000;background-color:#fff;font-family:bookman
old style, new york, times,
serif;font-size:16px;">
<div dir="ltr"><span>Pessoal boa tarde,</span></div>
<div dir="ltr"><span>estou com dois
bancos de dados com a estrutura do
USArrests que vem no R só que
repetido no tempo. USArrests2007,
USArrests2009, USArrests2010.
Gostaria de orientação de como
anexar estes dados para em um banco
de dados único e rodar PCA, Plotar,
dendograma e a análise de fatores.
Pelo que percebi em meus devo fazer
o seguinte. Montar um banco de
dados com os três bancos e tirar a
soma ou média. Não sei se fere
princípios da análise. <br
clear="none">
</span></div>
<div dir="ltr"><br clear="none">
</div>
<div dir="ltr"><span>dados=
data.frame(USArrests2007,
USArrests2009, USArrests2010</span>)
# Não sei se é data.frame aqui mesmo!
<br clear="none">
<span></span></div>
<div dir="ltr"><span><br clear="none">
</span></div>
<div dir="ltr"><span>V1=tapply(V1,grupos,mean)
<br clear="none">
</span></div>
<div dir="ltr"><span>.</span></div>
<div dir="ltr"><span>.</span></div>
<div dir="ltr"><span>.</span></div>
<div dir="ltr"><span>.</span></div>
<div dir="ltr"><span>.</span></div>
<div dir="ltr"><span>Vn=tapply(Vn,grupos,mean) </span></div>
<div dir="ltr"><span><br clear="none">
</span></div>
<div dir="ltr"><span>dadosmedio=cbind(V1,
........................,Vn)</span></div>
<div dir="ltr"><span>hclust(dist(dadosmedio)</span></div>
<div dir="ltr"><span>prcomp(dadosmédio)</span></div>
<div dir="ltr"><br clear="none">
</div>
<div dir="ltr">Ou usar a estrutura com
<span>Vn=tapply(Vn,grupos,sum) </span>em
vez de <span>Vn=tapply(Vn,grupos,mean).
</span>Ou usar o próprio banco formado
por <span>USArrests2007,
USArrests2009, USArrests2010. </span><br
clear="none">
<span></span></div>
<div dir="ltr"><br clear="none">
<span></span></div>
<div dir="ltr"><span>Caso falei
besteira, desconsidere. Ficarei
agradecido por sugestões. <br
clear="none">
</span></div>
<div dir="ltr"><span><br clear="none">
</span></div>
<div dir="ltr"><span><br clear="none">
</span></div>
<div dir="ltr"><span>obrigado pela
atenção<br clear="none">
</span></div>
<div dir="ltr"><br clear="none">
<span></span></div>
<div dir="ltr"><span><br clear="none">
</span></div>
<div dir="ltr"><br clear="none">
<span></span></div>
<div dir="ltr"><span><br clear="none">
</span></div>
<div dir="ltr"><br clear="none">
<span></span></div>
<div dir="ltr"><span></span></div>
<div> </div>
<div>
<div>
<br clear="none">
<br clear="none">
André Oliveira Souza. <br
clear="none">
Graduação em Matemática, mestrado em
estatística aplicada.Instituto
Federal de Educação, Ciência e
Tecnologia do Espirito Santo. IFES
<br clear="none">
<br clear="none">
</div>
</div>
</div>
</div>
</div>
<br clear="none">
_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a moz-do-not-send="true" rel="nofollow"
shape="rect"
ymailto="mailto:R-br@listas.c3sl.ufpr.br"
target="_blank"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br
clear="none">
<a moz-do-not-send="true" rel="nofollow"
shape="rect" target="_blank"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br
clear="none">
Leia o guia de postagem (<a
moz-do-not-send="true" rel="nofollow"
shape="rect" target="_blank"
href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>)
e forneça código mínimo reproduzível.<br
clear="none">
</blockquote>
</div>
<br clear="none">
</div>
</div>
</div>
<br>
<br>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>