# \donttest{
# the Harman (1967) correlation matrix
PCAoutput <- PCA(data_Harman, Nfactors = 2, Ncases=305, rotate='none', verbose=TRUE)
VARIMAX(PCAoutput$loadingsNOROT, verbose=TRUE)
# Rosenberg Self-Esteem scale items
PCAoutput <- PCA(data_RSE, corkind='polychoric', Nfactors = 2, rotate='none', verbose=TRUE)
VARIMAX(PCAoutput$loadingsNOROT, verbose=TRUE)
# NEO-PI-R scales
PCAoutput <- PCA(data_NEOPIR, Nfactors = 5, rotate='none', verbose=TRUE)
VARIMAX(PCAoutput$loadingsNOROT, verbose=TRUE)
# }
Run the code above in your browser using DataLab