1 Introduction
The “double descent” risk curve was proposed by *belkin2018reconciling to qualitatively describe the outofsample prediction performance of several variablyparameterized machine learning models. This risk curve reconciles the classical biasvariance tradeoff with the behavior of predictive models that interpolate training data, as observed for several model families (including neural networks) in a wide variety of applications bos1998dynamics,advani2017high,spigler2018jamming,belkin2018reconciling. In these studies, a predictive model with
parameters is fit to a training sample of size , and the test risk (i.e., outofsample error) is examined as a function of . When is below the sample size , the test risk is governed by the usual biasvariance decomposition. As is increased towards , the training risk (i.e., insample error) is driven to zero, but the test risk shoots up towards infinity. The classical biasvariance analysis identifies a “sweet spot” value of at which the bias and variance are balanced to achieve low test risk. However, as grows beyond , the test risk again decreases, provided that the model is fit using a suitable inductive bias (e.g., least norm solution). In many (but not all) cases, the limiting risk as is lower than what is achieved at the “sweet spot” value of .In this article, we study the key aspects of the “double descent” risk curve for the least squares/least norm predictor in two simple random features models. The first is a Gaussian model, which was studied by breiman1983many in the regime. The second is a Fourier series model for functions on the circle. In both cases, we prove that the risk is infinite around , and decreases to towards its minimum as increases beyond . Our results provide a precise mathematical analysis of the mechanism described by belkin2018reconciling. The transition from under to overparametrized regimes was also analyzed by *spigler2018jamming by drawing a connection to the physical phenomenon of “jamming” in particle systems.
We note that in both of the models, the features are selected randomly, which makes them useful for studying scenarios where features are plentiful but individually too “weak” to be selected in an informed manner. Such scenarios are commonplace in machine learning practice, but they should be contrasted with scenarios where features are carefully designed or curated, as is often the case in scientific applications. For comparison, we give an example of “prescient” feature selection, where the
most useful features are included in the model. In this case, the optimal test risk is achieved at some , which is consistent with the classical analysis of breiman1983many.2 Gaussian model
We consider a regression problem where the response is equal to a linear function of realvalued variables plus noise :
The learner observes iid copies of , but fits a linear model to the data only using a subset of variables.
Let be the design matrix, and let
be the vector of responses. For a subset
and a dimensional vector , we use to denote its dimensional subvector of entries from ; we also use to denote the design matrix with variables from . For , we denote its complement by . Finally, denotes the Euclidean norm.The learner fits regression coefficients with
Above, the symbol denotes the MoorePenrose pseudoinverse. In other words, the learner uses the solution to the normal equations of least norm for and forces to allzeros.
In the remainder of this section, we analyze the risk of in the case where the distribution of is the standard normal in , and then specialize the risk under particular selection models for . The Gaussian model was also studied by breiman1983many, although their analysis is restricted to the case where the number of variables used is always at most ; our analysis will also consider the regime. We show that the “interpolating” regime is preferred to the “classical” regime in our model.
We focus on the noisefree setting in Section 2.1 and Section 2.2; we consider noise in Section 2.3.
2.1 Risk analysis
In this section, we derive an expression for the risk of for an arbitrary choice of features .
We assume (i.e., the noisefree setting), so . Recall that we also assume
follows a standard normal distribution in
; since is isotropic (i.e., zero mean and identity covariance), the mean squared prediction error of any can be written asSince , it follows that the risk of is
Classical regime.
The risk of was computed by breiman1983many in the regime where :
Interpolating regime.
We consider the regime where . Recall that the pseudoinverse of can be written as . Thus,
On the right hand side, the first term is the orthogonal projection of onto the null space of , while the second term is a vector in the row space of . By the Pythagorean theorem, the squared norm of their sum is equal to the sum of their squared norms, so
We analyze the expected values of these two terms by exploiting properties of the standard normal distribution.
 First term.

Note that is the orthogonal projection matrix for the row space of . So, by the Pythagorean theorem, we have
By rotational symmetry of the standard normal distribution, it follows that
Therefore
 Second term.

