EN FR
EN FR


Section: New Results

Approximation of Martingale Optimal Transport problems

With J. Corbetta, A. Alfonsi and B. Jourdain study sampling methods preserving the convex order for two probability measures μ and ν on 𝐑d, with ν dominating μ. When (Xi)1iI (resp. (Yj)1jJ) are independent and identically distributed according μ (resp. ν), in general μI=1Ii=1IδXi and νJ=1Jj=1JδYj are not rankable for the convex order. They investigate modifications of μI (resp. νJ) smaller than νJ (resp. greater than μI) in the convex order and weakly converging to μ (resp. ν) as I,J. They first consider the one dimensional case d=1, where, according to Kertz and Rösler, the set of probability measures with a finite first order moment is a lattice for the increasing and the decreasing convex orders. Given μ and ν in this set, they define μν (resp. μν) as the supremum (resp. infimum) of μ and ν for the decreasing convex order when 𝐑xμ(dx)𝐑xν(dx) and for the increasing convex order otherwise. This way, μν (resp. μν) is greater than μ (resp. smaller than ν) in the convex order. They give efficient algorithms permitting to compute μν and μν (and therefore μIνJ and μIνJ) when μ and ν are convex combinations of Dirac masses. In general dimension, when μ and ν have finite moments of order ρ1, they define the projection μρν (resp. μρν) of μ (resp. ν) on the set of probability measures dominated by ν (resp. larger than μ) in the convex order for the Wasserstein distance with index ρ. When ρ=2, μI2νJ can be computed efficiently by solving a quadratic optimization problem with linear constraints. It turns out that, in dimension d=1, the projections do not depend on ρ and their quantile functions are explicit in terms of those of μ and ν, which leads to efficient algorithms for convex combinations of Dirac masses. Last, they illustrate by numerical experiments the resulting sampling methods that preserve the convex order and their application to approximate Martingale Optimal Transport problems.

With V. Ehrlacher, D. Lombardi and R. Coyaud, A. Alfonsi has started to develop and analyze numerical methods to approximate the optimal transport between two probability measures.