In this paper, a regularization of Wasserstein barycenters for random measures supported on @math is introduced via convex penalization. The existence and uniqueness of such barycenters is first proved for a large class of penalization functions. The Bregman divergence associated to the penalization term is then considered to obtain a stability result on penalized barycenters. This allows the comparison of data made of @math absolutely continuous probability measures, within the more realistic setting where one only has access to a dataset of random variables sampled from unknown distributions. The convergence of the penalized empirical barycenter of a set of @math iid random probability measures towards its population counterpart is finally analyzed. This approach is shown to be appropriate for the statistical analysis of either discrete or absolutely continuous random measures. It also allows to construct, from a set of discrete measures, consistent estimators of population Wasserstein barycenters that are absolutely continuous.
A detailed characterization of empirical Wasserstein barycenters in terms of existence, uniqueness and regularity for probability measures with support in @math is given in the seminal paper @cite_26 . The relation of such barycenters with the solution of the multi-marginal problem is also studied in both @cite_26 and @cite_31 . The notion of Wasserstein barycenter has been first generalized in @cite_16 , by establishing existence, uniqueness and consistency for random probability measures supported on a locally compact geodesic space. The more general case of probability measures supported on Riemannian manifolds has then been studied in @cite_24 . Subsequently, trimmed barycenters in the Wasserstein space @math have been introduced in @cite_20 for the purpose of combining informations from different experimental units in a parallelized or distributed estimation setting. The framework of optimal transport has been recently adapted for nonnegative measures supported on a compact of @math with different masses, independently by @cite_12 and @cite_33 . However, in all these papers, incorporating regularization into the computation of Wasserstein barycenters has not been considered, which is of interest when the data are irregular probability measures.
@cite_26 @cite_33 @cite_24 @cite_31 @cite_16 @cite_12 @cite_20
1479849574 2298250094 2949090854 2963587665
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. Wasserstein barycenters correspond to optimal solutions of transportation problems for several marginals, and as such have a wide range of applications ranging from economics to statistics and computer science. When the marginal probability measures are absolutely continuous (or vanish on small sets) the theory of Wasserstein barycenters is well-developed [see the seminal paper (Agueh and Carlier in SIAM J Math Anal 43(2):904–924, 2011)]. However, exact continuous computation of Wasserstein barycenters in this setting is tractable in only a small number of specialized cases. Moreover, in many applications data is given as a set of probability measures with finite support. In this paper, we develop theoretical results for Wasserstein barycenters in this discrete setting. Our results rely heavily on polyhedral theory which is possible due to the discrete structure of the marginals. The results closely mirror those in the continuous case with a few exceptions. In this discrete setting we establish that Wasserstein barycenters must also be discrete measures and there is always a barycenter which is provably sparse. Moreover, for each Wasserstein barycenter there exists a non-mass-splitting optimal transport to each of the discrete marginals. Such non-mass-splitting transports do not generally exist between two discrete measures unless special mass balance conditions hold. This makes Wasserstein barycenters in this discrete setting special in this regard. We illustrate the results of our discrete barycenter theory with a proof-of-concept computation for a hypothetical transportation problem with multiple marginals: distributing a fixed set of goods when the demand can take on different distributional shapes characterized by the discrete marginal distributions. A Wasserstein barycenter, in this case, represents an optimal distribution of inventory facilities which minimize the squared distance transportation cost totaled over all demands. 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. We study barycenters in the space of probability measures on a Riemannian manifold, equipped with the Wasserstein metric. Under reasonable assumptions, we establish absolute continuity of the barycenter of general measures Ω∈P(P(M))Ω∈P(P(M)) on Wasserstein space, extending on one hand, results in the Euclidean case (for barycenters between finitely many measures) of Agueh and Carlier [1] to the Riemannian setting, and on the other hand, results in the Riemannian case of Cordero-Erausquin, McCann, Schmuckenschlager [12] for barycenters between two measures to the multi-marginal setting. Our work also extends these results to the case where Ω is not finitely supported. As applications, we prove versions of Jensen's inequality on Wasserstein space and a generalized Brunn–Minkowski inequality for a random measurable set on a Riemannian manifold.