A Critical Comparison of Rejection-Based Algorithms for Simulation of Large Biochemical Reaction Networks

No Thumbnail Available

Access rights

openAccess
acceptedVersion

URL

Journal Title

Journal ISSN

Volume Title

A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä

Date

2019-08-01

Major/Subject

Mcode

Degree programme

Language

en

Pages

21

Series

Bulletin of Mathematical Biology, Volume 81, issue 8, pp. 3053-3073

Abstract

The rejection-based simulation technique has been applying to improve the computational efficiency of the stochastic simulation algorithm (SSA) in simulating large reaction networks, which are required for a thorough understanding of biological systems. We compare two recently proposed simulation methods, namely the composition–rejection algorithm (SSA-CR) and the rejection-based SSA (RSSA), aiming for this purpose. We discuss the right interpretation of the rejection-based technique used in these algorithms in order to make an informed choice when dealing with different aspects of biochemical networks. We provide the theoretical analysis as well as the detailed runtime comparison of these algorithms on concrete biological models. We highlight important factors that are omitted in previous analysis of these algorithms. The numerical comparison shows that for reaction networks where the search cost is expensive then SSA-CR is more efficient, and for reaction networks where the update cost is dominant, often the case in practice, then RSSA should be the choice.

Description

Keywords

Computational biology, Rejection-based simulation technique, Stochastic simulation

Other note

Citation

Thanh, V H 2019, ' A Critical Comparison of Rejection-Based Algorithms for Simulation of Large Biochemical Reaction Networks ', Bulletin of Mathematical Biology, vol. 81, no. 8, pp. 3053-3073 . https://doi.org/10.1007/s11538-018-0462-y