This paper is focused on the statistical analysis of probability measures @math on @math that can be viewed as independent realizations of an underlying stochastic process. We consider the situation of practical importance where the random measures @math are absolutely continuous with densities @math that are not directly observable. In this case, instead of the densities, we have access to datasets of real random variables @math organized in the form of @math experimental units, such that @math are iid observations sampled from a random measure @math for each @math . In this setting, we focus on first-order statistics methods for estimating, from such data, a meaningful structural mean measure. For the purpose of taking into account phase and amplitude variations in the observations, we argue that the notion of Wasserstein barycenter is a relevant tool. The main contribution of this paper is to characterize the rate of convergence of a (possibly smoothed) empirical Wasserstein barycenter towards its population counterpart in the asymptotic setting where both @math and @math may go to infinity. The optimality of this procedure is discussed from the minimax point of view with respect to the Wasserstein metric. We also highlight the connection between our approach and the curve registration problem in statistics. Some numerical experiments are used to illustrate the results of the paper on the convergence rate of empirical Wasserstein barycenters.
The construction of consistent estimators of a population Wasserstein barycenter for semi-parametric models of random measures can also be found in @cite_10 and @cite_6 , together with a discussion on their connection to the well known curve registration problem in statistics @cite_17 @cite_11 .
@cite_11 @cite_10 @cite_6 @cite_17
1479849574 2949090854 2043083193 2264931217
We define a notion of barycenter for random probability measures in the Wasserstein space. We study the population barycenter in terms of existence and uniqueness. Using a duality argument, we give a precise characterization of the population barycenter for compactly supported measures, and we make a connection between averaging in the Wasserstein space and taking the expectation of optimal transport maps. Then, the problem of estimating this barycenter from @math independent and identically distributed random probability measures is considered. To this end, we study the convergence of the empirical barycenter proposed in Agueh and Carlier (2011) to its population counterpart as the number of measures @math tends to infinity. To illustrate the benefits of this approach for data analysis and statistics, we show the usefulness of averaging in the Wasserstein space for curve and image warping. In this setting, we also study the rate of convergence of the empirical barycenter to its population counterpart for some semi-parametric models of random densities. This paper is concerned by the study of barycenters for random probability measures in the Wasserstein space. Using a duality argument, we give a precise characterization of the population barycenter for various parametric classes of random probability measures with compact support. In particular, we make a connection between averaging in the Wasserstein space as introduced in Agueh and Carlier (2011), and taking the expectation of optimal transport maps with respect to a fixed reference measure. We also discuss the usefulness of this approach in statistics for the analysis of deformable models in signal and image processing. In this setting, the problem of estimating a population barycenter from n independent and identically distributed random probability measures is also considered. In this paper we tackle the problem of comparing distributions of random variables and defining a mean pattern between a sample of random events. Using barycenters of measures in the Wasserstein space, we propose an iterative version as an estimation of the mean distribution. Moreover, when the distributions are a common measure warped by a centered random operator, then the barycenter enables to recover this distribution template. In this paper, based on the Frechet mean, we define a notion of barycenter corresponding to a usual notion of statistical mean. We prove the existence of Wasserstein barycenters of random distributions defined on a geodesic space (E, d). We also prove the consistency of this barycenter in a general setting, that includes taking barycenters of empirical versions of the distributions or of a growing set of distributions.