Variational Bayesian partially linear mean shift models for high-dimensional Alzheimer's disease neuroimaging data

Alzheimer's disease can be diagnosed by analyzing brain images (eg, magnetic resonance imaging, MRI) and neuropsychological tests (eg, mini-mental state examination, MMSE). A partially linear mean shift model (PLMSM) is here proposed to investigate the relationship between MMSE score and high-dimensional regions of interest in MRI, and detect the outliers. In the presence of high-dimensional data, existing Bayesian approaches (eg, Markov chain Monte Carlo) to analyze a PLMSM take intensive computational cost and require huge memory, and have low convergence rate. To address these issues, a variational Bayesian inference is developed to simultaneously estimate parameters and nonparametric functions and identify outliers in a PLMSM. A Bayesian P-splines method is presented to approximate nonparametric functions, a Bayesian adaptive Lasso approach is employed to select predictors, and outliers are detected by the classification variable. Two simulation studies are conducted to assess the finite sample performance of the proposed method. An MRI dataset with elderly cognitive ability is provided to corroborate the proposed method.