Hi all,
I have been wracking my brain and searching the literature for answers, but yet none seem clear to me so I feel I must ask. Are "fit statistics (e.g., RMSEA, TLI, CFI, etc..) important or even useful for mediation models?
So my issue is as follows:
When I specify the following model in Lavaan, using DWLS estimation (due to ordinal factors and mediators):
outcome ~ a*fact1+fact2+fact3+fact4+fact5+fact6+b*med1+c*med2 # DIRECT EFFECT
med1 ~ a2*fact1+fact2+fact3+fact4+fact5+fact6
med2 ~ a3*fact1+fact2+fact3+fact4+fact5+fact6
#indirect effect(s)
indirect_mod1:=a2*b
indirect_mod2:=a3*c
#total effect(s)
total_crp:=a +(b*c)
I fail to meet the thresholds for the fit statistics, but the modification indicies suggest allowing covariance between mod1 and mod2, or one predicting the other (which technically could make sense), however, this then renders my model saturated and the fit statistics become useless.
Any ideas what I can do please? The output of the model is interesting, but I am concerned that the fit is too far off, but can't think of any other way to assess mediators without SEM. I have tried looking at the model-implied (fitted) covariance matrix and the residuals of a fitted model, but I cant see anything drastically out of the ordinary there.
Below are Fit outputs for both the model before mediator covariance and after, for reference.
Unsaturated model:
lavaan 0.6-7 ended normally after 53 iterations
Estimator DWLS
Optimization method NLMINB
Number of free parameters 54
Number of observations 9619
Model Test User Model:
Test statistic 35.907
Degrees of freedom 1
P-value (Chi-square) 0.000
Model Test Baseline Model:
Test statistic 474.621
Degrees of freedom 24
P-value 0.000
User Model versus Baseline Model:
Comparative Fit Index (CFI) 0.923
Tucker-Lewis Index (TLI) -0.859
Root Mean Square Error of Approximation:
RMSEA 0.060
90 Percent confidence interval - lower 0.044
90 Percent confidence interval - upper 0.078
P-value RMSEA