Tutorials › Econometrics (Graduate) › Calibration and Validation of CATE Predictions

Treatment Effect Heterogeneity · Tutorial 842 of 952

Calibration and Validation of CATE Predictions

You will learn to evaluate CATE predictions on held-out data, interpret calibration summaries, and diagnose overlap-related instability.

Advanced 12 min read

What You'll Learn

  • Define calibration of predicted conditional average treatment effects.
  • State the assumptions that identify a CATE from observed data.
  • Use an orthogonal pseudo-outcome to check calibration on held-out observations.
  • Interpret calibration-in-the-large, calibration slopes, and grouped calibration plots.
  • Diagnose how weak overlap can make validation summaries unstable.
  • Distinguish evidence about average calibration from evidence about individual predictions.

Why Validate Conditional Effect Predictions?

A model can produce a conditional average treatment effect (CATE) prediction for every covariate profile and still be poorly calibrated. It may rank people reasonably well while systematically overstating effect sizes, or it may look accurate in the data used to fit it but fail on new observations. Validation asks whether predicted effects correspond to observed effect differences in a population not used to construct those predictions.

The target is a conditional mean effect, not an observable individual treatment effect. Write \(\tau(x)=E[Y(1)-Y(0)\mid X=x]\), where \(Y(1)\) and \(Y(0)\) are potential outcomes and \(X\) denotes observed pretreatment covariates. A fitted model, evaluated at \(X_i\), gives a prediction \(\hat{\tau}(X_i)\). Validation compares these predictions with an appropriately adjusted estimate of the conditional effect, not with the individual’s observed outcome or an unadjusted treated-versus-control difference.

Definition: A CATE prediction is calibrated over a target population if its predicted values agree, on average and across relevant prediction levels, with the corresponding conditional average treatment effects in that population. Calibration is a property of predictions relative to a specified population and validation procedure; it does not establish that every prediction is accurate or that the model identifies individual treatment effects.

Identification and a Validation Pseudo-Outcome

For observed data \((Y,W,X)\), let \(W\) be a binary treatment indicator and define \(e(x)=P(W=1\mid X=x)\) and \(m_w(x)=E[Y\mid W=w,X=x]\) for \(w\in\{0,1\}\). To interpret \(m_1(x)-m_0(x)\) as \(\tau(x)\), the usual identification conditions are needed: consistency, no interference between units, conditional exchangeability \((Y(1),Y(0))\perp W\mid X\), and positivity \(0<e(X)<1\) in the target population. These conditions connect observed conditional outcome means to potential-outcome means; they do not guarantee precise estimation where treatment probabilities are near zero or one.

A useful validation quantity is the doubly robust, or augmented inverse-probability-weighted, pseudo-outcome. If the nuisance functions \(e,m_1,m_0\) are known, define

$$ \Gamma = m_1(X)-m_0(X) +\frac{W}{e(X)}\{Y-m_1(X)\} -\frac{1-W}{1-e(X)}\{Y-m_0(X)\}. $$

In practice, nuisance functions are estimated. For validation, use estimates obtained without fitting to the validation observation—for example, fit nuisance models on training folds and predict them for held-out folds. This cross-fitting helps prevent overfitting from making the pseudo-outcome appear artificially aligned with the CATE prediction. The CATE prediction itself should also be produced without using the validation outcome.

Result: Under consistency, no interference, conditional exchangeability, positivity, and correctly defined conditional means and propensity, the pseudo-outcome satisfies \(E[\Gamma\mid X]=\tau(X)\). Consequently, \(E[\Gamma\mid\hat{\tau}(X)]\) describes the average conditional effect among observations with a given prediction, provided the predictions are evaluated in the target population.

The identity follows by conditioning on \(X\). The treated residual term has conditional expectation zero:

$$ E\left[\frac{W}{e(X)}\{Y-m_1(X)\}\mid X\right] = \frac{e(X)}{e(X)}E[Y-m_1(X)\mid W=1,X] =0. $$

The same argument gives zero for the control residual term, because \(E[Y-m_0(X)\mid W=0,X]=0\) and \(P(W=0\mid X)=1-e(X)\). The remaining term is \(m_1(X)-m_0(X)=\tau(X)\) under the identification conditions. Thus \(\Gamma\) is an unbiased, conditional-mean target for calibration. It is not a measured treatment effect for one person: it can be noisy, and its value depends on the observed outcome and treatment assignment.

Worked Example: Computing a Held-Out Pseudo-Outcome

