Abstract:Sensitivity analysis is a process of computing sensitivity indices, which are certain measures of importance of parameters in influencing the outputs of mathematical models. Sensitivity indices computed in variance-based sensitivity analysis yield quantitative answers to questions like how much on average the variance of model output, measuring its uncertainty, decreases, if exact values of certain unknown parameters are determined, e. g. in an experiment. We propose new schemes for estimation of variance-based sensitivity indices of outputs of stochastic models, their conditional expectations, and histograms given the parameters. Unbiased estimators obtained in these schemes can be used in a Monte Carlo (MC) procedure approximating sensitivity indices. We derive relations between variances of final estimators of MC procedures making the same number of evaluations of given function, but using different schemes, both for the newly introduced schemes and for some used before in the literature. Numerical experiment for a discrete state stochastic Markov model of a chemical reaction network (DM) shows that our method can lead to much lower error than method analogous to the one offered by Degasperi et al. Further numerical experiments demonstrate that the application of random time change (RTC) algorithm due to Rathinam et al. for simulation of DM can lead to over 30 times lower variance of estimators of certain sensitivity indices than when Gillespie's direct (GD) method is used, and that this variance may significantly depend on the order of reactions in GD method. We provide some intuitions explaining these effects. We generalize measures used for comparing dispersion of different distributions, such as coefficient of variation and Fano factor to the random parameters case, in a way that they can be computed along with variance-based sensitivity indices.

Sensitivity analysis for multiscale stochastic reaction networks using hybrid approximations

Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks

Reduction for stochastic biochemical reaction networks with multiscale conservations

Unbiased estimation of parameter sensitivities for stochastic chemical reaction networks

Stochastic filtering for multiscale stochastic reaction networks based on hybrid approximations

Stochastic filters based on hybrid approximations of multiscale stochastic reaction networks

Estimation of parameter sensitivities for stochastic reaction networks using tau-leap simulations

Efficiency of the Girsanov transformation approach for parametric sensitivity analysis of stochastic chemical kinetics

Multiscale Analysis of Reaction Networks

Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting

Variance-based sensitivity analysis for stochastic chemical kinetics

Unbiased estimation of second-order parameter sensitivities for stochastic reaction networks

Direct ab initio MD simulation of silver ion diffusion in chalcogenide glasses

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

Error Analysis of Diffusion Approximation Methods for Multiscale Systems in Reaction Kinetics

Adjoint Sensitivity Analysis on Multi-Scale Bioprocess Stochastic Reaction Network

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Efficient Sensitivity Analysis for Parametric Robust Markov Chains

Parametric Sensitivity Analysis for Biochemical Reaction Networks based on Pathwise Information Theory

Parametric Sensitivity Analysis for Models of Reaction Networks within Interacting Compartments