Variance decomposition in stochastic simulators

O. P. Le Maître, O. M. Knio, Alvaro Moraes

Research output: Contribution to journalArticlepeer-review

5 Scopus citations


This work aims at the development of a mathematical and computational approach that enables quantification of the inherent sources of stochasticity and of the corresponding sensitivities in stochastic simulations of chemical reaction networks. The approach is based on reformulating the system dynamics as being generated by independent standardized Poisson processes. This reformulation affords a straightforward identification of individual realizations for the stochastic dynamics of each reaction channel, and consequently a quantitative characterization of the inherent sources of stochasticity in the system. By relying on the Sobol-Hoeffding decomposition, the reformulation enables us to perform an orthogonal decomposition of the solution variance. Thus, by judiciously exploiting the inherent stochasticity of the system, one is able to quantify the variance-based sensitivities associated with individual reaction channels, as well as the importance of channel interactions. Implementation of the algorithms is illustrated in light of simulations of simplified systems, including the birth-death, Schlögl, and Michaelis-Menten models.
Original languageEnglish (US)
Pages (from-to)244115
JournalThe Journal of Chemical Physics
Issue number24
StatePublished - Jun 30 2015


Dive into the research topics of 'Variance decomposition in stochastic simulators'. Together they form a unique fingerprint.

Cite this