A proximal-gradient homotopy method for the sparse least-squares problem

From MaRDI portal



Abstract: We consider solving the ell1-regularized least-squares (ell1-LS) problem in the context of sparse recovery, for applications such as compressed sensing. The standard proximal gradient method, also known as iterative soft-thresholding when applied to this problem, has low computational cost per iteration but a rather slow convergence rate. Nevertheless, when the solution is sparse, it often exhibits fast linear convergence in the final stage. We exploit the local linear convergence using a homotopy continuation strategy, i.e., we solve the ell1-LS problem for a sequence of decreasing values of the regularization parameter, and use an approximate solution at the end of each stage to warm start the next stage. Although similar strategies have been studied in the literature, there have been no theoretical analysis of their global iteration complexity. This paper shows that under suitable assumptions for sparse recovery, the proposed homotopy strategy ensures that all iterates along the homotopy solution path are sparse. Therefore the objective function is effectively strongly convex along the solution path, and geometric convergence at each stage can be established. As a result, the overall iteration complexity of our method is O(log(1/epsilon)) for finding an epsilon-optimal solution, which can be interpreted as global geometric rate of convergence. We also present empirical results to support our theoretical analysis.


The authors deal with the \(l_1\)-regularized least squares problem, proposing and analyzing an efficient numerical method for solving it. An approximate homotopy continuation method is used, solving the \(l_1\)-regularized least squares problem with a large parameter \(\lambda\) first and gradually decreasing \(\lambda\) until the targeted regularization is obtained. For each \(\lambda\), the problem is solved by a proximal gradient method up to the desired precision, getting a solution, which is the initial point in the next iteration corresponding to the next value \(\lambda\). The resulting technique is called proximal gradient homotopy method in this paper. The method has a provable low iteration complexity in specific algorithmic conditions. The overall iteration complexity and the overall computational cost are theoretically evaluated. Empirical results, which support the authors' theoretical analysis, are finally presented.




Cited in
(45)


Describes a project that uses

Uses Software






This page was built for publication: A proximal-gradient homotopy method for the sparse least-squares problem

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2848186)