1 of 48

TOPOLOGICAL DATA ANALYSIS FOR MULTIPARAMETER DATA

Abigail Hickok

Columbia University

2 of 48

TOPOLOGICAL DATA ANALYSIS (TDA)

 

Persistent homology uses homology to captures the “shape” of a data set (e.g., a point cloud)

3 of 48

PERSISTENT HOMOLOGY (FORMALLY)

  •  

4 of 48

PERSISTENT HOMOLOGY (FORMALLY)

  •  

simplicial complexes

 

5 of 48

SUBLEVEL FILTRATIONS

  •  

 

6 of 48

MULTIPARAMETER TOPOLOGICAL DATA ANALYSIS

 

7 of 48

MULTIPARAMETER TDA

  • Previous examples involved a single parameter (e.g., the “scale parameter” for a point cloud filtration)
  • Often desirable/necessary to consider multiple parameters

 

1. Multiparameter persistent homology (G. Carlsson and A. Zomorodian, 2009)

Multiparameter persistent homology

8 of 48

MULTIPARAMETER TDA

  • Previous examples involved a single parameter (e.g., the “scale parameter” for a point cloud filtration)
  • Often desirable/necessary to consider multiple parameters

A vineyard is a 1-parameter continuously-varying set of PDs (along with “connecting information” )

 

2. Vineyards (Cohen-Steiner, Edelsbrunner, Morozov 2006)

 

 

 

9 of 48

MULTIPARAMETER TDA

  • Previous examples involved a single parameter (e.g., the “scale parameter” for a point cloud filtration)
  • Often desirable/necessary to consider multiple parameters

3. Persistence Diagram Bundles (AH, 2023+)

 

 

 

10 of 48

PERSISTENCE DIAGRAM BUNDLES

 

11 of 48

Specials Cases Of PDBs

 

 

Image (right): Turner, Mukherjee, Boyer. Persistent Homology Transform for Modeling Shapes and Surfaces (2014).

 

12 of 48

RELATIONSHIP TO FIBERED BARCODE OF MULTIPARAMETER PERSISTENCE MODULES

  •  

Fibered barcode of a bifiltration

13 of 48

MOTIVATING QUESTIONS

  • Are there discrete ways to encode the information in a persistence diagram bundle (PDB)?
  • What kinds of geometric structure do these objects have, and does the geometry of the PDB say anything about the geometry of the data?
  • How can we calculate PDBs so that they might be used for applications?

14 of 48

MOTIVATING QUESTIONS

  • Are there discrete ways to encode the information in a persistence diagram bundle (PDB)?
  • What kinds of geometric structure do these objects have, and does the geometry of the PDB say anything about the geometry of the data?
  • How can we calculate PDBs so that they might be used for applications?

15 of 48

BACKGROUND: BIRTH AND DEATH SIMPLICES OF SINGLE-PARAMETER FILTRATIONS

  •  

 

 

 

 

 

 

16 of 48

PERSISTENCE DIAGRAMS ARE DETERMINED BY BIRTH, DEATH SIMPLICES

  •  

 

 

 

 

 

 

17 of 48

Generic PDBs are Determined by Finitely Many Base Points

 

 

 

 

 

 

18 of 48

STRATIFICATION EXAMPLE

 

 

 

19 of 48

STRATIFICATION EXAMPLE

The PDB is determined by the stratification, the (birth, death) pairs in each stratum, and the filtration values of each simplex

 

20 of 48

CELLULAR SHEAVES

  •  

Image: Hansen and Ghrist, 2019

 

21 of 48

A compatible cellular sheaf for a PDB

  •  

 

 

22 of 48

CELLULAR SHEAF EXAMPLE

 

 

Morphisms in an associated cellular sheaf

23 of 48

MOTIVATING QUESTIONS

  • Are there discrete ways to encode the information in a persistence diagram bundle (PDB)?
  • What kinds of geometric structure do these objects have, and does the geometry of the PDB say anything about the geometry of the data?
  • How can we calculate PDBs so that they might be used for applications?

24 of 48

SECTIONS

  • In a vineyard, every point can be extended to a section
  • But in a PDB, sections may not exist

 

Each “vine” (curve) is a section of the vineyard

