Skip to content

Benign overfitting

When exact fitting of noisy training data can coexist with vanishing excess prediction risk—and why interpolation and minimum norm alone do not guarantee it.

Benign overfitting is the phenomenon in which a predictor fits noisy training observations exactly yet predicts nearly as well as the underlying signal rule. In an asymptotic formulation, its excess prediction risk tends to zero as sample size grows. For a quant researcher, zero training error therefore leaves a question unanswered: what happens on fresh data? Bartlett and colleagues characterise this possibility for minimum-norm linear regression.

“Benign” refers to excess risk

Take the correctly specified model y=x⊤β∗+εy=x^{\top}\beta^{*}+\varepsilon, with E[ε∣x]=0\mathbb{E}[\varepsilon\mid x]=0, constant conditional noise variance σ2>0\sigma^{2}>0, and finite second moments. For a fitted coefficient vector β^\hat\beta, evaluated on a fresh observation independent of training, expand the square:

E[(y−x⊤β^)2∣β^]=σ2+E(β^),E(β^)=(β^−β∗)⊤Σ(β^−β∗),Σ=E[xx⊤].\begin{aligned} \mathbb{E}[(y-x^{\top}\hat\beta)^2\mid\hat\beta] &=\sigma^{2}+\mathcal{E}(\hat\beta),\\ \mathcal{E}(\hat\beta) &=(\hat\beta-\beta^{*})^{\top}\Sigma(\hat\beta-\beta^{*}), \qquad \Sigma=\mathbb{E}[xx^{\top}]. \end{aligned}

The cross-term vanishes by the conditional mean-zero assumption. With centred features, Σ\Sigma is the covariance matrix. Benign overfitting asks for E(β^)→0\mathcal{E}(\hat\beta)\to0, not total test MSE tending to zero: irreducible noise remains. Nor does small prediction risk require every coefficient error to be small; the covariance weights errors by how much those directions vary on new observations.

The bias–variance derivation develops the fresh-noise distinction in detail (full lesson requires membership; public preview available).

A perfect fit that fails

Here is a deliberately finite, noiseless counterexample to the claim that interpolation plus minimum norm guarantees good prediction. Use one training row X=(1,0)X=(1,0) and response y=1y=1. Every vector

β=(1,t),t∈R,\beta=(1,t),\qquad t\in\mathbb{R},

fits exactly because Xβ=1X\beta=1. Minimising the squared Euclidean coefficient norm 1+t21+t^{2} picks β^=(1,0)\hat\beta=(1,0).

Suppose the true coefficient vector is β∗=(1,2)\beta^{*}=(1,2). It agrees with the training observation. At the test row xtest=(0,1)x_{\mathrm{test}}=(0,1), however,

y^=0,y∗=2,(y^−y∗)2=4.\hat y=0,\qquad y^{*}=2,\qquad (\hat y-y^{*})^{2}=4.
Constructed noiseless counterexample: training row (1, 0) with response 1 permits every coefficient vector (1, t). Minimum Euclidean norm chooses (1, 0), while the true vector is (1, 2). At test row (0, 1), prediction 0 misses truth 2 and gives squared error 4.
Constructed noiseless counterexample: training row (1, 0) with response 1 permits every coefficient vector (1, t). Minimum Euclidean norm chooses (1, 0), while the true vector is (1, 2). At test row (0, 1), prediction 0 misses truth 2 and gives squared error 4.
An original counterexample, not a benign-overfitting experiment. The blue line contains all exact training fits; minimum norm chooses a point that misses the unobserved signal direction.

The sample contains no information about the second coefficient. A norm preference resolves the ambiguity, but it cannot establish that zero is the true value. “Minimum norm” also refers to a specified feature representation and scaling; changing them changes the selection rule. This single test-point error is not a population-risk estimate or an asymptotic result.

What can make interpolation harmless?

In the paper’s linear setting, a sufficiently broad tail of small covariance eigenvalues can distribute fitted label noise across directions that contribute little to prediction risk. Signal distortion must also vanish: the noise-fitting condition alone is insufficient. The theorem uses independent training observations, distributional tail assumptions and rank conditions—not parameter count alone.

One concrete sufficient construction follows from its Theorem 4. For each nn, use pn=n2+1p_n=n^{2}+1 centred Gaussian features with covariance eigenvalues

1,n−3/2,…,n−3/2⏟n2 entries.1,\quad \underbrace{n^{-3/2},\ldots,n^{-3/2}}_{n^{2}\text{ entries}}.

Assume i.i.d. training pairs, independent Gaussian label noise with fixed positive variance, and true coefficient vectors whose Euclidean norms stay bounded. This is a sequence of changing feature spaces and covariances, not repeated fitting of one fixed dataset.

To check the construction, define the paper’s tail ranks for decreasing eigenvalues λj\lambda_j:

rk=∑j>kλjλk+1,Rk=(∑j>kλj)2∑j>kλj2.r_k=\frac{\sum_{j>k}\lambda_j}{\lambda_{k+1}},\qquad R_k=\frac{(\sum_{j>k}\lambda_j)^2}{\sum_{j>k}\lambda_j^2}.
Direct summation gives r0=1+nr_0=1+\sqrt n and r1=R1=n2r_1=R_1=n^{2}. For the theorem’s fixed constant b>1b>1, the first index with rk≥bnr_k\ge bn is k∗=1k^{*}=1 eventually, not necessarily at a particular finite sample size. Thus
r0n=1+nn→0,k∗n=1n→0,nRk∗=1n→0.\frac{r_0}{n}=\frac{1+\sqrt n}{n}\to0,\qquad \frac{k^{*}}{n}=\frac1n\to0,\qquad \frac{n}{R_{k^{*}}}=\frac1n\to0.

Under these assumptions, the theorem yields vanishing excess risk in probability. Here R1R_1 describes the small-eigenvalue tail; it is not R0R_0, which grows only on the order of nn. Neither the construction nor the counterexample proves success at a chosen finite nn.

Reading a financial modelling result

A low-variance feature direction is not automatically irrelevant: its coefficient and alignment with the signal matter, as the counterexample shows. Financial observations can also be dependent, and their feature covariance can change over time. Those features of a research problem need separate treatment before applying an i.i.d. theorem.

Out-of-sample evaluation still has to assess the fitted rule and the selection process. A double-descent curve describes how error changes across model sizes or training stages; benign overfitting describes the quality of interpolating prediction. Observing a second descent does not by itself establish vanishing excess risk.

Sources

Last updated
  • generalisation
  • interpolation
  • minimum norm
  • covariance spectrum
← Back to the glossaryBack to top ↑