Lavaan与Mplus的增长曲线建模：能否复刻Mplus中斜率与周得分的关联？

当然可以！Lavaan提供了和Mplus高度对齐的语法逻辑，完全能实现增长曲线模型里斜率因子与特定周观测得分的关联设定——不管是协方差还是直接效应，都能和Mplus的结果完美匹配。下面我结合常见场景拆解具体操作：

1. 核心逻辑对齐 #

Mplus里对增长因子（截距i、斜率s）与周得分（观测变量y1-y4这类）的关联，本质是通过WITH（协方差）或ON（回归效应）语句实现；Lavaan的语法逻辑完全对应，只是用=~替代Mplus的|来定义增长因子，其余关联设定的写法几乎一模一样。

2. 具体场景示例 #

场景1：设定斜率因子与某周得分的协方差 #

假设你的Mplus脚本是这样的（以4个周的重复测量为例）：

MODEL:
  ! 定义线性增长模型，时间编码为0,1,2,3
  i s | y1@0 y2@1 y3@2 y4@3;
  ! 设定斜率s与第2周得分y2的协方差
  s WITH y2;

对应的Lavaan脚本可以完美复刻：

library(lavaan)

# 模型语法
model <- '
  # 定义线性增长因子，和Mplus逻辑一致
  i =~ 1*y1 + 1*y2 + 1*y3 + 1*y4
  s =~ 0*y1 + 1*y2 + 2*y3 + 3*y4
  
  # 设定斜率因子s与第2周得分y2的协方差，对应Mplus的WITH语句
  s ~~ y2
'

# 拟合模型
fit <- cfa(model, data = your_data)
summary(fit, fit.measures = TRUE)

场景2：设定斜率因子对某周得分的直接效应 #

如果你的Mplus脚本是让斜率因子预测某周得分：

MODEL:
  i s | y1@0 y2@1 y3@2 y4@3;
  ! 让斜率s预测第3周的得分y3
  y3 ON s;

对应的Lavaan脚本：

model <- '
  i =~ 1*y1 + 1*y2 + 1*y3 + 1*y4
  s =~ 0*y1 + 1*y2 + 2*y3 + 3*y4
  
  ! 斜率s对y3的直接效应，对应Mplus的ON语句
  y3 ~ s
'
fit <- cfa(model, data = your_data)

3. 关键注意事项 #

• 增长因子负荷设定：Lavaan用=~定义因子负荷，固定负荷的写法和Mplus完全对齐（比如0*y1就是固定y1在斜率因子s上的负荷为0）。
• 协方差与残差：Lavaan中~~既可以表示因子间协方差，也能表示因子与观测变量的协方差，和Mplus的WITH完全等价。
• 结果验证：拟合后可以对比Lavaan和Mplus的拟合指数（如CFI、RMSEA）、参数估计值，确认两者结果完全一致。

只要按照这个对应逻辑转换，就能在Lavaan中实现和Mplus完全相同的斜率与周得分关联设定。

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