因子得分

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

因子分析中的因子得分如何在R中实现

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：因子得分 因子分析 得分

一个小例子，希望可以帮到你——  

data(USArrests)
fa <- factanal(~., factors = 1, data = USArrests, score = "Bartlett", rotation = "none") #不做旋转  
D <- diag (fa $ uniquenesses)   #特殊方差
  
A <- as.matrix(fa$loadings[,1])  #载荷矩阵

D1 <- solve(D)  #D的逆矩阵
x <- t(as.matrix(USArrests))

#Bartlett方法（最小二乘法）
  
因子得分为f <- solve(t(A) %*% D1 %*% A) %*% t(A) %*% D1 %*% x  
  
然后标准化就是fa$scores了
fa <- factanal(~., factors = 1, data = USArrests, score = "regression", rotation = "none")
r <- fa$correlation #x的相关矩阵
Thompson方法（回归方法）
f <- t(A) %*% solve(r) %*% x