Abstract:In this work we develop an effective Monte Carlo method for estimating sensitivities, or gradients of expectations of sufficiently smooth functionals, of a reflected diffusion in a convex polyhedral domain with respect to its defining parameters --- namely, its initial condition, drift and diffusion coefficients, and directions of reflection. Our method, which falls into the class of infinitesimal perturbation analysis (IPA) methods, uses a probabilistic representation for such sensitivities as the expectation of a functional of the reflected diffusion and its associated derivative process. The latter process is the unique solution to a constrained linear stochastic differential equation with jumps whose coefficients, domain and directions of reflection are modulated by the reflected diffusion. We propose an asymptotically unbiased estimator for such sensitivities using an Euler approximation of the reflected diffusion and its associated derivative process. Proving that the Euler approximation converges is challenging because the derivative process jumps whenever the reflected diffusion hits the boundary (of the domain). A key step in the proof is establishing a continuity property of the related derivative map, which is of independent interest. We compare the performance of our IPA estimator to a standard likelihood ratio estimator (whenever the latter is applicable), and provide numerical evidence that the variance of the former is substantially smaller than that of the latter. We illustrate our method with an example of a rank-based interacting diffusion model of equity markets. Interestingly, we show that estimating certain sensitivities of the rank-based interacting diffusion model using our method for a reflected Brownian motion description of the model outperforms a finite difference method for a stochastic differential equation description of the model.

Pathwise Estimation of Probability Sensitivities Through Terminating or Steady-State Simulations

Conditional Monte Carlo Estimation of Quantile Sensitivities

Estimating Sensitivities of Portfolio Credit Risk Using Monte Carlo

Conditional expectation method for estimating parameter sensitivity of unreliability

Estimating Quantile Sensitivities

Path weighting sensitivities

A new efficient simulation method based on Bayes' theorem and importance sampling Markov chain simulation to estimate the failure-probability-based global sensitivity measure

A probabilistic simulation method for sensitivity analysis of input epistemic uncertainties on failure probability

An efficient method based on Bayes’ theorem to estimate the failure-probability-based sensitivity measure

Meaningful Sensitivities: A New Family of Simulation Sensitivity Measures

On the Asymptotic Analysis of Quantile Sensitivity Estimation by Monte Carlo Simulation.

A Monte Carlo method for estimating sensitivities of reflected diffusions in convex polyhedral domains

Estimation of the Hurst Parameter in the Simultaneous Presence of Jumps and Noise

Uncertainty quantification for generalized Langevin dynamics

Simulating Sensitivities of Conditional Value at Risk

Nonparametric Inference for Sensitivity of Haezendonck–Goovaerts Risk Measure

Efficient metamodel-based importance sampling coupled with single-loop estimation method for parameter global reliability sensitivity analysis

Differential Quantile-Based Sensitivity in Discontinuous Models

Strong convergence of path sensitivities

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

A general framework for probabilistic sensitivity analysis with respect to distribution parameters