我正在尝试使用 Stata 手册中的示例 42g 作为指南,在 stata 中构建多级结构方程模型。我的代码正在运行,但我在计算间接影响方面有点偏离书本,可以只使用专业人士的健全性检查。
我正在测试一个有节制的中介模型,其中
Path A = fear <- i.theme + i.condition + i.theme#i.condition + easyseec
Path B = ppt_ave <- fear
Path C' = ppt_ave <- i.theme + i.condition + i.theme#i.condition + easyseec
该模型使用嵌套数据——easyseec 嵌套在主题下,因此我在变量 responseid 下添加了随机截距和随机斜率。
.gsem 输出与 R 中使用 lmer 构建的相同模型的输出大致跟踪,所以我对此感觉很好。但我真的可以在下面的 .nlcom 上使用眼球。这应该计算特定主题/条件对比的间接影响,对吧?
. gsem (ppt_ave <- fear i.theme i.condition i.theme#i.condition easyseec M1[responseid]) (fear <- i.theme i.condition i.theme#i.condition easyseec M2[responseid])
. nlcom _b[ppt_ave:fear]*(_b[fear:2.theme]+_b[fear:2.condition]+_b[fear:2.theme#2.condition]+_b[fear:easysee])