Suppose a hypothetical held-out patient has \(W=1\), an outcome of \(Y=10\) blood-pressure units, fitted conditional outcome means \(m_1(X)=8\) and \(m_0(X)=5\) units, and estimated treatment probability \(e(X)=0.05\). Assume the identification conditions are plausible for the target population and the nuisance predictions were obtained without fitting to this patient.

Substitute into the pseudo-outcome formula:

$$ \Gamma = (8-5)+\frac{1}{0.05}(10-8)-\frac{0}{0.95}(10-5) =3+20(2)-0 =43\text{ units}. $$

The value \(43\) is not a plausible claim that this patient’s treatment effect was 43 units. It is one noisy validation score; its large inverse-probability contribution arises because the patient was treated despite a fitted treatment probability of only \(0.05\). The score can still have the right conditional expectation under the assumptions, but a few such values can dominate a sample average. The calculation therefore motivates checking overlap and score variability before interpreting calibration summaries.

Calibration-in-the-Large and Calibration Slope

On held-out observations, let \(\hat{\tau}_i\) denote the predicted effect and \(\Gamma_i\) its pseudo-outcome. A simple calibration-in-the-large summary is the mean prediction error:

$$ \widehat{\mathrm{CIL}} = \frac{1}{n_{\mathrm{val}}}\sum_{i\in\mathrm{val}}(\Gamma_i-\hat{\tau}_i). $$

Here \(n_{\mathrm{val}}\) is the number of validation observations. A value near zero indicates little average offset in the validation sample, while a positive value indicates that the pseudo-outcomes average above predictions. This summary can conceal important variation: overprediction in one region may cancel underprediction in another.

A calibration slope summarizes whether differences in predictions track differences in pseudo-outcomes. Fit the validation regression

$$ \Gamma_i=a+b\hat{\tau}_i+\varepsilon_i. $$

The ideal linear calibration relationship has \(a=0\) and \(b=1\). A slope below one often indicates predictions that vary more than the corresponding validation effects, while a slope above one often indicates predictions that vary less. Those interpretations are conditional on the fitted intercept and on the linear summary being useful. A slope is not a test of all possible calibration failures: nonlinear departures can remain even when the slope is one. Since \(\Gamma_i\) is noisy, uncertainty in the slope can also be substantial, especially with limited validation data or weak overlap.

Worked Example: Estimating a Calibration Slope

Consider four hypothetical held-out observations. Their predicted effects and pseudo-outcomes are in blood-pressure units.

ObservationPrediction \(\hat{\tau}_i\)Pseudo-outcome \(\Gamma_i\)
111.0
221.5
332.5
443.0

The prediction mean is \(\bar{\hat{\tau}}=2.5\) units and the pseudo-outcome mean is \(\bar{\Gamma}=2.0\) units. The least-squares slope is the sum of cross-products of centered values divided by the sum of squared centered predictions:

$$ \hat{b} = \frac{(-1.5)(-1.0)+(-0.5)(-0.5)+(0.5)(0.5)+(1.5)(1.0)} {(-1.5)^2+(-0.5)^2+(0.5)^2+(1.5)^2} = \frac{3.5}{5} =0.70. $$

The fitted intercept is \(\hat{a}=\bar{\Gamma}-\hat{b}\bar{\hat{\tau}}=2.0-(0.70)(2.5)=0.25\) units. The sample mean prediction error is \[ \widehat{\mathrm{CIL}} =\frac{(1.0-1)+(1.5-2)+(2.5-3)+(3.0-4)}{4} =\frac{-2.0}{4} =-0.50\text{ units}. \]

These small illustrative data show predictions that average \(0.50\) units above their pseudo-outcomes and have a slope below one. The slope suggests that predicted variation is relatively large compared with the validation variation in this sample. Four observations do not support a stable general conclusion; the arithmetic illustrates the summaries, not a basis for claiming population miscalibration.

Grouped Calibration Plots

A grouped plot offers a visual check that does not force the entire calibration relationship to be linear. Sort held-out observations by \(\hat{\tau}_i\), divide them into groups (often prediction quantiles), and calculate within each group the mean prediction and mean pseudo-outcome. Plot the latter against the former, with the 45-degree line as the ideal reference. Group sizes and uncertainty intervals should be shown or reported where feasible; intervals should account for the validation design, including clustering when relevant.

