A comparison of computational efficiencies of stochastic algorithms in terms of two infection models

  • Received: 01 November 2011 Accepted: 29 June 2018 Published: 01 July 2012
  • MSC : 60J27, 60J22, 92D25.

  • In this paper, we investigate three particular algorithms: a stochastic simulation algorithm (SSA), and explicit and implicit tau-leaping algorithms. To compare these methods, we used them to analyze two infection models: a Vancomycin-resistant enterococcus (VRE) infection model at the population level, and a Human Immunodeficiency Virus (HIV) within host infection model. While the first has a low species count and few transitions, the second is more complex with a comparable number of species involved. The relative efficiency of each algorithm is determined based on computational time and degree of precision required. The numerical results suggest that all three algorithms have the similar computational efficiency for the simpler VRE model, and the SSA is the best choice due to its simplicity and accuracy. In addition, we have found that with the larger and more complex HIV model, implementation and modification of tau-Leaping methods are preferred.

    Citation: H. Thomas Banks, Shuhua Hu, Michele Joyner, Anna Broido, Brandi Canter, Kaitlyn Gayvert, Kathryn Link. A comparison of computational efficiencies of stochastic algorithms in terms of two infection models[J]. Mathematical Biosciences and Engineering, 2012, 9(3): 487-526. doi: 10.3934/mbe.2012.9.487

    Related Papers:

    [1] Damilola Olabode, Libin Rong, Xueying Wang . Stochastic investigation of HIV infection and the emergence of drug resistance. Mathematical Biosciences and Engineering, 2022, 19(2): 1174-1194. doi: 10.3934/mbe.2022054
    [2] Linda J. S. Allen, P. van den Driessche . Stochastic epidemic models with a backward bifurcation. Mathematical Biosciences and Engineering, 2006, 3(3): 445-458. doi: 10.3934/mbe.2006.3.445
    [3] Linda J. S. Allen, Vrushali A. Bokil . Stochastic models for competing species with a shared pathogen. Mathematical Biosciences and Engineering, 2012, 9(3): 461-485. doi: 10.3934/mbe.2012.9.461
    [4] ZongWang, Qimin Zhang, Xining Li . Markovian switching for near-optimal control of a stochastic SIV epidemic model. Mathematical Biosciences and Engineering, 2019, 16(3): 1348-1375. doi: 10.3934/mbe.2019066
    [5] Damilola Olabode, Jordan Culp, Allison Fisher, Angela Tower, Dylan Hull-Nye, Xueying Wang . Deterministic and stochastic models for the epidemic dynamics of COVID-19 in Wuhan, China. Mathematical Biosciences and Engineering, 2021, 18(1): 950-967. doi: 10.3934/mbe.2021050
    [6] Sukhitha W. Vidurupola, Linda J. S. Allen . Basic stochastic models for viral infection within a host. Mathematical Biosciences and Engineering, 2012, 9(4): 915-935. doi: 10.3934/mbe.2012.9.915
    [7] Ying He, Junlong Tao, Bo Bi . Stationary distribution for a three-dimensional stochastic viral infection model with general distributed delay. Mathematical Biosciences and Engineering, 2023, 20(10): 18018-18029. doi: 10.3934/mbe.2023800
    [8] Qi Zhou, Huaimin Yuan, Qimin Zhang . Dynamics and approximation of positive solution of the stochastic SIS model affected by air pollutants. Mathematical Biosciences and Engineering, 2022, 19(5): 4481-4505. doi: 10.3934/mbe.2022207
    [9] Gianni Gilioli, Sara Pasquali, Fabrizio Ruggeri . Nonlinear functional response parameter estimation in a stochastic predator-prey model. Mathematical Biosciences and Engineering, 2012, 9(1): 75-96. doi: 10.3934/mbe.2012.9.75
    [10] Rubayyi T. Alqahtani, M. A. Abdelkawy . An efficient numerical algorithm for solving fractional SIRC model with salmonella bacterial infection. Mathematical Biosciences and Engineering, 2020, 17(4): 3784-3793. doi: 10.3934/mbe.2020212
  • In this paper, we investigate three particular algorithms: a stochastic simulation algorithm (SSA), and explicit and implicit tau-leaping algorithms. To compare these methods, we used them to analyze two infection models: a Vancomycin-resistant enterococcus (VRE) infection model at the population level, and a Human Immunodeficiency Virus (HIV) within host infection model. While the first has a low species count and few transitions, the second is more complex with a comparable number of species involved. The relative efficiency of each algorithm is determined based on computational time and degree of precision required. The numerical results suggest that all three algorithms have the similar computational efficiency for the simpler VRE model, and the SSA is the best choice due to its simplicity and accuracy. In addition, we have found that with the larger and more complex HIV model, implementation and modification of tau-Leaping methods are preferred.


  • This article has been cited by:

    1. P.A. Maginnis, M. West, G.E. Dullerud, Variance-reduced simulation of lattice discrete-time Markov chains with applications in reaction networks, 2016, 322, 00219991, 400, 10.1016/j.jcp.2016.06.019
    2. P. A. Maginnis, M. West, G. E. Dullerud, Exact Variance-Reduced Simulation of Lattice Continuous-Time Markov Chains with Applications in Reaction Networks, 2019, 81, 0092-8240, 3159, 10.1007/s11538-019-00576-2
  • Reader Comments
  • © 2012 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3365) PDF downloads(504) Cited by(2)

Article outline

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog