The R Journal: accepted article

This article will be copy edited and may be changed before publication.

Benchmarking R packages for calculation of persistent homology PDF download
Eashwar V. Somasundaram, Shael E. Brown, Adam Litzler, Jacob G. Scott and Raoul R. Wadhwa

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 fastest. Ripser and TDAstats are also the most memory-efficient tools to calculate persistent homology.

Received: 2020-05-01; online 2021-06-07
CRAN packages: TDA, TDAstats, readr, ggplot2, scatterplot3d, recexcavAAR, deldir, magick, bench, pryr
CRAN Task Views implied by cited CRAN packages: Graphics, Multivariate, Phylogenetics, Spatial, TeachingStatistics

CC BY 4.0
This article is licensed under a Creative Commons Attribution 4.0 International license.

  author = {Eashwar V. Somasundaram and Shael E. Brown and Adam Litzler
          and Jacob G. Scott and Raoul R. Wadhwa},
  title = {{Benchmarking R packages for calculation of persistent
  year = {2021},
  journal = {{The R Journal}},
  doi = {10.32614/RJ-2021-033},
  url = {}