Detecting symmetry from data is a fundamental problem in signal analysis, providing insight into underlying structures and constraints. When data emerge as trajectories of dynamical systems, symmetries encode structural properties of the dynamics that enable model reduction, principled comparison across conditions, and detection of regime changes. While recent optimal transport methods provide practical tools for data-driven symmetry detection in this setting, they rely on deterministic thresholds and lack uncertainty quantification, limiting robustness to noise and ability to resolve hierarchical symmetry structures. We present a Bayesian framework that formulates symmetry detection as probabilistic model selection over a lattice of candidate subgroups, using a Gibbs posterior constructed from Wasserstein distances between observed data and group-transformed copies. We establish three theoretical guarantees: (i) a Bayesian Occam's razor favoring minimal symmetry consistent with data, (ii) conjugation equivariance ensuring frame-independence, and (iii) stability bounds under perturbations for robustness to noise. Posterior inference is performed via Metropolis-Hastings sampling and numerical experiments on equivariant dynamical systems and synthetic point clouds demonstrate accurate symmetry recovery under high noise and small sample sizes. An application to human gait dynamics reveals symmetry degradation induced by mechanical constraints, demonstrating the framework's utility for statistical inference in biomechanical and dynamical systems.
Ghanem et al. (Wed,) studied this question.