Hi He Zhang, >> Is the following right for extracting the scores? >> ... >> pca$loadings >> pca$score
Yes. But you should be aware that the function principal() in package psych is standardizing your data internally, which you might not want. That is, the analysis is being based on the correlation matrix rather than on the covariance matrix. The two analyses (and the "default" biplots that issue from them) have somewhat different properties. You should know about these. Regards, Mark. -- View this message in context: http://r.789695.n4.nabble.com/pca-analysis-extract-rotated-scores-tp3065061p3067370.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.