Search papers, labs, and topics across Lattice.
The paper introduces the Directional Kernel Mean Difference (DKMD), a novel signed statistic designed for univariate distribution comparison that retains directional information about distributional shifts. By integrating kernel mean embeddings with a fixed odd weighting function, DKMD exhibits properties such as antisymmetry and directional monotonicity, while also providing a computationally efficient $O(N \log N)$ algorithm for large datasets. Experimental results show that DKMD effectively distinguishes directional shifts from symmetric perturbations and is resilient to outliers, demonstrating its practical applicability in high-dimensional settings.
DKMD isolates directional shifts in distributions while being robust to outliers and efficient enough to handle millions of samples in seconds.
We introduce the Directional Kernel Mean Difference (DKMD), a signed statistic for univariate distribution comparison that preserves the direction of distributional shifts. Unlike the squared Maximum Mean Discrepancy (MMD), which discards directional information by squaring the RKHS distance, DKMD integrates the difference of kernel mean embeddings against a fixed odd weighting function. This construction yields three structural properties: antisymmetry, immunity to symmetric distributional differences, and directional monotonicity under stochastic dominance. We derive a data-driven Riemann estimator that ensures asymptotic consistency with the continuous formulation, strictly preserving the theoretical guarantees of the signed statistic in empirical evaluations. To overcome the quadratic computational cost of kernel methods, we develop an $O(N \log N)$ prefix--suffix scanning algorithm that exploits the total order of the real line while requiring only $O(N)$ memory. Experiments on synthetic benchmarks demonstrate that DKMD correctly isolates directional shifts from symmetric perturbations, remains robust to heavy-tailed outliers that can flip the sign of the mean difference, and scales to millions of samples in seconds.