title: Nested sampling for parameter inference in systems biology: application to an exemplar circadian model.
abstract: BACKGROUND
  Model selection and parameter inference are complex problems that have yet to be fully addressed in systems biology. In contrast with parameter optimisation, parameter inference computes both the parameter means and their standard deviations (or full posterior distributions), thus yielding important information on the extent to which the data and the model topology constrain the inferred parameter values.
  
  
  RESULTS
  We report on the application of nested sampling, a statistical approach to computing the Bayesian evidence Z, to the inference of parameters, and the estimation of log Z in an established model of circadian rhythms. A ten-fold difference in the coefficient of variation between degradation and transcription parameters is demonstrated. We further show that the uncertainty remaining in the parameter values is reduced by the analysis of increasing numbers of circadian cycles of data, up to 4 cycles, but is unaffected by sampling the data more frequently. Novel algorithms for calculating the likelihood of a model, and a characterisation of the performance of the nested sampling algorithm are also reported. The methods we develop considerably improve the computational efficiency of the likelihood calculation, and of the exploratory step within nested sampling.
  
  
  CONCLUSIONS
  We have demonstrated in an exemplar circadian model that the estimates of posterior parameter densities (as summarised by parameter means and standard deviations) are influenced predominately by the length of the time series, becoming more narrowly constrained as the number of circadian cycles considered increases. We have also shown the utility of the coefficient of variation for discriminating between highly-constrained and less-well constrained parameters.
authors.affiliation: MRC Human Genetics Unit, IGMM, University of Edinburgh, Edinburgh EH4 2XU, UK. S.Aitken@ed.ac.uk
authors.forename: Stuart
authors.initials: S
authors.lastname: Aitken
authors.forename: Ozgur E
authors.initials: OE
authors.lastname: Akman
mesh_terms: D000465:Algorithms; D002940:Circadian Rhythm; D008954:Models, Biological; D009010:Monte Carlo Method; D049490:Systems Biology / Q000379:methods*
references.citation: J Bioinform Comput Biol. 2005 Apr;3(2):415-36
references.pmid: 15852513
references.citation: Bioinformatics. 1998;14(10):869-83
references.pmid: 9927716
references.citation: Curr Opin Biotechnol. 2013 Aug;24(4):767-74
references.pmid: 23578462
references.citation: PLoS Comput Biol. 2011 Oct;7(10):e1002215
references.pmid: 22022255
references.citation: Bioinformatics. 2008 Mar 15;24(6):833-9
references.pmid: 18057018
references.citation: Mol Syst Biol. 2010 Sep 21;6:416
references.pmid: 20865009
references.citation: Mol Syst Biol. 2008;4:164
references.pmid: 18277380
references.citation: Nat Rev Mol Cell Biol. 2010 Nov;11(11):764-76
references.pmid: 20966970
references.citation: Ecol Lett. 2009 Jun;12(6):583-92
references.pmid: 19504722
references.citation: PLoS Comput Biol. 2007 Jan 19;3(1):e7
references.pmid: 17238282
references.citation: Bioinformatics. 2008 Sep 1;24(17):1933-4
references.pmid: 18632751
references.citation: Nat Rev Genet. 2009 Feb;10(2):122-33
references.pmid: 19139763
references.citation: Bioinformatics. 2010 Jan 1;26(1):104-10
references.pmid: 19880371
references.citation: J Biol Rhythms. 1999 Dec;14(6):433-48
references.pmid: 10643740
references.citation: Science. 2005 Jul 22;309(5734):630-3
references.pmid: 16040710
references.citation: Nature. 2002 Nov 14;420(6912):206-10
references.pmid: 12432404
references.citation: Bioinformatics. 2007 Oct 1;23(19):2596-603
references.pmid: 17660527