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

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@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.

Tiago, Obrigado pela dica. Consegui plotar os gráficos usando o 'FactoMineR', que na verdade usa o scatterplot3d: descobri isto ontem. Baseado no que vi ontem, e na tua dica. Vou explorar o pacote para customizar meus próximos gráficos de PCA. Cordialmente, Vinícius Em 30-09-2013 21:26, Tiago Souza Marçal escreveu:
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@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.
_______________________________________________ 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.
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
participantes (2)
-
Tiago Souza Marçal
-
Vinícius Lionel Mateus