Multi-index Importance Sampling for McKean-Vlasov Stochastic Differential Equation

07/11/2023
βˆ™
by   Nadhir Ben Rached, et al.
βˆ™
0
βˆ™

This work introduces a novel approach that combines the multi-index Monte Carlo (MC) method with importance sampling (IS) to estimate rare event quantities expressed as an expectation of a smooth observable of solutions to a broad class of McKean-Vlasov stochastic differential equations. We extend the double loop Monte Carlo (DLMC) estimator, previously introduced in our works (Ben Rached et al., 2022a,b), to the multi-index setting. We formulate a new multi-index DLMC estimator and conduct a comprehensive cost-error analysis, leading to improved complexity results. To address rare events, an importance sampling scheme is applied using stochastic optimal control of the single level DLMC estimator. This combination of IS and multi-index DLMC not only reduces computational complexity by two orders but also significantly decreases the associated constant compared to vanilla MC. The effectiveness of the proposed multi-index DLMC estimator is demonstrated using the Kuramoto model from statistical physics. The results confirm a reduced complexity from π’ͺ(TOL_r^-4) for the single level DLMC estimator (Ben Rached et al., 2022a) to π’ͺ(TOL_r^-2 (logTOL_r^-1)^2) for the considered example, while ensuring accurate estimation of rare event quantities within the prescribed relative error tolerance TOL_r.

READ FULL TEXT

page 26

page 28

page 29

research
βˆ™ 08/05/2022

Multilevel Importance Sampling for McKean-Vlasov Stochastic Differential Equation

This work combines multilevel Monte Carlo methods with importance sampli...
research
βˆ™ 07/14/2022

Single Level Importance Sampling for McKean-Vlasov Stochastic Differential Equation

This paper investigates Monte Carlo methods to estimate probabilities of...
research
βˆ™ 06/05/2023

Automated Importance Sampling via Optimal Control for Stochastic Reaction Networks: A Markovian Projection-based Approach

We propose a novel alternative approach to our previous work (Ben Hammou...
research
βˆ™ 10/27/2021

Efficient Importance Sampling via Stochastic Optimal Control for Stochastic Reaction Networks

We explore the efficient estimation of statistical quantities, particula...
research
βˆ™ 11/15/2021

Single-Index Importance Sampling with Stratification

In many stochastic problems, the output of interest depends on an input ...
research
βˆ™ 11/21/2021

Stochastic viscosity approximations of Hamilton-Jacobi equations and variance reduction

We consider the computation of free energy-like quantities for diffusion...
research
βˆ™ 03/31/2020

Moment-Preserving and Mesh-Adaptive Reweighting Method for Rare-Event Sampling in Monte-Carlo Algorithms

We present novel roulette schemes for rare-event sampling that are both ...

Please sign up or login with your details

Forgot password? Click here to reset