Stochastic algorithmic differentiation of (expectations of) discontinuous functions (indicator functions)

From MaRDI portal
Publication:5063446

DOI10.1080/00207160.2021.1883593zbMATH Open1499.65005arXiv1811.05741OpenAlexW3128652807MaRDI QIDQ5063446FDOQ5063446


Authors: Christian Fries Edit this on Wikidata


Publication date: 21 March 2022

Published in: International Journal of Computer Mathematics (Search for Journal in Brave)

Abstract: In this paper, we present a method for the accurate estimation of the derivative (aka.~sensitivity) of expectations of functions involving an indicator function by combining a stochastic algorithmic differentiation and a regression. The method is an improvement of the approach presented in [Risk Magazine April 2018]. The finite difference approximation of a partial derivative of a Monte-Carlo integral of a discontinuous function is known to exhibit a high Monte-Carlo error. The issue is evident since the Monte-Carlo approximation of a discontinuous function is just a finite sum of discontinuous functions and as such, not even differentiable. The algorithmic differentiation of a discontinuous function is problematic. A natural approach is to replace the discontinuity by continuous functions. This is equivalent to replacing a path-wise automatic differentiation by a (local) finite difference approximation. We present an improvement (in terms of variance reduction) by decoupling the integration of the Dirac delta and the remaining conditional expectation and estimating the two parts by separate regressions. For the algorithmic differentiation, we derive an operator that can be injected seamlessly - with minimal code changes - into the algorithm resulting in the exact result.


Full work available at URL: https://arxiv.org/abs/1811.05741




Recommendations




Cites Work


Cited In (1)





This page was built for publication: Stochastic algorithmic differentiation of (expectations of) discontinuous functions (indicator functions)

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5063446)