confusion about factor scores from princomp function

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

confusion about factor scores from princomp function

Larry D'Agostino
If I calculate factor scores with princomp() I don't get the same factor
scores using eigen() method.

fctscr1 <- princomp(USArrests, cor=T)$scores
fctscr2 <- scale(as.matrix(USArrests)) %*% eigen(cor(USArrests))$vectors

identical(fctscr1,fctscr2)       # results in FALSE

The values are close but not quite.  Is there something I'm doing wrong?
It's my understanding that princomp() uses eigen() to calculate the factor
loadings.

Larry

        [[alternative HTML version deleted]]

______________________________________________
[hidden email] 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.