The plot estimates average behavior within prediction ranges, not accuracy for each observation. Grouping smooths noisy individual pseudo-outcomes, but it also hides within-group variation and can make results sensitive to the number and boundaries of groups. A calibration curve that departs from the reference line may reveal systematic over- or underprediction in a region even when the overall mean error is near zero.

Worked Example: Reading a Grouped Calibration Plot

In a hypothetical held-out sample, predictions are divided into three equally sized groups of five observations. The table reports the mean prediction and mean pseudo-outcome in each group, in blood-pressure units.

Prediction groupMean predictionMean pseudo-outcomePseudo-outcome minus prediction
Lowest2.01.8\(1.8-2.0=-0.2\)
Middle5.04.9\(4.9-5.0=-0.1\)
Highest8.07.2\(7.2-8.0=-0.8\)

The plotted points are \((2.0,1.8)\), \((5.0,4.9)\), and \((8.0,7.2)\). Each lies below the 45-degree line, and the largest group shows the greatest gap. Since the groups have equal size, the overall average of the three group differences is \((-0.2-0.1-0.8)/3=-0.367\) units, matching the mean error across all 15 observations if the displayed values are exact group means. The pattern suggests overprediction, particularly at the high-prediction end. It does not establish that every observation in that group has an effect of 7.2 units or less.

Held-Out Design and Overlap Diagnostics

A credible assessment separates model development from validation. Fit the CATE prediction model using training data, then evaluate predictions on observations whose outcomes did not inform that model. In cross-validation, each observation’s prediction should be made by a model that excluded that observation. If nuisance models are also estimated, their predictions for the validation observations should likewise avoid using those observations’ outcomes. Repeatedly adjusting a model in response to the same validation results turns that set into development data; a final assessment then requires a fresh holdout or a clearly defined resampling procedure.

The validation sample should represent the population for which the predictions are intended. If the covariate distribution or treatment assignment mechanism differs, calibration summaries may change even if the fitted model is unchanged. Report the target population, the validation construction, and the distribution of predicted effects. For observational data, inspect treatment probabilities and the size of inverse-probability contributions; positivity in theory does not prevent extreme weights in finite samples.

Weak overlap can make pseudo-outcomes highly variable because residuals are divided by \(e(X)\) for treated observations and by \(1-e(X)\) for controls. In such regions, a few observations may drive a slope or group mean. Report where predictions are supported by data, inspect group counts and propensity distributions, and consider whether the estimand should be restricted to a population with adequate overlap. Trimming or restricting changes the target population and must be stated; it is not merely a numerical repair.

Common Pitfalls and Interpretation

  • Validating on training outcomes. A close fit on observations used to learn the CATE model can reflect overfitting. Use held-out or properly cross-fitted predictions.
  • Treating a pseudo-outcome as an individual effect. \(\Gamma_i\) is a noisy score with a conditional-mean property, not the unobserved \(Y_i(1)-Y_i(0)\).
  • Using raw outcome differences as validation. In observational data, differences by treatment status can reflect confounding. Use an adjustment that targets the conditional effect under stated assumptions.
  • Calling a slope of one proof of calibration. A linear slope summarizes one feature of the relationship and can miss nonlinear miscalibration or offsetting errors across groups.
  • Ignoring weak overlap. Large inverse-probability factors can make validation summaries unstable. Report this limitation rather than interpreting an extreme pseudo-outcome as a large treatment effect.
  • Overlooking the target population. Calibration in one population does not automatically transfer to a population with a different covariate or treatment distribution.

Calibration and validation complement model-building diagnostics: they ask whether predicted conditional effects correspond to effect differences in data not used to create the predictions. No single summary answers that question completely. A useful assessment combines held-out evaluation, a mean-error summary, a slope, grouped plots, and explicit checks of overlap and target-population relevance.

Check Your Understanding

Use the definitions and calculations in this tutorial to answer the following questions.

  1. What assumptions allow \(E[\Gamma\mid X]\) to be interpreted as a CATE?
  2. Why is a pseudo-outcome of 43 units not an estimate of one patient’s individual treatment effect?
  3. In the calibration regression \(\Gamma_i=a+b\hat{\tau}_i+\varepsilon_i\), what values of \(a\) and \(b\) represent ideal linear calibration?
  4. How could a grouped plot reveal a problem that a near-zero calibration-in-the-large summary would conceal?
  5. Why can a small estimated treatment probability make a held-out calibration summary unstable?
  6. What changes when validation results are used repeatedly to modify the prediction model?
Clinical Biostats