  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 three-dimensional 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 tree-space, 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\'echet mean of $k+1$ points in tree-space, 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
tree-space. 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 tree-space, revealing
structure within these complex data sets.
