A finite difference method for estimating second order parameter sensitivities of discrete stochastic chemical reaction networks
Abstract
We present an efficient finite difference method for the approximation of second derivatives, with respect to system parameters, of expectations for a class of discrete stochastic chemical reaction networks. The method uses a coupling of the perturbed processes that yields a much lower variance than existing methods, thereby drastically lowering the computational complexity required to solve a given problem. Further, the method is simple to implement and will also prove useful in any setting in which continuous time Markov chains are used to model dynamics, such as population processes. We expect the new method to be useful in the context of optimization algorithms that require knowledge of the Hessian.
- Publication:
-
arXiv e-prints
- Pub Date:
- August 2012
- DOI:
- 10.48550/arXiv.1208.0843
- arXiv:
- arXiv:1208.0843
- Bibcode:
- 2012arXiv1208.0843S
- Keywords:
-
- Quantitative Biology - Quantitative Methods;
- Mathematics - Numerical Analysis;
- 60H35;
- 65C99;
- 92C40
- E-Print:
- New format (two columns). 14 pages, 9 figures, 7 tables