REVIEW 3 major objections 4 minor 1 cited by
Robust Representation and Estimation of Barycenters and Modes of Probability Measures on Metric Spaces
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Classical Fréchet means and medians can flip under tiny perturbations, but the barycentric merge tree is Lipschitz-stable under Wasserstein perturbations and consistently estimable from samples.
desk verdict The barycentric merge tree is a real addition to non-Euclidean statistics and the main stability bound holds up, but the connectivity-constant hypothesis is under-advertised and Lemma 6.1 has a fixable constant error. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The object that carries the argument is the barycentric merge tree (BMT): the quotient space $\mathcal{T}_p(\mathcal{X})=X/{\sim}$, where $x\sim y$ when $x$ and $y$ lie in the same connected component of the sublevel set $\sigma_p^{-1}((-\infty,t])$ at level $t=\sigma_p(x)=\sigma_p(y)$. The quotient map $\alpha_p$ pushes $\mu$ to $\mu_p$, and the merge height defines the pseudo-metric $d_p$; the tree also carries the descending function $\kappa_p$. Three ingredients make stability work: Proposition 2.2 shows the $p$-deviation function is $\theta$-Lipschitz; Proposition 3.12 shows $\alpha_p$ is $KL$-Lipschitz whenever the space has connectivity constant $K$ and $\theta$ is $L$-admissible; and the Coupling Lemma (Lemma 4.4) converts a uniform bound $|\sigma_p-\sigma'_p|\le r$ into a metric coupling of the two trees with distortion at most $2r$. The final inequality follows by coupling the underlying measures with a Wasserstein-optimal transport plan and applying Minkowski's inequality to the structural and functional offsets.
What would settle it
On the Hawaiian earring, the closed union of circles tangent at the origin with radii $1/n$, the connectivity modulus is infinite: take $x_n$ and $y_n$ on adjacent circles at the leftmost points; their metric distance is $1/(n(n+1))$, while every path between them passes through the origin at distance about $1/n$ from both endpoints. Computing the BMT stability ratio on this space would show that no finite Lipschitz constant can hold without the $K$ assumption, so the theorem's scope is exactly the class of spaces with finite connectivity modulus.
Extended reading notes
Core claim
The discovery is that the right object to estimate is not the set of minimizers of the deviation function but the whole merging pattern of its sublevel sets, made into a functional metric measure space. The barycentric merge tree $\mathcal{T}_p(\mathcal{X})$ is the quotient of $X$ by the relation that identifies points lying in the same connected component of a sublevel set $\sigma_p^{-1}((-\infty,t])$; the quotient map pushes $\mu$ forward to $\mu_p$, and the merge height defines a pseudo-metric $d_p$ with the property that the induced function $\kappa_p$ is 1-Lipschitz. Theorem 4.5 proves $D_{\mathrm{KS},p}(\mathcal{F}_p,\mathcal{F}'_p)\le L(1+K)w_p(\mu,\mu')$ for connected, locally path-connected Polish spaces, where $\mathcal{F}_p=(\mathcal{T}_p,d_p,\mu_p,\kappa_p)$. Because $K=1$ for geodesic spaces, the bound becomes $D_{\mathrm{KS},p}\le 2w_p$ when $\theta=d_X$. Corollary 4.7 then uses empirical Wasserstein convergence rates to conclude almost-sure consistency and $\mathbb{E}[D_{\mathrm{KS},p}]\le C\,\mathrm{diam}(X)\,n^{-1/s}$ for any $s$ above the upper $p$-Wasserstein dimension of $X$.
Load-bearing premise
The stability bound requires a finite connectivity constant $K$: every pair of points must be joinable by a path whose maximum distance to the two endpoints is at most $K$ times their distance, and connected Polish spaces need not satisfy this.
Editorial extensions
If this is right
- On geodesic spaces with $\theta=d_X$, the stability bound is $D_{\mathrm{KS},p}\le 2w_p(\mu,\mu')$, so BMTs are Lipschitz-stable summaries rather than only continuous ones.
- Empirical BMTs converge almost surely to the population BMT, with expected error at most $C\,n^{-1/s}$ for every $s$ above the upper Wasserstein dimension; sample size directly controls the accuracy of the topological summary.
- Mode merge trees built from diffusion distances of uniformly Lipschitz kernels, including heat kernels on compact manifolds, inherit the same stability and consistency guarantees.
- For compact spaces, a $\delta$-covering graph with edges between points at distance at most $3\delta$ yields a combinatorial BMT that is $C\delta$-close in the Kantorovich–Sturm distance, and a binning step keeps the simplified tree within $\varepsilon$.
- The framework applies even when barycenters are highly non-unique, such as uniform antipodal masses on the circle, because the BMT records the full merging pattern instead of selecting one minimizer.
Reading between the lines
- A natural two-sample test for distributions on manifolds would compare empirical BMTs in $D_{\mathrm{KS},p}$; the $n^{-1/s}$ rate gives a control on the null distribution, though the authors do not develop such a test.
- The diffusion-kernel construction yields a one-parameter family of BMTs as the scale $t$ varies, and the paper notes phase transitions; choosing $t$ automatically from the largest gap in $D_{\mathrm{KS},p}$ across scales is a testable extension the authors leave open.
- Because the bound scales linearly with the admissibility constant $L$, using a diffusion distance can inflate constants; whether a $\theta$-adapted Wasserstein distance would absorb $L$ and sharpen the rates is a natural next step.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces the barycentric merge tree (BMT), a functional merge-tree summary of the p-deviation function of a probability measure on a metric space, and proposes it as a robust representation of barycenters and modes. The BMT is a quotient of the space by connected components of sublevel sets, equipped with a merge metric and pushforward measure; modes are treated as barycenters for diffusion pseudo-metrics. The main theoretical results are: a Lipschitz stability bound DKS,p(Fp,F'p) ≤ L(1+K) wp(μ,μ') (Theorem 4.5), empirical consistency with rates n^{-1/s} (Corollary 4.7), a discrete δ-approximation scheme with O(δ) error (Theorem 6.2 and Corollary 6.3), and a binning simplification step. Numerical examples on the circle, sphere, and polygon shape space illustrate the framework, and code is provided.
Significance. A stable, computable summary of barycenters and modes would be a meaningful advance, and the proof of Theorem 4.5 is clean and self-contained given standard optimal transport machinery, with explicit constants and a clear Coupling Lemma. The paper also gives a discrete approximation pipeline with publicly available code and derives explicit consistency rates from Weed–Bach; these are concrete strengths. However, the hypotheses needed for the main results are narrower than the abstract suggests: the finite connectivity constant K is not implied by the standing Polish-space assumptions, and the discrete approximation lemma contains a proof gap. The central stability claim itself appears sound, but the scope statements and the discrete approximation guarantee need correction.
major comments (3)
- [§2.2, Definition 2.3; §4, Theorem 4.5; §7] The standing hypotheses 'connected and locally path-connected Polish' do not imply the existence of a finite connectivity constant K. For the Hawaiian earring H = ⋃ C_n, where C_n is the circle of radius 1/n centered at (1/n,0), take x_n=(2/n,0) and y_n=(2/(n+1),0); then d_X(x_n,y_n)=2/(n(n+1)), but every path between x_n and y_n passes through (0,0), so r_X(x_n,y_n) ≥ 2/n and the ratio is at least n+1, giving K_H=∞. Since Proposition 3.12, Theorem 4.5, and Corollary 4.7 all require a finite K, the abstract's 'more general metric spaces' and Section 7's claim of results 'at the generality of all Borel probability measures on a Polish metric space' overstate the scope. The authors should add finite connectivity modulus to the standing assumptions in the abstract and Section 7, or prove the results under weaker assumptions.
- [§6.2, Lemma 6.1, Eq. (58)] The proof of the second inequality in Lemma 6.1 applies the connectivity condition to the wrong endpoints. Equation (58) bounds d_X(v,γ_i(t))∨d_X(γ_i(t),w), where v,w are the endpoints of the full path, but Definition 2.3 applied to the subpath from z_{i-1} to z_i would only justify a bound on d_X(z_{i-1},γ_i(t))∨d_X(γ_i(t),z_i). Consequently, the displayed conclusion 'd_X(γ(t), z_{i(t)}) ≤ ε+Kδ' does not follow from (58). Re-deriving the estimate with per-edge connectivity yields an additional factor, at least 3 since d_V(z_{i-1},z_i)≤3δ, so Lemma 6.1 and therefore Theorem 6.2 and Corollary 6.3 need a corrected proof or modified constants. The main stability theorem, Theorem 4.5, is unaffected.
- [§4, Corollary 4.7(ii)] The statement of Corollary 4.7(ii) omits the bounded-diameter hypothesis that its proof invokes. The proof cites Weed–Bach 'under the hypothesis that diam(X)≤1', but the corollary as stated asserts E[DKS,p] ≤ C diam(X) n^{-1/s} without any finiteness or boundedness condition on diam(X). If X is unbounded, the bound is not meaningful, and if diam(X) is only finite, the constant in the Weed–Bach estimate generally depends on the diameter. The bounded-diameter assumption (or an equivalent moment/support condition) should be added to the statement.
minor comments (4)
- [§4.2, heading] The heading of Section 4 refers to 'FMTs' while the rest of the paper uses 'BMTs'; the notation should be made consistent.
- [§3, Definition 3.5] The phrase 'cf. [13, 19]' after the definition is helpful, but the merge-height distance is essentially the cophenetic distance for merge trees; a brief explanatory sentence connecting these notions would improve readability.
- [§6.2, Theorem 6.2 proof] The proof of Theorem 6.2 refers to 'Proposition 3.12 and (6.2)' immediately before the structural-offset computation; the equation reference should be made explicit, since the displayed equation defining δ_r is not labeled in the text.
- [§7, last paragraph] The GitHub link is a useful reproducibility statement, but the corresponding repository should be cited in a stable archival form, such as a DOI or Zenodo record, for the published version.
Circularity Check
No circularity: BMT stability is a parameter-free derivation under explicit connectivity assumptions.
full rationale
The central derivation is self-contained. Theorem 4.5 is proved directly: equations (36)-(38) bound the pointwise deviation difference by L times the Wasserstein distance using Minkowski's inequality and L-admissibility; the Coupling Lemma (Lemma 4.4) converts that bound into a metric coupling on the merge trees; Proposition 3.12 uses the finite connectivity constant K to make the quotient map alpha_p KL-Lipschitz; and equations (39)-(42) assemble the final Lipschitz bound. No fitted parameter is relabeled as a prediction: L and K are explicit standing hypotheses, and the Wasserstein distance is an independent external quantity, not a function of the barycentric merge trees. Corollary 4.7 imports the empirical Wasserstein convergence estimate of Weed and Bach [38], which is external to the authors' prior work. Citations to [13] and [14] concern standard merge-tree metrics and Reeb-space arguments, but the key proof of Proposition 3.8 is reproduced rather than assumed, and no ‘uniqueness theorem' or unverified self-citation is load-bearing. Citation [15] supplies a routine diffusion-distance formula, but the admissibility of the heat-kernel diffusion distance is also directly established in Example 5.7. The finite connectivity constant K is a genuine substantive hypothesis: some connected, locally path-connected Polish spaces have infinite connectivity modulus, so the abstract's and Section 7's broad phrasing overstates the scope of Theorem 4.5 and Corollary 4.7. That is a scope/correctness caveat, not circularity, because the theorem states the assumption rather than hiding the conclusion in it. Similarly, the alleged proof issue around inequality (58) in Lemma 6.1 would be a constant-factor correctness matter, not a reduction of a conclusion to its input. Overall, the claimed stability and consistency results are derived, not defined, into existence.
Assumptions & free parameters
free parameters (2)
- Diffusion scale t =
user-selected; not fitted
- Grid spacing δ and binning level ε =
user-selected accuracy controls
assumptions (9)
- standard math Minkowski inequality and standard Lp norm properties
- standard math Properties of Wasserstein distances and Kantorovich-Sturm distances (Sturm 2006, Mémoli 2011)
- standard math Burago-Burago-Ivanov proposition on metric couplings (Prop 4.3)
- standard math Weed-Bach empirical Wasserstein convergence rate E[w_p(µ,µ_n)] ≤ C n^{-1/s}
- domain assumption The space (X,d_X) is connected, locally path-connected and Polish, with finite connectivity modulus K
- domain assumption The pseudo-metric θ is L-admissible (Lipschitz with respect to d_X)
- domain assumption Heat kernel Lipschitz estimates on closed manifolds (Kasue-Kumura) and finiteness of reference measure ν
- domain assumption Compactness and finite δ-covering for the discrete model
- domain assumption For the shape-space example, the moduli space of n-gons is identified with Grassmannian Gr2(R^n) with its canonical metric
invented entities (1)
-
Barycentric merge tree (BMT)
independent evidence
Cite this review
Pith. "Pith review of Robust Representation and Estimation of Barycenters and Modes of Probability Measures on Metric Spaces." pith.science (2026). https://pith.science/paper/PTSZUTKW
@misc{pith2026250509609,
author = {Pith},
title = {Pith review of: Robust Representation and Estimation of Barycenters and Modes of Probability Measures on Metric Spaces},
year = {2026},
howpublished = {\url{https://pith.science/paper/PTSZUTKW}},
note = {Machine review of arXiv:2505.09609}
}
read the original abstract
This paper is concerned with the problem of defining and estimating statistics for distributions on spaces such as Riemannian manifolds and more general metric spaces. The challenge comes, in part, from the fact that statistics such as means and modes may be unstable: for example, a small perturbation to a distribution can lead to a large change in Fr\'echet means on spaces as simple as a circle. We address this issue by introducing a new merge tree representation of barycenters called the barycentric merge tree (BMT), which takes the form of a measured metric graph and summarizes features of the distribution in a multiscale manner. Modes are treated as special cases of barycenters through diffusion distances. In contrast to the properties of classical means and modes, we prove that BMTs are stable -- this is quantified as a Lipschitz estimate involving optimal transport metrics. This stability allows us to derive a consistency result for approximating BMTs from empirical measures, with explicit convergence rates. We also give a provably accurate method for discretely approximating the BMT construction and use this to provide numerical examples for distributions on spheres and shape spaces.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 1 Pith paper
-
Observable Covariance and Principal Observable Analysis for Data on Metric Spaces
Principal Observable Analysis defines stable mean and covariance statistics for metric measure spaces using 1-Lipschitz functions and proves Wasserstein and Kantorovich-Sturm stability bounds for them.
Reference graph
Works this paper leans on
-
[1]
B. Afsari. Riemannian Lp center of mass: existence, uniqueness, and convexity. Proc. Amer. Math. Soc., 139:655–673, 2011
work page 2011
-
[2]
M. Agueh and G. Carlier. Barycenters in the Wasserstein space. SIAM J. Math. Anal. , 43(2):904–924, 2011
work page 2011
-
[3]
S. Anbouhi, W. Mio, and O. B. Okutan. On metrics for analysis of functional data on geometric domains. Found. Data Sci., 7(3):671–704, 2025
work page 2025
-
[4]
M. Arnaudon and L. Miclo. Means in complete manifolds: uniqueness and approximation. ESAIM: Probability and Statistics , 18:185–206, 2014
work page 2014
-
[5]
K. Beketayev, D. Yeliussizov, D. Morozov, G. Weber, and B. Hamann. Measuring the distance between merge trees. In Topological Methods in Data Analysis and Visualization III , pages 151–166, 2014
work page 2014
-
[6]
M. Berger. A Panoramic View of Riemannian Geometry . Springer-Verlag, 2003
work page 2003
-
[7]
R. Bhattacharya and V. Patrangenaru. Large sample theory of intrinsic and extrinsic sample means on manifolds. I. Ann. Statist., 31:1–29, 2003
work page 2003
-
[8]
R. Bhattacharya and V. Patrangenaru. Large sample theory of intrinsic and extrinsic sample means on manifolds. II. Ann. Statist., 33:1225–1259, 2005. 25
work page 2005
Show all 38 references
-
[9]
Burago, Y
D. Burago, Y. Burago, and S. Ivanov. A Course in Metric Geometry . American Mathematical Society, 2001
2001
-
[10]
Cantarella, T
J. Cantarella, T. Needham, C. Shonkwiler, and G. Stewart. Random triangles and polygons in the plane. The American Mathematical Monthly , 126(2):113–134, 2019
2019
-
[11]
Chaudhuri and J
P. Chaudhuri and J. S. Marron. Scale space view of curve estimation.Ann. Statist., 28:408–428, 2000
2000
-
[12]
R. R. Coifman and S. Lafon. Diffusion maps. Appl. Comput. Harmon. Anal., 21(1):5–30, 2006
2006
-
[13]
Curry, H
J. Curry, H. Hang, W. Mio, T. Needham, and O. B. Okutan. Decorated merge trees for persistent topology. J Appl. and Comput. Topology , 6(3):371–428, 2022
2022
-
[14]
Curry, W
J. Curry, W. Mio, T. Needham, O. B. Okutan, and F. Russold. Stability and approximations for decorated Reeb spaces. In Symposium on Computational Geometry , Athens, Greece, June 2024
2024
-
[15]
D. H. Diaz Martinez, C. H. Lee, P. T. Kim, and W. Mio. Probing the geometry of data with diffusion Fr´ echet functions.Appl. Comput. Harmon. Anal. , 47(3):935–947, 2019
2019
-
[16]
Edelman, T
A. Edelman, T. A. Arias, and S. T. Smith. The geometry of algorithms with orthogonality constraints. SIAM journal on Matrix Analysis and Applications , 20(2):303–353, 1998
1998
-
[17]
Flamary, N
R. Flamary, N. Courty, A. Gramfort, M. Z. Alaya, A. Boisbunon, S. Chambon, L. Chapel, A. Corenflos, K. Fatras, N. Fournier, et al. Pot: Python optimal transport. Journal of Machine Learning Research, 22(78):1–8, 2021
2021
-
[18]
Fr´ echet
M. Fr´ echet. Les ´ el´ ements al´ eatoires de nature quelconque dans un espace distanci´ e.Ann. Inst. H. Poincar´ e, 10:215–310, 1948
1948
-
[19]
Gasparovic, E
E. Gasparovic, E. Munch, S. Oudot, K. Turner, B. Wang, and Y. Wang. Intrinsic interleaving distance for merge trees. La Matematica, 4(1):40–65, 2025
2025
-
[20]
Grove and H
K. Grove and H. Karcher. How to conjugate C1-close group actions. Math. Z. , 132:11–20, 1973
1973
-
[21]
Grove, H
K. Grove, H. Karcher, and E. A. Ruh. Group actions and curvature. Invent. Math., 23:31–48, 1974
1974
-
[22]
Grove, H
K. Grove, H. Karcher, and E. A. Ruh. Jacobi fields and Finsler metrics on compact Lie groups with an application to differentiable pinching problems. Math. Ann., 211:7–21, 1974
1974
-
[23]
H. Hang, F. M´ emoli, and W. Mio. A topological study of functional data and Fr´ echet functions of metric measure spaces. J Appl. and Comput. Topology , 3(4):359–380, 2019
2019
-
[24]
Hausmann and A
J.-C. Hausmann and A. Knutson. Polygon spaces and grassmannians. L’Enseignement Math´ ematique, 43, 1997
1997
-
[25]
Hundrieser, B
S. Hundrieser, B. Eltzner, and S. Huckemann. A lower bound for estimating Fr´ echet means. arXiv:2402.12290, 2024
2024 arXiv
-
[26]
H. Karcher. Riemannian center of mass and mollifier smoothing. Comm. Pure Appl. Math. , 30:509–541, 1977. 26
1977
-
[27]
Kasue and H
A. Kasue and H. Kumura. Spectral convergence of Riemannian manifolds. Tohoku Math. J. , 46(2):147–179, 1994
1994
-
[28]
Kim and B
Y.-H. Kim and B. Pass. Wasserstein barycenters over Riemannian manifolds. Adv. Math , 307:640–683, 2017
2017
-
[29]
H. Le. Locating Fr´ echet means with application to shape spaces.Adv. Appl. Prob., 33:324–338, 2001
2001
-
[30]
Lindeberg
T. Lindeberg. Scale Space Theory in Computer Vision . Kluwer, Boston, 1994
1994
-
[31]
M´ emoli
F. M´ emoli. On the use of Gromov-Hausdorff distances for shape comparison. In Proceedings Point Based Graphics, pages 81–90, 01 2007
2007
-
[32]
M´ emoli
F. M´ emoli. Gromov-Wasserstein distances and the metric approach to object matching.Found. Comput. Math., 11(4):417–487, 2011
2011
-
[33]
Morozov, K
D. Morozov, K. Beketayev, and G. Weber. Interleaving distance between merge trees. Discrete Comput. Geom., 49:22–45, 2013
2013
-
[34]
R. R. Sokal and F. J. Rohlf. The comparison of dendrograms by objective methods. Taxon, 11:33–40, 1962
1962
-
[35]
K.-T. Sturm. On the geometry of metric measure spaces. Acta Mathematica, 196(1):65–131, 2006
2006
-
[36]
Vayer, L
T. Vayer, L. Chapel, R. Flamary, R. Tavenard, and N. Courty. Fused Gromov-Wasserstein distance for structured objects. Algorithms, 13(9):212, 2020
2020
-
[37]
C. Villani. Optimal Transport: Old and New , volume 338. Springer Science & Business Media, 2008
2008
-
[38]
Weed and F
J. Weed and F. Bach. Sharp asymptotic and finite-sample rates of convergence of empirical measures in Wasserstein distance. Bernoulli, 25(4A):2620–2648, 2019. 27
2019
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.