Open Access
December 2011 Error analysis of tau-leap simulation methods
David F. Anderson, Arnab Ganguly, Thomas G. Kurtz
Ann. Appl. Probab. 21(6): 2226-2262 (December 2011). DOI: 10.1214/10-AAP756

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.

Citation

Download Citation

David F. Anderson. Arnab Ganguly. Thomas G. Kurtz. "Error analysis of tau-leap simulation methods." Ann. Appl. Probab. 21 (6) 2226 - 2262, December 2011. https://doi.org/10.1214/10-AAP756

Information

Published: December 2011
First available in Project Euclid: 23 November 2011

zbMATH: 1234.60066
MathSciNet: MR2895415
Digital Object Identifier: 10.1214/10-AAP756

Subjects:
Primary: 60H35 , 65C99
Secondary: 92C40

Keywords: chemical master equation , error analysis , Markov chain , Reaction networks , simulation , Tau-leaping

Rights: Copyright © 2011 Institute of Mathematical Statistics

Vol.21 • No. 6 • December 2011
Back to Top