Principal component analysis and the locus of the Fréchet mean in the space of phylogenetic trees
MetadataShow full item record
Evolutionary relationships are represented by phylogenetic trees, and a phylogenetic analysis of gene sequences typically produces a collection of these trees, one for each gene in the analysis. Analysis of samples of trees is difficult due to the multi-dimensionality of the space of possible trees. In Euclidean spaces, principal component analysis is a popular method of reducing high-dimensional data to a low-dimensional representation that preserves much of the sample’s structure. However, the space of all phylogenetic trees on a fixed set of species does not form a Euclidean vector space, and methods adapted to tree space are needed. Previous work introduced the notion of a principal geodesic in this space, analogous to the first principal component. Here we propose a geometric object for tree space similar to the kth principal component in Euclidean space: the locus of the weighted Fréchet mean of k + 1 vertex trees when the weights vary over the k-simplex. We establish some basic properties of these objects, in particular showing that they have dimension k, and propose algorithms for projection onto these surfaces and for finding the principal locus associated with a sample of trees. Simulation studies demonstrate that these algorithms perform well, and analyses of two datasets, containing Apicomplexa and African coelacanth genomes respectively, reveal important structure from the second principal components.
The article of record as published may be found at http://dx.doi.org/10.1093/biomet/asx047
Showing items related by title, author, creator and subject.
Page, Robert L. (Monterey, CA; Naval Postgraduate School, 2019-06);A known challenge in the rapidly growing area of phylogenomics is the lack of tools to analyze the large volume of genome data. Genomic data includes information on the evolution, structure and mapping of genomes. Phylogenetic ...
Lin, Bo; Monod, Anthea; Yoshida, Ruriko (2018-06-04);We introduce a novel framework for the statistical analysis of phylogenetic trees: Palm tree space is con- structed on principles of tropical algebraic geometry, and represents phylogenetic trees as a point in a space ...
Forecasting tropical cyclone recurvature using an empirical othogonal [sic] function representation of vorticity fields Ford, Debra M. (Monterey, California: Naval Postgraduate School, 1990-09);An empirical orthogonal function (EOF) representation of relative vorticity is used to forecast recurvature (change in storm heading from west to east of 000 deg N) of western North Pacific tropical cyclones. The time-dependent ...