Lavaan多组结构方程模型标准误计算失败及协方差矩阵非正定求助

核心问题诊断 #

1. 多组比较的操作逻辑错误：你用于比较的mod.mg1a和mod.mg1b语法完全一致，没有设置任何跨组参数约束的差异，二者自由度完全相等，所以运行anova()时会直接报"部分模型自由度相同"的警告，无法完成似然比检验。
2. 模型未识别：你收到的"信息矩阵无法逆、vcov矩阵非正定"警告本质都是模型未识别，诱因有三点：

• lavaan默认fixed.x = TRUE，会将外显自变量warm、crit的方差、截距固定为样本观测值，跨组分析时该默认设置会导致自动生成的截距被强制设为0，触发对应警告。
• 你没有遵循测量不变性的检验逻辑，一次性将所有跨组参数全部约束相等，超过了模型可识别的自由度阈值，直接导致信息矩阵不可逆。
• 未做数据前置检查，若自变量共线性过高、题项存在完美相关、样本量远小于待估参数数量，都会触发非正定警告。

具体修正方案 #

1. 调整多组比较逻辑，按不变性检验顺序逐步添加约束：
   首先构建无跨组约束的构型不变性基准模型，再依次添加载荷、截距、残差、结构参数的约束，每次仅比较相邻两个嵌套模型：

# 通用模型语法，无需手动声明方差，lavaan会默认估计
mod <- "
# 结构路径
trust ~ warm + crit
# 测量模型
trust =~ never + moral + han + good
# 显式指定估计外显自变量的截距，避免自动置0
warm ~ 1
crit ~ 1
"

# 1. 构型不变性模型：无跨组约束
fit.config <- lavaan(mod, data = mach, group = "gen", fixed.x = FALSE)
# 先检查基准模型的适配度，适配合格再往下操作
summary(fit.config, fit.measures = T, standardized = T)

# 2. 弱不变性：约束因子载荷跨组相等
fit.metric <- lavaan(mod, data = mach, group = "gen", 
                     group.equal = "loadings", fixed.x = FALSE)
# 似然比检验，若p>0.05说明载荷约束可接受
anova(fit.config, fit.metric)

# 3. 强不变性：在载荷约束基础上，加测量截距跨组相等
fit.scalar <- lavaan(mod, data = mach, group = "gen", 
                     group.equal = c("loadings", "intercepts"), fixed.x = FALSE)
anova(fit.metric, fit.scalar)

# 4. 后续再依次检验严格不变性、结构参数约束，每一步验证通过再向下进行

2. 前置检查数据质量：

• 核查warm、crit的相关系数，若大于0.8需处理共线性问题
• 核查所有观测变量的相关矩阵，是否存在相关系数大于0.9的题项
• 确认样本量至少达到待估参数数量的10倍以上
• 核查缺失值是否已做合理插补或删除

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