In the three - path mediation
analysis, we used the
multivariate delta estimator using the following equation from [109]: This variance estimate was used in Empirical Bayes estimation
procedure for second - level bootstrapping of the path coefficients [8].