Re: [R] plotting other axes for PCA

From: Mark Difford <>
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:
Sent from the R help mailing list archive at

______________________________________________ mailing list
PLEASE do read the posting guide
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 Please read the posting guide before posting to the list.

list of date sections of archive