Benchmarking R Packages for Calculation of Persistent Homology

Authors

Eashwar V. Somasundaram, Case Western Reserve University
Shael E. Brown, McGill University
Adam Litzler, Cleveland Clinic Foundation
Jacob G. Scott, Cleveland Clinic Foundation
Raoul R. Wadhwa, Case Western Reserve University

Date of this Version

6-2021

Document Type

Article

Citation

The R Journal (June 2021) 13(1); Editor: Dianne Cook

Comments

Copyright 2021, The R Foundation. Open access material. License: CC BY 4.0 International

Abstract

Several persistent homology software libraries have been implemented in R. Specifically, the Dionysus, GUDHI, and Ripser libraries have been wrapped by the TDA and TDAstats CRAN packages. These software represent powerful analysis tools that are computationally expensive and, to our knowledge, have not been formally benchmarked. Here, we analyze runtime and memory growth for the 2 R packages and the 3 underlying libraries. We find that datasets with less than 3 dimensions can be evaluated with persistent homology fastest by the GUDHI library in the TDA package. For higher-dimensional datasets, the Ripser library in the TDAstats package is the fastest. Ripser and TDAstats are also the most memory-efficient tools to calculate persistent homology.