Error analysis of tau-leap simulation methods

Abstract

We perform an error analysis for numerical approximation methods of continuous time Markov chain models commonly found in the chemistry and biochemistry literature. The motivation for the analysis is to be able to compare the accuracy of different approximation methods and, specifically, Euler tau-leaping and midpoint tau-leaping. We perform our analysis under a scaling in which the size of the time discretization is inversely proportional to some (bounded) power of the norm of the state of the system. We argue that this is a more appropriate scaling than that found in previous error analyses in which the size of the time discretization goes to zero independent of the rest of the model. Under the present scaling, we show that midpoint tau-leaping achieves a higher order of accuracy, in both a weak and a strong sense, than Euler tau-leaping; a result that is in contrast to previous analyses. We present examples that demonstrate our findings.

Published in at http://dx.doi.org/10.1214/10-AAP756 the Annals of Applied Probability (http://www.imstat.org/aap/) by the Institute of Mathematical Statistics (http://www.imstat.org)

References in corpus (2)

• A modified Next Reaction Method for simulating chemical systems with time dependent propensities and delays
• Incorporating postleap checks in tau-leaping

Cited by in corpus (13)

• Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
• Stochastic Representations of Ion Channel Kinetics and Exact Stochastic Simulation of Neuronal Dynamics
• Extending the multi-level method for the simulation of stochastic biological systems
• State and parameter estimation from exact partial state observation in stochastic reaction networks
• Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes
• Multilevel Hybrid Split-Step Implicit Tau-Leap
• Quasi-Monte Carlo methods applied to tau-leaping in stochastic biological systems
• Tensor product approach to modelling epidemics on networks
• Pathwise error bounds in Multiscale variable splitting methods for spatial stochastic kinetics
• Strong convergence for split-step methods in stochastic jump kinetics
• Slow-scale split-step tau-leap method for stiff stochastic chemical systems
• Importance sampling for a robust and efficient multilevel Monte Carlo estimator for stochastic reaction networks
• Fast approximate simulation of finite long-range spin systems