Annals of Statistics
- Ann. Statist.
- Volume 39, Number 5 (2011), 2716-2739.
Principal components analysis in the space of phylogenetic trees
Phylogenetic analysis of DNA or other data commonly gives rise to a collection or sample of inferred evolutionary trees. Principal Components Analysis (PCA) cannot be applied directly to collections of trees since the space of evolutionary trees on a fixed set of taxa is not a vector space. This paper describes a novel geometrical approach to PCA in tree-space that constructs the first principal path in an analogous way to standard linear Euclidean PCA. Given a data set of phylogenetic trees, a geodesic principal path is sought that maximizes the variance of the data under a form of projection onto the path. Due to the high dimensionality of tree-space and the nonlinear nature of this problem, the computational complexity is potentially very high, so approximate optimization algorithms are used to search for the optimal path. Principal paths identified in this way reveal and quantify the main sources of variation in the original collection of trees in terms of both topology and branch lengths. The approach is illustrated by application to simulated sets of trees and to a set of gene trees from metazoan (animal) species.
Ann. Statist., Volume 39, Number 5 (2011), 2716-2739.
First available in Project Euclid: 22 December 2011
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Primary: 92D15: Problems related to evolution
Secondary: 62H25: Factor analysis and principal components; correspondence analysis
Nye, Tom M. W. Principal components analysis in the space of phylogenetic trees. Ann. Statist. 39 (2011), no. 5, 2716--2739. doi:10.1214/11-AOS915. https://projecteuclid.org/euclid.aos/1324563353
- Supplementary material: Principal components analysis in the space of phylogenetic trees: Supplementary information. This contains further information about the simulation studies in Section 4.1.