When Can Exact Fit Still Predict Well?
An interpolating estimator fits every observed response, including its noise. That fact alone says little about its prediction error: the fitted noise may be spread across many directions that matter little for new observations, or it may be amplified along directions with substantial population variance. For minimum-norm linear interpolation, this distinction can be made precise by separating two sources of error: the part of the true signal the design cannot identify, and the noise the estimator uses to achieve exact fit.
The result below is a sufficient condition, not a claim that interpolation is always safe. It makes the roles of the covariance spectrum and signal alignment explicit. One important qualification is that population covariance conditions do not by themselves guarantee a well-conditioned sample: the stated sample-design conditions must also hold.
Model and Spectral Quantities
Consider a sequence of linear regression problems with \(n\) observations and \(p\) features, where \(p>n\). Write the design matrix as \(X\in\mathbb{R}^{n\times p}\), with rows \(x_i^T\), and suppose the responses satisfy \(y=X\beta^*+\varepsilon\). The independent test feature \(x\) has mean zero and covariance \(\Sigma=E[xx^T]\). For the risk calculation, assume the training errors are independent of \(X\), have mean zero, and have covariance \(\sigma^2 I_n\). No Gaussian assumption is needed for the conditional calculation below.
Diagonalize the population covariance as \(\Sigma=U\operatorname{diag}(\lambda_1,\ldots,\lambda_p)U^T\), where \(\lambda_1\ge\cdots\ge\lambda_p\ge0\). Choose an integer \(k<n\), and call the coordinates after the first \(k\) the covariance tail. Its mass, squared spectral mass, and effective rank are
These quantities are defined when \(q_k>0\). Tail mass measures total test-feature variance in those directions; \(q_k\) records how concentrated that variance is in large eigenvalues. Effective rank is large when the mass is distributed across many directions rather than dominated by a few. Since \(\lambda_j\le\lambda_{k+1}\) in the tail, \(q_k\le\lambda_{k+1}\tau_k\), and therefore \(r_k\ge \tau_k/\lambda_{k+1}\).
A Conditional Benign-Overfitting Result
Let \(X_T\) denote the \(n\times(p-k)\) design matrix formed from the tail coordinates in the eigenbasis of \(\Sigma\), and let \(P=X^\dagger X\), the orthogonal projection in parameter space onto the row space of \(X\). Two conditions describe the sample geometry. First, assume the tail alone supplies enough independent variation to interpolate stably. Second, assume the observed design does not have an unusually large covariance-weighted squared size.
In particular, along a sequence of problems, if \(c,C\) remain bounded away from zero and infinity, the signal-alignment term tends to zero, and \(r_k/n\to\infty\), then the conditional prediction risk tends to zero. This conclusion follows because \(nq_k/\tau_k^2=n/r_k\).
Proof. Since \(X^\dagger X=P\), the coefficient error is \(\hat\beta-\beta^*=-(I-P)\beta^*+X^\dagger\varepsilon\). The first term is fixed conditional on \(X\), while the second has conditional mean zero and covariance \(\sigma^2X^\dagger(X^\dagger)^T\). Thus their cross term has conditional expectation zero. For any fixed vector \(v\), \(E[(x^Tv)^2]=v^T\Sigma v=\|\Sigma^{1/2}v\|_2^2\). Applying this identity to the coefficient error gives the displayed equality, including the trace expression for the variance contribution.
For the bound, put \(A=XX^T\). The tail condition implies \(A\succeq X_TX_T^T\succeq c\tau_k I_n\), because the head contribution to \(XX^T\) is positive semidefinite. Hence \(\|A^{-1}\|_{\mathrm{op}}\le(c\tau_k)^{-1}\). Using \(X^\dagger=X^TA^{-1}\), cyclicity of trace, and positive semidefiniteness,
The inequality uses \(\operatorname{tr}(A^{-1}MA^{-1})=\operatorname{tr}(MA^{-2})\le\|A^{-1}\|_{\mathrm{op}}^2\operatorname{tr}(M)\) for positive semidefinite \(M\). This proves the finite-sample bound. Finally \(q_k\le\lambda_{k+1}\tau_k\) also gives \(nq_k/\tau_k^2\le n/r_k\), so the stated effective-rank condition makes the noise bound vanish. \(\square\)
Worked Examples
Worked Example: Isotropic Features
Suppose \(\Sigma=I_p\), \(p>n\), and choose \(k=0\). Then \(\tau_0=p\), \(q_0=p\), and \(r_0=p\). If the realized design satisfies \(\lambda_{\min}(XX^T)\ge cp\) and \(\operatorname{tr}(XX^T)\le Cnp\), the theorem gives a noise contribution at most
For the zero signal \(\beta^*=0\), the alignment term is exactly zero. Thus, when \(p/n\to\infty\) and the two sample conditions hold with fixed constants, the conditional prediction risk vanishes even though the estimator interpolates all \(n\) noisy responses. The conclusion depends on the sample conditions as well as on \(p>n\); overparameterization alone is not the guarantee.
Worked Example: Many Small Tail Eigenvalues
Suppose the tail has \(m\) eigenvalues, each equal to \(1/m\). Then \(\tau_k=m(1/m)=1\), \(q_k=m(1/m^2)=1/m\), and \(r_k=1/(1/m)=m\). If \(m\gg n\), the theorem’s noise factor is \(nq_k/\tau_k^2=n/m\), which tends to zero. A tail with fixed total variance can therefore suppress noise amplification when that variance is spread over sufficiently many directions. This is a covariance calculation; the theorem still requires the realized tail design to have \(\lambda_{\min}(X_TX_T^T)\ge c\) and the stated trace bound.
Worked Example: Why Signal Alignment Matters
Take \(\Sigma=\operatorname{diag}(10,1,0.01)\). If the row space of a design omits the first coordinate, then \(I-P\) retains \(e_1\). For \(\beta^*=5e_1\), the squared alignment term is
If instead the row space omits the third coordinate and \(\beta^*=5e_3\), the same calculation gives \(5^2(0.01)=0.25\). In both cases the coefficient component is unobserved, but its prediction cost differs sharply because the corresponding population variances differ. The alignment condition concerns the unobserved signal in the \(\Sigma\)-weighted prediction metric, not merely the Euclidean size of the signal.
What Each Condition Rules Out
- Tail sample conditioning. The lower eigenvalue bound ensures that the tail features alone span the observation space with enough strength. Without it, \(XX^T\) can be nearly singular and the exact-fit coefficients can become very large.
- Tail effective rank. The requirement \(r_k/n\to\infty\) is a sufficient way to make the noise term small. It says that the tail’s variance is not concentrated in too few directions relative to the sample size.
- Signal alignment. The term \(\|\Sigma^{1/2}(I-P)\beta^*\|_2^2\) measures prediction error from the part of the signal outside the observed row space. It must be small; a favorable tail spectrum cannot remove this bias by itself.
- Trace control. The bound on \(\operatorname{tr}(X\Sigma X^T)\) rules out an atypically large covariance-weighted design. It is a sample-level condition, not a consequence of the effective-rank definition alone.
The theorem is deliberately conditional: it separates a transparent risk calculation from the separate problem of proving that particular random designs satisfy the sample conditions with high probability. It also describes prediction risk, not training error. Training error is zero by construction, regardless of whether prediction risk is small.
Common Pitfall
A tempting but incorrect argument is: “The model has more parameters than observations, so the minimum-norm interpolator can distribute the noise over many coordinates; therefore it must generalize.” The number of parameters does not say whether the tail covariance has large effective rank, whether the sample tail is well-conditioned, or whether the unobserved signal lies in low-variance directions. The correction is to check all three: tail sample conditioning, a small noise factor \(nq_k/\tau_k^2\), and small covariance-weighted signal misalignment. Benign overfitting is a joint spectral, geometric, and signal condition—not a consequence of exact interpolation alone.
Check Your Understanding
Use the definitions and theorem to answer the following.
- For tail eigenvalues \(1,1,1,1\), calculate \(\tau_k\), \(q_k\), and \(r_k\). How do they change if the same total mass is spread evenly over 100 eigenvalues?
- Why does a large tail effective rank not, by itself, guarantee a small conditional prediction risk?
- In the theorem, what does the projection \(P\) represent, and what prediction error is measured by \(\|\Sigma^{1/2}(I-P)\beta^*\|_2^2\)?
- If \(nq_k/\tau_k^2\) tends to zero but the alignment term stays bounded away from zero, what prevents the theorem from implying vanishing risk?
- Explain why zero training error is compatible with nonzero prediction risk.