Basically, I managed to write the following:
Code: Select all
efa = fa(data,
fm = "minres",
nfactors = fa.parallel(data, fm = "minres")$nfact,
rotate = "oblimin")
efasort = fa.sort(efa)$loadings
print(efasort, cutoff = 0.3)
However, what's being given in R with the psych package is ever so slightly different from what jamovi (both the app and jmvcore on R) is giving me. Usually, the two results' item loadings would be off by around 0.005 to 0.030.
Is there an option that I'm missing in the psych fa function that's causing the differences?