Best subset selection, persistence in high-dimensional statistical learning and optimization under l1 constraint

Abstract

Let (Y, X₁, …, X_m) be a random vector. It is desired to predict Y based on (X₁, …, X_m). Examples of prediction methods are regression, classification using logistic regression or separating hyperplanes, and so on.

We consider the problem of best subset selection, and study it in the context m=n^α, α>1, where n is the number of observations. We investigate procedures that are based on empirical risk minimization. It is shown, that in common cases, we should aim to find the best subset among those of size which is of order o(n / log(n)). It is also shown, that in some “asymptotic sense,” when assuming a certain sparsity condition, there is no loss in letting m be much larger than n, for example, m=n^α, α>1. This is in comparison to starting with the “best” subset of size smaller than n and regardless of the value of α.

We then study conditions under which empirical risk minimization subject to l₁ constraint yields nearly the best subset. These results extend some recent results obtained by Greenshtein and Ritov.

Finally we present a high-dimensional simulation study of a “boosting type” classification procedure.