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.

Unbiased estimation of parameter sensitivities for stochastic chemical reaction networks

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

Variance-based sensitivity analysis for stochastic chemical kinetics

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

Sensitivity analysis for multiscale stochastic reaction networks using hybrid approximations

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

Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks

Efficient Finite-Difference Estimation of Second-Order Parametric Sensitivities for Stochastic Discrete Biochemical Systems

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

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

Steady-state parameter sensitivity in stochastic modeling via trajectory reweighting

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

A New Framework of Global Sensitivity Analysis for the Chemical Kinetic Model Using PSO-BPNN.

Parametric Sensitivity Analysis for Models of Reaction Networks within Interacting Compartments

CERENA: ChEmical REaction Network Analyzer—A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics

Bayesian Parameter Estimation for Biochemical Reaction Networks Using Region-Based Adaptive Parallel Tempering

Transitory Sensitivity in Automatic Chemical Kinetic Mechanism Analysis

Global parameter identification of stochastic reaction networks from single trajectories

Sensitivity of chemical reaction networks: a structural approach 3. Regular multimolecular systems

Cramer-Rao bound and absolute sensitivity in chemical reaction networks

Efficient Sensitivity Analysis for Parametric Robust Markov Chains