Tutorials › Causal Inference › Common Errors in Semiparametric Causal Estimation

Doubly Robust and Semiparametric Methods · Tutorial 489 of 972

Common Errors in Semiparametric Causal Estimation

Diagnose common failures in AIPW and related semiparametric causal analyses, and identify the conditions needed for reliable estimation and inference.

Advanced 10 min read

What You'll Learn

  • Explain why double robustness does not repair a positivity violation or guarantee stable estimates.
  • Check how nuisance-estimation errors combine in the product-rate condition for root-n inference.
  • Distinguish out-of-fold predictions from predictions evaluated on the same observations used to fit a model.
  • Calculate an influence-function standard error from estimated score contributions.
  • Identify when a standard error does not match the estimator or the dependence structure of the data.

Why Semiparametric Estimates Can Still Go Wrong

An estimator can have an attractive theoretical property and still be used incorrectly. In particular, double robustness is not a general guarantee of validity: it does not fix a lack of treatment overlap, remove every bias from flexible nuisance models, or make an inappropriate standard error trustworthy. Reliable inference requires checking the conditions that support identification, estimation, and uncertainty quantification separately.

This tutorial focuses on the average treatment effect (ATE) estimated with the augmented inverse probability weighted (AIPW) score. The same diagnostic habits apply to related semiparametric estimators. Let \(A\) be a binary treatment, \(Y\) an outcome, and \(X\) measured covariates. Write \(e_0(X)=P(A=1\mid X)\) for the true propensity score and \(m_{a,0}(X)=E(Y\mid A=a,X)\) for the true outcome regression under treatment level \(a\). Under consistency, conditional exchangeability given \(X\), and positivity, the ATE is identified by \(\psi_0=E[m_{1,0}(X)-m_{0,0}(X)]\).

Definition: Given estimates \(\hat e\), \(\hat m_1\), and \(\hat m_0\), an AIPW estimate averages the score contributions $$ \hat q_i=\hat m_1(X_i)-\hat m_0(X_i) +\frac{A_i}{\hat e(X_i)}\{Y_i-\hat m_1(X_i)\} -\frac{1-A_i}{1-\hat e(X_i)}\{Y_i-\hat m_0(X_i)\}, \qquad \hat\psi=\frac{1}{n}\sum_{i=1}^n\hat q_i. $$ The estimated influence-function contribution is \(\hat\phi_i=\hat q_i-\hat\psi\). The score has units of the outcome, as does the ATE.

The usual double-robustness claim is conditional: under the identification and regularity conditions, the AIPW estimator can be consistent if either the propensity model or the outcome regression is correctly specified (and the relevant fitted nuisance estimates behave appropriately). That claim is not the same as saying every finite-sample estimate is stable, or that a confidence interval is valid regardless of how the nuisance estimates were obtained.

Error 1: Treating Double Robustness as Protection Against Poor Overlap

Positivity requires that both treatment levels be possible at covariate values in the target population: \(0<e_0(X)<1\) almost surely. A stronger uniform-positivity condition requires some \(\epsilon>0\) such that \(\epsilon\le e_0(X)\le 1-\epsilon\) almost surely. When the true propensity is close to zero or one, the inverse-propensity factors in the AIPW score can be large. This can make estimates highly variable even when the propensity model is correct.

If a treatment level has probability zero in a covariate stratum, the observed data provide no treated-versus-untreated comparison there. Double robustness cannot supply the missing information or identify the ATE for that stratum without additional assumptions. If probabilities are positive but very small, identification may hold while estimation remains unstable.

Worked Example: A Large Score Contribution Under Near-Violation

Suppose, hypothetically, a person has an estimated propensity \(\hat e(X)=0.02\), is treated, and has outcome \(Y=10\). The fitted outcome regressions are \(\hat m_1(X)=8\) and \(\hat m_0(X)=6\), in outcome units. This person's AIPW contribution is

$$ \hat q=8-6+\frac{1}{0.02}(10-8)-0=2+50(2)=102. $$

The inverse-propensity correction is 100 outcome units, so the total contribution is 102. The arithmetic is correct, but one observation can now have substantial influence on the sample average. For comparison, if the same fitted values and residual occurred at \(\hat e(X)=0.5\), the contribution would be \(2+(1/0.5)(2)=6\). The contrast illustrates instability from a small propensity; it does not imply that replacing the propensity by 0.5 would be justified.

A common mistake is to say, “The estimator is doubly robust, so extreme propensities are harmless.” The correction is that double robustness concerns consistency under specified nuisance-model conditions. It does not bound inverse weights, guarantee small variance, or overcome structural nonpositivity.

Clipping estimated propensities to a chosen range may limit extreme weights, but it changes the nuisance estimates used in the score and does not create observations from unsupported treatment groups. Its effects should be treated as a sensitivity or modeling decision, not as proof that positivity holds.

Error 2: Checking Each Nuisance Rate but Ignoring Their Product

For root-\(n\) inference, it is not enough to say that each nuisance estimate is “reasonably accurate.” A standard sufficient condition, alongside appropriate consistency, moment, and positivity conditions, is that the product of the propensity and outcome-regression errors is small enough. In \(L_2\) notation, a common version is

$$ \|\hat e-e_0\|_2\,\|\hat m_a-m_{a,0}\|_2=o_p(n^{-1/2}),\qquad a\in\{0,1\}. $$

The notation \(o_p(n^{-1/2})\) means that the product is negligible relative to \(n^{-1/2}\) in probability as the sample size grows. If both errors are of order \(o_p(n^{-1/4})\), their product is \(o_p(n^{-1/2})\). The condition also allows one nuisance to converge faster than the other. The relevant error norms and inverse-propensity factors must be controlled; a rate statement alone does not replace the other regularity conditions.

Worked Example: Comparing Products of Nuisance Rates

Consider hypothetical rate descriptions for a propensity estimate and an outcome-regression estimate. In the first case, their errors are \(0.8n^{-0.22}\) and \(0.7n^{-0.30}\). The product is

$$ (0.8n^{-0.22})(0.7n^{-0.30})=0.56n^{-0.52}. $$

Because \(0.52>0.50\), this product is \(o(n^{-1/2})\): its ratio to \(n^{-1/2}\) is \(0.56n^{-0.02}\), which tends to zero. At \(n=10{,}000\), the product is approximately \(0.56(10{,}000^{-0.52})=0.00466\), while \(n^{-1/2}=0.0100\).

Now suppose both errors instead have rate \(n^{-0.20}\), with the same constants. Their product is \(0.56n^{-0.40}\), whose ratio to \(n^{-1/2}\) is \(0.56n^{0.10}\). That ratio grows rather than tends to zero. At \(n=10{,}000\), the product is approximately \(0.0141\), greater than \(0.0100\). These calculations compare rate behavior; a single finite-sample comparison does not establish that an asymptotic condition holds in an applied analysis.

The pitfall is to report that both nuisance models were estimated with flexible methods, or that each appeared accurate, and infer that root-\(n\) inference follows. The correction is to assess their joint error behavior under suitable conditions. If one nuisance is correctly specified, the usual double-robust consistency result may apply, but that alone does not certify every claimed rate, limiting distribution, or standard error.

A Short Justification: Why the Product Appears

The product-rate condition is connected to the bias remaining in the expected AIPW score when nuisance functions are estimated. Conditional on \(X\), take the expectation of the score using the true data distribution but fixed fitted functions \(\hat e,\hat m_1,\hat m_0\). Subtract the target conditional contrast \(m_{1,0}(X)-m_{0,0}(X)\). Directly collecting the treated and untreated terms gives

$$ \{\hat e(X)-e_0(X)\} \left\{ \frac{\hat m_1(X)-m_{1,0}(X)}{\hat e(X)} + \frac{\hat m_0(X)-m_{0,0}(X)}{1-\hat e(X)} \right\}. $$

This expression is zero if either the propensity equals its truth or both outcome regressions equal their truths. That cancellation explains the double-robust structure. When neither is exact, the remaining bias depends on products of nuisance errors, with denominators that can amplify errors near the boundaries. Under suitable bounds, taking expectations and applying standard product bounds leads to the product-rate requirement. This is a bias calculation, not a proof that all other sources of error or all inferential conditions are automatically controlled.

Error 3: Evaluating Training Predictions on Their Own Observations

Flexible nuisance learners can adapt closely to the observations used to fit them. If the AIPW score is then evaluated on those same observations, the fitted residuals and predictions may be unusually favorable in-sample. The resulting empirical score can behave differently from the score that would be obtained on new observations, complicating control of overfitting bias and the usual asymptotic argument.

Cross-fitting addresses this concern by splitting observations into folds, fitting nuisance models without a fold, and using those fitted models to predict for observations in the held-out fold. Each observation's score is computed from nuisance predictions made without using that observation to train the nuisance fits. Cross-fitting is not a remedy for positivity violations, inadequate nuisance rates, or incorrect data splitting for dependent observations.

Worked Example: How Reusing a Training Prediction Changes a Score

Suppose a hypothetical treated observation has \(Y=12\), \(\hat e(X)=0.5\), and \(\hat m_0(X)=8\). If an out-of-fold prediction for the treated outcome is \(\hat m_1(X)=10\), the score is

$$ \hat q_{\text{out}}=10-8+\frac{12-10}{0.5}=2+4=6. $$

Imagine instead that a highly adaptive in-sample fit memorizes this observation and returns \(\hat m_1(X)=12\). The score using that training prediction becomes

$$ \hat q_{\text{in}}=12-8+\frac{12-12}{0.5}=4+0=4. $$

