<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><a class="t_atc ICName" id="rmic1_senderName">Vinícius, já usei o pacote scatterplot3d para plotar gráficos de variáveis canônicas em 3D.</a><BR> <BR>Segue um CMR:<BR> <BR>dat <- data.frame(x=1:10, y=1:10, z=1:10, label=letters[1:10])<br>library(scatterplot3d)<br>s3d <- scatterplot3d(dat$x, dat$y, dat$z)<br>text(s3d$xyz.convert(dat$x, dat$y, dat$z), labels=dat$label, pos=1)<BR> <BR>Att.<BR> <BR>Tiago.<BR> <br><br><BR><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div><div> </div><div>################################################################# </div><br> <BR><div>> Date: Thu, 26 Sep 2013 09:24:17 -0300<br>> From: vinynegrelli@gmail.com<br>> To: listas.c3sl.ufpr.br: R-br@listas.c3sl.ufpr.br<br>> Subject: [R-br] Produção de Gráficos: PCA<br>> <br>> Prezados,<br>> <br>> Estou começando a aplicar PCA usando R, mas anteriormente o fazia no <br>> Statistica.<br>> Tenho utilizado os pacotes 'psych' e 'GPArotation' pois prefiro utilizar <br>> rotações em Varimax, porém os gráficos não ficaram ao meu gosto.<br>> <br>> Depois de algum esforço, eu consegui trocar as labels do meu gráfico, <br>> mas além disto gostaria de retirar o histograma, e plotar gráficos da <br>> dispersão das amostras(labels) e das variáveis entres os componentes, em <br>> 2D e 3D.<br>> <br>> Segue meu CMR.<br>> <br>> Desde já, obrigado.<br>> <br>> <br>> # instalando os pacotes<br>> install.packages("psych")<br>> install.packages("GPArotation")<br>> #carregando<br>> require (psych)<br>> install.packages (GPArotation)<br>> #entrada de dados (https://www.dropbox.com/s/l8v9g6qub27oig0/enqalog.xls)<br>> read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> d<br>> #entrada de dados1<br>> read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> enqa<br>> enqa -> d<br>> attach(d)<br>> names(d)<br>> summary(d)<br>> #PCA(Varimax)/psych<br>> dpca<-principal(d[,2:16], nfactors=3, rotate="varimax", scores=T)<br>> #Scree plot<br>> plot (dpca$values, type="b", ylab="Eigenvalue", xlab="Component", <br>> lab=c(5,5,5))<br>> #biplot<br>> biplot(dpca, labels=d[,1], cex=c(0.75,1))<br>> <br>> Vinícius<br>> --<br>> -- <br>> Atenciosamente,<br>> <br>> VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665)<br>> Bacharel em Química - Doutorando em Química Analítica<br>> Laboratório de Química Atmosférica - Departamento de Química<br>> Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio)<br>> Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900<br>> Telefone: (+55) (21) 3527-1327<br>> (+55) (21) 9358-8051<br>> www.puc-rio.br <br>> <br>> _______________________________________________<br>> R-br mailing list<br>> R-br@listas.c3sl.ufpr.br<br>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<br>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.<br></div> </div></body>
</html>