R语言-因子分析中因子综合得分求解

原创

已于 2023-11-25 10:57:31 修改 · 1w 阅读

58 ·

CC 4.0 BY-SA版权

文章标签：

#r语言

于 2021-06-14 12:00:31 首次发布

综合得分求解原理：公因子各自的方差贡献率占累计贡献率的比重为权重加权计算。公式：公因子×权重的和。

由于本文的主要目的是计算综合因子得分。所以对于数据的处理和代码的解读不做过多描述。

#载入数据，求出相关系数矩阵
library(xlsx)
a0 <- read.xlsx("D:/RData/sjhz.xlsx",1,encoding = "UTF-8")
a0.cor <- cov2cor(cov(a0))

library(psych)
fa.parallel(a0.cor,n.obs = 70,fa="fa",n.iter = 100)#碎石图检验得到最佳公共因子个数为2
a0.varimax <- fa(a0.cor,nfactors = 2,n.obs = 70,rotate = "varimax",score=TRUE,fm="ml")#直接进行正交旋转，使用方法为极大似然估计。
a0.varimax
Factor Analysis using method =  ml
Call: fa(r = a0.cor, nfactors = 2, n.obs = 7