Convergence of Probability Measures. By P. Billingsley. Chichester, Sussex, Wiley, 1968. xii, 253 p. 9 1/4“. 117s.

Convergence of Probability Measures

The author's preface gives an outline: "This book is about weakconvergence methods in metric spaces, with applications sufficient to show their power and utility. The Introduction motivates the definitions and indicates how the theory will yield solutions to problems arising outside it. Chapter 1 sets out the basic general theorems, which are then specialized in Chapter 2 to the space C[0, l ] of continuous functions on the unit interval and in Chapter 3 to the space D [0, 1 ] of functions with discontinuities of the first kind. The results of the first three chapters are used in Chapter 4 to derive a variety of limit theorems for dependent sequences of random variables. " The book develops and expands on Donsker's 1951 and 1952 papers on the invariance principle and empirical distributions. The basic random variables remain real-valued although, of course, measures on C[0, l ] and D[0, l ] are vitally used. Within this framework, there are various possibilities for a different and apparently better treatment of the material. More of the general theory of weak convergence of probabilities on separable metric spaces would be useful. Metrizability of the convergence is not brought up until late in the Appendix. The close relation of the Prokhorov metric and a metric for convergence in probability is (hence) not mentioned (see V. Strassen, Ann. Math. Statist. 36 (1965), 423-439; the reviewer, ibid. 39 (1968), 1563-1572). This relation would illuminate and organize such results as Theorems 4.1, 4.2 and 4.4 which give isolated, ad hoc connections between weak convergence of measures and nearness in probability. In the middle of p. 16, it should be noted that C*(S) consists of signed measures which need only be finitely additive if 5 is not compact. On p. 239, where the author twice speaks of separable subsets having nonmeasurable cardinal, he means "discrete" rather than "separable." Theorem 1.4 is Ulam's theorem that a Borel probability on a complete separable metric space is tight. Theorem 1 of Appendix 3 weakens completeness to topological completeness. After mentioning that probabilities on the rationals are tight, the author says it is an

Convergence of probability measures

A criterion is given for the convergence of numerical solutions of the Navier-Stokes equations in two dimensions under steady conditions. The criterion applies to all cases, of steady viscous flow in two dimensions and shows that if the local ' mesh Reynolds number ', based on the size of the mesh used in the solution, exceeds a certain fixed value, the numerical solution will not converge.

On the Navier-Stokes equations

Stochastic equations in infinite dimensions

The flow of homogeneous fluids through porous media: analogies with other physical problems

In this article, we give a brief review of some recent results concerning the study of the Euler-Maruyama scheme and its high-order extensions. These numerical schemes are used to approximate solutions of stochastic differential equations, which enables to approximate various important quantities including solutions of partial differential equations. Some have been implemented in Premia [56]. In this article we mainly consider results about weak approximation, the most important for financial applications.

A Review of Recent Results on Approximation of Solutions of Stochastic Differential Equations

This paper investigates the use of stratified sampling as a variance reduction technique for approximating integrals over large dimensional spaces. The accuracy of this method critically depends on the choice of the space partition, the strata, which should be ideally fitted to thesubsets where the functions to integrate is nearly constant, and on the allocation of the number of samples within each strata. When the dimension is large and the function to integrate is complex, finding such partitions and allocating the sample is a highly non-trivial problem. In this work, we investigate a novel method to improve the efficiency of the estimator "on the fly", by jointly sampling and adapting the strata and the allocation within the strata. The accuracy of estimators when this method is used is examined in detail, in the so-called asymptotic regime (i.e. when both the number of samples and the number of strata are large). We illustrate the use of the method for the computation of the price of path-dependent options in models with both constant and stochastic volatility. The use of this adaptive technique yields variance reduction by factors sometimes larger than 1000 compared to classical Monte Carlo estimators.

On adaptive stratification

