Skip to content
Prev 13700 / 398502 Next

Request for Help: Rotation of PCA Solution or Eigenvectors

On Fri, 21 Sep 2001, Jonathan Baron wrote:

            
You can do it, but whether it is meaningful is another question. Christian
Hoffmann has sent me examples.  He just calls princomp, and then calls
varimax on the loadings.  As an example (not his)

library(mva)
data(USArrests)
pc.cr <- princomp(USArrests, cor=TRUE)
loadings(pc.cr)
varimax(loadings(pc.cr)[, 1:2])  # rotate the first two components

I don't believe that is the best way to use rotation in PCA (although it is
commonly used).  There is a discussion in the on-line Statistical
Complements to MASS3 (via http://www.stats.ox.ac.uk/pub/MASS3)
with recommendations from various sources and S-PLUS examples (that
would work with small changes in R).  Caveat: I am not an expert and the
authorities disagree.   The issues seem much more contentious than
using rotation in factor analysis.