Abstract:In multivariate statistics, the question of finding direct interactions can be formulated as a problem of network inference - or network reconstruction - for which the Gaussian graphical model (GGM) provides a canonical framework. Unfortunately, the Gaussian assumption does not apply to count data which are encountered in domains such as genomics, social sciences or ecology. To circumvent this limitation, state-of-the-art approaches use two-step strategies that first transform counts to pseudo Gaussian observations and then apply a (partial) correlation-based approach from the abundant literature of GGM inference. We adopt a different stance by relying on a latent model where we directly model counts by means of Poisson distributions that are conditional to latent (hidden) Gaussian correlated variables. In this multivariate Poisson lognormal-model, the dependency structure is completely captured by the latent layer. This parametric model enables to account for the effects of covariates on the counts. To perform network inference, we add a sparsity inducing constraint on the inverse covariance matrix of the latent Gaussian vector. Unlike the usual Gaussian setting, the penalized likelihood is generally not tractable, and we resort instead to a variational approach for approximate likelihood maximization. The corresponding optimization problem is solved by alternating a gradient ascent on the variational parameters and a graphical-Lasso step on the covariance matrix. We show that our approach is highly competitive with the existing methods on simulation inspired from microbiological data. We then illustrate on three various data sets how accounting for sampling efforts via offsets and integrating external covariates (which is mostly never done in the existing literature) drastically changes the topology of the inferred network.

Topology Inference of Directed Graphs by Gaussian Processes with Sparsity Constraints

Network Topology Inference with Sparsity and Laplacian Constraints

Learning directed acyclic graphs based on sparsest permutations

High dimensional sparse covariance estimation via directed acyclic graphs

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

Sparse Graph Learning from Spatiotemporal Time Series

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

Tensor Graphical Model: Non-Convex Optimization and Statistical Inference

Statistical Inference on High Dimensional Gaussian Graphical Regression Models

Spectral Inference Methods on Sparse Graphs: Theory and Applications

Online Graph Topology Learning from Matrix-valued Time Series

Interpretable and Scalable Graphical Models for Complex Spatio-temporal Processes

Efficient Learning Algorithms for Gaussian Processes

Dynamic Chain Graph Models for Ordinal Time Series Data

On Topology Inference for Networked Dynamical Systems: Principles and Performances.

Topological Signal Processing and Learning: Recent Advances and Future Challenges

Variational inference for sparse network reconstruction from count data

Bayesian Inference for Relational Graph in a Causal Vector Autoregressive Time Series

Learning Cartesian Product Graphs with Laplacian Constraints

Bayesian inference of transition matrices from incomplete graph data with a topological prior

A Hierarchical Graph Signal Processing Approach to Inference from Spatiotemporal Signals