Abstract:The Wasserstein barycenter (WB) is an important tool for summarizing sets of probabilities. It finds applications in applied probability, clustering, image processing, etc. When the probability supports are finite and fixed, the problem of computing a WB is formulated as a linear optimization problem whose dimensions generally exceed standard solvers' capabilities. For this reason, the WB problem is often replaced with a simpler nonlinear optimization model constructed via an entropic regularization function so that specialized algorithms can be employed to compute an approximate WB efficiently. Contrary to such a widespread inexact scheme, we propose an exact approach based on the Douglas-Rachford splitting method applied directly to the WB linear optimization problem for applications requiring accurate WB. Our algorithm, which has the interesting interpretation of being built upon averaging marginals, operates series of simple (and exact) projections that can be parallelized and even randomized, making it suitable for large-scale datasets. As a result, our method achieves good performance in terms of speed while still attaining accuracy. Furthermore, the same algorithm can be applied to compute generalized barycenters of sets of measures with different total masses by allowing for mass creation and destruction upon setting an additional parameter. Our contribution to the field lies in the development of an exact and efficient algorithm for computing barycenters, enabling its wider use in practical applications. The approach's mathematical properties are examined, and the method is benchmarked against the state-of-the-art methods on several data sets from the literature.

Multiplier bootstrap for Bures-Wasserstein barycenters

Distributed Bootstrap Simultaneous Inference for High-Dimensional Quantile Regression

Multiplier subsample bootstrap for statistics of time series

Cluster-Robust Bootstrap Inference in Quantile Regression Models

High-dimensional bootstrap and asymptotic expansion

Continuous Wasserstein-2 Barycenter Estimation without Minimax Optimization

Stochastic gradient descent for barycenters in Wasserstein space

Bootstrap quantification of estimation uncertainties in network degree distributions

Simultaneous likelihood-based bootstrap confidence sets for a large number of models

Stratified sampling and bootstrapping for approximate Bayesian computation

Bootstrapping Estimators based on the Block Maxima Method

Gaussian Approximation and Multiplier Bootstrap for Polyak-Ruppert Averaged Linear Stochastic Approximation with Applications to TD Learning

Computing Wasserstein Barycenter via operator splitting: the method of averaged marginals

Confidence surfaces for the mean of locally stationary functional time series

The Big Data Bootstrap

High-dimensional Simultaneous Inference on Non-Gaussian VAR Model via De-biased Estimator

On the validity of bootstrap uncertainty estimates in the Mallows-Binomial model

A Subsampled Double Bootstrap for Massive Data