Sparse Directed Network Recovery From Multivariate Count Time Series
Abstract
Count time series arise naturally in many scientific applications. Examples include species abundances collected over time, and gene expression counts over pseudotime. Network inference from such data can yield important insight into underlying biological mechanisms, but is challenging because count observations are discrete, overdispersed, and temporally dependent. We propose a multivariate Poisson log-normal model in which observed counts arise from a Poisson distribution with rates determined by a latent vector autoregressive (VAR) process. The latent VAR transition matrix induces directed temporal dependencies among components and serves as the target network parameter. We establish identifiability of the model under suitable conditions and develop two procedures for sparse estimation of the VAR transition matrix. Our first estimator is a sparse method-of-moments estimator, for which we establish consistent network recovery and parameter estimation in high-dimensional settings. Our second estimator is an L1-penalized variational inference estimator that is motivated by recent approaches for count data. We evaluate the proposed estimators through simulations under a range of sampling regimes. We further demonstrate the method on a macroinvertebrate abundance time series collected from European freshwater ecosystems, obtaining an interpretable ecological network.
Related
- Graphon Hawkes Processes: Sparse Estimation and Network Recovery
- Penalized Likelihood Inference in Highly Multivariate Sparse Spectral Graphical Gaussian Processes
- The spectral edge of sparse directed Erdös-Rényi graphs
- Low-Rank and Structured Sparse Tensor Decomposition for Anomaly Detection in Multivariate Functional Data
- Multi-Index Recovery for Total-Degree Multivariate Vandermonde Matrices
BibTeX
@misc{mathur2026sparse,
title = {{Sparse Directed Network Recovery From Multivariate Count Time Series}},
author = {Shirley Mathur and Nir Keret and Ali Shojaie},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.39236},
eprint = {2610.39236},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.39236},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Sparse Directed Network Recovery From Multivariate Count Time Series
AU - Shirley Mathur
AU - Nir Keret
AU - Ali Shojaie
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2610.39236
UR - https://arxiv.org/abs/2610.39236
AB - Count time series arise naturally in many scientific applications. Examples include species abundances collected over time, and gene expression counts over pseudotime. Network inference from such data can yield important insight into underlying biological mechanisms, but is challenging because count observations are discrete, overdispersed, and temporally dependent. We propose a multivariate Poisson log-normal model in which observed counts arise from a Poisson distribution with rates determined by a latent vector autoregressive (VAR) process. The latent VAR transition matrix induces directed temporal dependencies among components and serves as the target network parameter. We establish identifiability of the model under suitable conditions and develop two procedures for sparse estimation of the VAR transition matrix. Our first estimator is a sparse method-of-moments estimator, for which we establish consistent network recovery and parameter estimation in high-dimensional settings. Our second estimator is an L1-penalized variational inference estimator that is motivated by recent approaches for count data. We evaluate the proposed estimators through simulations under a range of sampling regimes. We further demonstrate the method on a macroinvertebrate abundance time series collected from European freshwater ecosystems, obtaining an interpretable ecological network.
ER -