{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,1]],"date-time":"2026-01-01T13:48:00Z","timestamp":1767275280961},"reference-count":32,"publisher":"Oxford University Press (OUP)","issue":"20","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2015,10,15]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>Motivation: The field of phylodynamics focuses on the problem of reconstructing population size dynamics over time using current genetic samples taken from the population of interest. This technique has been extensively used in many areas of biology but is particularly useful for studying the spread of quickly evolving infectious diseases agents, e.g. influenza virus. Phylodynamic inference uses a coalescent model that defines a probability density for the genealogy of randomly sampled individuals from the population. When we assume that such a genealogy is known, the coalescent model, equipped with a Gaussian process prior on population size trajectory, allows for nonparametric Bayesian estimation of population size dynamics. Although this approach is quite powerful, large datasets collected during infectious disease surveillance challenge the state-of-the-art of Bayesian phylodynamics and demand inferential methods with relatively low computational cost.<\/jats:p>\n               <jats:p>Results: To satisfy this demand, we provide a computationally efficient Bayesian inference framework based on Hamiltonian Monte Carlo for coalescent process models. Moreover, we show that by splitting the Hamiltonian function, we can further improve the efficiency of this approach. Using several simulated and real datasets, we show that our method provides accurate estimates of population size dynamics and is substantially faster than alternative methods based on elliptical slice sampler and Metropolis-adjusted Langevin algorithm.<\/jats:p>\n               <jats:p>Availability and implementation: The R code for all simulation studies and real data analysis conducted in this article are publicly available at http:\/\/www.ics.uci.edu\/\u223cslan\/lanzi\/CODES.html and in the R package phylodyn available at https:\/\/github.com\/mdkarcher\/phylodyn.<\/jats:p>\n               <jats:p>Contact: \u00a0S.Lan@warwick.ac.uk or babaks@uci.edu<\/jats:p>\n               <jats:p>Supplementary information: Supplementary data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btv378","type":"journal-article","created":{"date-parts":[[2015,6,21]],"date-time":"2015-06-21T00:00:20Z","timestamp":1434844820000},"page":"3282-3289","source":"Crossref","is-referenced-by-count":37,"title":["An efficient Bayesian inference framework for coalescent-based nonparametric phylodynamics"],"prefix":"10.1093","volume":"31","author":[{"given":"Shiwei","family":"Lan","sequence":"first","affiliation":[{"name":"1 Department of Statistics, University of Warwick, Coventry CV4 7AL, UK,"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Julia A.","family":"Palacios","sequence":"additional","affiliation":[{"name":"2 Department of Organismic and Evolutionary Biology, Harvard University, MA 02138, US,"},{"name":"3 Department of Ecology and Evolutionary Biology, Brown University, RI 02912, US,"},{"name":"4 Center for Computational Molecular Biology, Brown University,"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Michael","family":"Karcher","sequence":"additional","affiliation":[{"name":"5 Department of Statistics, University of Washington, WA 98195, US,"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Vladimir N.","family":"Minin","sequence":"additional","affiliation":[{"name":"5 Department of Statistics, University of Washington, WA 98195, US,"},{"name":"6 Department of Biology, University of Washington and"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Babak","family":"Shahbaba","sequence":"additional","affiliation":[{"name":"7 Department of Statistics, University of California, Irvine, CA 92697, US"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2015,6,20]]},"reference":[{"key":"2023020202314601100_btv378-B1","first-page":"733","article-title":"On conditional and intrinsic autoregressions","volume":"82","author":"Besag","year":"1995","journal-title":"Biometrika"},{"key":"2023020202314601100_btv378-B2","doi-asserted-by":"crossref","first-page":"1969","DOI":"10.1093\/molbev\/mss075","article-title":"Bayesian phylogenetics with BEAUti and the BEAST 1.7","volume":"29","author":"Drummond","year":"2012","journal-title":"Mol. Biol. Evol."},{"key":"2023020202314601100_btv378-B3","doi-asserted-by":"crossref","first-page":"1307","DOI":"10.1093\/genetics\/161.3.1307","article-title":"Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data","volume":"161","author":"Drummond","year":"2002","journal-title":"Genetics"},{"key":"2023020202314601100_btv378-B4","doi-asserted-by":"crossref","first-page":"1185","DOI":"10.1093\/molbev\/msi103","article-title":"Bayesian coalescent inference of past population dynamics from molecular sequences","volume":"22","author":"Drummond","year":"2005","journal-title":"Mol. Biol. Evol."},{"key":"2023020202314601100_btv378-B5","doi-asserted-by":"crossref","first-page":"216","DOI":"10.1016\/0370-2693(87)91197-X","article-title":"Hybrid Monte Carlo","volume":"195","author":"Duane","year":"1987","journal-title":"Phys. Lett. B"},{"key":"2023020202314601100_btv378-B6","first-page":"473","article-title":"Practical Markov chain Monte Carlo","volume":"7","author":"Geyer","year":"1992","journal-title":"Stat. Sci."},{"key":"2023020202314601100_btv378-B7","doi-asserted-by":"crossref","first-page":"713","DOI":"10.1093\/molbev\/mss265","article-title":"Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci","volume":"30","author":"Gill","year":"2013","journal-title":"Mol. Biol. Evol."},{"key":"2023020202314601100_btv378-B8","doi-asserted-by":"crossref","first-page":"123","DOI":"10.1111\/j.1467-9868.2010.00765.x","article-title":"Riemann manifold Langevin and Hamiltonian Monte Carlo methods","volume":"73","author":"Girolami","year":"2011","journal-title":"J. R. Stat. Soc. B"},{"key":"2023020202314601100_btv378-B9","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1098\/rstb.1994.0079","article-title":"Sampling theory for neutral alleles in a varying environment","volume":"344","author":"Griffiths","year":"1994","journal-title":"Philos. Trans. R. Soc. Lond. B Biol. Sci."},{"key":"2023020202314601100_btv378-B10","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1186\/1471-2148-8-289","article-title":"Bayesian inference of population size history from multiple loci","volume":"8","author":"Heled","year":"2008","journal-title":"BMC Evol. Biol."},{"key":"2023020202314601100_btv378-B11","doi-asserted-by":"crossref","first-page":"35","DOI":"10.1115\/1.3662552","article-title":"A new approach to linear filtering and prediction problems","volume":"82(Series D)","author":"Kalman","year":"1960","journal-title":"Trans. ASME J. Basic Eng."},{"key":"2023020202314601100_btv378-B12","doi-asserted-by":"crossref","first-page":"235","DOI":"10.1016\/0304-4149(82)90011-4","article-title":"The coalescent","volume":"13","author":"Kingman","year":"1982","journal-title":"Stochastic Processes Appl."},{"key":"2023020202314601100_btv378-B13","doi-asserted-by":"crossref","first-page":"597","DOI":"10.1111\/1467-9469.00308","article-title":"On block updating in Markov random field models for disease mapping","volume":"29","author":"Knorr-Held","year":"2002","journal-title":"Scand. J. Stat."},{"key":"2023020202314601100_btv378-B14","doi-asserted-by":"crossref","first-page":"429","DOI":"10.1093\/genetics\/149.1.429","article-title":"Maximum likelihood estimation of population growth rates based on the coalescent","volume":"149","author":"Kuhner","year":"1998","journal-title":"Genetics"},{"key":"2023020202314601100_btv378-B15","author":"Lan","year":"2013"},{"key":"2023020202314601100_btv378-B16","volume-title":"Simulating Hamiltonian Dynamics","author":"Leimkuhler","year":"2004"},{"key":"2023020202314601100_btv378-B17","doi-asserted-by":"crossref","first-page":"1459","DOI":"10.1093\/molbev\/msn090","article-title":"Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics","volume":"25","author":"Minin","year":"2008","journal-title":"Mol. Biol. Evol."},{"key":"2023020202314601100_btv378-B18","first-page":"25: 451","article-title":"Log Gaussian Cox processes","author":"M\u00f8ller","year":"1998","journal-title":"Scand. J. Stat."},{"key":"2023020202314601100_btv378-B19","first-page":"541","article-title":"Elliptical slice sampling","volume":"9","author":"Murray","year":"2010","journal-title":"J. Machine Learn. Res. Workshop Conf. Proc."},{"key":"2023020202314601100_btv378-B20","first-page":"113","article-title":"MCMC using Hamiltonian dynamics","volume-title":"Handbook of Markov Chain Monte Carlo","author":"Neal","year":"2010"},{"key":"2023020202314601100_btv378-B21","doi-asserted-by":"crossref","first-page":"6","DOI":"10.1186\/1471-2148-5-6","article-title":"Inference of demographic history from genealogical trees using reversible jump Markov chain Monte Carlo","volume":"5","author":"Opgen-Rhein","year":"2005","journal-title":"BMC Evol. Biol."},{"key":"2023020202314601100_btv378-B22","doi-asserted-by":"crossref","first-page":"518","DOI":"10.1080\/10618600.2013.788448","article-title":"Exact Hamiltonian Monte Carlo for truncated multivariate Gaussians","volume":"23","author":"Pakman","year":"2014","journal-title":"J. Comput. Graphical Stat."},{"key":"2023020202314601100_btv378-B23","first-page":"726","article-title":"Integrated nested Laplace approximation for Bayesian nonparametric phylodynamics","volume-title":"UAI","author":"Palacios","year":"2012"},{"key":"2023020202314601100_btv378-B24","doi-asserted-by":"crossref","first-page":"8","DOI":"10.1111\/biom.12003","article-title":"Gaussian process-based Bayesian nonparametric inference of population size trajectories from gene genealogies","volume":"69","author":"Palacios","year":"2013","journal-title":"Biometrics"},{"key":"2023020202314601100_btv378-B25","doi-asserted-by":"crossref","first-page":"615","DOI":"10.1038\/nature06945","article-title":"The genomic and epidemiological dynamics of human influenza A virus","volume":"453","author":"Rambaut","year":"2008","journal-title":"Nature"},{"key":"2023020202314601100_btv378-B26","doi-asserted-by":"crossref","first-page":"337","DOI":"10.1023\/A:1023562417138","article-title":"Langevin diffusions and Metropolis-Hastings algorithms","volume":"4","author":"Roberts","year":"2002","journal-title":"Methodol. Comput. Appl. Probability"},{"key":"2023020202314601100_btv378-B27","doi-asserted-by":"crossref","first-page":"341","DOI":"10.2307\/3318418","article-title":"Exponential convergence of Langevin distributions and their discrete approximations","volume":"2","author":"Roberts","year":"1996","journal-title":"Bernoulli"},{"key":"2023020202314601100_btv378-B28","first-page":"233","article-title":"Coalescent approaches to HIV population genetics","volume-title":"The Evolution of HIV","author":"Rodrigo","year":"1999"},{"key":"2023020202314601100_btv378-B29","doi-asserted-by":"crossref","DOI":"10.1201\/9780203492024","volume-title":"Gaussian Markov Random Fields: Theory and Applications, volume 104 of Monographs on Statistics and Applied Probability","author":"Rue","year":"2005"},{"key":"2023020202314601100_btv378-B30","doi-asserted-by":"crossref","first-page":"319","DOI":"10.1111\/j.1467-9868.2008.00700.x","article-title":"Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations","volume":"71","author":"Rue","year":"2009","journal-title":"J. R. Stat. Soc. B"},{"key":"2023020202314601100_btv378-B31","first-page":"339","article-title":"Split Hamiltonian Monte Carlo","volume-title":"Statistics and Computing","author":"Shahbaba","year":"2013"},{"key":"2023020202314601100_btv378-B32","doi-asserted-by":"crossref","first-page":"2298","DOI":"10.1093\/oxfordjournals.molbev.a003776","article-title":"Exploring the demographic history of DNA sequences using the generalized skyline plot","volume":"18","author":"Strimmer","year":"2001","journal-title":"Mol. Biol. Evol."}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/31\/20\/3282\/49035918\/bioinformatics_31_20_3282.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/31\/20\/3282\/49035918\/bioinformatics_31_20_3282.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,2,2]],"date-time":"2023-02-02T03:50:42Z","timestamp":1675309842000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/31\/20\/3282\/196012"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2015,6,20]]},"references-count":32,"journal-issue":{"issue":"20","published-print":{"date-parts":[[2015,10,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btv378","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2015,10,15]]},"published":{"date-parts":[[2015,6,20]]}}}