data(WansbeekMeijer, package = "GPArotation")
# Correlation matrix
round(cov2cor(NetherlandsTV$cov), 2)
# factanal picks up n.obs automatically from the list
factanal(factors = 2, covmat = NetherlandsTV, rotation = "none")
# Two-step oblique rotation
fa.unrotated <- factanal(factors = 3, covmat = NetherlandsTV,
rotation = "none")
oblimin(loadings(fa.unrotated), randomStarts = 100)
Run the code above in your browser using DataLab