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.

Sensitivities via Rough Paths

Pathwise Itô Calculus for Rough Paths and Rough PDEs with Path Dependent Coefficients

Rough volatility, path-dependent PDEs and weak rates of convergence

Strong convergence of path sensitivities

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

Rough Path Renormalization from Stratonovich to Itô for Fractional Brownian Motion

On the rough-paths approach to non-commutative stochastic calculus

On small time asymptotics for rough differential equations driven by fractional Brownian motions

Error calculus and path sensitivity in financial models

Lipschitz continuity in the Hurst parameter of functionals of stochastic differential equations driven by a fractional Brownian motion

Pathwise convergence of the Euler scheme for rough and stochastic differential equations

Itô stochastic differentials

Rough path theory

Path weighting sensitivities

Stochastic analysis, rough path analysis and fractional Brownian motions

Regularization by noise for rough differential equations driven by Gaussian rough paths

Rough Differential Equations Containing Path-Dependent Bounded Variation Terms

An Introduction to (Stochastic) Calculus with Respect to Fractional Brownian Motion

Rough paths and symmetric-Stratonovich integrals driven by singular covariance gaussian processes

Variability and the existence of rough integrals with irregular coefficients