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})\).
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
For \(t=K-1,\ldots,0\), define
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:
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,
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.
| Approach | What it does | Key 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 weighting | Reweights observed data to represent a treatment regime. | Uses treatment probabilities in the weights rather than the same backward outcome-and-covariate averaging. |
| Ordinary outcome regression | Models 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.
- For \(K=0\), write the identifying formula for \(E[Y^{a_0}]\).
- 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\)?
- Why does the induction require a conditional distribution for each covariate transition, rather than only the final outcome regression?
- 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.