Tutorials › Propensity Score Matching and Instrumental Variables › Treatment Prediction Versus Confounder Control

Time-Varying and High-Dimensional Confounding · Tutorial 619 of 699

Treatment Prediction Versus Confounder Control

Distinguish treatment prediction from confounder control, and assess how predictive instruments and weak confounders can affect propensity-score adjustment.

Advanced 12 min read

What You'll Learn

  • Explain why treatment-prediction accuracy does not establish conditional exchangeability.
  • Derive the balancing property of the propensity score for measured covariates.
  • Show how a strong treatment predictor can leave an important confounder uncontrolled.
  • Calculate how extreme propensity scores reduce inverse-weighted effective sample size.
  • Recognize why a weak treatment predictor can still be an important confounder.
  • Compare propensity adjustment with outcome adjustment and instrumental-variable approaches.

Why Treatment Prediction Is Not the Causal Objective

A propensity model estimates how treatment assignment varies with measured characteristics. It is tempting to judge such a model by how accurately it predicts who was treated. But accurate classification answers a prediction question, not the causal question: have treated and untreated people been made comparable with respect to causes of the outcome? A model can predict treatment very well because it captures a powerful treatment-assignment mechanism, yet leave an important confounder uncontrolled. Conversely, a measured confounder may add little to treatment prediction while substantially affecting the outcome.

Let \(A\) be a binary treatment, \(Y\) the observed outcome, and \(Y(1)\) and \(Y(0)\) the potential outcomes under treatment and control. Let \(X\) denote measured pretreatment covariates, and define the propensity score \(e(X)=P(A=1\mid X)\). For an average treatment effect, a central target is \(E[Y(1)-Y(0)]\). Estimating it from observed outcomes requires assumptions about treatment assignment and the covariates used to adjust for it.

Definition: The propensity score \(e(X)=P(A=1\mid X)\) is the conditional probability of receiving treatment given measured covariates \(X\). It is a balancing score: conditional on the true score, treated and untreated groups have the same distribution of the measured \(X\). This property alone does not guarantee balance on omitted variables or identify a causal effect.

For identification by adjustment, the covariates must be measured before treatment, treatment assignment must be conditionally exchangeable given those covariates, and both treatment options must be possible at covariate values relevant to the target population. In potential-outcome notation, these conditions include consistency, \(Y=AY(1)+(1-A)Y(0)\); conditional exchangeability, \((Y(1),Y(0))\perp A\mid X\); and positivity, \(0<P(A=1\mid X)<1\). Consistency also requires a sufficiently well-defined treatment, and the usual individual-level interpretation assumes no interference between units. The propensity model must estimate the relevant assignment probabilities adequately for the chosen adjustment procedure. Prediction accuracy is not a substitute for any of these assumptions.

The Balancing Property—and Its Boundary

The balancing property has a short derivation. Fix a score value \(e(X)=s\), where \(0<s<1\), and any set of covariate values \(B\). The probability of both being treated and having \(X\) in \(B\), conditional on the score, is \(sP(X\in B\mid e(X)=s)\). Dividing by the probability of treatment at that score gives the conditional distribution among treated people:

$$ P(X\in B\mid A=1,e(X)=s) =\frac{P(A=1,X\in B\mid e(X)=s)}{P(A=1\mid e(X)=s)} =\frac{sP(X\in B\mid e(X)=s)}{s} =P(X\in B\mid e(X)=s). $$

For controls, the same argument uses \(1-s\), so their distribution of \(X\) also equals \(P(X\in B\mid e(X)=s)\). Thus treatment status is independent of the measured \(X\) within an exact propensity-score stratum. This is a distributional statement, not just a statement about means.

The boundary is crucial: the result concerns the covariates used to define the score. If \(X\) contains only a strong treatment predictor \(Z\) and omits a confounder \(U\), balancing \(Z\) does not balance \(U\). A high area under the receiver operating characteristic curve, classification accuracy, or treatment-model likelihood does not establish conditional exchangeability. The introductory material on propensity scores establishes their role in adjustment; here the issue is why optimizing treatment prediction alone may undermine that role.

Worked Example: A Highly Predictive Instrument Leaves Confounding

Consider a hypothetical large population with two independent, equally prevalent binary variables: \(Z\), a treatment-predictive instrument, and \(U\), an unmeasured confounder. The table gives treatment probabilities in each covariate cell.

\(Z\)\(U\)\(P(A=1\mid Z,U)\)
000.01
010.19
100.81
110.99

A rule that predicts treatment for \(Z=1\) and control for \(Z=0\) is correct with probability \(0.90\), since the treatment rates are \(0.90\) for \(Z=1\) and \(0.10\) for \(Z=0\). Yet \(Z\) is independent of \(U\), and \(Z\) has no direct effect on the outcome in this example. Suppose \(Y=4A+12U+\varepsilon\), where \(\varepsilon\) has mean zero conditional on \(A,Z,U\). The treatment effect for every unit is 4 outcome units.

Because \(P(A=1)=0.50\), the treated-group prevalence of \(U=1\) is \(\{0.5[0.5(0.19)+0.5(0.99)]\}/0.50=0.59\). Among controls it is \(\{0.5[0.5(0.81)+0.5(0.01)]\}/0.50=0.41\). The unadjusted difference in mean outcomes is therefore \(4+12(0.59-0.41)=6.16\) outcome units, rather than the causal effect of 4.

Adjusting for \(Z\) alone does not fix the problem. Within \(Z=0\), the prevalence of \(U=1\) is \(0.19/(0.01+0.19)=0.95\) among treated people and \(0.81/(0.19+0.81)=0.81\) is not the control prevalence: the control probability mass for \(U=1\) is \(0.5(1-0.19)\), while total control mass is \(0.5(1-0.01)+0.5(1-0.19)=0.90\), giving \(0.405/0.90=0.45\). Thus the within-\(Z=0\) outcome contrast is \(4+12(0.95-0.45)=10\). Within \(Z=1\), the corresponding treated and control prevalences are \(0.99/(0.81+0.99)=0.55\) and \(0.01/(0.19+0.01)=0.05\), again yielding \(4+12(0.50)=10\). Adjustment using \(Z\) alone can therefore produce a contrast of 10, despite the 90% treatment-classification accuracy of \(Z\). The omitted \(U\), not prediction accuracy, determines the remaining confounding.

Strong Predictors, Poor Overlap, and Weight Instability

Even when a strong treatment predictor is not a confounder, including it in a propensity model can have a practical cost. If it drives estimated scores close to zero or one, few comparable treated and untreated people remain in parts of the covariate space. In inverse probability weighting for the average treatment effect, treated observations receive weight \(1/e(X)\) and controls receive weight \(1/[1-e(X)]\). Scores near the boundaries can consequently produce large weights, high variance, and a target estimate dominated by a few observations. Under ideal population conditions, including a valid instrument alongside all confounders does not by itself invalidate identification; the concern is often loss of precision, finite-sample instability, or amplified bias when confounding remains unresolved.

One useful summary is the weighted effective sample size, \(ESS=(\sum_i w_i)^2/\sum_i w_i^2\). It describes how much information the weights retain relative to an equally weighted sample; it is not a replacement for checking balance or overlap.

Worked Example: A Predictive Covariate Shrinks Effective Sample Size

Suppose a hypothetical cohort has 1,000 people, 500 treated and 500 controls. The estimated propensity scores are 0.90 for 500 people and 0.10 for the other 500. Suppose the realized treated group contains 450 people with score 0.90 and 50 with score 0.10; the controls contain 50 and 450 people in those respective score groups. These counts are consistent with strong treatment prediction: the score-0.90 group is mostly treated and the score-0.10 group mostly control.

For the treated group, the total of the ATE weights is \(450/0.90+50/0.10=500+500=1{,}000\). The sum of squared weights is \(450(1/0.90)^2+50(1/0.10)^2=555.6+5{,}000=5{,}555.6\). Therefore \(ESS_{\text{treated}}=1{,}000^2/5{,}555.6\approx180\), far below the 500 treated observations counted without weights.

For controls, the weights are \(1/(1-e)\). Their total weight is \(50/0.10+450/0.90=1{,}000\), and their squared-weight sum is \(50(1/0.10)^2+450(1/0.90)^2=5{,}555.6\). Thus \(ESS_{\text{control}}\approx180\) as well. The calculation does not prove that the estimate is biased, but it signals that weighted precision may be poor and that results may be sensitive to a small number of observations.

In practice, inspect the estimated-score distributions by treatment group, the number and characteristics of observations near the boundaries, weight summaries, and effective sample sizes. After matching, weighting, or subclassification, assess covariate balance directly, commonly with standardized mean differences and distributional comparisons. Check important confounders and prognostic covariates rather than relying on a single overall balance summary. Large residual imbalances, extreme weights, or very small effective sample sizes are warnings that the design may not support stable adjustment for the intended target population.

Weak Treatment Predictors Can Be Important Confounders

Treatment prediction also provides a poor filter for deciding whether a covariate matters causally. A variable can have a weak association with treatment but a strong association with the outcome. If it is a common cause of treatment and outcome, even a modest treatment association can generate meaningful confounding. The direction and importance of confounding depend on the joint relationship with treatment and outcome, not on treatment-model performance alone.

