Lavaan多组SEM分析：单组协变量无方差的解决方法问询

报错原因 #

组2的协变量c取值完全相同，方差为0，lavaan无法估计该变量的回归系数（回归系数计算依赖自变量的变异），因此抛出报错。

实现建议 #

利用lavaan的组特定模型语法，为两组分别定义模型：组1保留c的控制路径，组2移除c的路径（因无变异，控制无意义）。同时可根据需求对组间核心系数（如b对a的效应）设置约束，实现组间比较。

1. 基础拟合代码（无系数约束） #

library(lavaan)

N<-1000
a<-rnorm(N,0,1)
b<-a*0.5+rnorm(N,0,0.2)
c<-c(rnorm((N/2),0,0.2),rep(0,times=(N/2))) 
group<-rep(c("1","2"),each=(N/2))
dtest<-cbind.data.frame(a,b,c,group)

# 组特定模型定义
m <- '
  # 组1：控制c的效应
  group: 1
    a ~ b + c
  # 组2：无c的路径（因c无变异）
  group: 2
    a ~ b
'

fit <- sem(model=m, data=dtest, group="group")
summary(fit, fit.measures=TRUE)

2. 组间系数比较（检验b的效应是否组间有差异） #

若要检验b对a的效应在两组是否相同，可通过模型比较实现：

• 先拟合约束b系数组间相等的模型
• 再拟合无约束的模型
• 用卡方检验比较两个模型的拟合度差异

# 约束模型：组1和组2的b系数相等
m_constrained <- '
  group: 1
    a ~ b + c
  group: 2
    a ~ b
  # 约束组1与组2的b系数一致
  a~b == a~b.g2
'

fit_constrained <- sem(model=m_constrained, data=dtest, group="group")

# 模型比较（卡方差异检验）
anova(fit, fit_constrained)

若卡方检验结果显著，则说明两组b的效应存在差异；反之则无显著差异。

3. 结果解释 #

• 组1的c系数：代表组1年龄相对于组2均值的偏差对a的效应，符合你“控制组1年龄”的需求。
• 组2因年龄无变异，无需控制c，模型仅估计b对a的效应，避免了方差为0的报错。

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