Why Augment Inverse-Probability Weighting?
Inverse-probability weighting estimates an average treatment effect by reweighting observed outcomes according to treatment probabilities. Outcome regression takes a different route: it predicts each unit’s outcome under treatment and under control, then averages the predicted differences. The augmented inverse-probability-weighted (AIPW) estimator combines these ideas. It starts with predicted treatment effects and adds weighted residual corrections that account for the difference between observed outcomes and those predictions.
The augmentation is not simply an extra adjustment to an already-computed effect. Its form follows from an identity for the average treatment effect under standard causal identification assumptions. This tutorial derives that identity, turns it into a sample estimator, and works through the calculation. The focus is constructing and interpreting the estimator, not conducting inference about it.
Setup and Identification Assumptions
For unit \(i\), let \(A_i\) equal 1 if the unit receives treatment and 0 otherwise, let \(X_i\) denote measured pre-treatment covariates, and let \(Y_i\) be the observed outcome. Write \(Y_i(1)\) and \(Y_i(0)\) for the potential outcomes under treatment and control. The target is the average treatment effect (ATE), \(E[Y(1)-Y(0)]\), for the population represented by the data.
The causal interpretation of the estimator relies on assumptions about how treatment is assigned and outcomes are observed. Consistency means the observed outcome is the potential outcome corresponding to the treatment actually received: \(Y=A Y(1)+(1-A)Y(0)\). The usual no-interference condition means one unit’s treatment does not change another unit’s potential outcomes. Conditional exchangeability means \((Y(1),Y(0))\) is independent of \(A\) given \(X\): after conditioning on the measured covariates, treatment assignment carries no additional information about the potential outcomes. Positivity requires \(0<e(X)<1\) throughout the target population, so units with each covariate pattern have a chance of receiving either treatment.
Together, these assumptions identify the ATE from observed data. They do not guarantee that a particular finite-sample estimate is accurate: in practice, propensity scores and outcome regressions must be estimated, and limited overlap can make weighted residual corrections unstable.
Deriving the Augmentation
Begin with the regression-based contrast \(m_1(X)-m_0(X)\). It predicts the difference between treatment and control outcomes for a unit with covariates \(X\). The observed treated outcome \(Y\), when \(A=1\), can provide a correction to the prediction \(m_1(X)\); the observed control outcome, when \(A=0\), can correct \(m_0(X)\). To use those residuals to represent the corresponding population differences, weight them by the inverse probability of observing that treatment status.
With the true propensity score \(e(X)\), conditional exchangeability and consistency imply the two identities below. For the treated residual, conditioning on \(X\) gives \(E[A\{Y-m_1(X)\}/e(X)\mid X]=E[Y(1)-m_1(X)\mid X]\). For the control residual, it gives \(E[(1-A)\{Y-m_0(X)\}/\{1-e(X)\}\mid X]=E[Y(0)-m_0(X)\mid X]\). These equalities follow because the factor \(A/e(X)\) selects treated units and corrects for their conditional treatment probability; the control factor works analogously.
Add the first residual correction to the predicted treatment outcome and subtract the second correction from the predicted control outcome. The outcome-regression predictions then cancel with the corresponding parts of the residual expectations, leaving the conditional mean treatment effect. Averaging over \(X\) gives the ATE.
This is the construction’s central justification: a predicted contrast is paired with weighted residual corrections, and their conditional expectation recovers the conditional causal contrast. Averaging that expression over the covariate distribution targets the population ATE.
Each unit contributes one estimated effect score: a predicted difference plus a treated residual correction or a control residual correction, but not both. The estimator is the ordinary sample average of those scores. In particular, this expression does not divide the treated and control weights by their separate sums; that would be a different normalization.
Calculating the Estimator
Worked Example: Four Individual Contributions
Suppose a hypothetical sample of four patients has a binary outcome measured in points on a hypothetical symptom scale. The table gives the observed treatment, outcome, estimated propensity score, and estimated outcome predictions. All propensity scores lie strictly between 0 and 1, as required for the displayed weights.
| Unit | \(A_i\) | \(Y_i\) (points) | \(\hat e(X_i)\) | \(\hat m_1(X_i)\) | \(\hat m_0(X_i)\) |
|---|---|---|---|---|---|
| 1 | 1 | 14 | 0.60 | 12 | 9 |
| 2 | 0 | 8 | 0.40 | 11 | 8 |
| 3 | 1 | 10 | 0.50 | 13 | 7 |
| 4 | 0 | 5 | 0.25 | 9 | 6 |
For unit 1, the predicted contrast is \(12-9=3\) points. The unit is treated, so its treated residual is weighted by \(1/0.60\): \(3+(14-12)/0.60=6.333\) points.
For unit 2, the contrast is \(11-8=3\) points. The unit is a control, so the control residual is subtracted with weight \(1/(1-0.40)\): \(3-(8-8)/0.60=3\) points.
For unit 3, the contrast is \(13-7=6\) points. Its treated residual correction is \((10-13)/0.50=-6\) points, giving \(6-6=0\) points.
For unit 4, the contrast is \(9-6=3\) points. Its control residual is \((5-6)/(1-0.25)=-1.333\) points, and subtracting that negative correction gives \(3-(-1.333)=4.333\) points.
Thus the sample estimate is \((6.333+3+0+4.333)/4=3.417\) points. Under the identification assumptions and for the population represented by this hypothetical sample, the estimate is an AIPW estimate of the average treatment effect. It is an estimated contrast in symptom-scale points, not a statement that treatment caused a 3.417-point change for every individual.
Seeing the Treated and Control Corrections Separately
The formula’s signs are easiest to check by writing the two cases separately. A treated unit has score \(\hat m_1-\hat m_0+(Y-\hat m_1)/\hat e\). A control unit has score \(\hat m_1-\hat m_0-(Y-\hat m_0)/(1-\hat e)\). The minus sign on the control residual is essential because the target contrast subtracts the control outcome.
Worked Example: A Treated and a Control Unit
Consider two hypothetical patients with outcomes measured in minutes of recovery time. For the treated patient, \(Y=20\), \(\hat e=0.50\), \(\hat m_1=18\), and \(\hat m_0=12\). The predicted contrast is \(18-12=6\) minutes and the treated correction is \((20-18)/0.50=4\) minutes. This patient’s score is \(6+4=10\) minutes.
For the control patient, \(Y=9\), \(\hat e=0.25\), \(\hat m_1=15\), and \(\hat m_0=10\). The predicted contrast is \(15-10=5\) minutes. The control residual correction before applying the sign is \((9-10)/(1-0.25)=-1/0.75=-1.333\) minutes. The score is therefore \(5-(-1.333)=6.333\) minutes.
If these were the only two units, the AIPW estimate would be \((10+6.333)/2=8.167\) minutes. This small calculation illustrates the score arithmetic; it does not establish that the estimate is precise or representative of a larger target population.
How AIPW Relates to Inverse-Probability Weighting
The AIPW expression can also be rearranged to show the inverse-probability-weighted outcomes explicitly. For a treated unit, the outcome part becomes \(\hat m_1+(Y-\hat m_1)/\hat e\); for a control unit it becomes \(\hat m_0+(Y-\hat m_0)/(1-\hat e)\), with the control quantity subtracted. The predictions thus serve as baselines, while the weighted residuals adjust those baselines using observed outcomes.
Worked Example: Comparing the Scores Algebraically
Suppose a hypothetical treated unit has \(Y=14\), \(\hat e=0.50\), \(\hat m_1=12\), and \(\hat m_0=9\), with outcome measured in points. The AIPW score is \(12-9+(14-12)/0.50=3+4=7\) points.
The same score written as a baseline prediction plus a residual correction is \(9+\{12+(14-12)/0.50-9\}=9+(12+4-9)=16\) only if the control baseline is added incorrectly. The correct rearrangement keeps the contrast baseline and adds the treated correction: \((12-9)+(14-12)/0.50=7\). Equivalently, its treated-outcome component is \(12+(14-12)/0.50=16\), from which the control prediction \(9\) is subtracted, again giving \(16-9=7\) points.
For a control unit, the parallel rearrangement is to retain the treated prediction and subtract the control component \(\hat m_0+(Y-\hat m_0)/(1-\hat e)\). This algebra is useful for checking that the regression baseline and the residual correction have not been counted twice.
Practical Checks and a Common Pitfall
Before interpreting an AIPW estimate, verify that the inputs and calculations match the target and data structure. The outcome regressions must predict the outcome under each treatment status using pre-treatment covariates; the propensity score must model treatment probability given those covariates. Predictions and observed outcomes must use the same outcome scale. Propensity estimates close to zero or one produce large inverse weights, so overlap and weight behavior matter even when every estimated probability is technically within the interval \((0,1)\).
A common mistake is to add the treated and control residual corrections with the same sign, or to use \(1/\hat e\) for both groups. The treated correction uses \(A_i/\hat e(X_i)\), while the control correction uses \((1-A_i)/\{1-\hat e(X_i)\}\) and is subtracted. Another error is to treat the average predicted contrast alone as the AIPW estimate; that omits the augmentation. Conversely, an estimate of the propensity score does not make the estimator automatically causal: the identification assumptions, covariate choice, and overlap still require substantive justification.
Check Your Understanding
Work through these questions using the definitions and formula above.
- Write the AIPW score for a treated unit and for a control unit. Which residual correction is subtracted?
- A treated unit has \(Y=11\), \(\hat e=0.50\), \(\hat m_1=9\), and \(\hat m_0=6\). Calculate its score in the outcome’s units.
- For a control unit with \(Y=7\), \(\hat e=0.30\), \(\hat m_1=10\), and \(\hat m_0=8\), calculate its score and explain the sign of its residual correction.
- Why does the treated residual use the inverse propensity score, while the control residual uses the inverse probability of control?
- Name the causal identification assumptions needed to interpret the estimator as an ATE, and explain why estimated propensities near zero or one deserve attention.