Large-Scale Precision Matrix Estimation With SQUIC
15 Pages Posted: 16 Aug 2021
Date Written: August 12, 2021
Abstract
High-dimensional sparse precision matrix estimation is a ubiquitous task in multivariate analysis with applications that cross many disciplines. In this paper, we introduce the SQUIC package, which benefits from superior runtime performance and scalability, significantly exceeding the available state-of-the-art packages. This package is a second-order method that solves the L1--regularized maximum likelihood problem using highly optimized linear algebra subroutines, which leverage the underlying sparsity and the intrinsic parallelism in the computation. We provide two sets of numerical tests; the first one consists of didactic examples using synthetic datasets highlighting the performance and accuracy of the package, and the second one is a real-world classification problem of high dimensional medical datasets. The base algorithm is implemented in C++ with interfaces for R and Python.
Keywords: SQUIC, precision matrix, LASSO, maximum likelihood,R,Python
JEL Classification: C02, C10, C13
Suggested Citation: Suggested Citation