Worked Example: A Weak Predictor Creates Meaningful Bias

Suppose a hypothetical population has a binary pretreatment confounder \(U\), with \(P(U=1)=0.50\). Treatment probabilities are \(P(A=1\mid U=1)=0.55\) and \(P(A=1\mid U=0)=0.45\). Classifying treatment as \(A=1\) when \(U=1\) gives only 55% accuracy. Let the outcome be \(Y=5A+20U+\varepsilon\), with mean-zero error conditional on \(A,U\); treatment has a constant causal effect of 5 outcome units, while \(U\) shifts the outcome by 20 units.

The overall treatment probability is \(0.50\). The fraction with \(U=1\) among treated people is \(0.5(0.55)/0.50=0.55\); among controls it is \(0.5(0.45)/0.50=0.45\). The unadjusted mean difference is \(5+20(0.55-0.45)=7\) outcome units. Within either level of \(U\), the treatment contrast is 5. Standardizing those within-\(U\) contrasts to the same \(U\) distribution therefore recovers 5 under the stated model and assumptions. Despite its weak treatment prediction, \(U\) is important for confounding control because it strongly predicts the outcome.

This distinction informs model assessment: treatment discrimination and causal balance are different objectives. A treatment model can be compared using predictive measures, but the design should be assessed through overlap and covariate balance. Outcome prognostic information can help prioritize which imbalances are consequential, while a variable's causal role still depends on the assignment and outcome structure. Adding variables indiscriminately is not a solution: covariates should be pretreatment, and adjustment must respect the causal structure.

Choosing and Diagnosing an Adjustment Strategy

Propensity-score matching, weighting, and subclassification use different implementations of the same broad idea: compare treatment groups after adjustment for measured pretreatment covariates. None repairs an omitted confounder merely by achieving excellent treatment prediction. Outcome regression models the outcome conditional on covariates instead; it may be more efficient in some settings, but it also relies on adequate covariate measurement and model specification. A carefully designed analysis can use propensity adjustment, outcome adjustment, or both, while explicitly stating the target population and assumptions.

Instrumental-variable analysis addresses a different identification strategy. A variable that predicts treatment but is independent of outcome causes and affects the outcome only through treatment may serve as an instrument under additional assumptions. Its treatment-prediction strength is relevant there, but it does not make it a confounder-control variable. Simply adding an instrument to a propensity model does not convert propensity adjustment into an instrumental-variable analysis. When residual confounding is present, conditioning on a strong instrument can, in some settings, amplify bias; it can also worsen overlap. Neither effect is automatic in every design, so the causal structure and observed diagnostics matter.

Common mistake: “This propensity model has high treatment-classification accuracy, so it must control confounding.” Accuracy describes prediction of \(A\), not comparability in the causes of \(Y\). Instead, justify the pretreatment covariate set, assess exchangeability as a substantive assumption, and inspect balance, overlap, and weight stability. A low-accuracy confounder may still be essential; a high-accuracy instrument may add little to confounding control and may reduce precision.

The limits of diagnostics should also be explicit. Balance checks evaluate measured covariates, not unmeasured confounding. Overlap plots cannot prove that positivity holds for all relevant covariate patterns, and acceptable effective sample size does not establish correct model specification. When estimates rely on trimming or restricting the sample to improve overlap, the target population may change and should be described. Robust causal interpretation requires substantive knowledge alongside statistical checks.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. A propensity model has excellent treatment discrimination but omits a measured pretreatment cause of both treatment and outcome. Does discrimination establish exchangeability?
  2. In the weak-confounder example, what is the unadjusted mean difference, and what is the within-\(U\) treatment contrast?
  3. For 450 treated observations with \(e=0.90\) and 50 with \(e=0.10\), calculate the treated-group effective sample size using ATE weights.
  4. Why can including a strong valid instrument in propensity adjustment be problematic even if it is not itself a confounder?

Solutions. (1) No. Discrimination concerns prediction of treatment; the omitted common cause can leave the groups nonexchangeable. (2) The unadjusted difference is \(5+20(0.55-0.45)=7\) outcome units, while the within-\(U\) contrast is 5 units. (3) The weight sum is \(450/0.90+50/0.10=1{,}000\), and the squared-weight sum is \(450/0.90^2+50/0.10^2=5{,}555.6\), so \(ESS=1{,}000^2/5{,}555.6\approx180\). (4) It can push scores toward zero or one, reduce overlap and effective sample size, and increase variance; with residual confounding, it can amplify bias in some settings. It does not automatically do so in every design.

Clinical Biostats