我使用psych package的fa命令进行因子分析,因此有一个fa类的对象。我可以使用fac $ load来查询载荷,但是如何捕获载荷的下半部分,这会给出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
答案 0 :(得分:2)
运行str(fac)
后,您似乎可以使用
df <- fac$Vaccounted
> df
ML1 ML2 ML3 ML4 ML5
SS loadings 4.4291169 2.42278894 1.5623383 1.3308373 0.96561620
Proportion Var 0.1581827 0.08652818 0.0557978 0.0475299 0.03448629
Cumulative Var 0.1581827 0.24471092 0.3005087 0.3480386 0.38252492
Proportion Explained 0.4135227 0.22620272 0.1458671 0.1242531 0.09015437
Cumulative Proportion 0.4135227 0.63972545 0.7855925 0.9098456 1.00000000