Tutorials › Propensity Score Matching and Instrumental Variables › Deriving Covariate Standardization from Exchangeability

Causal Effects, Ignorability and Positivity · Tutorial 16 of 699

Deriving Covariate Standardization from Exchangeability

Learn why treatment-specific outcome means must be averaged over the target population’s covariate distribution to identify population mean potential outcomes.

Beginner 10 min read

What You'll Learn

  • Derive the covariate-standardization identity from iterated expectation, conditional exchangeability, and consistency.
  • Distinguish a target-population covariate distribution from the covariate distribution within a treatment group.
  • Calculate standardized mean potential outcomes and their contrast in discrete examples.
  • Explain how crude outcome differences can disagree with a standardized causal contrast.
  • Identify assumptions and data-support problems that limit interpretation.
  • Compare standardization with crude comparisons and model-based adjustment.

Why Average Over the Covariate Distribution?

In an observational study, the people who receive treatment may differ from those who do not. If a covariate predicts both treatment and outcome, a raw comparison of observed outcome means can mix treatment differences with differences in who received treatment. Conditional exchangeability gives a way to address this: compare outcomes within covariate values, then combine those conditional means using the covariate distribution of the population whose effect is of interest.

The key is the choice of weights. To identify a population mean potential outcome, the weights must describe the covariates in that target population, not merely the covariates among treated or untreated people. This procedure is called covariate standardization.

Definition: Let \(A\) be a binary treatment, with treatment states \(a=0\) and \(a=1\); let \(Y(a)\) be the potential outcome under treatment state \(a\); and let \(X\) denote pretreatment covariates. The standardized mean potential outcome for state \(a\) is the average of the conditional mean \(E[Y\mid A=a,X=x]\) over the covariate distribution in the target population. For discrete \(X\), it is \(\sum_x E[Y\mid A=a,X=x]P(X=x)\).

Here \(P(X=x)\) refers to the target population’s distribution. If the study sample represents that population, its covariate proportions may estimate these weights. If the target population is different, its covariate distribution is needed instead. The result can be calculated separately for \(a=1\) and \(a=0\); subtracting the resulting means identifies the population average causal effect, under the assumptions below.

Assumptions and the Identification Argument

The argument uses potential outcomes and conditional exchangeability. Consistency links an individual’s observed outcome to the potential outcome under the treatment actually received: when \(A=a\), \(Y=Y(a)\). Conditional exchangeability states that, within levels of the measured pretreatment covariates, treatment assignment is independent of the potential outcomes. In symbols, \(Y(a)\perp A\mid X\), for each treatment state \(a\).

The conditional outcome mean \(E[Y\mid A=a,X=x]\) must also be defined for covariate values included in the target average. In practical terms, the data need to contain relevant observations at the treatment state being used to estimate that conditional mean. Consistency and exchangeability alone cannot reveal an outcome mean for a covariate group with no observations receiving that treatment.

1
Average within covariate values.
By the law of total expectation, \(E[Y(a)]=\sum_x E[Y(a)\mid X=x]P(X=x)\) for discrete \(X\). This expresses the population mean as a weighted average of conditional potential-outcome means.
2
Use conditional exchangeability.
Because \(Y(a)\) is independent of \(A\) given \(X\), its conditional mean at \(X=x\) equals its conditional mean among those with \(A=a\): \(E[Y(a)\mid X=x]=E[Y(a)\mid A=a,X=x]\), where that treatment group is observed.
3
Use consistency.
Among units with \(A=a\), consistency gives \(Y=Y(a)\). Therefore, \(E[Y(a)\mid A=a,X=x]=E[Y\mid A=a,X=x]\).
4
Substitute into the population average.
Combining the equalities yields \(E[Y(a)]=\sum_x E[Y\mid A=a,X=x]P(X=x)\). The weights remain the target population’s \(P(X=x)\) throughout.

Thus, for a binary treatment and discrete covariates, the two identified population means are

$$ E[Y(1)]=\sum_x E[Y\mid A=1,X=x]P(X=x), \qquad E[Y(0)]=\sum_x E[Y\mid A=0,X=x]P(X=x). $$

The population average causal effect is their difference:

$$ E[Y(1)-Y(0)] = \sum_x \left\{E[Y\mid A=1,X=x]-E[Y\mid A=0,X=x]\right\}P(X=x). $$

The final expression shows two equivalent ways to calculate the contrast: standardize each treatment-specific mean separately and subtract, or first take the conditional treatment contrast in each covariate group and average those contrasts using the target covariate weights. For continuous covariates, the same reasoning uses an integral over the target distribution of \(X\) rather than a sum.

Worked Applications

Worked Example: Standardizing a Treated Mean

Suppose a hypothetical target population has two covariate groups: \(X=0\) comprises 60% of the population and \(X=1\) comprises 40%. In a hypothetical study, the observed mean outcome under treatment is 12 units for \(X=0\) and 18 units for \(X=1\). Assume consistency, conditional exchangeability, and that the relevant treated observations exist in both groups. Find the mean potential outcome under treatment for the target population.

The required weights are the target group proportions, so

$$ E[Y(1)] = 12(0.60)+18(0.40)=7.20+7.20=14.40\text{ units}. $$

The identified mean potential outcome under treatment is 14.40 units. This is not necessarily the observed mean among treated people: if treated people have a different mix of \(X=0\) and \(X=1\) than the target population, their raw mean uses different weights. Standardization answers the target-population question by applying the target weights to the within-group means.

