Efficient parameter sensitivity computation for spatially extended reaction networks.
Reaction-diffusion models are widely used to study spatially extended chemical reaction systems. In order to understand how the dynamics of a reaction-diffusion model are affected by changes in its input parameters, efficient methods for computing parametric sensitivities are required. In this work, we focus on the stochastic models of spatially extended chemical reaction systems that involve partitioning the computational domain into voxels. Parametric sensitivities are often calculated using Monte Carlo techniques that are typically computationally expensive; however, variance reduction techniques can decrease the number of Monte Carlo simulations required. By exploiting the characteristic dynamics of spatially extended reaction networks, we are able to adapt existing finite difference schemes to robustly estimate parametric sensitivities in a spatially extended network. We show that algorithmic performance depends on the dynamics of the given network and the choice of summary statistics. We then describe a hybrid technique that dynamically chooses the most appropriate simulation method for the network of interest. Our method is tested for functionality and accuracy in a range of different scenarios.
Mathematical Institute, University of Oxford, Woodstock Road, Oxford OX2 6GG, United Kingdom.
C
C
Lester
Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, United Kingdom.
C A
0000000304617297
CA
Yates
Mathematical Institute, University of Oxford, Woodstock Road, Oxford OX2 6GG, United Kingdom.
R E
0000000263049333
RE
Baker