학술논문

Prediction uncertainty assessment of a systems biology model requires a sample of the full probability distribution of its parameters
Document Type
Academic Journal
Source
PeerJ. June 17, 2014, Vol. 2 e433
Subject
Monte Carlo methods -- Models
Markov processes -- Models
Algorithms -- Models
Probability distributions -- Models
Biological sciences
Algorithm
Models
Language
English
ISSN
2167-8359
Abstract
Multi-parameter models in systems biology are typically 'sloppy': some parameters or combinations of parameters may be hard to estimate from data, whereas others are not. One might expect that parameter uncertainty automatically leads to uncertain predictions, but this is not the case. We illustrate this by showing that the prediction uncertainty of each of six sloppy models varies enormously among different predictions. Statistical approximations of parameter uncertainty may lead to dramatic errors in prediction uncertainty estimation. We argue that prediction uncertainty assessment must therefore be performed on a per-prediction basis using a full computational uncertainty analysis. In practice this is feasible by providing a model with a sample or ensemble representing the distribution of its parameters. Within a Bayesian framework, such a sample may be generated by a Markov Chain Monte Carlo (MCMC) algorithm that infers the parameter distribution based on experimental data. Matlab code for generating the sample (with the Differential Evolution Markov Chain sampler) and the subsequent uncertainty analysis using such a sample, is supplied as Supplemental Information.
Author(s): Simon van Mourik (1,2), Cajo ter Braak (1), Hans Stigter (1), Jaap Molenaar (1,2) Introduction By combining experiments and mathematical model analysis, systems biology tries to unravel the key [...]