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]] ______________________________________________ R-help@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.