Interpretable discriminant analysis for functional data supported on random nonlinear domains with an application to Alzheimer's disease

Abstract

We introduce a novel framework for the classification of functional data supported on nonlinear, and possibly random, manifold domains. The motivating application is the identification of subjects with Alzheimer's disease from their cortical surface geometry and associated cortical thickness map. The proposed model is based upon a reformulation of the classification problem as a regularized multivariate functional linear regression model. This allows us to adopt a direct approach to the estimation of the most discriminant direction while controlling for its complexity with appropriate differential regularization. Our approach does not require prior estimation of the covariance structure of the functional predictors, which is computationally prohibitive in our application setting. We provide a theoretical analysis of the out-of-sample prediction error of the proposed model and explore the finite sample performance in a simulation setting. We apply the proposed method to a pooled dataset from the Alzheimer's Disease Neuroimaging Initiative and the Parkinson's Progression Markers Initiative. Through this application, we identify discriminant directions that capture both cortical geometric and thickness predictive features of Alzheimer's disease that are consistent with the existing neuroscience literature.

References in corpus (7)

• A reproducing kernel Hilbert space approach to functional linear regression
• Kernel Methods on Riemannian Manifolds with Gaussian RBF Kernels
• Statistical analysis of trajectories on Riemannian manifolds: Bird migration, hurricane tracking and video surveillance
• Functional single index models for longitudinal data
• Optimal Penalized Function-on-Function Regression under a Reproducing Kernel Hilbert Space Framework
• Three-dimensional Cardiovascular Imaging-Genetics: A Mass Univariate Framework
• Multivariate Spline Estimation and Inference for Image-On-Scalar Regression