English

Global multivariate model learning from hierarchically correlated data

Disordered Systems and Neural Networks 2021-02-12 v1 Statistical Mechanics Quantitative Methods

Abstract

Inverse statistical physics aims at inferring models compatible with a set of empirical averages estimated from a high-dimensional dataset of independently distributed equilibrium configurations of a given system. However, in several applications such as biology, data result from stochastic evolutionary processes, and configurations are related through a hierarchical structure, typically represented by a tree, and therefore not independent. In turn, empirical averages of observables superpose intrinsic signal related to the equilibrium distribution of the studied system and spurious historical (or phylogenetic) signal resulting from the structure underlying the data-generating process. The naive application of inverse statistical physics techniques therefore leads to systematic biases and an effective reduction of the sample size. To advance on the currently open task of extracting intrinsic signals from correlated data, we study a system described by a multivariate Ornstein-Uhlenbeck process defined on a finite tree. Using a Bayesian framework, we can disentangle covariances in the data corresponding to their multivariate Gaussian equilibrium distribution from those resulting from the historical correlations. Our approach leads to a clear gain in accuracy in the inferred equilibrium distribution, which corresponds to an effective two- to fourfold increase in sample size.

Keywords

Cite

@article{arxiv.2102.06036,
  title  = {Global multivariate model learning from hierarchically correlated data},
  author = {Edwin Rodriguez Horta and Alejandro Lage and Martin Weigt and Pierre Barrat-Charlaix},
  journal= {arXiv preprint arXiv:2102.06036},
  year   = {2021}
}

Comments

34 pages 10 figures

R2 v1 2026-06-23T23:04:16.786Z