我正在使用 psych 包的 fa 命令进行因子分析,因此有一个 fa 类的对象。我可以使用 fac$loadings 查询负载,但是如何捕获负载的下半部分,它给出了 SS 负载以及比例和累积方差。
示例代码:
library(psych)
data(bfi)
fac <- fa(r=cor(bfi, use="complete.obs"), nfactors=5, fm="ml", rotate="none")
fac$loadings
以及如何提取以下信息?
ML1 ML2 ML3 ML4 ML5
SS loadings 4.429 2.423 1.562 1.331 0.966
Proportion Var 0.158 0.087 0.056 0.048 0.034
Cumulative Var 0.158 0.245 0.301 0.348 0.383