EMS Coreset: An Efficient Expectation-Maximization Algorithm for Sinkhorn Coreset
Organizations: Department of Statistics, Purdue University West Lafayette, IN 47906, USA
Abstract
Coresets distill large datasets into small, representative subsets for efficient downstream learning. Yet Optimal Transport (OT)-based selection typically requires intensive computation of transport plans, limiting scalability. We introduce a scalable Sinkhorn coreset method that permits closed-form updates of the entropically regularized OT coupling by allowing non-uniform coreset weights. This produces centroids that generalize k-means via soft assignments. We establish asymptotic consistency of the selected measure and Lipschitz stability to data perturbations, providing accuracy and robustness guarantees. Across synthetic and real-world benchmarks, the proposed method achieves competitive or improved approximation quality while substantially reducing runtime compared to Wasserstein- and standard Sinkhorn-based coreset selection, especially at large scale.
Figures & tables
| Dataset | Coreset size | Random | K-means | WCSL | EMS | |||
| dist / time | dist / time | dist / time | dist / time | dist / time | dist / time | dist / time | ||
| Gaussian | 18.28 / - | 10.87 / 0.09 | 10.85 / 14.9 | 10.87 / 0.08 | 10.86 / 0.15 | 10.87 / 0.15 | 10.87 / 0.12 | |
| 16.05 / - | 10.14 / 0.14 | 10.11 / 24.4 | 10.14 / 0.09 | 10.14 / 0.10 | 10.14 / 0.17 | 10.12 / 0.22 | ||
| 14.31 / - | 9.25 / 0.08 | 9.19 / 56.9 | 9.24 / 0.07 | 9.24 / 0.11 | 9.23 / 0.18 | 9.20 / 0.39 | ||
| 12.92 / - | 8.69 / 0.06 | 8.53 / 92.5 | 8.63 / 0.06 | 8.63 / 0.08 | 8.61 / 0.14 | 8.56 / 0.54 | ||
| Varying | Setup | WCSL | EMS | |||
| distance / time | distance / time | distance / time | distance / time | distance / time | ||
| Data size | ||||||
| Dimension | ||||||
| Census ( ) | |||
| Method | MeanAE | MaxAE | BVE |
| EMS( ) | |||
| EMS( ) | |||
| K-means | |||
| WCSL | |||
| Diabetes ( ) | |||
Appendix figures & tables1 asset
Supplementary material from the paper’s appendix.
Appendix
| Dataset | Random | K-means | WCSL | EMS | ||||
| distance / time | distance / time | distance / time | distance / time | distance / time | distance / time | distance / time | ||
| Transactions | 20 | 74.47 / – | 46.47 / 1.43 | 46.45 / 76.57 | 46.47 / 1.37 | 46.48 / 1.66 | 46.48 / 1.45 | 46.47 / 1.85 |
| 50 | 68.85 / – | 44.88 / 1.63 | 44.82 / 148.02 | 44.88 / 2.08 | 44.87 / 1.95 | 44.87 / 2.26 | 44.84 / 2.20 | |
| 100 | 65.65 / – | 43.67 / 1.76 | 43.56 / 203.11 | 43.66 / 2.24 | 43.65 / 2.81 | 43.65 / 3.00 | 43.61 / 3.00 | |
| 200 | 63.13 / – | 42.48 / 3.15 | 42.26 / 307.59 | 42.42 / 2.86 | 42.42 / 3.65 | 42.40 / 6.36 | 42.34 / 6.46 | |
Coreset Selection
Coreset selection chooses compact subsets that preserve a dataset’s learning, evaluation, or optimization value while reducing cost.
See all 57 papers