Coupling sample paths to the partial thermodynamic limit in stochastic chemical reaction networks

Abstract

Many biochemical systems appearing in applications have a multiscale structure so that they converge to piecewise deterministic Markov processes in a thermodynamic limit. The statistics of the piecewise deterministic process can be obtained much more efficiently than those of the exact process. We explore the possibility of coupling sample paths of the exact model to the piecewise deterministic process in order to reduce the variance of their difference. We then apply this coupling to reduce the computational complexity of a Monte Carlo estimator. In addition to rigorous results concerning the asymptotic computational complexity of the Monte Carlo estimator, numerical simulations are performed on some simple biological models confirming that computational gains are made.

arXiv commentSome significant organizational changes have been made. The discussion of models with feedback is now correct and rigorous. More convincing numerical experiments have been performed

References in corpus (6)

• Intrinsic noise in systems with switching environments
• Stochastic Representations of Ion Channel Kinetics and Exact Stochastic Simulation of Neuronal Dynamics
• Hybrid framework for the simulation of stochastic chemical kinetics
• Extending the multi-level method for the simulation of stochastic biological systems
• Coupling Control Variates for Markov Chain Monte Carlo
• Pathwise error bounds in Multiscale variable splitting methods for spatial stochastic kinetics