Sparse High-Dimensional Linear Regression. Algorithmic Barriers and a Local Search Algorithm
Abstract
We consider a sparse high dimensional regression model where the goal is to recover a k-sparse unknown vector β* from n noisy linear observations of the form Y=Xβ*+W ∈ Rn where X ∈ Rn × p has iid N(0,1) entries and W ∈ Rn has iid N(0,σ2) entries. Under certain assumptions on the parameters, an intriguing assymptotic gap appears between the minimum value of n, call it n*, for which the recovery is information theoretically possible, and the minimum value of n, call it nalg, for which an efficient algorithm is known to provably recover β*. In gamarnikzadik it was conjectured that the gap is not artificial, in the sense that for sample sizes n ∈ [n*,nalg] the problem is algorithmically hard. We support this conjecture in two ways. Firstly, we show that the optimal solution of the LASSO provably fails to 2-stably recover the unknown vector β* when n ∈ [n*,c nalg], for some sufficiently small constant c>0. Secondly, we establish that nalg, up to a multiplicative constant factor, is a phase transition point for the appearance of a certain Overlap Gap Property (OGP) over the space of k-sparse vectors. The presence of such an Overlap Gap Property phase transition, which originates in statistical physics, is known to provide evidence of an algorithmic hardness. Finally we show that if n>C nalg for some large enough constant C>0, a very simple algorithm based on a local search improvement rule is able both to 2-stably recover the unknown vector β* and to infer correctly its support, adding it to the list of provably successful algorithms for the high dimensional linear regression problem.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.