> Model15 <- '
MED ~ a1*IV
Dep ~ c1*IV + c2*MOD+ c3* MOD:IV + b1*MED + b2* MOD:MED
#Indirect
a1b1 := a1*b1 #Will I use b2 in here as well?
#direct
c1c3 := c1*c3 #Will the operator be “+”?
#Total
total := a1b1 + c1c3
#Index of moderated mediation
index.mod.med := b1*c3 #Should I add “b2” as well?
#MOD~ MOD.mean*1
MOD~ MOD.mean* 1
#MOD~~ MOD.var* MOD
MOD~~MOD.var* MOD
#Are the following codes about indirect, direct and total effects conditional on moderator correct?
# Indirect effects conditional on moderator
indirect.SDbelow := (a1 + b2*( MOD.mean-sqrt(MOD.var))) * b1
indirect.SDabove := (a1 + b2*( MOD.mean+sqrt(MOD.var))) * b1
# Direct effects conditional on moderator
direct.SDbelow := c1 + c3*( MOD.mean-sqrt(MOD.var))
direct.SDabove := c1 + c3*( MOD.mean+sqrt(MOD.var))
# Total effects conditional on moderator
total.SDbelow := direct.SDbelow + indirect.SDbelow
total.SDabove := direct.SDabove + indirect.SDabove
# Proportion mediated conditional on moderator
prop.mediated.SDbelow =~ indirect.SDbelow / total.SDbelow
prop.mediated.SDabove =~ indirect.SDabove / total.SDabove
## Ensure country variable is character
IntOrnek$country <- as.character(IntOrnek$country) # Add this line'
# Fit the multigroup model
Model15Sem <- sem(Model15, data = IntOrnek, estimator = "MLR", missing = "ML", group = "country")
Model15Sem.constrained <- sem(Model15, data = IntOrnek, estimator = "MLR", missing = "ML", group = "country", group.equal = c("regressions"))
summary(Model15Sem, modindices = TRUE, fit.measures = TRUE, standardized = TRUE)