Why Specify an Order for Mediator Interventions?
When two mediators have a meaningful temporal or causal order, an intervention can be described in stages: first assign \(M_1\), then assign \(M_2\) conditional on the value of \(M_1\) that was assigned. This is more specific than saying only that the pair \((M_1,M_2)\) should follow some joint distribution. It makes explicit how the second assignment responds to the first mediator and allows a contrast to change one stage while holding the other stage's policy fixed.
For example, suppose a hypothetical intervention first assigns a measure of early treatment engagement, \(M_1\), and then assigns a later adherence measure, \(M_2\), based on the assigned engagement level. The outcome exposure can be held at a chosen value while these mediator policies are varied. The resulting mean is an interventional estimand: it averages outcomes under a specified assignment procedure, not under each person's naturally occurring mediator values.
The distributions \(g_1\) and \(g_2\) can be chosen policies or distributions taken from specified reference regimes. Superscripts such as \(r\) and \(s\) will label the policies used at the first and second stages; they do not, by themselves, imply that the policies are observed exposure levels. Let \(\mu(a_y,m_1,m_2,c)\) denote the conditional mean outcome under exposure \(a_y\) and mediator values \(m_1,m_2\), given \(C=c\).
The Sequential Mediational G-Formula
Under the identification conditions stated below, the mean outcome under first-stage policy \(g_1^r\) and second-stage policy \(g_2^s\) is
The outer expectation averages over the population distribution of baseline covariates. It applies whether \(C\) is discrete or continuous; when \(C\) is continuous, it is an integral with respect to the distribution of \(C\), not a sum over covariate values. The displayed sums over mediators apply to discrete mediators. For continuous mediators, replace the corresponding sums by integrals with respect to the intervention densities.
This formula follows the intervention's order. For each \(C\), average the outcome mean over \(M_2\) under the second-stage policy conditional on \(M_1\); then average that result over \(M_1\) under the first-stage policy; finally average over \(C\). The conditioning of \(g_2^s\) on \(m_1\) is essential: it specifies which second-stage distribution is used after each possible first-stage assignment.
Worked Example: Calculating a Sequential Intervention Mean
Suppose there are no baseline covariates, both mediators are binary, and the outcome is a hypothetical indicator with conditional mean \( \mu(a_y,m_1,m_2)=0.10+0.20m_1+0.30m_2 \). Use a first-stage policy with \(P(M_1=1)=0.40\). Under the second-stage policy, let \(P(M_2=1\mid M_1=0)=0.20\) and \(P(M_2=1\mid M_1=1)=0.70\). Assume these policy assignments are supported by the data and the identification conditions hold.
For \(M_1=0\), the outcome mean averaged over the second-stage assignment is \(0.80(0.10)+0.20(0.40)=0.16\), since the outcome mean is \(0.10\) when \(M_2=0\) and \(0.40\) when \(M_2=1\).
For \(M_1=1\), it is \(0.30(0.30)+0.70(0.60)=0.51\). The first-stage policy assigns \(M_1=0\) with probability \(0.60\) and \(M_1=1\) with probability \(0.40\). Therefore,
The estimated target under these hypothetical policies is a mean outcome of \(0.300\), or 30.0% if the outcome is binary. This is the mean under the specified sequence of stochastic assignments and exposure value \(a_y\); it is not a claim that changing the mediators in an observed population would necessarily cause this mean without the stated assumptions.
Sequential Pathway Contrasts
A sequential pathway contrast changes a policy at one stage while holding the other stage's policy fixed. For example, holding the second-stage policy at \(g_2^0\), a first-stage contrast is
It compares two intervention means that differ in the distribution used to assign \(M_1\), while both use the same rule for assigning \(M_2\) conditional on the assigned \(M_1\). A second-stage contrast can instead hold \(g_1^0\) fixed and compare \(g_2^1\) with \(g_2^0\). These are contrasts between well-defined intervention regimes. Their interpretation depends on which policy is held fixed, and they should not automatically be treated as unique, context-free shares of an overall effect.
Worked Example: Changing One Stage at a Time
Continue with binary mediators and suppose \( \mu(a_y,m_1,m_2)=0.05+0.15m_1+0.25m_2+0.10m_1m_2 \). The common second-stage policy assigns \(M_2=1\) with probability \(0.10\) if \(M_1=0\), and with probability \(0.80\) if \(M_1=1\). Consider first-stage policies with \(P(M_1=1)=0.20\) and \(P(M_1=1)=0.60\). Assume both policies are supported and the identification conditions hold.
If \(M_1=0\), the conditional mean averaged over \(M_2\) is \(0.90(0.05)+0.10(0.30)=0.075\). If \(M_1=1\), the mean is \(0.20(0.20)+0.80(0.55)=0.480\).
Thus the two intervention means are \(0.80(0.075)+0.20(0.480)=0.156\) and \(0.40(0.075)+0.60(0.480)=0.318\). The first-stage contrast is \(0.318-0.156=0.162\), a difference of 16.2 percentage points for a binary outcome.
Now hold the first-stage policy at \(P(M_1=1)=0.40\), and change only the second-stage policy. Let policy \(g_2^0\) use probabilities \(0.10\) and \(0.80\) for \(M_1=0\) and \(M_1=1\), respectively; let \(g_2^1\) use \(0.30\) and \(0.90\). Under \(g_2^0\), the stage-specific means are \(0.075\) and \(0.480\), so the overall mean is \(0.60(0.075)+0.40(0.480)=0.237\). Under \(g_2^1\), they are \(0.125\) and \(0.515\), giving \(0.60(0.125)+0.40(0.515)=0.281\). The second-stage contrast is \(0.281-0.237=0.044\), or 4.4 percentage points. Both contrasts are tied to their stated fixed policies; neither is an unconditional measure of the importance of a mediator.
Identification Conditions
The g-formula is an identification result: it expresses the target intervention mean using observed-data distributions under assumptions. In the basic setting here, \(C\) contains the relevant baseline confounders, and there are no additional exposure-induced mediator-outcome confounders outside the ordered mediator history. The main conditions are:
- Consistency and well-defined interventions: observed outcomes agree with the corresponding potential outcomes for the exposure and mediator values actually received, and the stochastic assignment rules are sufficiently specified to define the intervention.
- No interference: one unit's outcome is not affected by other units' exposure or mediator assignments.
- Exposure exchangeability: conditional on \(C\), exposure assignment is independent of the potential outcomes relevant to the estimand. In a randomized exposure study, randomization supports this condition.
- Sequential mediator-outcome exchangeability: after conditioning on \(A,C,M_1\), there is no unmeasured confounding of the relation between \(M_2\) and the outcome; and, after conditioning on \(A,C\), there is no unmeasured confounding of the relation between \(M_1\) and the outcome that is needed for the intervention mean.
- Positivity: the exposure and mediator values needed by the target policies have support in the relevant observed-data strata. In particular, outcome regression must be identified where \(g_1\) and \(g_2\) place positive probability.
These conditions permit the outcome regression and mediator distributions in the formula to be learned from observed data. If important mediator-outcome confounders are omitted, or a policy assigns mediator values outside observed support, the formula alone does not guarantee identification. As with other interventional mediation estimands, the definition does not require assuming that a person's potential outcome under one exposure is independent of that same person's mediator value under a different exposure.
Continuous Baseline Covariates and the Outer Average
The expectation over \(C\) is part of the estimand, not a notational shortcut that can be replaced by a sum in every setting. If the baseline covariate is continuous, the inner mediator averages are first computed at each covariate value and then integrated over the covariate distribution.
Worked Example: Averaging Over a Continuous Covariate
Let \(C\) be uniformly distributed on \([0,1]\). Suppose \(M_1\) and \(M_2\) are binary, with \(P(M_1=1\mid C=c)=c\), and the sequential second-stage policy has \(P(M_2=1\mid M_1=m_1,C=c)=0.20+0.40m_1\). Let \( \mu(a_y,m_1,m_2,c)=0.10+0.20c+0.10m_1+0.20m_2 \). Assume all required mediator values have support and the identification conditions hold.
At a fixed \(c\), the first-stage policy gives \(E(M_1\mid c)=c\). The second mediator's conditional mean, averaged over the first-stage assignment, is \(E(M_2\mid c)=0.20(1-c)+0.60c=0.20+0.40c\). Consequently, the outcome mean averaged over both stages, conditional on \(C=c\), is \(0.10+0.20c+0.10c+0.20(0.20+0.40c)=0.14+0.38c\).
Now average this expression over the continuous distribution of \(C\):
The target mean is \(0.330\), or 33.0% for a binary outcome. Replacing the outer expectation by a sum over a few arbitrary covariate values would not, in general, yield this population average.
How Sequential Policies Relate to Joint Interventions
A joint stochastic intervention specifies a distribution \(q(m_1,m_2\mid C)\) for the mediator pair. Its intervention mean is
Again, the outer average is an expectation over \(C\), so it applies to continuous as well as discrete baseline covariates. For continuous mediators, replace the sums over mediator values with integrals. A sequential policy induces the joint distribution \(q(m_1,m_2\mid C)=g_1(m_1\mid C)g_2(m_2\mid m_1,C)\). When the joint intervention uses exactly this distribution, both definitions give the same intervention mean. A joint distribution can also be expressed in this sequential factorization wherever the relevant conditional distributions are defined.
The distinction is therefore not that sequential and joint interventions must produce different mediator distributions. It is that the sequential specification makes the assignment order and the conditional second-stage rule explicit. That specification is useful for contrasts that change one stage's policy while retaining the other. A joint policy instead states the target distribution of the pair as a whole; it need not be presented as a procedure that assigns one mediator before the other.
Worked Example: A Joint Policy With a Sequential Representation
Suppose \(C\) is absent and a joint policy for binary mediators assigns probabilities \(q(0,0)=0.50\), \(q(0,1)=0.10\), \(q(1,0)=0.20\), and \(q(1,1)=0.20\). Let the outcome mean be \( \mu(a_y,m_1,m_2)=0.10+0.20m_1+0.30m_2 \). The marginal probability of \(M_1=1\) is \(0.20+0.20=0.40\). Conditional on \(M_1=0\), the probability of \(M_2=1\) is \(0.10/(0.50+0.10)=1/6\); conditional on \(M_1=1\), it is \(0.20/(0.20+0.20)=0.50\). Thus the same joint policy can be implemented sequentially with these first- and second-stage probabilities.
Under the joint policy, \(E(M_1)=0.40\) and \(E(M_2)=0.10+0.20=0.30\). Its outcome mean is \(0.10+0.20(0.40)+0.30(0.30)=0.270\). Under the sequential representation, the implied joint probabilities are \(0.60(5/6)=0.50\), \(0.60(1/6)=0.10\), \(0.40(0.50)=0.20\), and \(0.40(0.50)=0.20\), in the order \((0,0),(0,1),(1,0),(1,1)\). They match \(q\), so the sequential intervention mean is also \(0.270\).
Derivation and a Common Pitfall
For fixed \(C=c\), the sequential policy first gives \(M_1\) probability \(g_1^r(m_1\mid c)\), then gives \(M_2\) probability \(g_2^s(m_2\mid m_1,c)\). By the multiplication rule, the probability of the ordered pair under the intervention is their product. Averaging the outcome mean over this pair distribution yields the inner sums; averaging over baseline covariates yields the outer expectation. Explicitly, for discrete mediators:
The second line is the iterated average: first over \(M_2\) conditional on \(M_1,C\), then over \(M_1\) conditional on \(C\). The derivation also shows why the \(M_2\) policy must be evaluated at each value of \(M_1\), rather than averaged independently of it.
Common mistake: use the marginal distribution of \(M_2\) in place of \(g_2(m_2\mid m_1,C)\), then multiply separate marginal probabilities for the two mediators. That generally defines a different intervention because it removes the specified dependence of the second-stage assignment on the first-stage value. The correction is to apply the second-stage rule conditional on the realized \(M_1\), and only then average over the first-stage assignment. Another error is to write a sum over \(C\) without restricting \(C\) to be discrete; the population average must be written as \(E_C[\cdot]\), or as an integral when \(C\) is continuous.
Check Your Understanding
Use the definitions and formula above to answer the following questions.
- In a sequential intervention, what does \(g_2(m_2\mid m_1,C)\) specify, and why is conditioning on \(m_1\) important?
- How should the outer average in the g-formula be written when \(C\) is continuous?
- Suppose a first-stage contrast changes \(g_1^0\) to \(g_1^1\). What must be held fixed for this to be a contrast that changes only the first-stage policy?
- Under what condition will a joint intervention \(q(m_1,m_2\mid C)\) and a sequential intervention have the same intervention mean?
- Why can multiplying the marginal distributions of \(M_1\) and \(M_2\) fail to represent a specified sequential policy?