title: Bayesian parameter inference and model selection by population annealing in systems biology.
abstract: Parameter inference and model selection are very important for mathematical modeling in systems biology. Bayesian statistics can be used to conduct both parameter inference and model selection. Especially, the framework named approximate Bayesian computation is often used for parameter inference and model selection in systems biology. However, Monte Carlo methods needs to be used to compute Bayesian posterior distributions. In addition, the posterior distributions of parameters are sometimes almost uniform or very similar to their prior distributions. In such cases, it is difficult to choose one specific value of parameter with high credibility as the representative value of the distribution. To overcome the problems, we introduced one of the population Monte Carlo algorithms, population annealing. Although population annealing is usually used in statistical mechanics, we showed that population annealing can be used to compute Bayesian posterior distributions in the approximate Bayesian computation framework. To deal with un-identifiability of the representative values of parameters, we proposed to run the simulations with the parameter ensemble sampled from the posterior distribution, named "posterior parameter ensemble". We showed that population annealing is an efficient and convenient algorithm to generate posterior parameter ensemble. We also showed that the simulations with the posterior parameter ensemble can, not only reproduce the data used for parameter inference, but also capture and predict the data which was not used for parameter inference. Lastly, we introduced the marginal likelihood in the approximate Bayesian computation framework for Bayesian model selection. We showed that population annealing enables us to compute the marginal likelihood in the approximate Bayesian computation framework and conduct model selection depending on the Bayes factor.
authors.affiliation: Department of Biophysics, Division of Biology, Graduate School of Science, Kyoto University, Kyoto, Japan.
authors.forename: Yohei
authors.initials: Y
authors.lastname: Murakami
mesh_terms: D001499:Bayes Theorem; D003198:Computer Simulation; D007328:Insulin / Q000235:genetics; D008962:Models, Theoretical; D009010:Monte Carlo Method; D051057:Proto-Oncogene Proteins c-akt / Q000235:genetics; D015398:Signal Transduction / Q000235:genetics; D049490:Systems Biology
chemical_list: D007328:Insulin; D051057:Proto-Oncogene Proteins c-akt
references.citation: Bioinformatics. 1998;14(10):869-83
references.pmid: 9927716
references.citation: Trends Endocrinol Metab. 2002 Dec;13(10):444-51
references.pmid: 12431841
references.citation: Bioinformatics. 2008 Mar 15;24(6):833-9
references.pmid: 18057018
references.citation: Bioinformatics. 2010 Jan 1;26(1):104-10
references.pmid: 19880371
references.citation: Nat Genet. 2002 May;31(1):64-8
references.pmid: 11967538
references.citation: Proc Natl Acad Sci U S A. 2003 Dec 23;100(26):15324-8
references.pmid: 14663152
references.citation: PLoS Comput Biol. 2013;9(1):e1002803
references.pmid: 23341757
references.citation: Proc Natl Acad Sci U S A. 2003 Oct 14;100(21):11980-5
references.pmid: 14530388
references.citation: Syst Biol. 2006 Apr;55(2):195-207
references.pmid: 16522570
references.citation: Phys Rev E Stat Nonlin Soft Matter Phys. 2013 Apr;87(4):042716
references.pmid: 23679458
references.citation: Proc Natl Acad Sci U S A. 2009 Jul 7;106(27):11107-12
references.pmid: 19549830
references.citation: Cell. 2007 Jun 29;129(7):1261-74
references.pmid: 17604717
references.citation: PLoS One. 2013 Sep 27;8(9):e74178
references.pmid: 24086320
references.citation: Mol Biosyst. 2009 Dec;5(12):1925-35
references.pmid: 19798456
references.citation: Proc Natl Acad Sci U S A. 2007 Feb 6;104(6):1760-5
references.pmid: 17264216
references.citation: Proc Natl Acad Sci U S A. 2011 Aug 2;108(31):12675-80
references.pmid: 21768338
references.citation: Mol Cell. 2012 Jun 29;46(6):820-32
references.pmid: 22633957
references.citation: Proc Natl Acad Sci U S A. 2011 Sep 13;108(37):15190-5
references.pmid: 21876136
references.citation: J Theor Biol. 2005 May 7;234(1):133-43
references.pmid: 15721042
references.citation: Science. 2002 Mar 1;295(5560):1662-4
references.pmid: 11872829
references.citation: J R Soc Interface. 2009 Feb 6;6(31):187-202
references.pmid: 19205079
references.citation: FEBS J. 2009 Feb;276(4):886-902
references.pmid: 19215296
references.citation: Syst Biol (Stevenage). 2005 Dec;152(4):221-8
references.pmid: 16986264
references.citation: Curr Opin Biotechnol. 2013 Aug;24(4):767-74
references.pmid: 23578462
references.citation: BMC Genomics. 2006 May 05;7:108
references.pmid: 16677373
references.citation: Mol Biosyst. 2012 Jul 6;8(7):1921-9
references.pmid: 22555461
references.citation: Genome Res. 2003 Nov;13(11):2467-74
references.pmid: 14559783
references.citation: Mol Biol Evol. 1999 Dec;16(12):1791-8
references.pmid: 10605120
references.citation: Biophys J. 2008 Jul;95(2):540-9
references.pmid: 18456830
references.citation: Nat Rev Genet. 2007 Jun;8(6):450-61
references.pmid: 17510665
references.citation: Mol Syst Biol. 2013;9:644
references.pmid: 23385484
references.citation: Mol Syst Biol. 2013 May 14;9:664
references.pmid: 23670537
references.citation: PLoS One. 2014 Feb 11;9(2):e88419
references.pmid: 24523891
references.citation: Genetics. 2002 Dec;162(4):2025-35
references.pmid: 12524368
references.citation: Nature. 2002 Nov 14;420(6912):206-10
references.pmid: 12432404
references.citation: Phys Rev E Stat Nonlin Soft Matter Phys. 2010 Aug;82(2 Pt 2):026704
references.pmid: 20866937
references.citation: Mol Syst Biol. 2006;2:58
references.pmid: 17102803
references.citation: Science. 2002 Oct 25;298(5594):824-7
references.pmid: 12399590