我试图找出如何用下面模型的R用Lavaan编写语法。我正在关注博让的书,但是他所举的例子非常简单,我不确定我的说法是否正确。
https://ibb.co/YfYYJYX#此处是模型
我的目标是利用简单调查中的数据进行结构方程建模。我有3个介体,4个独立变量和一个结果。所有变量(结果除外)均用李克特量表(1-7)进行度量。 X1是例如item1,item2,item3,item4的结果。结果(Y)是消耗的百分比,范围在10到100之间。因此,结果变量不会由其他变量显示。
我希望模型的语法正确吗? (例如,我应该像以前那样编写间接效果吗?是否应该添加协方差?) 如何在semPaths中添加错误? 我应该注意以下警告吗:“ lavaan警告:某些估计的ov方差为负”?
谢谢!
#variables specification
mod <-
'
X1 =~ v1 + va2 + v3 + v4
M1 =~ v5+ v6+ v7+ v8 +v9
M2 =~ v10 + v11 + v13
X2 =~ v14 + v15 + v16
X3 =~ v17 + v18 + v19
X4 =~ v21 + v20
M3 =~ v22 + v23 + v24 + v25 + v26
Y =~ y
y ~~ 0*y
#model
M1 ~ a1*X1
M2 ~ a2*X2 + a3*X3
M3 ~ b2*M2
Y ~ c1*M3 + c2*X4
#indirect effect
Ind1 := a1*b1*b2*c1
Ind2 := a2*b2* c1
Ind3 := a3 *b2*c1
#total effect
Tot := (a1*b1*b2*c1) + (a2*b2* c1) + (a3 *b2*c1)
'
fit <-lavaan::sem(mod, data=data, std.lv=TRUE, orthogonal=TRUE)
summary(fit, fit.measures= TRUE, standardized = TRUE)