A framework to efficiently smooth L1 penalties for linear regression
摘要
Penalized linear regression operators that include an L1 penalty, such as the least absolute shrinkage and selection operator (Lasso), have become an important tool in statistical data analysis. Although the L1 penalty makes their objective function convex, it is not differentiable everywhere, motivating the development of proximal gradient algorithms. In this article, we introduce a unified framework to compute closed-form smooth surrogates of a whole class of L1-penalized regression operators using Nesterov smoothing in a fixed-parameter setting (the problem size n and number of parameters p are fixed). Our contributions are: (1) Our smooth surrogates preserve the convexity of the original (unsmoothed) objective functions, are uniformly close to them, and have closed-form derivatives everywhere for efficient minimization via gradient descent; (2) We prove that the estimates obtained with the smooth surrogates can be made arbitrarily close to the ones of the original (unsmoothed) objective functions and provide explicitly computable a priori error bounds on the accuracy of our estimates; (3) We propose an iterative algorithm to progressively smooth the L1 penalty that increases accuracy and is virtually free of tuning parameters. Since the resulting estimates are typically dense, we explore the use of thresholding to enforce sparsity. We compare our framework to the current gold standards on simulated and real data, demonstrating that our smooth surrogates provide equal or higher accuracy than the gold standards while keeping the aforementioned theoretical guarantees and having roughly the same asymptotic runtime scaling.