  Persistent homology barcodes and diagrams are a cornerstone of topological
data analysis. Widely used in many real data settings, they relate variation in
topological information (as measured by cellular homology) with variation in
data, however, they are challenging to use in statistical settings due to their
complex geometric structure. In this paper, we revisit the persistent homology
rank function -- an invariant measure of ``shape" that was introduced before
barcodes and persistence diagrams and captures the same information in a form
that is more amenable to data and computation. In particular, since they are
functions, techniques from functional data analysis -- a domain of statistics
adapted for functions -- apply directly to persistent homology when represented
by rank functions. Rank functions, however, have been less popular than
barcodes because they face the challenge that stability -- a property that is
crucial to validate their use in data analysis -- is difficult to guarantee,
mainly due to metric concerns on rank function space. However, rank functions
extend more naturally to the increasingly popular and important case of
multiparameter persistent homology. In this paper, we study the performance of
rank functions in functional inferential statistics and machine learning on
both simulated and real data, and in both single and multiparameter persistent
homology. We find that the use of persistent homology captured by rank
functions offers a clear improvement over existing approaches. We then provide
theoretical justification for our numerical experiments and applications to
data by deriving several stability results for single- and multiparameter
persistence rank functions under various metrics with the underlying aim of
computational feasibility and interpretability.
