{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,29]],"date-time":"2026-04-29T18:10:43Z","timestamp":1777486243200,"version":"3.51.4"},"reference-count":18,"publisher":"Oxford University Press (OUP)","issue":"6","license":[{"start":{"date-parts":[[2022,1,5]],"date-time":"2022-01-05T00:00:00Z","timestamp":1641340800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/academic.oup.com\/journals\/pages\/open_access\/funder_policies\/chorus\/standard_publication_model"}],"funder":[{"DOI":"10.13039\/100000002","name":"NIH","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000057","name":"NIGMS","doi-asserted-by":"publisher","award":["R01GM111510"],"award-info":[{"award-number":["R01GM111510"]}],"id":[{"id":"10.13039\/100000057","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100007000","name":"LDRD","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100007000","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100008902","name":"Los Alamos National Laboratory","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100008902","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2022,3,4]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Summary<\/jats:title><jats:p>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.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded\/cloned from GitHub and used and modified under terms of the BSD-3 license (https:\/\/github.com\/lanl\/pybnf). Online documentation covering installation\/usage is available (https:\/\/pybnf.readthedocs.io\/en\/latest\/). A tutorial video is available on YouTube (https:\/\/www.youtube.com\/watch?v=2aRqpqFOiS4&amp;t=63s).<\/jats:p><\/jats:sec><jats:sec><jats:title>Supplementary information<\/jats:title><jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btac004","type":"journal-article","created":{"date-parts":[[2022,1,3]],"date-time":"2022-01-03T20:09:57Z","timestamp":1641240597000},"page":"1770-1772","source":"Crossref","is-referenced-by-count":12,"title":["Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit"],"prefix":"10.1093","volume":"38","author":[{"given":"Jacob","family":"Neumann","sequence":"first","affiliation":[{"name":"Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ 86011, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6893-8423","authenticated-orcid":false,"given":"Yen Ting","family":"Lin","sequence":"additional","affiliation":[{"name":"Information Sciences Group, Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory , Los Alamos, NM 87545, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-3588-5358","authenticated-orcid":false,"given":"Abhishek","family":"Mallela","sequence":"additional","affiliation":[{"name":"Department of Mathematics, University of California , Davis, CA 95616, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Ely F","family":"Miller","sequence":"additional","affiliation":[{"name":"Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ 86011, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Joshua","family":"Colvin","sequence":"additional","affiliation":[{"name":"Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ 86011, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Abell T","family":"Duprat","sequence":"additional","affiliation":[{"name":"Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ 86011, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Ye","family":"Chen","sequence":"additional","affiliation":[{"name":"Department of Mathematics and Statistics, Northern Arizona University, Flagstaff, AZ 86011, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-4383-8711","authenticated-orcid":false,"given":"William S","family":"Hlavacek","sequence":"additional","affiliation":[{"name":"Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory , Los Alamos, NM 87545, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Richard G","family":"Posner","sequence":"additional","affiliation":[{"name":"Department of Biological Sciences, Northern Arizona University , Flagstaff, AZ 86011, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2022,1,5]]},"reference":[{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"5","DOI":"10.1023\/A:1020281327116","article-title":"An introduction to MCMC for machine learning","volume":"50","author":"Andrieu","year":"2003","journal-title":"Mach. Learn"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"343","DOI":"10.1007\/s11222-008-9110-y","article-title":"A tutorial on adaptive MCMC","volume":"18","author":"Andrieu","year":"2008","journal-title":"Stat. Comput"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"74","DOI":"10.1016\/j.biosystems.2018.07.006","article-title":"Tellurium: an extensible python-based modeling environment for systems and synthetic biology","volume":"171","author":"Choi","year":"2018","journal-title":"Biosystems"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"644","DOI":"10.1038\/msb.2012.69","article-title":"Properties of cell death models calibrated and compared using Bayesian approaches","volume":"9","author":"Eydgahi","year":"2013","journal-title":"Mol. Syst. Biol"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"113","DOI":"10.1007\/978-1-59745-525-1_5","article-title":"Rule-based modeling of biochemical systems with BioNetGen","volume":"500","author":"Faeder","year":"2009","journal-title":"Methods Mol. Biol"},{"key":"2023033101030311500_","volume-title":"Bayesian Data Analysis","author":"Gelman","year":"2014","edition":"3rd edn"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"690","DOI":"10.1109\/PDP2018.2018.00114","article-title":"Evaluation of parallel tempering to accelerate Bayesian parameter estimation in systems biology","author":"Gupta","year":"2018","journal-title":"2018 26th Euromicro International Conference on Parallel, Distributed and Network-based Processing (PDP)"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"15586","DOI":"10.1038\/s41598-017-15568-2","article-title":"Timescale separation of positive and negative signaling creates history-dependent responses to IgE receptor stimulation","volume":"7","author":"Harmon","year":"2017","journal-title":"Sci. Rep"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"3366","DOI":"10.1093\/bioinformatics\/btw469","article-title":"BioNetGen 2.2: advances in rule-based modeling","volume":"32","author":"Harris","year":"2016","journal-title":"Bioinformatics"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"123","DOI":"10.1038\/373123a0","article-title":"Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection","volume":"373","author":"Ho","year":"1995","journal-title":"Nature"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"e9110","DOI":"10.15252\/msb.20199110","article-title":"SBML Level 3: an extensible format for the exchange and reuse of biological models","volume":"16","author":"Keating","year":"2020","journal-title":"Mol. Syst. Biol"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"767","DOI":"10.3201\/eid2703.203364","article-title":"Daily forecasting of regional epidemics of Coronavirus Disease with Bayesian uncertainty quantification, United States","volume":"27","author":"Lin","year":"2021","journal-title":"Emerg. Inf. Dis"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"646","DOI":"10.1038\/msb.2013.1","article-title":"Programming biological models in Python using PySB","volume":"9","author":"Lopez","year":"2013","journal-title":"Mol. Syst. Biol"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"1012","DOI":"10.1016\/j.isci.2019.08.045","article-title":"PyBioNetFit and the biological property specification language","volume":"19","author":"Mitra","year":"2019","journal-title":"iScience"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"1582","DOI":"10.1126\/science.271.5255.1582","article-title":"HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time","volume":"271","author":"Perelson","year":"1996","journal-title":"Science"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"695","DOI":"10.1093\/bioinformatics\/btx626","article-title":"PyDREAM: high-dimensional parameter inference for biological models in python","volume":"34","author":"Shockley","year":"2018","journal-title":"Bioinformatics"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"3315","DOI":"10.1093\/bioinformatics\/btv363","article-title":"LibRoadRunner: a high performance SBML simulation and analysis library","volume":"31","author":"Somogyi","year":"2015","journal-title":"Bioinformatics"},{"key":"2023033101030311500_","doi-asserted-by":"crossref","first-page":"798","DOI":"10.1093\/bioinformatics\/btv655","article-title":"BioNetFit: a fitting tool compatible with BioNetGen, NFsim and distributed computing environments","volume":"32","author":"Thomas","year":"2016","journal-title":"Bioinformatics"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btac004\/42169834\/btac004.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/38\/6\/1770\/49692951\/btac004.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/38\/6\/1770\/49692951\/btac004.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,11,15]],"date-time":"2023-11-15T07:04:31Z","timestamp":1700031871000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/38\/6\/1770\/6497784"}},"subtitle":[],"editor":[{"given":"Christina","family":"Kendziorski","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2022,1,5]]},"references-count":18,"journal-issue":{"issue":"6","published-print":{"date-parts":[[2022,3,4]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btac004","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2022,3,15]]},"published":{"date-parts":[[2022,1,5]]}}}