{"title":"Sampling of probability measures in the convex order by Wasserstein projection","authors":"A. Alfonsi, Jacopo Corbetta, B. Jourdain","doi":"10.1214/19-AIHP1014","DOIUrl":null,"url":null,"abstract":"In this paper, for $\\mu$ and $\\nu$ two probability measures on $\\mathbb{R}^d$ with finite moments of order $\\rho\\ge 1$, we define the respective projections for the $W_\\rho$-Wasserstein distance of $\\mu$ and $\\nu$ on the sets of probability measures dominated by $\\nu$ and of probability measures larger than $\\mu$ in the convex order. The $W_2$-projection of $\\mu$ can be easily computed when $\\mu$ and $\\nu$ have finite support by solving a quadratic optimization problem with linear constraints. In dimension $d=1$, Gozlan et al.~(2018) have shown that the projections do not depend on $\\rho$. We explicit their quantile functions in terms of those of $\\mu$ and $\\nu$. The motivation is the design of sampling techniques preserving the convex order in order to approximate Martingale Optimal Transport problems by using linear programming solvers. We prove convergence of the Wasserstein projection based sampling methods as the sample sizes tend to infinity and illustrate them by numerical experiments.","PeriodicalId":7902,"journal":{"name":"Annales De L Institut Henri Poincare-probabilites Et Statistiques","volume":"62 1","pages":"1706-1729"},"PeriodicalIF":1.2000,"publicationDate":"2017-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"30","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Annales De L Institut Henri Poincare-probabilites Et Statistiques","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1214/19-AIHP1014","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"STATISTICS & PROBABILITY","Score":null,"Total":0}
引用次数: 30
Abstract
In this paper, for $\mu$ and $\nu$ two probability measures on $\mathbb{R}^d$ with finite moments of order $\rho\ge 1$, we define the respective projections for the $W_\rho$-Wasserstein distance of $\mu$ and $\nu$ on the sets of probability measures dominated by $\nu$ and of probability measures larger than $\mu$ in the convex order. The $W_2$-projection of $\mu$ can be easily computed when $\mu$ and $\nu$ have finite support by solving a quadratic optimization problem with linear constraints. In dimension $d=1$, Gozlan et al.~(2018) have shown that the projections do not depend on $\rho$. We explicit their quantile functions in terms of those of $\mu$ and $\nu$. The motivation is the design of sampling techniques preserving the convex order in order to approximate Martingale Optimal Transport problems by using linear programming solvers. We prove convergence of the Wasserstein projection based sampling methods as the sample sizes tend to infinity and illustrate them by numerical experiments.
期刊介绍:
The Probability and Statistics section of the Annales de l’Institut Henri Poincaré is an international journal which publishes high quality research papers. The journal deals with all aspects of modern probability theory and mathematical statistics, as well as with their applications.