Approximate Algorithms For $k$-Sparse Wasserstein Barycenter With Outliers
Qingyuan Yang, Hu Ding
TL;DR
This work tackles the computation of a $k$-sparse Wasserstein barycenter in the presence of outliers, a problem that blends optimal transport with robust clustering. The authors connect $k$-sparse WB with outliers to $k$-means clustering and develop two main approaches: a clustering-based LP method that yields constant-factor guarantees and a low-dimensional coreset technique that achieves a $(1+ ext{ε})$-approximation when the ambient dimension is manageable. They provide rigorous approximation bounds, including a $(2+ ext{√α})^2$-type guarantee when using an $ ext{α}$-approximate clustering method, and extend the results to doubling metrics via anchor-based coresets. Empirical results on synthetic and real data, including MNIST, demonstrate practical effectiveness and robustness to outliers. Overall, the paper advances robust WB with provable guarantees and scalable techniques, enabling practical use in noisy settings.
Abstract
Wasserstein Barycenter (WB) is one of the most fundamental optimization problems in optimal transportation. Given a set of distributions, the goal of WB is to find a new distribution that minimizes the average Wasserstein distance to them. The problem becomes even harder if we restrict the solution to be ``$k$-sparse''. In this paper, we study the $k$-sparse WB problem in the presence of outliers, which is a more practical setting since real-world data often contains noise. Existing WB algorithms cannot be directly extended to handle the case with outliers, and thus it is urgently needed to develop some novel ideas. First, we investigate the relation between $k$-sparse WB with outliers and the clustering (with outliers) problems. In particular, we propose a clustering based LP method that yields constant approximation factor for the $k$-sparse WB with outliers problem. Further, we utilize the coreset technique to achieve the $(1+ε)$-approximation factor for any $ε>0$, if the dimensionality is not high. Finally, we conduct the experiments for our proposed algorithms and illustrate their efficiencies in practice.
