Abstract:We consider the problem of learning a sparse graph under the Laplacian constrained Gaussian graphical models. This problem can be formulated as a penalized maximum likelihood estimation of the Laplacian constrained precision matrix. Like in the classical graphical lasso problem, recent works made use of the $\ell_1$-norm regularization with the goal of promoting sparsity in Laplacian constrained precision matrix estimation. However, we find that the widely used $\ell_1$-norm is not effective in imposing a sparse solution in this problem. Through empirical evidence, we observe that the number of nonzero graph weights grows with the increase of the regularization parameter. From a theoretical perspective, we prove that a large regularization parameter will surprisingly lead to a complete graph, i.e., every pair of vertices is connected by an edge. To address this issue, we introduce the nonconvex sparsity penalty, and propose a new estimator by solving a sequence of weighted $\ell_1$-norm penalized sub-problems. We establish the non-asymptotic optimization performance guarantees on both optimization error and statistical error, and prove that the proposed estimator can recover the edges correctly with a high probability. To solve each sub-problem, we develop a projected gradient descent algorithm which enjoys a linear convergence rate. Finally, an extension to learn disconnected graphs is proposed by imposing additional rank constraint. We propose a numerical algorithm based on based on the alternating direction method of multipliers, and establish its theoretical sequence convergence. Numerical experiments involving synthetic and real-world data sets demonstrate the effectiveness of the proposed method.

Estimation of the L2-Norm and Testing in Sparse Linear Regression with Unknown Variance

Minimax rate of testing in sparse linear regression

Sparse Estimation Via ℓ_q Optimization Method in High-Dimensional Linear Regression

Adaptive robust estimation in sparse vector model

A Study of Error Variance Estimation in Lasso Regression

Regularity Properties for Sparse Regression

Sparse Estimation Via Lower-Order Penalty Optimization Methods in High-Dimensional Linear Regression.

Minimax estimation of functionals in sparse vector model with correlated observations

Distributed Sparse Total Least-Squares over Networks

Uniform Post Selection Inference for LAD Regression and Other Z-estimation problems

Adaptive posterior concentration rates for sparse high-dimensional linear regression with random design and unknown error variance

Scaled sparse linear regression

Conditional Sparse Linear Regression

Linear Programming-Based Sparse Kernel Regression with L1-Norm Minimization for Nonlinear System Modeling

Non-uniform Nyström approximation for sparse kernel regression: Theoretical analysis and experimental evaluation

Does the $\ell_1$-norm Learn a Sparse Graph under Laplacian Constrained Graphical Models?

A note on the minimax risk of sparse linear regression

Variance estimation for sparse ultra-high dimensional varying coefficient models

The Influence Function of Penalized Regression Estimators

Minimax Instrumental Variable Regression and L_2 Convergence Guarantees Without Identification or Closedness

Sparse estimation of a covariance matrix