English

Compressive Recovery of Sparse Precision Matrices

Machine Learning 2023-12-13 v3 Machine Learning

Abstract

We consider the problem of learning a graph modeling the statistical relations of the dd variables from a dataset with nn samples XRn×dX \in \mathbb{R}^{n \times d}. Standard approaches amount to searching for a precision matrix Θ\Theta representative of a Gaussian graphical model that adequately explains the data. However, most maximum likelihood-based estimators usually require storing the d2d^{2} values of the empirical covariance matrix, which can become prohibitive in a high-dimensional setting. In this work, we adopt a compressive viewpoint and aim to estimate a sparse Θ\Theta from a \emph{sketch} of the data, i.e. a low-dimensional vector of size md2m \ll d^{2} carefully designed from XX using non-linear random features. Under certain assumptions on the spectrum of Θ\Theta (or its condition number), we show that it is possible to estimate it from a sketch of size m=Ω((d+2k)log(d))m=\Omega\left((d+2k)\log(d)\right) where kk is the maximal number of edges of the underlying graph. These information-theoretic guarantees are inspired by compressed sensing theory and involve restricted isometry properties and instance optimal decoders. We investigate the possibility of achieving practical recovery with an iterative algorithm based on the graphical lasso, viewed as a specific denoiser. We compare our approach and graphical lasso on synthetic datasets, demonstrating its favorable performance even when the dataset is compressed.

Keywords

Cite

@article{arxiv.2311.04673,
  title  = {Compressive Recovery of Sparse Precision Matrices},
  author = {Titouan Vayer and Etienne Lasalle and Rémi Gribonval and Paulo Gonçalves},
  journal= {arXiv preprint arXiv:2311.04673},
  year   = {2023}
}
R2 v1 2026-06-28T13:15:06.774Z