Я пытаюсь выяснить, как написать синтаксис в Lavaan с R из приведенной ниже модели. Я слежу за книгой Божана, но примеры, которые он приводит, очень просты, и я очень не уверен, правильно ли это.
https://ibb.co/YfYYJYX # здесь модель
Моя цель - выполнить моделирование структурного уравнения с использованием данных простого опроса. У меня есть 3 посредника, 4 независимых переменных и один результат. Все переменные (кроме результата) измеряются по шкале Лайкерта (1-7). X1 - это результаты, например, item1, item2, item3, item4. Результат (Y) - это процент потребления, в диапазоне от 10 до 100. Таким образом, переменная результата не проявляется другими переменными.
Синтаксис модели, которую я хотел бы иметь, правильный? (например, я должен написать косвенные эффекты как я? Я должен добавить ковариацию?)
Как я могу добавить ошибки в semPaths?
Мне нужно позаботиться об этом предупреждении: "Лаваан ПРЕДУПРЕЖДЕНИЕ: некоторые предполагаемые отклонения 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)