PhyloMeasures: a package for computing phylogenetic biodiversity measures and their statistical moments

We present PhyloMeasures, a new software package including both a C++ and R version, that provides very fast computation of popular phylogenetic diversity measures. PhyloMeasures introduces two major advances over existing methods. First, it uses efficient algorithms for calculating basic phylogenetic metrics (such as Faith's PD and the mean pairwise distance, MPD) and two-sample measures (such as common branch length, CBL, and the unique fraction) that are designed to perform well even on very large trees. Second, it computes exact richness-standardised versions of these measures (such as the widely used net relatedness index, NRI) by efficiently evaluating analytical expressions for the mean and variance of the basic measures, rather than by the slow and inexact randomization techniques that are the current standard. Together, these lead to massive improvements in performance compared to the current state of the art. For example, running on a standard laptop, PhyloMeasures functions can provide the NRI for 20 samples from a tree of 100 000 tips in about 1.5 s, compared to an estimated 37 d using standard resampling approaches. This will allow analyses on larger data sets than were previously possible.