data(lipset)
library(psych)
loa <- as.data.frame(unclass(principal(lipset,
nfactors=3, rotate="varimax")$loadings))
flagged <- qflag(nqsorts=9, nstat=33, loa=loa)
qmzsc <- qzscores(lipset, nstat=33, nfactors=3, nqsorts=9,
flagged=flagged, loa=loa,
rotation="varimax")
qmzsc
Run the code above in your browser using DataLab