Hello everyone,
I am trying to fit a cross-lagged moderation model (two waves), with one moderator and two dependent variables. I would really appreciate and be grateful for any help in understanding whether I am running a correct model, or any corrections.
Thank you
Emanuele
--
YA1 = First dependent.
YA2 = Second dependent.
X1 = Independent
W1 = Moderator.
INT = Interaction.
All variables were previously scaled.
model <- "
# Regressions
YA2 ~ YA1 + a1*X1 + W1 + c1*INT
YB2 ~ YB1 + a2*X1 + W1 + c2*INT
# Residual covariances at Time 1
X1 ~~ W1 + YA1 + YB1 + INT
W1 ~~ YA1 + YB1 + INT
YA1 ~~ YB1 + INT
YB1 ~~ INT
# Residual covariances at Time 2
X2 ~~ W1_2 + YA2 + YB2
W1_2 ~~ YA2 + YB2
YA2 ~~ YB2
# Simple slopes for condition effect
MOD_Y1_LOW := a1 + c1*(W1.mean - sqrt(W1.var))
MOD_Y1_MID := a1 + c1*(mean.mean)
MOD_Y1_HIGH := a1 + c1*(mean.mean + sqrt(W1.var))
MOD_Y2_LOW := a2 + c2*(mean.mean - sqrt(W1.var))
MOD_Y2_MID := a2 + c2*(mean.mean)
MOD_Y2_HIGH := a2 + c2*(mean.mean + sqrt(W1.var))
"