Error analysis of tau-leap simulation methods. (English) Zbl 1234.60066
Summary: 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.
MSC:
| 60H35 | Computational methods for stochastic equations (aspects of stochastic analysis) |
| 65C40 | Numerical analysis or methods applied to Markov chains |
| 92C40 | Biochemistry, molecular biology |
Keywords:
tau-leaping; simulation; error analysis; reaction networks; Markov chain; chemical master equationReferences:
| [1] | Anderson, D. F. (2007). A modified next reaction method for simulating chemical systems with time dependent propensities and delays. J. Chem. Phys. 127 214107. |
| [2] | Anderson, D. F. (2008). Incorporating postleap checks in tau-leaping. J. Chem. Phys. 128 054103. |
| [3] | Cao, Y., Gillespie, D. T. and Petzold, L. R. (2005). Avoiding negative populations in explicit Poisson tau-leaping. J. Chem. Phys. 123 054104. |
| [4] | Cao, Y., Gillespie, D. T. and Petzold, L. R. (2006). Efficient step size selection for the tau-leaping simulation method. J. Chem. Phys. 124 044109. |
| [5] | Chatterjee, A., Vlachos, D. G. and Katsoulakis, M. A. (2005). Binomial distribution based tau-leap accelerated stochastic simulation. J. Chem. Phys. 122 024112. |
| [6] | Ethier, S. N. and Kurtz, T. G. (1986). Markov Processes : Characterization and Convergence . Wiley, New York. · Zbl 0592.60049 |
| [7] | Ethier, S. N. and Kurtz, T. G. (2005). Markov Processes : Characterization and Convergence , 2nd ed. Wiley, New York. · Zbl 1089.60005 |
| [8] | Gibson, M. A. and Bruck, J. (2000). Efficient exact stochastic simulation of chemical systems with many species and many channels. J. Phys. Chem. A 105 1876-1889. |
| [9] | Gillespie, D. T. (1976). A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J. Comput. Phys. 22 403-434. · doi:10.1016/0021-9991(76)90041-3 |
| [10] | Gillespie, D. T. (1977). Exact stochastic simulation of coupled chemical reactions. J. Phys. Chem. 81 2340-2361. |
| [11] | Gillespie, D. T. (2001). Approximate accelerated simulation of chemically reaction systems. J. Chem. Phys. 115 1716-1733. |
| [12] | Gillespie, D. T. and Petzold, L. R. (2003). Improved leap-size selection for accelerated stochastic simulation. J. Chem. Phys. 119 8229-8234. |
| [13] | Kloeden, P. E. and Platen, E. (1992). Numerical Solution of Stochastic Differential Equations. Applications of Mathematics ( New York ) 23 . Springer, Berlin. · Zbl 0752.60043 |
| [14] | Kurtz, T. G. (1972). The relationship between stochastic and deterministic models for chemical reactions. J. Chem. Phys. 57 2976-2978. |
| [15] | Kurtz, T. G. (1977/78). Strong approximation theorems for density dependent Markov chains. Stochastic Processes Appl. 6 223-240. · Zbl 0373.60085 · doi:10.1016/0304-4149(78)90020-0 |
| [16] | Kurtz, T. G. (1981). Approximation of Population Processes. CBMS-NSF Regional Conference Series in Applied Mathematics 36 . SIAM, Philadelphia, PA. · Zbl 0465.60078 |
| [17] | Li, T. (2007). Analysis of explicit tau-leaping schemes for simulating chemically reacting systems. Multiscale Model. Simul. 6 417-436 (electronic). · Zbl 1145.92041 · doi:10.1137/06066792X |
| [18] | Rathinam, M., Petzold, L. R., Cao, Y. and Gillespie, D. T. (2005). Consistency and stability of tau-leaping schemes for chemical reaction systems. Multiscale Model. Simul. 4 867-895 (electronic). · Zbl 1089.60516 · doi:10.1137/040603206 |
| [19] | Wilkinson, D. J. (2006). Stochastic Modelling for Systems Biology . Chapman and Hall/CRC, Boca Raton, FL. · Zbl 1099.92004 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
