text(s3d$xyz.convert(dat$x, dat$y, dat$z), labels=dat$label,
pos=1)
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@gmail.com
> To: listas.c3sl.ufpr.br:
R-br@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@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.