Geometrically Convergent Simulation of the Extrema of Lévy Processes

by   Jorge González Cázares, et al.

We develop a novel Monte Carlo algorithm for the simulation from the joint law of the position, the running supremum and the time of the supremum of a general Lévy process at an arbitrary finite time. We prove that the bias decays geometrically, in contrast to the power law for the random walk approximation (RWA). We identify the law of the error and, inspired by the recent work of Ivanovs [Iva18] on RWA, characterise its asymptotic behaviour. We establish a central limit theorem, construct non-asymptotic and asymptotic confidence intervals and prove that the multilevel Monte Carlo (MLMC) estimator has optimal computational complexity (i.e. of order ϵ^-2 if the L^2-norm of the error is at most ϵ) for locally Lipschitz and barrier-type functionals of the triplet. If the increments of the Lévy process cannot be sampled directly, we combine our algorithm with the Asmussen-Rosiński approximation [AR01] by choosing the rate of decay of the cutoff level for small jumps so that the corresponding MC and MLMC estimators have minimal computational complexity. Moreover, we give an unbiased version of our estimator using ideas from Rhee-Glynn [RG15] and Vihola [Vih18].



page 1

page 2

page 3

page 4


Monte Carlo algorithm for the extrema of tempered stable processes

We develop a novel Monte Carlo algorithm for the vector consisting of th...

Unbiased Multilevel Monte Carlo methods for intractable distributions: MLMC meets MCMC

Constructing unbiased estimators from Markov chain Monte Carlo (MCMC) ou...

A Note on Monte Carlo Integration in High Dimensions

Monte Carlo integration is a commonly used technique to compute intracta...

A Scrambled Method of Moments

Quasi-Monte Carlo (qMC) methods are a powerful alternative to classical ...

The Langevin Monte Carlo algorithm in the non-smooth log-concave case

We prove non-asymptotic polynomial bounds on the convergence of the Lang...

Full version: An evaluation of estimation techniques for probabilistic reachability

We evaluate numerically-precise Monte Carlo (MC), Quasi-Monte Carlo (QMC...

Edgeworth corrections for spot volatility estimator

We develop Edgeworth expansion theory for spot volatility estimator unde...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.