Tutorials › Propensity Score Matching and Instrumental Variables › Combining Stratum-Specific Effects for an ATE

Weighting and Stratification · Tutorial 212 of 699

Combining Stratum-Specific Effects for an ATE

Learn why ATE estimation combines stratum-specific treatment effects using each stratum’s share of the target population, and how to calculate and assess that estimate.

Intermediate 12 min read

What You'll Learn

  • Derive the ATE decomposition across pre-treatment strata
  • Distinguish target-population shares from equal weights or sample shares
  • Calculate the stratum-weighted estimator from treatment and control means
  • State the identification conditions needed for a causal interpretation
  • Diagnose empty treatment groups, poor overlap, and mismatched target weights
  • Compare stratum weighting with pooled and inverse-probability-weighted estimates

Why Stratum Effects Need to Be Combined

Propensity-score subclassification divides a population into strata, then compares treated and untreated outcomes within each stratum. Those comparisons answer a stratum-specific question. To obtain an average treatment effect (ATE) for a target population, the stratum effects must then be combined in a way that reflects how much of that target population belongs to each stratum.

The key is the target population’s composition—not an automatic equal average across strata. If one stratum contains 70% of the target population and another contains 30%, their effects contribute in those proportions. This remains true when the stratum-specific effects differ: the combined ATE averages the differences, not the treatment and control means across strata without regard to their membership.

Definition: Let \(A\) be a binary treatment, \(Y(1)\) and \(Y(0)\) the potential outcomes under treatment and control, and \(S\) a pre-treatment stratum taking values \(1,\ldots,K\). The stratum-specific average treatment effect is \(\tau_k=E[Y(1)-Y(0)\mid S=k]\). For a target population in which stratum \(k\) has share \(q_k=P(S=k)\), the ATE is \(\tau_{\mathrm{ATE}}=\sum_{k=1}^{K}q_k\tau_k\), where \(q_k\geq 0\) and \(\sum_{k=1}^{K}q_k=1\).

When the target population is the study sample and that sample represents the target, the usual estimated share is \(\hat q_k=n_k/n\), where \(n_k\) is the number of sample members in stratum \(k\), and \(n\) is the total sample size. If target-population shares are known from an external source, use those shares instead. The choice of \(q_k\) is part of defining the estimand; it is not merely a computational detail.

Identification and the Estimator

A stratum-weighted difference in observed means has a causal interpretation only under conditions that connect observed outcomes to potential outcomes. Let \(Y\) be the observed outcome. The consistency condition says \(Y=Y(A)\): a person’s observed outcome equals the potential outcome under the treatment actually received. Assume also that treatment assignment is conditionally exchangeable within each stratum, \(\{Y(1),Y(0)\}\perp A\mid S\), and that treatment positivity holds: \(0<P(A=1\mid S=k)<1\) for every stratum with positive target share. These conditions, together with an appropriate target-population distribution of \(S\), identify the ATE by stratum-specific comparisons.

For propensity-score strata, conditional exchangeability within a stratum is a substantive requirement, not something guaranteed simply by dividing the score into bins. Stratification may make covariate distributions more comparable, but a coarse grouping can leave residual differences within strata. The earlier material on subclassification provides the setting for forming such groups; here, the central task is to combine their effects correctly.

Estimator: In stratum \(k\), let \(\bar Y_{1k}\) be the observed mean outcome among treated units and \(\bar Y_{0k}\) the observed mean among control units. Define \(\hat\tau_k=\bar Y_{1k}-\bar Y_{0k}\). Given target shares \(q_k\), the plug-in estimator is \(\hat\tau_{\mathrm{ATE}}=\sum_{k=1}^{K}q_k\hat\tau_k\). If the sample estimates its own target shares, use \(\hat\tau_{\mathrm{ATE}}=\sum_{k=1}^{K}(n_k/n)\hat\tau_k\).

The arithmetic has two stages: estimate a treated-versus-control contrast separately within each stratum, then average those contrasts using the target shares. A stratum’s share is not its treatment-group proportion, nor the proportion treated within it. The same share \(q_k\) weights that stratum’s estimated treatment effect.

Why the Decomposition Holds

