Module 5 — Modern Causal Analysis

Double Machine Learning

Estimate low-dimensional causal parameters with orthogonal scores, cross-fitting and flexible nuisance models

Double Machine Learning

Prediction enters as a nuisance task

Consider the partially linear model

Y=θD+g(X)+ε,D=m(X)+v,Y=\theta D+g(X)+\varepsilon,\qquad D=m(X)+v,

with E[εX,D]=0E[\varepsilon\mid X,D]=0 and E[vX]=0E[v\mid X]=0. The target is the scalar θ\theta; gg and mm may be complex nuisance functions.

Residualise both sides:

Y~=YE^[YX],D~=DE^[DX],qquadθ^=iD~iY~iiD~i2.\tilde Y=Y-\hat E[Y\mid X],\qquad \tilde D=D-\hat E[D\mid X],qquad \hat\theta=\frac{\sum_i\tilde D_i\tilde Y_i}{\sum_i\tilde D_i^2}.

This resembles FWL, but flexible learners create overfitting bias if a record helps train the function used to residualise itself.

Cross-fitting gives observations separate jobs

For two folds:

  1. fit both nuisance functions on fold A;
  2. predict and residualise fold B;
  3. reverse A and B;
  4. combine out-of-fold residuals;
  5. estimate θ\theta and design-aligned uncertainty.

All tuning must occur inside the training data for each fold. Preprocessing on the full sample leaks information.

Py

Cross-fit nonlinear nuisance functions

Idle

The polynomial basis is only a transparent learner for teaching. In practice, choose learners through out-of-fold performance and substantive structure, then inspect whether the causal result is stable across reasonable learner libraries and splits.

Orthogonality buys local robustness, not immunity

The score

ψ(W;θ,η)=[Dm(X)]{Y(X)θ[Dm(X)]}\psi(W;\theta,\eta)= [D-m(X)]\{Y-\ell(X)-\theta[D-m(X)]\}

is constructed so small first-order nuisance errors have limited effect on θ\theta. Chernozhukov et al. (2018) formalise this orthogonality and cross-fitting framework.

It still needs:

  • an identifying model such as conditional exchangeability;
  • overlap and residual treatment variation;
  • suitable nuisance convergence and moments;
  • clustering/splitting that respects grouped or time-dependent data;
  • a low-dimensional target defined before model search.

Practical audit

Report the learner library, tuning process, fold construction, random seed, overlap, out-of-fold nuisance performance, repeated-split stability and final score/inference formula. The DoubleML documentation is useful implementation guidance, not a replacement for this record.

Quick check

A powerful learner predicts scholarship receipt almost perfectly. Is this ideal for DML?

Answer
Not necessarily. Near-perfect treatment prediction can reveal poor overlap: residual treatment variation becomes tiny and the effect is weakly supported. Inspect propensity distributions and the denominator D~i2\sum\tilde D_i^2.

Next: Heterogeneity and Policy Learning

Copyright © 2026