Abstract:Phylogenetic trees are believed to contain a wealth of information on diversification processes. Comparing phylogenetic trees is not straightforward due to their high dimensionality. Researchers have therefore defined a wide range of one-dimensional summary statistics. However, it remains unexplored to what extent these summary statistics cover the same underlying information and what summary statistics best explain observed variation across phylogenies. Furthermore, a large subset of available summary statistics focusses on measuring the topological features of a phylogenetic tree, but are often only explored at the extreme edge cases of the fully balanced or unbalanced tree and not for trees of intermediate balance. Here, we introduce a new R package that provides speed optimized code to compute 54 summary statistics. We study correlations between summary statistics on empirical trees and on trees simulated using several diversification models. Furthermore, we introduce an algorithm to create intermediately balanced trees in a well-defined manner, in order to explore variation in summary statistics across a balance gradient. We find that almost all summary statistics are correlated with tree size, and it is difficult if not impossible to correct for tree size, unless the tree generating model is known. Furthermore, we find that across empirical and simulated trees, at least two large clusters of correlated summary statistics can be found, where statistics group together based on information used (topology or branching times). However, the finer grained correlation structure appears to depend strongly on either the taxonomic group studied (in empirical studies) or the diversification model (in simulation studies). Nevertheless, we can identify multiple groups of summary statistics that are strongly and consistently correlated, indicating that these statistics measure the same underlying property of a tree. Lastly, we find that almost all topological summary statistics vary non-linearly and sometimes even non-monotonically with our intuitive balance gradient. Therefore, in order to avoid introducing biases and missing underlying information, we advocate for selecting as many summary statistics as possible in phylogenetic analyses. With the introduction of the treestats package, which provides fast and reliable calculations, such an approach is now routinely possible.

ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R

DendroPy 5: a mature Python library for phylogenetic computing

phytools: an R package for phylogenetic comparative biology (and other things)

EvoPhylo: An r package for pre‐ and postprocessing of morphological data from relaxed clock Bayesian phylogenetics

pegas: an R package for population genetics with an integrated-modular approach

BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis

'treesliceR': a package for slicing phylogenies and inferring phylogenetic patterns over evolutionary time

DAMBE5: A Comprehensive Software Package for Data Analysis in Molecular Biology and Evolution

A dendrochronology program library in R (dplR)

BEAGLE: An Application Programming Interface and High-Performance Computing Library for Statistical Phylogenetics

Efficient comparative phylogenetics on large trees

PAReTT: A Python Package for the Automated Retrieval and Management of Divergence Time Data from the TimeTree Resource for Downstream Analyses

POY version 5: phylogenetic analysis using dynamic homologies under multiple optimality criteria

BEAGLE 3: Improved Performance, Scaling, and Usability for a High-Performance Computing Library for Statistical Phylogenetics

Phylogenetic tree statistics: a systematic overview using the new R package ‘treestats’

Arothron: An R package for geometric morphometric methods and virtual anthropology applications

PEAS V1.0: a Package for Elementary Analysis of SNP Data

phangorn: phylogenetic analysis in R

TreeExp1.0: R Package for Analyzing Expression Evolution Based on RNA-Seq Data.

DAMBE7: New and Improved Tools for Data Analysis in Molecular Biology and Evolution

Phylogeny.fr: robust phylogenetic analysis for the non-specialist