Principal component analysis and the locus of the Frechet mean in the space of phylogenetic trees
Abstract
Most biological data are multidimensional, posing a major challenge to human comprehension and computational analysis. Principal component analysis is the most popular approach to rendering two or threedimensional representations of the major trends in such multidimensional data. The problem of multidimensionality is acute in the rapidly growing area of phylogenomics. Evolutionary relationships are represented by phylogenetic trees, and very typically a phylogenomic analysis results in a collection of such trees, one for each gene in the analysis. Principal component analysis offers a means of quantifying variation and summarizing a collection of phylogenies by dimensional reduction. However, the space of all possible phylogenies on a fixed set of species does not form a Euclidean vector space, so principal component analysis must be reformulated in the geometry of treespace, which is a CAT(0) geodesic metric space. Previous work has focused on construction of the first principal component, or principal geodesic. Here we propose a geometric object which represents a $k$th order principal component: the locus of the weighted Fréchet mean of $k+1$ points in treespace, where the weights vary over the standard $k$dimensional simplex. We establish basic properties of these objects, in particular that locally they generically have dimension $k$, and we propose an efficient algorithm for projection onto these surfaces. Combined with a stochastic optimization algorithm, this projection algorithm gives a procedure for constructing a principal component of arbitrary order in treespace. Simulation studies confirm these algorithms perform well, and they are applied to data sets of Apicomplexa gene trees and the African coelacanth genome. The results enable visualizations of slices of treespace, revealing structure within these complex data sets.
 Publication:

arXiv eprints
 Pub Date:
 September 2016
 arXiv:
 arXiv:1609.03045
 Bibcode:
 2016arXiv160903045N
 Keywords:

 Statistics  Methodology;
 60D05 (Primary) 62H25;
 92D15 (Secondary)
 EPrint:
 26 pages, 5 figures