R语言CFA单因子拟合出现极端指标问题咨询

问题背景 #

在R中做多因子CFA时，包含A1、A2、A3三个因子的模型拟合指标正常：

A <- 'A1 =~ A_1_1 + A_1_2 + A_1_3 
A2 =~ A_2_1 + A_2_2
A3 =~ A_3_1 + A_3_2 + A_3_3
A_1_1 ~~ A_1_1
A_1_2 ~~ A_1_2
A_1_3 ~~ A_1_3
A_2_1 ~~ A_2_1
A_2_2 ~~ A_2_2
A_3_1 ~~ A_3_1
A_3_2 ~~ A_3_2
A_3_3 ~~ A_3_3'

A_fit <- cfa(A, fixed.x=F, data=data)
summary(A_fit, standardized=T, fit.measures=T)

拟合结果：CMIN=112.793、DF=41.00、P-Value=0.000、RMSEA=0.053、SRMR=0.021、CFI=0.986、TLI=0.981。

但单独对A1因子做CFA时：

A <- 'A1 =~ A_1_1 + A_1_2 + A_1_3 
A_1_1 ~~ A_1_1
A_1_2 ~~ A_1_2
A_1_3 ~~ A_1_3'

A_fit <- cfa(A, fixed.x=F, data=data)
summary(A_fit, standardized=T, fit.measures=T)

得到结果：CMIN=0、DF=0、P-Value=NA、RMSEA=0、SRMR=0、CFI=1、TLI=1，无报错但结果看似异常。

核心问题解答 #

1. 本次操作是否有误？ #

操作完全正确，结果看似异常是模型特性导致的，并非代码错误。

2. 能否获取单个因子的拟合度？ #

无法通过常规拟合优度指标评估这类单因子模型，原因在于：
单因子+3个指标的CFA模型属于恰好识别模型——模型待估参数的数量与观测数据协方差矩阵的元素数量完全相等，自由度（DF）为0。这类模型会完全拟合观测数据，所有拟合优度指标都会显示“完美值”，但这些指标没有实际参考意义，因为拟合优度的本质是对比假设模型与饱和模型的差异，而恰好识别模型本身就是饱和模型的一种，不存在差异可比较。

3. 后续该如何处理？ #

无需纠结拟合优度指标，转而通过信度与效度指标评估单因子结构的合理性：

• 信度指标：
  • Cronbach's α：使用psych包计算，代码示例：

    library(psych)
    alpha(data[, c("A_1_1", "A_1_2", "A_1_3")])
    

  • 组合信度（CR）与平均方差提取（AVE）：使用semTools包直接计算模型信度：

    library(semTools)
    reliability(A_fit)
    

• 效度指标：
  • 查看标准化因子载荷：要求载荷显著（p<0.05）且绝对值通常建议≥0.7，反映观测指标对潜因子的解释力
  • 可搭配探索性因子分析（EFA），确认3个指标是否仅提取出1个核心因子
• 若一定要做模型对比，可将单因子模型与虚无模型（所有指标两两不相关）比较，但这种对比的参考价值有限，核心仍以信度、载荷指标为主。

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