Bayesian causal mediation analysis through linear mixed-effect models

In mediational settings, the main focus is on the estimation of the indirect effect of an exposure on an outcome through a third variable called mediator. The traditional maximum likelihood estimation method presents several problems in the estimation of the standard error and the confidence interval of the indirect effect. In this paper, we propose a Bayesian approach to obtain the posterior distribution of the indirect effect through MCMC, in the context of mediational mixed models for longitudinal data. A simulation study shows that our method outperforms the traditional maximum likelihood approach in terms of bias and coverage rates.