Principal Component Analysis of Persistent Homology Rank Functions with case studies of Spatial Point Patterns, Sphere Packing and Colloids

Abstract

Persistent homology, while ostensibly measuring changes in topology, captures multiscale geometrical information. It is a natural tool for the analysis of point patterns. In this paper we explore the statistical power of the (persistent homology) rank functions. For a point pattern $X$ we construct a filtration of spaces by taking the union of balls of radius $a$ centered on points in $X$, $X_a={\cup }_{x\in X}B(x,a)$. The rank function ${\beta }_k(X):\{(a,b)\in {\mathbb{R}}^2:a\le b\}\to \mathbb{R}$ is then defined by ${\beta }_k(X)(a,b)=rank({\iota }_{\ast }:H_k(X_a)\to H_k(X_b))$ where ${\iota }_{\ast }$ is the induced map on homology from the inclusion map on spaces. We consider the rank functions as lying in a Hilbert space and show that under reasonable conditions the rank functions from multiple simulations or experiments will lie in an affine subspace. This enables us to perform functional principal component analysis which we apply to experimental data from colloids at different effective temperatures and of sphere packings with different volume fractions. We also investigate the potential of rank functions in providing a test of complete spatial randomness of 2D point patterns using the distances to an empirically computed mean rank function of binomial point patterns in the unit square.

25 pages, 15 figures

References in corpus (4)

• The Statistical Physics of Athermal Materials
• An invariant distribution in static granular media
• Betti number signatures of homogeneous Poisson point processes
• Fréchet Means for Distributions of Persistence diagrams

Cited by in corpus (17)

• A persistence landscapes toolbox for topological statistics
• Persistence weighted Gaussian kernel for topological data analysis
• Sliced Wasserstein Kernel for Persistence Diagrams
• The persistence landscape and some of its properties
• Persistent homology detects curvature
• Topological data analysis distinguishes parameter regimes in the Anderson-Chaplain model of angiogenesis
• A Short Survey of Topological Data Analysis in Time Series and Systems Analysis
• Machine Learning Approaches to Solar-Flare Forecasting: Is Complex Better?
• Topological Learning for Brain Networks
• Wasserstein Stability for Persistence Diagrams
• Data structures for real multiparameter persistence modules
• Kernel method for persistence diagrams via kernel embedding and weight factor
• Computational Topology Techniques for Characterizing Time-Series Data
• Persistent homology and topological statistics of hyperuniform point clouds
• Distributing Persistent Homology via Spectral Sequences
• Fast Topological Clustering with Wasserstein Distance
• Stability for Inference with Persistent Homology Rank Functions