0

我正在以 ocd 症状作为预测因子和抑郁 (mdd) 症状作为结果进行调解分析。我在 R 中使用 lavaan。我有来自 3 个时间点的数据，我想控制时间 1 mdd 和时间 2 mdd 症状。这看起来对吗？我不确定我是否正确地包含了两个协变量。另外，我使用 ocd 和 mdd 量表中的每个单独项目作为指标是否正确？如果我这样做，我应该使用 t1 和 t2 mdd 的单个项目作为协变量，而不是使用总数的项目吗？

med.ocd.sr.mdd.cov6 <- '# define predictor
        t1.pre =~ t1.ocd1 + t1.ocd2 + t1.ocd3 + t1.ocd4 + t1.ocd5 + t1.ocd6 + t1.ocd7 + t1.ocd8 + t1.ocd9 + t1.ocd10
        # define outcome
        t3.out =~ t3.mdd1+t3.mdd2+t3.mdd3+t3.mdd4+t3.mdd5+t3.mdd6
        # define mediator
        t2.med =~ t2.med1 + t2.med2 + t2.med3
        # define covariate 
        cov1 =~ t1.mdd.total
        # define covariate 
        cov2 =~ t2.mdd.total
        # direct effect (c path)
        t3.out ~ c*t1.pre + cov1 + cov2
        # mediator (a*b path)
        t2.med ~ a*t1.pre
        t3.out ~ b*t2.med
        # indirect effect (a*b path)
        ab := a*b
        # total effect
        abc := c + a*b’'
f.med.ocd.sr.mdd.cov6 <- lavaan::sem(med.ocd.sr.mdd.cov6, data=mds, estimator="MLR", bootstrap=10000) # structural equation model 
summary(f.med.ocd.sr.mdd.cov6, ci=T, fit=T, std = T)

谢谢！！

4

0 回答 0