Dear Forum
I have a question about calculating indirect (and total) effects in a SEM which relates to an environmental psychology study about the effects of greenspace on restoration.
My SEM has two 1st-order latent variables and one second order latent variable, constructed from the two 1st-order latent variables and a third (non-latent) variable. I am mainly interested in the effects of determinants 1 and 2 on "restoration" (the outcome), mediated by the second-order latent variable. My question is about the indirect effects (which I calculate at the end of the model). The indirect effects of both determinant1 and determinant2 are not significant and very small. This is not the case in simpler models I have run, but the effects disappear in this model.
model <-
"
# latent variables
2nd_order =~ a*1st_order1+1st_order2+f*R5
1st_order1 =~ R1 + R2
1st_order2 =~ R3 + R4
# regressions
1st_order1 ~ b*determinant1+control1
R5 ~ g*determinant2+control2+control3+control4+control5
determinant3~ h*determinant2+control2+control3+control4+control5
2nd_order ~ i*determinant2+e*determinant1+control1+control2+control3+control4+control5
restoration ~ c*2nd_order+d*determinant1+control1+j*determinant3+k*determinant2+control2+control3+control4+control5
# control1-5 are (socio-economic) control variables
# indirect effect determinant1
ind1:=a*b*c*e
# total effect determinant1
total1:=ind+d
# indirect effect determinant2
ind2:=c*f*g*h*j*i
# total effect determinant2
total2:=ind2+k
"
model.fit <- sem(model, data=data, missing="fiml", estimator="MLR")
summary(model.fit, fit.measures=T, standardized=T)
What puzzles me is that while the indirect effects are not significant, the individual paths are all significant. So focusing on determinant1, the effect of determinant1 on on 2nd-order is significant and substantial, as is the effect of 2nd-order on restoration, as is the effect of determinant1 on 1st-order1 etc. The diagram attached gives a visual illustration, with all the paths involved in the indirect effect of determinant1 on restoration marked in red (coefficients given are standardized betas).
Of course, given the way indirect effects are calculated, I'm not that surprised that I get non-significant and small effects, as I am multiplying 4 values between 0 and 1 with each other, which results in a very small number. Indeed, in simpler model specifications with fewer values to multiply, I get a significant indirect effect of determinant1.
Am I actually calculating the indirect effects correctly in my R script? Is there a limit to this method of multiplying coefficients to calculate the indirect effect? Is there a better way of calculating indirect effects? Or is my problem related to my sample size, which is around 300 cases.
Thanks for any help
Chris