The decomposition follows by conditioning the population average treatment effect on the stratum. First partition the population according to \(S\); the expectation of a variable equals the sum of its conditional expectations multiplied by the probabilities of the corresponding groups. Apply that identity to the individual treatment effect \(Y(1)-Y(0)\):

$$ \begin{aligned} \tau_{\mathrm{ATE}} &=E[Y(1)-Y(0)]\\ &=\sum_{k=1}^{K}P(S=k)E[Y(1)-Y(0)\mid S=k]\\ &=\sum_{k=1}^{K}q_k\tau_k. \end{aligned} $$

To see how the observed-data contrast identifies each \(\tau_k\), consider either treatment level \(a\in\{0,1\}\). Conditional exchangeability implies \(E[Y(a)\mid A=a,S=k]=E[Y(a)\mid S=k]\). By consistency, \(Y=Y(a)\) for units with \(A=a\), so \(E[Y\mid A=a,S=k]=E[Y(a)\mid S=k]\). Subtracting the result for \(a=0\) from that for \(a=1\) gives

$$ E[Y\mid A=1,S=k]-E[Y\mid A=0,S=k] =E[Y(1)-Y(0)\mid S=k] =\tau_k. $$

Positivity ensures both conditional observed means are defined in every stratum that contributes to the target ATE. In a finite sample, \(\hat\tau_k\) estimates that contrast, and the weighted estimator substitutes estimated or known target shares and sample stratum contrasts into the population decomposition.

Worked Example: ATE From Two Strata

Suppose a hypothetical study’s target population is represented by 100 people: 60 in a lower-score stratum and 40 in a higher-score stratum. In the lower-score stratum, the 30 treated people have a mean outcome of 8 units and the 30 controls have a mean of 6 units. In the higher-score stratum, 20 treated people have a mean of 12 units and 20 controls have a mean of 9 units. Assume consistency, conditional exchangeability within each stratum, and treatment positivity.

The estimated stratum effects are

$$ \hat\tau_1=8-6=2\text{ units},\qquad \hat\tau_2=12-9=3\text{ units}. $$

The target shares are \(q_1=60/100=0.60\) and \(q_2=40/100=0.40\). Therefore,

$$ \hat\tau_{\mathrm{ATE}} =(0.60)(2)+(0.40)(3) =1.20+1.20 =2.40\text{ units}. $$

Under the stated assumptions, the estimate is an average treatment effect of 2.40 outcome units for this target population. It is not the equal-weight average \((2+3)/2=2.50\), because the two strata do not have equal target shares.

Target Shares Are Not Always Sample Shares

Using \(n_k/n\) is appropriate when the study sample represents the target population’s stratum distribution. But a sample can deliberately include different numbers of people from different strata—for example, to ensure enough observations in a relatively uncommon group. In that case, the sample shares may not describe the target population. If the target shares are known, retain them when combining the stratum effects. Otherwise, the result describes a different mixture of strata.

Worked Example: Correcting for a Different Target Mix

Consider a hypothetical study with 80% of its participants in stratum 1 and 20% in stratum 2. Suppose the estimated effects are 1 unit in stratum 1 and 5 units in stratum 2. The target population, however, is known to have equal shares in the two strata. Assume the identification conditions hold in both strata and both treatment groups are observed.

If the study-sample shares were used, the calculation would be

$$ (0.80)(1)+(0.20)(5)=0.80+1.00=1.80\text{ units}. $$

That is the effect averaged over the sample’s 80%/20% stratum mix. For the specified target population, use \(q_1=q_2=0.50\):

$$ \hat\tau_{\mathrm{ATE}} =(0.50)(1)+(0.50)(5) =0.50+2.50 =3.00\text{ units}. $$

The target-population ATE estimate is 3.00 units. The difference from 1.80 units is not a contradiction: the two calculations average the same stratum effects over different populations.

Diagnostics and Failure Modes

Before combining estimates, check whether each stratum supplies a credible treated-versus-control comparison and whether the weights correspond to the intended target. These checks do not prove the causal assumptions, but they can reveal when the calculation is unsupported or answers the wrong question.

  • Check both treatment groups within every contributing stratum. If a stratum has no treated or no control observations, its observed contrast cannot be estimated. A very small group can also make its estimated mean highly unstable.
  • Inspect overlap and covariate balance within strata. Compare the distributions of important pre-treatment covariates between treated and control units within each stratum. Persistent imbalances warn that the stratum may not adequately control confounding. Positivity requires support for both treatment levels; a score bin alone does not ensure it.
  • Verify the target shares. Confirm that the \(q_k\) values sum to one and describe the intended population. If they are estimated from the study sample, consider whether the sample composition represents that population.
  • Review the stratum-level estimates before aggregation. A large overall estimate may be driven by a stratum with a large effect, a large target share, or both. Report the component effects and shares so the aggregate can be understood.

