Larry D'Agostino
2011-Dec-07 20:40 UTC
[R] confusion about factor scores from princomp function
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]]