English (United Kingdom)

https://curated-unify.zendy.io/wp-json/zendy-region/v1/featured_content/oa?rat=en

https://curated-unify.zendy.io/wp-json/zendy-region/v1/highlighted_journal/

Global: Zendy Tools annual

Presents the keyphrase highlighting as premium feature

Keyphrase Highlighting

Presents the summarisation as premium feature

Summarisation

Insights

Presents the pdf analysis as premium feature

PDF Analysis

Presents the zaia usage as premium feature

ZAIA

Global: Zendy Tools monthly

Global: Zendy Plus monthly

Presents the access of premium content as premium feature

Premium Content

Global: Zendy Plus annual

Global: Zendy Free Plan

Most common Monte Carlo methods for sensitivity analysis of stochastic reaction networks are the finite difference (FD), the Girsanov transformation (GT) and the regularized pathwise derivative (RPD) methods. It has been numerically observed in the literature, that the biased FD and RPD methods tend to have lower variance than the unbiased GT method and that centering the GT method (CGT) reduces its variance. We provide a theoretical justification for these observations in terms of system size asymptotic analysis under what is known as the classical scaling. Our analysis applies to GT, CGT and FD, and shows that the standard deviations of their estimators when normalized by the actual sensitivity, scale as $\mathcal{O}(N^{1/2}), \mathcal{O}(1)$ and $\mathcal{O}(N^{-1/2})$ respectively, as system size $N \to \infty$. In the case of the FD methods, the $N \to \infty$ asymptotics are obtained keeping the finite difference perturbation $h$ fixed. Our numerical examples verify that our order estimates are sharp and that the variance of the RPD method scales similarly to the FD methods. We combine our large $N$ asymptotics with previously known small $h$ asymptotics to obtain the best choice of $h$ in terms of $N$, and estimate the number $N_s$ of simulations required to achieve a prescribed relative $\mathcal{L}_2$ error $\delta$. This shows that $N_s$ depends on $\delta$ and $N$ as $\delta^{-2 - \frac{\gamma_2}{\gamma_1}} N^{-1}, \delta^{-2}$ and $N \delta^{-2}$, for FD, CGT and GT respectively. Here $\gamma_1 >0, \gamma_2>0$ depend on the type of FD method used.

Efficiency of the Girsanov Transformation Approach for Parametric Sensitivity Analysis of Stochastic Chemical Kinetics