Stata factor主因子分析在R中的等价实现及语法报错问题

问题场景 #

你需要在R中复现Stata里factor+score命令的主因子分析功能，Stata示例代码如下：

factor v1 v2 v3 v4 v5 
score newvar1

你尝试用R的psych包fa()函数实现，但语法理解有误，用mtcars数据集测试的代码如下：

mtcars <- mtcars
library('psych')
mtcars$newvar1 <- fa(mtcars, values=c(mpg, cyl, hp))

运行后出现报错：

Error in $<-.data.frame(*tmp*, fa, value = list(residual = c(0.124759140589603, :
replacement has 49 rows, data has 32
In addition: Warning message:
In fa.stats(r = r, f = f, phi = phi, n.obs = n.obs, np.obs = np.obs, :
The estimated weights for the factor scores are probably incorrect. Try a different factor score estimation method.

问题原因与解决方法 #

你的代码存在3个核心问题：

• 参数误用：fa()的values参数是用来传入初始因子载荷矩阵的，不是指定分析变量的；指定分析变量应直接传入变量子集。
• 未提取得分：fa()返回的是包含分析结果的列表，不是直接的得分向量，需提取列表中的得分元素才能赋值给数据框列。
• 方法与因子数不匹配：Stata的factor默认是主因子分析，且score newvar1默认生成单因子得分，你需要在fa()中明确指定对应参数。

正确的R代码示例 #

library(psych)

# 提取要做因子分析的变量子集
analysis_vars <- mtcars[, c("mpg", "cyl", "hp")]

# 执行主因子分析：指定1个因子，用最小残差法（匹配Stata默认逻辑），并计算得分
fa_result <- fa(analysis_vars, nfactors = 1, fm = "minres", score = TRUE)

# 将因子得分添加到原数据框
mtcars$newvar1 <- fa_result$scores

关键参数说明 #

• nfactors=1：指定提取1个因子，和Statascore newvar1生成单因子得分的行为一致。
• fm="minres"：指定主因子分析方法，匹配Statafactor的默认算法。
• score=TRUE：触发因子得分计算，结果会存储在返回列表的$scores元素中。

内容的提问来源于stack exchange，提问作者tchoup