2021 article
Local Convergence of an AMP Variant to the LASSO Solution in Finite Dimensions
2021 IEEE INTERNATIONAL SYMPOSIUM ON INFORMATION THEORY (ISIT), pp. 2256–2261.
TL;DR:
It is shown that whenever the AMP variant converges, it converges to the LASSO solution for arbitrary finite dimensional regression matrices and that the original AMP, which is a special case of the proposed AMP variants, is locally stable around the LassO solution.
(via Semantic Scholar)
open access via arxiv.org (repository)
Source: Web Of Science
Added: October 26, 2021
A common sparse linear regression formulation is $\ell_{1}$ regularized least squares, which is also known as least absolute shrinkage and selection operator (LASSO). Approximate message passing (AMP) has been proved to asymptotically achieve the LASSO solution when the regression matrix has independent and identically distributed (i.i.d.) Gaussian entries in the sense that the averaged per-coordinate $\ell_{2}$ distance between the AMP iterates and LASSO solution vanishes as the signal dimension goes to infinity before the iteration number. However, in finite dimensional settings, characterization of AMP iterates in the limit of large iteration number has not been established. In this work, we propose an AMP variant by including a parameter that depends on the largest singular value of the regression matrix. The proposed algorithm can also be considered as a primal dual hybrid gradient algorithm with adaptive stepsizes. We show that whenever the AMP variant converges, it converges to the LASSO solution for arbitrary finite dimensional regression matrices. Moreover, we show that our AMP variant is locally stable around the LASSO solution under the condition that the LASSO solution is unique and that the regression matrix is drawn from a continuous distribution. Our local stability result implies that when the regression matrix is large and has i.i.d. random entries, the original AMP, which is a special case of the proposed AMP variant, is locally stable around the LASSO solution.