data(ability.cov)
x <- factanal(factors = 3, covmat = ability.cov, rotation="none")
fssT.df(x$loadings, kij = 2)
fssQ.df(x$loadings, kij = 4)
# 3 different methods
data("WansbeekMeijer", package="GPArotation")
fa.unrotated <- factanal(factors = 3, covmat=NetherlandsTV, rotation="none")
#
fa.varimax <- GPForth.df(loadings(fa.unrotated), method = "varimax", normalize = TRUE)
fa.cubimax <- cubimax.df(loadings(fa.unrotated), normalize = TRUE)
fa.quartimax <- GPForth.df(loadings(fa.unrotated), method = "quartimax", normalize = TRUE)
print(cbind(loadings(fa.varimax), loadings(fa.cubimax), loadings(fa.quartimax)), digits = 2)
Run the code above in your browser using DataLab