I can't find my PAF Factor loadings < .3

May 19, 2010 8:14 PM

Statistics Filter: After performing factor analysis (PAF) using R, how do I display all the factor loadings, not just the ones above .3?

This is my first excursion with R, so I'm pretty much a novice. I'm no biostatistician either, please be kind.

Using the 'psych' package, I did the following factor analysis.

fa(mydata, nfactors=2, residuals = TRUE, rotate = "promax",n.obs = NA, scores = TRUE,SMC=TRUE, missing=TRUE,impute="median",min.err = 0.001,max.iter = 50,symmetric=TRUE,warnings=TRUE,fm="pa")

It worked great, but now I need to see all the factor loadings, not just those above .3.

So I've tried all manner of things and was thinking that the 'cut' argument in the print.psych command would be useful. But alas the fa command cut defaults to .3.

Am I doing something wrong? Is there anything I can do?

posted by bingoes to Science & Nature (2 answers total)

Best answer: There's probably a more elegant way to do this, but the following should show you all the loadings.

fa(mydata, nfactors=2, residuals = TRUE, rotate = "promax",n.obs = NA, scores = TRUE,SMC=TRUE, missing=TRUE,impute="median",min.err = 0.001,max.iter = 50,symmetric=TRUE,warnings=TRUE,fm="pa") $ loadings[,]

Let me know if there are any problems with the syntax (uploading R code via MeFi can sometimes get mangled).
posted by eisenkr at 9:31 PM on May 19, 2010 [1 favorite]

Response by poster: You've made me so happy. Thank you.
posted by bingoes at 9:38 PM on May 19, 2010