The two scores differ by 2 outcome units because the residual correction disappears in the memorized fit. This one-observation calculation does not establish the direction or size of bias in a whole analysis. It shows why using training predictions can change the score in a way that is tied to overfitting, and why the inferential justification must match how predictions were generated. With cross-fitting, the prediction used for this observation would come from a model trained on other folds.

Error 4: Reporting a Standard Error for the Wrong Quantity

For an asymptotically linear AIPW estimator under its required conditions, a standard influence-function standard error is the empirical standard deviation of estimated influence-function contributions divided by \(\sqrt n\). Equivalently, it is the square root of their sample variance divided by \(n\). The contributions should correspond to the actual estimator and, when cross-fitting is used, to the cross-fitted predictions. A standard error calculated from the raw outcome variation alone generally does not represent the sampling variability of an AIPW estimate, which also depends on treatment assignment, weighting, outcome predictions, and covariates.

Worked Example: Computing an Influence-Function Standard Error

Suppose, for arithmetic illustration, the cross-fitted AIPW score contributions for four hypothetical independent observations are \(2,4,6,8\) outcome units. Their average is

$$ \hat\psi=\frac{2+4+6+8}{4}=5. $$

The estimated influence-function contributions are the scores centered at their average: \(-3,-1,1,3\). Their sample variance is

$$ \widehat{\operatorname{Var}}(\hat\phi) =\frac{(-3)^2+(-1)^2+1^2+3^2}{4-1} =\frac{20}{3}=6.667. $$

Thus the plug-in standard error is

$$ \widehat{SE}(\hat\psi)=\sqrt{\frac{6.667}{4}}=1.291 $$

outcome units. This is a calculation of the standard-error formula, not a claim that four observations suffice for reliable normal-approximation inference. The method requires the influence-function approximation and the relevant regularity conditions; it does not validate them merely by producing a number.

A common invalid shortcut is to use the standard deviation of \(Y\) divided by \(\sqrt n\) as though the estimate were a simple sample mean of outcomes. Another is to calculate influence-function values from predictions that do not correspond to the predictions used in the reported estimate. Both can produce numbers, but neither automatically estimates the sampling variation of the stated causal estimator. If observations are clustered or otherwise dependent, an independent-observation variance calculation may also be inappropriate; the variance estimator must reflect the sampling and dependence structure.

A Practical Diagnostic Sequence

Before interpreting a semiparametric estimate or confidence interval, separate the questions rather than treating “doubly robust” as a single pass-fail label.

1
Check identification and overlap.
State the causal assumptions and inspect whether both treatment levels are represented across relevant covariate patterns. Look for extreme fitted propensities and assess their consequences.
2
Check nuisance behavior jointly.
Consider the required nuisance-rate and moment conditions together, including whether small propensity denominators amplify estimation errors.
3
Check how predictions were produced.
For cross-fitted estimation, verify that each score uses predictions from fits that excluded that observation and that the reported estimate uses those same predictions.
4
Match the standard error to the estimator.
Use the appropriate estimated influence-function contributions and account for the sampling design or dependence when required.

Common Mistake

Wrong reasoning: “The AIPW estimate is doubly robust, so it remains reliable even with extreme propensities, slow nuisance fits, in-sample predictions, and a standard error based only on outcome variation.” Each part of this claim confuses a specific theoretical property with a different requirement. Double robustness does not repair nonidentification or weak overlap; it does not eliminate product-rate and regularity conditions for root-\(n\) inference; it does not make training predictions equivalent to out-of-fold predictions; and it does not justify a variance formula unrelated to the estimator's influence function.

Correction: Treat identification, nuisance estimation, prediction construction, and variance estimation as separate checks. A defensible conclusion is only as strong as the conditions supporting each step. Sensitivity to extreme propensities, nuisance fitting choices, or variance procedures should be reported as uncertainty about the analysis, not hidden by the phrase “doubly robust.”

Key takeaway: Double robustness is a useful protection against certain forms of nuisance-model misspecification under stated conditions. It is not protection against absent overlap, inadequate joint nuisance rates, uncontrolled reuse of training data, or an invalid standard error.

Check Your Understanding

Use the distinctions in this tutorial to assess each situation.

  1. Why can a correct propensity model still produce an unstable AIPW estimate when some fitted propensities are close to zero?
  2. If both nuisance errors are of order \(n^{-0.20}\), what is the order of their product, and how does it compare with \(n^{-1/2}\)?
  3. For an observation with \(A=1\), \(Y=9\), \(\hat e(X)=0.25\), \(\hat m_1(X)=7\), and \(\hat m_0(X)=5\), calculate its AIPW score contribution.
  4. What feature of cross-fitting distinguishes the prediction used for an observation from a prediction obtained by fitting on that observation?
  5. Why is the sample standard deviation of the raw outcome generally not enough to obtain an AIPW standard error?
Clinical Biostats