Harmful Overfitting in Sobolev Spaces
Kedar Karhadkar, Alexander Sietsema, Deanna Needell, Guido Montufar
TL;DR
This work analyzes interpolation with noisy data in Sobolev spaces $W^{k,p}(\mathbb{R}^d)$ in fixed dimension. It proves that approximately norm-minimizing interpolants (ANM) cannot yield benign overfitting: under mild noise and regularity assumptions, any interpolant with Sobolev norm within a constant factor of the minimum incurs a constant excess risk, independent of the sample size $n$, for $kp>d$ and $k\in(d/p, 1.5d/p)$. The results extend prior $p=2$ Hilbert-space analyses to general $p\in[1,\infty)$ by constructing explicit bump-function interpolants, locating a noisy, well-separated subset, and applying Sobolev/Morrey-type inequalities to relate local oscillation to global generalization. In the Gaussian heteroskedastic noise setting with squared loss, the paper provides a concrete $L^2(\mu)$ lower bound on the error of ANM interpolants. Overall, the findings show that norm minimization combined with smoothness does not guarantee benign interpolation in fixed-dimensional Sobolev spaces, highlighting the need for alternative inductive biases or undertraining in such regimes.
Abstract
Motivated by recent work on benign overfitting in overparameterized machine learning, we study the generalization behavior of functions in Sobolev spaces $W^{k, p}(\mathbb{R}^d)$ that perfectly fit a noisy training data set. Under assumptions of label noise and sufficient regularity in the data distribution, we show that approximately norm-minimizing interpolators, which are canonical solutions selected by smoothness bias, exhibit harmful overfitting: even as the training sample size $n \to \infty$, the generalization error remains bounded below by a positive constant with high probability. Our results hold for arbitrary values of $p \in [1, \infty)$, in contrast to prior results studying the Hilbert space case ($p = 2$) using kernel methods. Our proof uses a geometric argument which identifies harmful neighborhoods of the training data using Sobolev inequalities.