When treatment positivity fails, do not insert a zero effect for the missing contrast. Nor does a weighted sum repair an unobserved stratum effect. The target ATE is not identified by these within-stratum comparisons unless additional defensible assumptions or information are introduced. Restricting the target population to strata with adequate support changes the estimand and should be stated explicitly.

Worked Example: A Stratum Without Controls

Suppose a hypothetical target population has two strata with shares \(q_1=0.70\) and \(q_2=0.30\). In stratum 1, the estimated treated mean is 7 units and the control mean is 5 units, so \(\hat\tau_1=2\) units. In stratum 2, all observed people are treated; there are no controls.

The first stratum’s contribution is calculable: \((0.70)(2)=1.40\) units. But \(\bar Y_{0,2}\) does not exist in the observed data, so \(\hat\tau_2=\bar Y_{1,2}-\bar Y_{0,2}\) cannot be calculated. Because \(q_2=0.30>0\), the target ATE formula requires that missing stratum effect:

$$ \hat\tau_{\mathrm{ATE}}=(0.70)(2)+(0.30)\hat\tau_2. $$

The data therefore do not yield a complete stratum-weighted ATE estimate. Reporting 1.40 units as the ATE would silently omit stratum 2 and fail to normalize the remaining weight. The appropriate conclusion is that the requested ATE is not identified by these observed within-stratum comparisons under the stated conditions.

How This Aggregation Relates to Other Approaches

The stratum-weighted estimator is a form of standardization: it uses within-stratum outcome contrasts and averages them over a specified distribution of strata. A pooled comparison that ignores strata generally does not make the same adjustment, and it need not target the same quantity. Inverse-probability weighting instead uses treatment-probability weights to construct a weighted outcome contrast; it is a different estimator, although under appropriate identification conditions it can target the same ATE. Agreement between methods is not guaranteed in finite samples, especially when overlap, model specification, or target-population representation differs.

Equal weighting of strata is not inherently invalid, but it answers a different question unless the target population gives every stratum equal share. Similarly, using sample shares when the sample was drawn with a different stratum mix changes the population being averaged over. Always identify the target population first, then use its stratum shares.

Common Mistake

Common mistake: “There are three strata, so I should average the three stratum effects equally.” This is correct only if the target population assigns one-third of its members to each stratum. Otherwise, use \(\sum_k q_k\hat\tau_k\), with each \(q_k\) representing the target population’s share in stratum \(k\). A related error is to use the study’s sample shares automatically even when that sample does not represent the target mix.
Key takeaway: An ATE across strata is the sum of each stratum-specific effect multiplied by that stratum’s share of the target population. The weighting rule, the within-stratum causal assumptions, and adequate representation of both treatment groups are all necessary for a defensible interpretation.

Check Your Understanding

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

  1. Two strata have target shares 0.25 and 0.75 and estimated effects of 4 and 2 outcome units. What is the combined ATE estimate?
  2. Why might an equal-weight average of stratum effects fail to estimate the target-population ATE?
  3. A study has sample stratum shares 0.90 and 0.10, while the target shares are 0.50 and 0.50. Which shares should be used for the target ATE, and why?
  4. What does an empty control group in a stratum with a positive target share imply for the usual stratum-weighted estimator?

Solutions. (1) \((0.25)(4)+(0.75)(2)=1+1.5=2.50\) outcome units. (2) Equal weighting gives each stratum one-third, or another equal share depending on the number of strata, whether or not that matches the target population’s composition; it therefore generally averages over the wrong distribution. (3) Use 0.50 and 0.50, because the requested estimand is for the target population, not the differently composed study sample. (4) The control mean and hence that stratum’s effect cannot be estimated from the observed comparison; the full target ATE is not identified by this estimator without additional defensible assumptions or information.

Clinical Biostats