English

Simulation of stochastic network dynamics via entropic matching

Quantitative Methods 2013-03-14 v2 Data Analysis, Statistics and Probability

Abstract

The simulation of complex stochastic network dynamics arising, for instance, from models of coupled biomolecular processes remains computationally challenging. Often, the necessity to scan a models' dynamics over a large parameter space renders full-fledged stochastic simulations impractical, motivating approximation schemes. Here we propose an approximation scheme which improves upon the standard linear noise approximation while retaining similar computational complexity. The underlying idea is to minimize, at each time step, the Kullback-Leibler divergence between the true time evolved probability distribution and a Gaussian approximation (entropic matching). This condition leads to ordinary differential equations for the mean and the covariance matrix of the Gaussian. For cases of weak nonlinearity, the method is more accurate than the linear method when both are compared to stochastic simulations.

Keywords

Cite

@article{arxiv.1209.3700,
  title  = {Simulation of stochastic network dynamics via entropic matching},
  author = {Tiago Ramalho and Marco Selig and Ulrich Gerland and Torsten A. Enßlin},
  journal= {arXiv preprint arXiv:1209.3700},
  year   = {2013}
}

Comments

23 pages, 6 figures; significantly revised version

R2 v1 2026-06-21T22:06:36.451Z