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.

Penalized estimation of the Gaussian graphical model from data with replicates

A sparse conditional Gaussian graphical model for analysis of genetical genomics data

Penalized Maximum Likelihood Estimation of Multi-layered Gaussian Graphical Models

Joint estimation of multiple dependent Gaussian graphical models with applications to mouse genomics

Gaussian Graphical Model Estimation and Selection for High-Dimensional Incomplete Data Using Multiple Imputation and Horseshoe Estimators

L1-Penalized Censored Gaussian Graphical Model

Semiblind subgraph reconstruction in Gaussian graphical models

Joint estimation of multiple high-dimensional Gaussian copula graphical models

The conditional censored graphical lasso estimator

Frequentist model averaging for undirected gaussian graphical models

Algorithm xxx: A Covariate-Dependent Approach to Gaussian Graphical Modeling in R

Estimation of multiple networks in Gaussian mixture models.

Generating Independent Replicates Directly from the Posterior Distribution for a Class of Spatial Hierarchical Models

Variational inference for sparse network reconstruction from count data

A Log-Linear Graphical Model for Inferring Genetic Networks from High-Throughput Sequencing Data

Estimating graphical models for count data with applications to single-cell gene network

Joint Mean and Covariance Estimation with Unreplicated Matrix-Variate Data

High-Dimensional Gaussian Graphical Regression Models with Covariates

Enhancing the Power of Gaussian Graphical Model Inference by Modeling the Graph Structure

Statistical Inference on High Dimensional Gaussian Graphical Regression Models

Nonparametric Bayesian Knockoff Generators for Feature Selection Under Complex Data Structure