Worked Example: A Standardized Causal Contrast

Consider a hypothetical population in which 70% of people are in a lower-risk covariate group and 30% are in a higher-risk group. The outcome is recovery time in days, with lower values being better. Suppose the observed conditional means are 8 and 12 days under treatment, and 10 and 15 days without treatment, in the lower-risk and higher-risk groups respectively. Assume consistency, conditional exchangeability, and observations in each treatment state within both groups.

First standardize the treatment mean:

$$ E[Y(1)] = 8(0.70)+12(0.30)=5.60+3.60=9.20\text{ days}. $$

Then standardize the no-treatment mean using the same population weights:

$$ E[Y(0)] = 10(0.70)+15(0.30)=7.00+4.50=11.50\text{ days}. $$

The population average causal effect, defined as treatment minus no treatment, is

$$ E[Y(1)-Y(0)] = 9.20-11.50=-2.30\text{ days}. $$

Under the stated assumptions, treatment reduces mean recovery time by 2.30 days in this target population. The units are days, and the negative sign indicates a lower mean under treatment. This conclusion concerns the population average; it does not imply that every individual’s recovery time would fall by 2.30 days.

Worked Example: Why the Target Weights Matter

In another hypothetical population, \(X=0\) makes up 80% of the target population and \(X=1\) makes up 20%. The conditional outcome means under treatment are 5 and 15 units for the two groups. Suppose, however, that only 25% of the treated group is in \(X=0\), while 75% is in \(X=1\). Assume the conditional means are identified under consistency and conditional exchangeability, and observations exist in both treated covariate groups.

The treated group’s raw mean is weighted by its own covariate distribution:

$$ 5(0.25)+15(0.75)=1.25+11.25=12.50\text{ units}. $$

But the target-population mean potential outcome under treatment uses the target proportions:

$$ E[Y(1)] = 5(0.80)+15(0.20)=4.00+3.00=7.00\text{ units}. $$

The two means differ because they answer different questions: 12.50 units describes the observed treated group, whereas 7.00 units is the standardized mean for the specified target population. The difference is not evidence that the arithmetic is inconsistent; it reflects different weights on groups with different conditional means.

Interpretation, Diagnostics, and Limitations

Standardization is an identification argument before it is a computational recipe. Consistency requires a well-defined treatment and a meaningful link between assigned treatment and observed outcome. Conditional exchangeability requires that, after conditioning on the chosen pretreatment covariates, treatment assignment does not retain an association with the potential outcomes. This is a substantive assumption: observed data alone cannot establish that all relevant confounding has been removed.

Useful diagnostics include checking that covariates are measured before treatment, examining whether important prognostic characteristics differ across treatment groups within covariate groups, and checking whether both treatment states are represented where the target population has support. Sparse or empty cells make conditional means unstable or unavailable. Many covariates can create this problem even when each covariate appears well represented separately. Diagnostics can reveal warning signs, but they do not prove exchangeability.

Standardization also depends on how the conditional means are obtained. With a small number of discrete covariate groups, they may be estimated directly within groups. With many or continuous covariates, they may be estimated using a statistical model and then averaged over the target covariate distribution. A model can make estimation practical, but its predictions may be inaccurate if its functional form is inadequate. Model-based adjustment should not be treated as a substitute for the identification assumptions.

The main alternatives answer related but not always identical questions. A crude difference in observed outcome means uses the covariate distributions of the treatment groups and is not generally a population causal contrast when assignment is confounded. A regression-adjusted contrast may account for covariates through a fitted model, while standardization makes the target averaging distribution explicit. The introductory material on conditional exchangeability explains why adjustment must address treatment-related differences; the distinctive step here is deriving which weights recover the target population mean.

Common pitfall: “The treated group’s mean outcome is the mean outcome everyone would have had under treatment.” This reasoning uses the treated group’s covariate mix as if it were the target population’s mix. The correction is to estimate the treatment-specific conditional mean within covariate values and average those means with the target population’s covariate weights.

Standardization does not recover individual potential outcomes, and it does not guarantee that an estimate is precise. Its population interpretation can fail if consistency or conditional exchangeability is implausible, if required treatment groups are absent at relevant covariate values, or if estimated conditional means are seriously inaccurate. These limitations should be considered alongside the reported standardized contrast.

Check Your Understanding

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

  1. A target population is 50% in each of two covariate groups. The conditional mean under treatment is 6 units in the first group and 10 in the second. What is the standardized treatment mean?
  2. State the two assumptions that let the conditional observed mean \(E[Y\mid A=a,X=x]\) stand in for \(E[Y(a)\mid X=x]\).
  3. Why can the raw mean among treated people differ from the standardized treatment mean even if the conditional means are accurately estimated?
  4. In one sentence, interpret a standardized contrast of \(-1.5\) hours when the outcome is time to an event.

Solutions. (1) Apply the target weights: \(6(0.50)+10(0.50)=3+5=8\) units. (2) Conditional exchangeability gives \(E[Y(a)\mid X=x]=E[Y(a)\mid A=a,X=x]\), and consistency gives \(E[Y(a)\mid A=a,X=x]=E[Y\mid A=a,X=x]\). (3) The raw treated mean uses the treated group’s covariate proportions, while the standardized mean uses the target population’s proportions; the weights can differ. (4) Under the identification assumptions, the treatment state has a mean event time 1.5 hours lower than the comparison state in the target population.

Clinical Biostats