We use the “trace trick” to write
where the second equality holds almost surely because is almost surely invertible. Since and are uncorrelated, it follows that
The distribution of is normal with mean zero and covariance , so
The distribution of is inverseWishart with identity scale matrix and degreesoffreedom, so
Combining the last two displayed equations yields a simple expression for .
Thus, we obtain expressions for and hence also .
We summarize the risk of in the following theorem.
Theorem 1.
Assume the distribution of is the standard normal in and for some . Pick any and of cardinality . The risk of , where and , is
2.2 Feature selection model
We now study the risk of under a random selection model for . Again, we assume .
Let be a uniformly random subset of of cardinality , so
We analyze the risk of , taking expectation with respect to the random choice of . First, consider . By Theorem 1, the risk of is
which increases with as long as . Now consider . By Theorem 1, the risk of is
which decreases with as long as .
Thus, we observe that the risk first increases with up to the “interpolation threshold” (), after which the risk decreases with . Moreover, the risk is smallest at . This is the “double descent” risk curve observed by belkin2018reconciling where the first “descent” is degenerate (i.e., the “sweet spot” that balances bias and variance is at ). See Figure 1 for an illustration. For a scenario where the first “descent” is nondegenerate, see Appendix A.
It is worth pointing out that the behavior under the random selection model of can be very different from that under a deterministic model of . Consider including variables in by decreasing order of —a kind of “prescient” selection model studied by breiman1983many. For simplicity, assume and , so
The behavior of the risk as a function of , illustrated in Figure 2, reveals a striking difference between the random selection model and the “prescient” selection model.
2.3 Noise
We now briefly discuss the effect of additive noise in the responses. The following is a straightforward generalization of Theorem 1.
Theorem 2.
Assume the distribution of is the standard normal in ,
is a standard normal random variable independent of
, and for some and . Pick any and of cardinality . The risk of , where and , isA similar analysis with a random selection model for can be obtained from Theorem 2.
3 Fourier series model
Let denote the
discrete Fourier transform matrix: its
th entry iswhere is a primitive root of unity. Let for some . Consider the following observation model:

and are independent and uniformly random subsets of of cardinalities and , respectively.

We observe the design matrix and dimensional vector of responses . Here, is the submatrix of with rows from and columns from , and is the subvector of of entries from .
The learner fits regression coefficients with
This can be regarded as a onedimensional version of the random Fourier features model studied by rahimi2008random for functions defined on the unit circle.
One important property of the discrete Fourier transform matrix that we use is that the matrix has rank for any . This is a consequence of the fact that is Vandermonde. Thus, for , we have
In the remainder of this section, we analyze the risk of under a random model for , where
(which implies ). The random choice of is independent of and . Considering the risk under this random model for is a form of averagecase analysis. For simplicity, we only consider the regime where , as it suffices to reveal some key aspects of the risk of .
Following the arguments from Section 2.1, we have
Now we take (conditional) expectations with respect to , given and :
(1) 
Since has rank , the first trace expression is equal to
For the second trace expression, we use the explicit formula for and the fact that to obtain
where
are the eigenvalues of
. Therefore, from Equation 1, we haveA precise characterization of is difficult to obtain. Under a slightly different model, in which membership in (respectively, ) is determined by independent Bernoulli variables with mean (respectively, ), we can use asymptotic arguments to characterize the empirical eigenvalue distribution for .^{1}^{1}1We can also derive essentially the same formula for the risk under this Bernoulli model, although the derivation is somewhat more cumbersome since we do not always have even if . Hence, we have opted to present the derivation only for the simpler model. Assuming the asymptotic equivalence of these random models for and , we find that the quantity approaches
as , where and are held fixed and farrell2011limiting. So, in this limit, we have
This quantity diverges to as , and decreases as . This is the same behavior as in the Gaussian model from Section 2 with random selection; we depict it empirically in Figure 3.
4 Discussion
Our analyses show that when features are chosen in an uninformed manner, it may be optimal to choose as many as possible—even more than the number of data—rather than limit the number to that which balances bias and variance. This choice is conceptually and algorithmically simple and avoids the need for precise control of regularization parameters. It is reminiscent of the practice in machine learning applications like image and speech recognition, where signal processingbased features are individually weak but in great abundance, and models that use all of the features are highly successful. This stands in contrast to scenarios where informed selection of features is possible; for example, in many science and medical applications, features are handcrafted and purposefully selected. As illustrated by the “prescient” selection model, choosing the number of features to balance bias and variance can be better than incurring the costs that come with using all of the features. The best practices for model and feature selection thus crucially depend on which regime the application is operating under.
Appendix A Nondegenerate double descent in a Fourier series model
To observe the “double descent” risk curve where the first “descent” is nondegenerate, we consider a model in which the distribution of the feature vector is nonisotropic; instead, has a diagonal covariance matrix with decaying eigenvalues (e.g., ). In such a scenario, it is natural to select features in decreasing order of the eigenvalues, as is done in principal components regression. Like random feature selection, this form of feature selection is also generally uninformed by the responses.
Formally, we let for some and positive sequence , and use the same observation model as in Section 3, except that is deterministically set to . We fit in the same manner as in Section 3. We are interested in the risk
which can be regarded as the mean squared error when is drawn uniformly at random from the rows of the design matrix .
We carried out the same simulation as from Figure 3 under the modified model, with and . We chose uniformly at random (once) from the unit sphere in for . Then, for each , we computed from independent random choices of (with ), and plotted the average value of . The plot is shown below. (The vertical axis is truncated for clarity; the curve peaks around with value on the order of .)
The plot shows the usual “U”shaped curve arising from the biasvariance tradeoff when , and a second “descent” towards the overall minimum for . This risk curve is qualitatively the same as those observed by belkin2018reconciling for neural networks and other predictive models.