Multiple mediator path models were estimated in MPlus 7.11 (Muthén & Muthén, 1998-2012) using full information maximum likelihood (FIML). Several options exist to assess mediated effects in Mplus. The first is the product-of-coefficients approach using the IND command. The second approach involves calculating bias-corrected bootstrap confidence intervals (BBCIs), which is more robust (Preacher & Hayes, 2008). In Mplus it is not possible to account for clustering while using resampling approaches. Therefore, indirect effects when controlling for cluster effects were compared to BBCIs, a procedure used in prior work (Trucco et al., 2014). Although results were largely the same across procedures, we present path estimates accounting for clustering effects, and report BBCIs for indirect effects. Study variables were all normally distributed (range of skew across endogenous variables: −0.48 to 1.65, range of kurtosis = −0.11 to 1.24).
