Maximum Likelihood Analysis of Linear Mediation Models with Treatment–Mediator Interaction

This research concerns a mediation model, where the mediator model is linear and the outcome model is also linear but with a treatment–mediator interaction term and a residual correlated with the residual of the mediator model. Assuming the treatment is randomly assigned, parameters in this mediation model are shown to be partially identifiable. Under the normality assumption on the residual of the mediator and the residual of the outcome, explicit full-information maximum likelihood estimates of model parameters are introduced given the correlation between the residual for the mediator and the residual for the outcome. A consistent variance matrix of these estimates is derived. Currently, the coefficients of this mediation model are estimated using the iterative feasible generalized least squares (IFGLS) method that is originally developed for seemingly unrelated regressions (SURs). We argue that this mediation model is not a system of SURs. While the IFGLS estimates are consistent, their variance matrix is not. Theoretical comparisons of the FIMLE variance matrix and the IFGLS variance matrix are conducted. Our results are demonstrated by simulation studies and an empirical study. The FIMLE method has been implemented in a freely available R package iMediate.

[1]  Peter Schmidt,et al.  ON THE ESTIMATION OF TRIANGULAR STRUCTURAL SYSTEMS , 1978 .

[2]  Kosuke Imai,et al.  Causal Mediation Analysis Using R , 2010 .

[3]  T. VanderWeele,et al.  Sharp sensitivity bounds for mediation under unmeasured mediator-outcome confounding , 2016, Biometrika.

[4]  J. Pearl,et al.  A New Identification Condition for Recursive Models With Correlated Errors , 2002 .

[5]  T. VanderWeele,et al.  Sensitivity analysis for direct and indirect effects in the presence of exposure-induced mediator-outcome confounders. , 2014, Epidemiology, biostatistics and public health.

[6]  L. Keele,et al.  A General Approach to Causal Mediation Analysis , 2010, Psychological methods.

[7]  S. Cessie.  Bias Formulas for Estimating Direct and Indirect Effects When Unmeasured Confounding Is Present. , 2016 .

[8]  M. Arellano.  An efficient GLS estimator of triangular models with covariance restrictions , 1989 .

[9]  L. Keele,et al.  Identification, Inference and Sensitivity Analysis for Causal Mediation Effects , 2010, 1011.1079.

[10]  Kristopher J Preacher,et al.  SPSS and SAS procedures for estimating indirect effects in simple mediation models , 2004, Behavior research methods, instruments, & computers : a journal of the Psychonomic Society, Inc.

[11]  David R. Cox.  Planning of Experiments , 1958 .

[12]  Tyler J. VanderWeele,et al.  Conceptual issues concerning mediation, interventions and composition , 2009 .

[13]  A. Zellner.  An Efficient Method of Estimating Seemingly Unrelated Regressions and Tests for Aggregation Bias , 1962 .

[14]  D. A. Kenny,et al.  The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. , 1986, Journal of personality and social psychology.

[15]  JEFFREY M. ALBERT,et al.  Sensitivity analyses for parametric causal mediation effect estimation. , 2015, Biostatistics.

[16]  Michael Eichler,et al.  Computing Maximum Likelihood Estimates in Recursive Linear Models with Correlated Errors , 2009, J. Mach. Learn. Res..

[17]  Tyler J. VanderWeele,et al.  Explanation in Causal Inference: Methods for Mediation and Interaction , 2015 .

[18]  D. Mackinnon.  Introduction to Statistical Mediation Analysis , 2008 .