Abstract:The Gaussian kernel and its traditional normalizations (e.g., row-stochastic) are popular approaches for assessing similarities between data points. Yet, they can be inaccurate under high-dimensional noise, especially if the noise magnitude varies considerably across the data, e.g., under heteroskedasticity or outliers. In this work, we investigate a more robust alternative -- the doubly stochastic normalization of the Gaussian kernel. We consider a setting where points are sampled from an unknown density on a low-dimensional manifold embedded in high-dimensional space and corrupted by possibly strong, non-identically distributed, sub-Gaussian noise. We establish that the doubly stochastic affinity matrix and its scaling factors concentrate around certain population forms, and provide corresponding finite-sample probabilistic error bounds. We then utilize these results to develop several tools for robust inference under general high-dimensional noise. First, we derive a robust density estimator that reliably infers the underlying sampling density and can substantially outperform the standard kernel density estimator under heteroskedasticity and outliers. Second, we obtain estimators for the pointwise noise magnitudes, the pointwise signal magnitudes, and the pairwise Euclidean distances between clean data points. Lastly, we derive robust graph Laplacian normalizations that accurately approximate various manifold Laplacians, including the Laplace Beltrami operator, improving over traditional normalizations in noisy settings. We exemplify our results in simulations and on real single-cell RNA-sequencing data. For the latter, we show that in contrast to traditional methods, our approach is robust to variability in technical noise levels across cell types.

Doubly Stochastic Graph Matrices

Doubly stochastic and combined matrices

Laplacian matrices of graphs: a survey

Structure Properties of a Doubly-Stochastic Process on a Network

On a question of Erdős on doubly stochastic matrices

Distributed strategies for generating weight-balanced and doubly stochastic digraphs

Convex Set of Doubly Substochastic Matrices

On the image of graph distance matrices

On the Eigenvalues of Graphs with Mixed Algebraic Structure

On Hoffman polynomials of $λ$-doubly stochastic irreducible matrices and commutative association schemes

The geometry of Gaussian double Markovian distributions

Positive semidefiniteness of -matrix on some families of graphs

Generalized unistochastic matrices

Ideal lattices of semigroup of doubly stochastic matrices

Structured BFGS Method for Optimal Doubly Stochastic Matrix Approximation.

Random matrices with row constraints and eigenvalue distributions of graph Laplacians

Generating Sets of Stochastic Matrices

On Steinerberger Curvature and Graph Distance Matrices

Robust Inference of Manifold Density and Geometry by Doubly Stochastic Scaling

Quotient graphs and stochastic matrices