[R-br] Produção de Gráficos: PCA

Tiago Souza Marçal tiagosouzamarcal em hotmail.com
Segunda Setembro 30 21:26:56 BRT 2013


Vinícius, já usei o pacote scatterplot3d para plotar gráficos de variáveis canônicas em 3D.
 
Segue um CMR:
 
dat <- data.frame(x=1:10, y=1:10, z=1:10, label=letters[1:10])
library(scatterplot3d)
s3d <- scatterplot3d(dat$x, dat$y, dat$z)
text(s3d$xyz.convert(dat$x, dat$y, dat$z), labels=dat$label, pos=1)
 
Att.
 
Tiago.
  


################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ################################################################# 
 
> Date: Thu, 26 Sep 2013 09:24:17 -0300
> From: vinynegrelli em gmail.com
> To: listas.c3sl.ufpr.br: R-br em listas.c3sl.ufpr.br
> Subject: [R-br] Produção de Gráficos: PCA
> 
> Prezados,
> 
> Estou começando a aplicar PCA usando R, mas anteriormente o fazia no 
> Statistica.
> Tenho utilizado os pacotes 'psych' e 'GPArotation' pois prefiro utilizar 
> rotações em Varimax, porém os gráficos não ficaram ao meu gosto.
> 
> Depois de algum esforço, eu consegui trocar as labels do meu gráfico, 
> mas além disto gostaria de retirar o histograma, e plotar gráficos da 
> dispersão das amostras(labels) e das variáveis entres os componentes, em 
> 2D e 3D.
> 
> Segue meu CMR.
> 
> Desde já, obrigado.
> 
> 
> # instalando os pacotes
> install.packages("psych")
> install.packages("GPArotation")
> #carregando
> require (psych)
> install.packages (GPArotation)
> #entrada de dados (https://www.dropbox.com/s/l8v9g6qub27oig0/enqalog.xls)
> read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> d
> #entrada de dados1
> read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> enqa
> enqa -> d
> attach(d)
> names(d)
> summary(d)
> #PCA(Varimax)/psych
> dpca<-principal(d[,2:16], nfactors=3, rotate="varimax", scores=T)
> #Scree plot
> plot (dpca$values, type="b", ylab="Eigenvalue", xlab="Component", 
> lab=c(5,5,5))
> #biplot
> biplot(dpca, labels=d[,1], cex=c(0.75,1))
> 
> Vinícius
> --
> -- 
> Atenciosamente,
> 
> VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665)
> Bacharel em Química - Doutorando em Química Analítica
> Laboratório de Química Atmosférica - Departamento de Química
> Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio)
> Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900
> Telefone: (+55) (21) 3527-1327
>                (+55) (21) 9358-8051
> www.puc-rio.br   
> 
> _______________________________________________
> R-br mailing list
> R-br em 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.
 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131001/531aeabe/attachment.html>


Mais detalhes sobre a lista de discussão R-br