Principal Component Analysis for Functional Data on Riemannian Manifolds and Spheres October 25, 2017 Short title: Functional Data on Riemannian Manifolds Xiongtao Dai Department of Statistics University of California, Davis Davis, CA 95616 U.S.A. Email:
[email protected] Hans-Georg M¨uller Department of Statistics University of California, Davis Davis, CA 95616 U.S.A. Email:
[email protected] arXiv:1705.06226v2 [math.ST] 24 Oct 2017 Research supported by NSF grants DMS-1407852 and DMS-1712864. We thank FlightAware for the permission to use the flight data in Section 5. ABSTRACT Functional data analysis on nonlinear manifolds has drawn recent interest. Sphere-valued functional data, which are encountered for example as movement trajectories on the surface of the earth, are an important special case. We consider an intrinsic principal component analysis for smooth Riemannian manifold-valued functional data and study its asymptotic properties. Riemannian functional principal component analysis (RFPCA) is carried out by first mapping the manifold-valued data through Riemannian logarithm maps to tangent spaces around the time-varying Fr´echet mean function, and then performing a classical multi- variate functional principal component analysis on the linear tangent spaces. Representations of the Riemannian manifold-valued functions and the eigenfunctions on the original manifold are then obtained with exponential maps. The tangent-space approximation through func- tional principal component analysis is shown to be well-behaved in terms of controlling the residual variation if the Riemannian manifold has nonnegative curvature. Specifically, we derive a central limit theorem for the mean function, as well as root-n uniform convergence rates for other model components, including the covariance function, eigenfunctions, and functional principal component scores.