Re: [R] plotting other axes for PCA

From: Mark Difford <mark_difford_at_yahoo.co.uk>
Date: Fri, 18 Jan 2008 10:11:45 -0800 (PST)

Hi Silvia,

>> What I need is exactly what I get using biplot (pca.object) but for other >> axes.

You need to look at ?biplot.prcomp (argument: choices=)

## Try

biplot(prcomp(USArrests), choices=c(1,2))     ## plots ax1 and ax2
biplot(prcomp(USArrests), choices=c(1,3))     ## plots ax1 and ax3

Also look at package ade4 (which in my view still sets the standard in R)

## Example
T.pca <- dudi.pca(USArrests, scannf=F, nf=4)

scatter(T.pca, clab.r=0.5, clab.c=.8)
scatter(T.pca, clab.r=0.5, clab.c=.8, xax=1, yax=3, posie="none")
add.scatter.eig(T.pca$eig, T.pca$nf, "bottomright", xax=1, yax=3)

## Exclude row labels
scatter(T.pca, clab.r=0, clab.c=.8)

HTH, Mark.

Silvia Lomascolo wrote:
>
> Hi R-community,
> I am doing a PCA and I need plots for different combinations of axes
> (e.g., PC1 vs PC3, and PC2 vs PC3) with the arrows indicating the loadings
> of each variables. What I need is exactly what I get using biplot
> (pca.object) but for other axes.
>
> I have plotted PC2 and 3 using the scores of the cases, but I don't get
> the arrows proportional to the loadings of each variables on each axis.
>
> Any hints on this?
>
> Thanks, Silvia.
>

-- 
View this message in context: http://www.nabble.com/plotting-other-axes-for-PCA-tp14952955p14955292.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Fri 18 Jan 2008 - 18:22:15 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 18 Jan 2008 - 18:30:07 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive