Tutorials › Marginal Structural Models › Inductive Identification Across Multiple Decision Times

Sequential Exchangeability, Consistency and Positivity · Tutorial 45 of 501

Inductive Identification Across Multiple Decision Times

Learn to derive a regime-specific mean by iterating conditional outcome and covariate distributions across multiple treatment decisions.

Intermediate 12 min read

What You'll Learn

  • Define observed and counterfactual histories for a fixed treatment regime.
  • State the assumptions needed for identification across multiple decision times.
  • Derive the recursive conditional-mean formula by backward induction.
  • Calculate regime-specific means for discrete and continuous covariates.
  • Distinguish sequential standardization from conditioning on people who followed the regime.
  • Recognize support and model failures that can undermine identification.

From Two Decisions to a Sequence

With several treatment decisions, a regime-specific mean is not obtained by adjusting for one collection of covariates in a single step. At each decision time, the covariate distribution depends on the treatment history that preceded it. Identification therefore proceeds through a sequence: identify the final conditional outcome mean, average it over the covariate distribution after the preceding treatment, and continue backward until only baseline covariates remain.

The two-time identification formula is the first instance of this pattern. Here the aim is to establish the pattern for any finite number of decisions, and to show why its recursive form follows from consistency and sequential exchangeability.

Histories, Regimes, and Assumptions

Suppose treatment is assigned at times \(t=0,\ldots,K\), where \(K\) is finite. Let \(L_0\) denote baseline covariates. For \(t=0,\ldots,K\), let \(A_t\) be treatment at time \(t\), and let \(L_{t+1}\) be the covariates measured after \(A_t\) when \(t<K\). The outcome \(Y\) is measured after the final treatment \(A_K\). Write \(\bar A_t=(A_0,\ldots,A_t)\) and \(\bar L_t=(L_0,\ldots,L_t)\). The history just before treatment at time \(t\) is \(H_t=(\bar L_t,\bar A_{t-1})\), with \(H_0=L_0\).

Fix a treatment regime \(\bar a=(a_0,\ldots,a_K)\). Write \(L_{t+1}^{\bar a_t}\) for the covariate that would follow treatment history \(\bar a_t=(a_0,\ldots,a_t)\), and \(Y^{\bar a}\) for the outcome under the full regime. The counterfactual history before decision \(t\) is \(H_t^{\bar a}=(L_0,L_1^{a_0},\ldots,L_t^{\bar a_{t-1}},\bar a_{t-1})\).

Definition: Consistency requires that, when a person's observed treatment history through time \(t\) is \(\bar a_t\), the observed covariate after that history equals \(L_{t+1}^{\bar a_t}\); when the observed full treatment history is \(\bar a\), the observed outcome equals \(Y^{\bar a}\). Sequential exchangeability requires, at every time \(t\), that the treatment \(A_t\) be conditionally independent of the relevant future counterfactual covariates and outcome under the regime, given the observed history \(H_t\), for histories compatible with the regime. In particular, the final-time condition includes \(Y^{\bar a}\mathrel{\perp\!\!\!\perp} A_K\mid H_K\). Identification also requires treatment support: the conditional probability of receiving \(a_t\) given \(H_t\) must be positive on histories that can occur under \(\bar a\).

These assumptions have distinct jobs. Consistency connects observed quantities to counterfactuals when the observed treatment matches the regime. Exchangeability lets the observed treatment-specific conditional distributions stand in for counterfactual distributions. Support ensures those observed conditional distributions exist where the regime needs them. The treatment versions represented by each \(a_t\) must also be well defined, and the outcome must have a finite mean for the mean formula below to be finite.

The Backward Recursion

The formula is easiest to state recursively. At the final decision, define the conditional outcome mean

$$ \mu_K(h_K)=E[Y\mid H_K=h_K,\ A_K=a_K]. $$

For \(t=K-1,\ldots,0\), define

$$ \mu_t(h_t)=\int \mu_{t+1}(h_t,a_t,l_{t+1})\, dF(l_{t+1}\mid H_t=h_t,A_t=a_t). $$

The argument \((h_t,a_t,l_{t+1})\) is the next history with treatment set to its regime value. When covariates are discrete, the integral is a sum. The final identified mean averages the resulting baseline-specific quantity over the baseline distribution:

$$ E[Y^{\bar a}]=\int \mu_0(l_0)\,dF(l_0). $$
Theorem: Under consistency, sequential exchangeability, treatment support along the regime, and integrability of \(Y^{\bar a}\), the backward recursion above identifies \(E[Y^{\bar a}]\). Equivalently, the mean is obtained by integrating the observed final-time conditional outcome mean over the sequence of observed covariate distributions conditional on the regime treatment at each time, beginning with the baseline distribution.

Proof. Begin at the final decision. Given a counterfactual history \(H_K^{\bar a}\), sequential exchangeability allows the conditional mean of \(Y^{\bar a}\) to be represented by the mean among people with that observed history who receive \(A_K=a_K\). Consistency identifies their observed outcome with \(Y^{\bar a}\). Thus the conditional counterfactual outcome mean is \(\mu_K(H_K^{\bar a})\).

Now suppose the conditional counterfactual mean from the next decision onward is represented by \(\mu_{t+1}\). By the law of total expectation,

$$ E[Y^{\bar a}\mid H_t^{\bar a}] = E\!\left[ E[Y^{\bar a}\mid H_{t+1}^{\bar a}] \mid H_t^{\bar a} \right]. $$

Under the induction hypothesis, the inner conditional mean is \(\mu_{t+1}(H_{t+1}^{\bar a})\). Sequential exchangeability at time \(t\), together with consistency, identifies the conditional distribution of \(L_{t+1}^{\bar a_t}\) given the regime history with the observed distribution \(F(L_{t+1}\mid H_t,A_t=a_t)\). Therefore the outer expectation is exactly the recursion defining \(\mu_t(H_t^{\bar a})\). This proves the recursion backward through all decision times. Finally, \(H_0^{\bar a}=L_0\), whose distribution is the observed baseline distribution; averaging \(\mu_0(L_0)\) over that distribution gives \(E[Y^{\bar a}]\). The finite number of steps and the stated integrability condition make each conditional expectation well defined. \(\square\)

The same argument also explains the equivalent integral form: each recursive averaging step inserts one conditional covariate distribution. No single joint distribution of future covariates has to be guessed; the regime mean is built from the baseline law and the successive treatment-specific transition laws.

Worked Examples

Worked Example: Two Decisions With Binary Covariates

Let \(K=1\), so treatment is assigned at times 0 and 1. Let \(L_0\) be binary with \(P(L_0=1)=0.4\). Under the fixed regime \((a_0,a_1)=(1,0)\), suppose \(P(L_1=1\mid L_0=0,A_0=1)=0.2\) and \(P(L_1=1\mid L_0=1,A_0=1)=0.7\). The final conditional outcome means under \(A_1=0\) are \(0.1,0.5,0.4,0.8\) for \((L_0,L_1)=(0,0),(0,1),(1,0),(1,1)\), respectively. Calculate the regime-specific mean.

For \(L_0=0\), the recursive mean is \(0.8(0.1)+0.2(0.5)=0.08+0.10=0.18\). For \(L_0=1\), it is \(0.3(0.4)+0.7(0.8)=0.12+0.56=0.68\). Averaging over baseline covariates gives \(0.6(0.18)+0.4(0.68)=0.108+0.272=0.38\). Thus \(E[Y^{(1,0)}]=0.38\). The calculation averages the final outcome mean over the covariate distribution following \(A_0=1\), then averages over baseline \(L_0\).

Worked Example: Three Decisions and Two Covariate Transitions

Let \(K=2\), with no baseline covariates, and consider regime \((1,1,1)\). Suppose \(P(L_1=1)=0.5\) under \(A_0=1\); after \(A_1=1\), let \(P(L_2=1\mid L_1=0)=0.2\) and \(P(L_2=1\mid L_1=1)=0.8\). Under \(A_2=1\), suppose the outcome means for \((L_1,L_2)=(0,0),(0,1),(1,0),(1,1)\) are \(0.1,0.4,0.5,0.9\).

Start at the final decision. For \(L_1=0\), averaging over \(L_2\) gives \(\mu_1(0)=0.8(0.1)+0.2(0.4)=0.08+0.08=0.16\). For \(L_1=1\), \(\mu_1(1)=0.2(0.5)+0.8(0.9)=0.10+0.72=0.82\). Now average over \(L_1\): \(E[Y^{(1,1,1)}]=0.5(0.16)+0.5(0.82)=0.08+0.41=0.49\). Both covariate transitions enter, in order, before the baseline average (which is trivial here because there are no baseline covariates).

Worked Example: A Continuous Intermediate Covariate

Consider two decisions, \(K=1\), and a binary baseline covariate \(L_0\) with \(P(L_0=1)=0.3\). Under \(A_0=1\), let \(L_1\mid L_0=0\sim N(2,1)\) and \(L_1\mid L_0=1\sim N(4,1)\). Under \(A_1=0\), suppose \(E[Y\mid L_0,L_1,A_1=0]=1+2L_0+0.5L_1\).

For \(L_0=0\), the conditional mean after integrating out \(L_1\) is \(1+2(0)+0.5E[L_1\mid L_0=0,A_0=1]=1+0+0.5(2)=2\). For \(L_0=1\), it is \(1+2(1)+0.5(4)=5\). The variance of \(L_1\) does not enter this calculation because the outcome regression is linear in \(L_1\). Finally, \(E[Y^{(1,0)}]=0.7(2)+0.3(5)=1.4+1.5=2.9\).

Why the Induction Matters

A tempting shortcut is to calculate the outcome mean among everyone whose observed treatment history exactly equals the regime. That generally does not reproduce the regime-specific mean. Treatment assignment can depend on past covariates, so selecting people who happened to follow the entire regime can change the distribution of their histories. The recursive formula instead averages each conditional distribution using the treatment-specific transition law at that decision, then integrates over the preceding history.

This formula is an identification result, not a guarantee that an estimate will be reliable. In practice, estimating it often requires models for the conditional outcome mean and for each covariate transition. Poorly specified models can yield poor estimates even when the identifying assumptions are plausible. Support failures also matter: if the regime treatment is never observed for some histories that arise under the regime, the relevant conditional distribution cannot be learned directly from the observed data.

ApproachWhat it doesKey distinction
Sequential standardization (the recursion here)Models conditional outcome and covariate distributions, then averages in sequence.Makes the regime-specific history distribution explicit.
Inverse-probability weightingReweights observed data to represent a treatment regime.Uses treatment probabilities in the weights rather than the same backward outcome-and-covariate averaging.
Ordinary outcome regressionModels outcomes conditional on selected observed covariates and treatments.A conditional regression coefficient alone need not equal the marginal regime-specific mean.

These approaches are alternatives for estimation, not different definitions of the causal target. The recursion is especially useful for seeing exactly which conditional distributions identify that target and where a required distribution is unsupported.

Common Mistake and Diagnostics

  • Wrong reasoning: “The people observed to follow the regime are the people who would follow it under intervention, so their outcome average is the regime mean.” This ignores selection into observed treatment histories.
  • Correction: Use the regime-specific conditional outcome and covariate distributions in the backward recursion, and integrate over the resulting histories rather than simply restricting to observed adherers.
  • Warning sign: The estimated regime mean changes sharply when a conditional model is modestly altered. This can indicate model dependence, sparse histories, or extrapolation.
  • Warning sign: Some treatment-history strata needed by the regime have no observations. The corresponding conditional law is not empirically identified without additional assumptions.
  • Limitation: The derivation cannot repair unmeasured confounding, ill-defined treatment versions, or unsupported treatment choices. Those are failures of the assumptions, not algebraic problems with the recursion.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. For \(K=0\), write the identifying formula for \(E[Y^{a_0}]\).
  2. Suppose a two-decision regime has \(P(L_1=1\mid L_0=0,A_0=a_0)=0.25\). The final outcome means are 2 when \(L_1=0\) and 6 when \(L_1=1\), at \(L_0=0,A_1=a_1\). What is the conditional regime mean at \(L_0=0\)?
  3. Why does the induction require a conditional distribution for each covariate transition, rather than only the final outcome regression?
  4. What is wrong with treating the mean among observed regime-adherers as automatically equal to the counterfactual regime mean?

Solutions. (1) Define \(\mu_0(l_0)=E[Y\mid L_0=l_0,A_0=a_0]\); then \(E[Y^{a_0}]=\int\mu_0(l_0)dF(l_0)\). (2) The mean is \(0.75(2)+0.25(6)=1.5+1.5=3\). (3) The covariate distributions determine which histories occur under the regime, so the final conditional outcome mean must be averaged over each successive regime-specific transition. (4) Observed adherence can select people according to their covariate history and treatment assignment; the selected histories need not have the counterfactual regime distribution. The recursive integrations account for that distribution explicitly.

Clinical Biostats