This paper investigates the use of stratified sampling as a variance reduction technique for approximating integrals over large dimensional spaces. The accuracy of this method critically depends on the choice of the space partition, the strata, which should be ideally fitted to the subsets where the functions to integrate is nearly constant, and on the allocation of the number of samples within each strata. When the dimension is large and the function to integrate is complex, finding such partitions and allocating the sample is a highly non-trivial problem. In this work, we investigate a novel method to improve the efficiency of the estimator “on the fly”, by jointly sampling and adapting the strata which are hyperrectangles and the allocation within the strata. The accuracy of estimators when this method is used is examined in detail, in the so-called asymptotic regime (i.e. when both the number of samples and the number of strata are large). It turns out that the limiting variance depends on the directions defining the hyperrectangles but not on the precise abscissa of their boundaries along these directions, which gives a mathematical justification to the common choice of equiprobable strata. So, only the directions are adaptively modified by our algorithm. We illustrate the use of the method for the computation of the price of path-dependent options in models with both constant and stochastic volatility. The use of this adaptive technique yields variance reduction by factors sometimes larger than 1000 compared to classical Monte Carlo estimators.

https://hal.archives-ouvertes.fr/hal-00319157/document

Motivated by the approximation of Martingale Optimal Transport problems, we study sampling methods preserving the convex order for two probability measures $\mu$ and $\nu$ on $\mathbb{R}^d$, with $\nu$ dominating $\mu$. When $(X_i)_{1\le i\le I}$ (resp. $(Y_j)_{1\le j\le J}$) are i.i.d. according $\mu$ (resp. $\nu$), the empirical measures $\mu_I$ and $\nu_J$ are not in the convex order. We investigate modifications of $\mu_I$ (resp. $\nu_J$) smaller than $\nu_J$ (resp. greater than $\mu_I$) in the convex order and weakly converging to $\mu$ (resp. $\nu$) as $I,J\to\infty$. In dimension 1, according to Kertz and R\"osler (1992), the set of probability measures with a finite first order moment is a lattice for the increasing and the decreasing convex orders. From this result, we can define $\mu\vee\nu$ (resp. $\mu\wedge\nu$) that is greater than $\mu$ (resp. smaller than $\nu$) in the convex order. We give efficient algorithms permitting to compute $\mu\vee\nu$ and $\mu\wedge\nu$ when $\mu$ and $\nu$ are convex combinations of Dirac masses. In general dimension, when $\mu$ and $\nu$ have finite moments of order $\rho\ge 1$, we define the projection $\mu\curlywedge_\rho \nu$ (resp. $\mu\curlyvee_\rho\nu$) of $\mu$ (resp. $\nu$) on the set of probability measures dominated by $\nu$ (resp. larger than $\mu$) in the convex order for the Wasserstein distance with index $\rho$. When $\rho=2$, $\mu_I\curlywedge_2 \nu_J$ can be computed efficiently by solving a quadratic optimization problem with linear constraints. It turns out that, in dimension 1, the projections do not depend on $\rho$ and their quantile functions are explicit, which leads to efficient algorithms for convex combinations of Dirac masses. Last, we illustrate by numerical experiments the resulting sampling methods that preserve the convex order and their application to approximate Martingale Optimal Transport problems.

Sampling of probability measures in the convex order and approximation of Martingale Optimal Transport problems

In the particular case of a concave flux function, we are interested in the long time behaviour of the nonlinear process associated to the one-dimensional viscous scalar conservation law. We also consider the particle system obtained by remplacing the cumulative distribution function in the drift coefficient of this nonlinear process by the empirical cdf. We first obtain trajectorial propagation of chaos result. Then, Poincare inequalities are used to get explicit estimates concerning the long time behaviour of both the nonlinear process and the particle system.

Benjamin Jourdain

Papers

A Review of Recent Results on Approximation of Solutions of Stochastic Differential Equations

On adaptive stratification

On adaptive stratification

Sampling of probability measures in the convex order and approximation of Martingale Optimal Transport problems

Propagation of chaos and Poincar\'e inequalities for a system of particles interacting through their cdf