data(lipset)
library(psych)
loa <- unclass(principal(lipset[[1]], nfactors = 3,
rotate = "varimax")$loadings)
flagged <- qflag(loa = loa, nstat = nrow(lipset[[1]]))
summary(flagged)
# Remember to manually inspect the automatic pre-flagging:
results=list(loa=loa, flagged=flagged, brief=list(nfactors = ncol(loa)))
loa.and.flags(results)
Run the code above in your browser using DataLab