25 of 48

PDB WITH NO (NONTRIVIAL) SECTIONS

 

 

 

26 of 48

PDB WITH NO (NONTRIVIAL) SECTIONS

“Monodromy” in the PDB

 

27 of 48

CASE STUDY: MONODROMY IN THE �PERSISTENT HOMOLOGY TRANSFORM

Image: Turner, Mukherjee, Boyer. Persistent Homology Transform for Modeling Shapes and Surfaces (2014).

 

 

 

Question:

What does monodromy in the PHT of a shape means for the geometry of the shape itself?

 

(Subset of ongoing work with S. Arya, B. Giunti, AH, L. Kanari, S. McGuire, K. Turner)

28 of 48

PHT EXAMPLE WITH MONODROMY

29 of 48

PHT EXAMPLE WITHOUT MONODROMY

30 of 48

WHAT’S DIFFERENT ABOUT THESE SHAPES?

PHT of the spiral has monodromy

PHT of the five-arm star does not have monodromy

31 of 48

A SHAPE’S GEOMETRY IS REFLECTED �IN ITS PHT’S GEOMETRY

A convex shape

 

32 of 48

A SHAPE’S GEOMETRY IS REFLECTED �IN ITS PHT’S GEOMETRY

 

Example of a star shape.

 

33 of 48

PHT EXAMPLE WITHOUT MONODROMY

34 of 48

A SHAPE’S GEOMETRY IS REFLECTED �IN ITS PHT’S GEOMETRY

 

Example of a star shape.

 

35 of 48

MOTIVATING QUESTIONS

  • Are there discrete ways to encode the information in a persistence diagram bundle (PDB)?
  • What kinds of geometric structure do these objects have, and does the geometry of the PDB say anything about the geometry of the data?
  • How do we calculate PDBs so that they might be used for applications?

36 of 48

COMPUTING A PDB

  • The existence of a stratification makes it possible to compute PDBs
  • We restrict to piecewise-linear fibered filtrations

 

37 of 48

COMPUTING A PDB

  • The existence of a stratification makes it possible to compute PDBs
  • We restrict to piecewise-linear fibered filtrations

Step 1: Compute the stratification (in this case, line arrangement)

38 of 48

COMPUTING A PDB�STEP 1: LINE ARRANGEMENT

39 of 48

COMPUTING A PDB

  • The existence of a stratification makes it possible to compute PDBs
  • We restrict to piecewise-linear fibered filtrations

Step 1: Compute the stratification (in this case, line arrangement)

Step 2: Calculate the (birth, death) pairs in each stratum via updating procedure

Step 3: Query for PDs

 

40 of 48

EXAMPLE: PARAMETERIZED SET OF POINT CLOUDS

 

 

Distance between cluster centers

 

 

“Width” of clusters

41 of 48

EXAMPLE: PARAMETERIZED SET OF POINT CLOUDS

 

TP = Total persistence

 

42 of 48

AN APPLICATION TO GRAPH DATA

 

Curvature filtrations: Filtering a graph by ORC encodes structural information (e.g., communities) that we can use for graph ML tasks (Southern et al, 2023)

ORC on edges of a graph with two “communities”

 

43 of 48

GRAPH CURVATURE PDB

 

 

 

A graph

Mean total persistence function, averaged over set of graphs

 

44 of 48

DISCUSSION AND CONCLUSIONS

 

45 of 48

STABILITY

PDB’s are “pointwise” stable, but the global structure isn’t guaranteed to be stable

Two vines intersect

A perturbation of the vineyard

A different perturbation

Question: Is the global structure “generically” stable?

46 of 48

ALGEBRAIC STRUCTURE

 

 

47 of 48

CONCLUSIONS

Summary

  • In persistent homology, we often have multiple parameters to account for
  • PDBs are useful for identifying topological features of a parameterized set of point clouds or a parameterized set of filtrations of any form
  • For generic PDBs, the information can be encoded discretely
  • Monodromy within the PDB can reflect something geometrically about the data

Open questions

  • What can we say about the algebraic structure?
  • Are (generic) PDBs stable?
  • When is it possible to get a “vine decomposition”?

48 of 48

THANK YOU

Questions?