Phase transition and higher order analysis of Lq regularization under dependence

Abstract

We study the problem of estimating a k-sparse signal β0∈ Rp from a set of noisy observations y∈ Rn under the model y= Xβ+ w, where X∈ Rn× p is the measurement matrix the row of which is drawn from distribution N(0,). We consider the class of Lq-regularized least squares (LQLS) given by the formulation β(λ,q)=argminβ∈ Rp12\| y- Xβ\|22+λ\|β\|qq, where \|·\|q (0 q 2) denotes the Lq-norm. In the setting p,n,k→∞ with fixed k/p=ε and n/p=δ, we derive the asymptotic risk of β(λ,q) for arbitrary covariance matrix which generalizes the existing results for standard Gaussian design, i.e. Xiji.i.dN(0,1). We perform a higher-order analysis for LQLS in the small-error regime in which the first dominant term can be used to determine the phase transition behavior of LQLS. Our results show that the first dominant term does not depend on the covariance structure of in the cases 0 q< 1 and 1< q 2 which indicates that the correlations among predictors only affect the phase transition curve in the case q=1 a.k.a. LASSO. To study the influence of the covariance structure of on the performance of LQLS in the cases 0 q< 1 and 1<q 2, we derive the explicit formulas for the second dominant term in the expansion of the asymptotic risk in terms of small error. Extensive computational experiments confirm that our analytical predictions are consistent with numerical results.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…