Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Jacob Neumann, Yen Ting Lin, Abhishek Mallela, Ely F. Miller, Joshua Colvin, Abell T. Duprat, Ye Chen, William S. Hlavacek, Richard G. Posner

Research output: Contribution to journalArticlepeer-review

5 Scopus citations

Abstract

Summary: Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty.

Original languageEnglish (US)
Pages (from-to)1770-1772
Number of pages3
JournalBioinformatics
Volume38
Issue number6
DOIs
StatePublished - Mar 15 2022

ASJC Scopus subject areas

  • Statistics and Probability
  • Biochemistry
  • Molecular Biology
  • Computer Science Applications
  • Computational Theory and Mathematics
  • Computational Mathematics

Fingerprint

Dive into the research topics of 'Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit'. Together they form a unique fingerprint.

Cite this