Preconditioned Bayesian Regression for Stochastic Chemical Kinetics

  • Authors:
  • Alen Alexanderian;Francesco Rizzi;Muruhan Rathinam;Olivier P. Le Maître;Omar M. Knio

  • Affiliations:
  • Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, USA 78712;Department of Mechanical Engineering and Materials Science, Duke University, Durham, USA 27708;Department of Mathematics and Statistics, University of Maryland, Baltimore County, Baltimore, USA 21250;LIMSI-CNRS, Orsay Cedex, France 91403;Department of Mechanical Engineering and Materials Science, Duke University, Durham, USA 27708

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2014

Quantified Score

Hi-index 0.00

Visualization

Abstract

We develop a preconditioned Bayesian regression method that enables sparse polynomial chaos representations of noisy outputs for stochastic chemical systems with uncertain reaction rates. The approach is based on the definition of an appropriate multiscale transformation of the state variables coupled with a Bayesian regression formalism. This enables efficient and robust recovery of both the transient dynamics and the corresponding noise levels. Implementation of the present approach is illustrated through applications to a stochastic Michaelis---Menten dynamics and a higher dimensional example involving a genetic positive feedback loop. In all cases, a stochastic simulation algorithm (SSA) is used to compute the system dynamics. Numerical experiments show that Bayesian preconditioning algorithms can simultaneously accommodate large noise levels and large variability with uncertain parameters, and that robust estimates can be obtained with a small number of SSA realizations.