Improved estimation of controlled direct effects in the presence of unmeasured confounding of intermediate variables

Adjusting for a causal intermediate is a common analytic strategy for estimating an average causal direct effect (ACDE). The ACDE is the component of the total exposure effect that is not relayed through the specified intermediate. Even if the total effect is unconfounded, the usual ACDE estimate may be biased when an unmeasured variable affects the intermediate and outcome variables. Using linear programming optimization to compute non-parametric bounds, we develop new ACDE estimators for binary measured variables in this causal structure, and use root mean square confounding bias (RMSB) to compare their performance with the usual stratified estimator in simulated distributions of target populations comprised of the 64 possible potential response types as well as distributions of target populations restricted to subsets of 18 or 12 potential response types defined by monotonicity or no-interactions assumptions of unit-level causal effects. We also consider target population distributions conditioned on fixed outcome risk among the unexposed, or fixed true ACDE in one stratum of the intermediate. Results show that a midpoint estimator constructed from the optimization bounds has consistently lower RMSB than the usual stratified estimator both unconditionally and conditioned on any risk in the unexposed. When conditioning on true ACDE, this midpoint estimator performs more poorly only when conditioned on an extreme true ACDE in one stratum of the intermediate, yet outperforms the stratified estimator in the other stratum when interaction is permitted. An alternate 'limit-modified crude' estimator can never perform less favourably than the stratified estimator, and often has lower RMSB.