Adaptive Multilevel Stochastic Approximation of the Value-at-Risk

Abstract

Crépey, Frikha and Louzi (2023) introduced a multilevel stochastic approximation scheme to compute the value-at-risk of a financial loss that is only simulatable by Monte Carlo. The optimal complexity of the scheme is in epsilon^(-5/2), epsilon>0 being a prescribed accuracy, which is suboptimal when compared to the canonical multilevel Monte Carlo performance. This suboptimality stems from the discontinuity of the Heaviside function involved in the biased stochastic gradient that is recursively evaluated to derive the value-at-risk. To mitigate this issue, this paper proposes and analyzes a multilevel stochastic approximation algorithm that adaptively selects the number of inner samples at each level, and proves that its optimal complexity is in epsilon^{-2} |ln{epsilon}|^(5/2). Our theoretical analysis is exemplified through numerical experiments.