Abstract:Suppose we observe samples of a subset of a collection of random variables. No additional information is provided about the number of latent variables, nor of the relationship between the latent and observed variables. Is it possible to discover the number of latent components, and to learn a statistical model over the entire collection of variables? We address this question in the setting in which the latent and observed variables are jointly Gaussian, with the conditional statistics of the observed variables conditioned on the latent variables being specified by a graphical model. As a first step we give natural conditions under which such latent-variable Gaussian graphical models are identifiable given marginal statistics of only the observed variables. Essentially these conditions require that the conditional graphical model among the observed variables is sparse, while the effect of the latent variables is "spread out" over most of the observed variables. Next we propose a tractable convex program based on regularized maximum-likelihood for model selection in this latent-variable setting; the regularizer uses both the $\ell_1$ norm and the nuclear norm. Our modeling framework can be viewed as a combination of dimensionality reduction (to identify latent variables) and graphical modeling (to capture remaining statistical structure not attributable to the latent variables), and it consistently estimates both the number of latent components and the conditional graphical model structure among the observed variables. These results are applicable in the high-dimensional setting in which the number of latent/observed variables grows with the number of samples of the observed variables. The geometric properties of the algebraic varieties of sparse matrices and of low-rank matrices play an important role in our analysis.

Lévy graphical models

Graphical models for infinite measures with applications to extremes

Graphical models, exponential families, and variational inference

Learning extremal graphical structures in high dimensions

Latent variable graphical model selection via convex optimization

Learning Latent Variable Gaussian Graphical Models

Learning Sparse High-Dimensional Matrix-Valued Graphical Models From Dependent Data

On Graphical Models via Univariate Exponential Family Distributions

An introduction to Lévy processes with applications in finance

Learning graphical models for stationary time series

Probability Based Independence Sampler for Bayesian Quantitative Learning in Graphical Log-Linear Marginal Models

Efficient Learning of Discrete Graphical Models

Bayesian Models of Graphs, Arrays and Other Exchangeable Random Structures

Multivariate Bernoulli distribution

Graphical model selection with latent variables

Sparse Graphical Linear Dynamical Systems

Application of Lévy Processes in Modelling (geodetic) Time Series with Mixed Spectra

On Sparse High-Dimensional Graphical Model Learning For Dependent Time Series

Multivariate Sparse Laplacian Shrinkage for Joint Estimation of Two Graphical Structures

Reliable and Efficient Inference of Bayesian Networks from Sparse Data by Statistical Learning Theory