Bayesian Semiparametric Orthogonal Tucker Factorized Mixed Models for Multi-dimensional Longitudinal Functional Data
Arkaprava Roy, Abhra Sarkar
TL;DR
This work tackles the challenge of analyzing ultra-high-resolution, multi-dimensional longitudinal functional data by introducing a Bayesian semiparametric orthogonal Tucker factorized mixed model. It integrates basis-function expansions with Tucker tensor decomposition under semi-orthogonal mode matrices enforced by CSC-PING priors and a graph-Laplacian smoothness structure to capture spatio-temporal dependencies while enabling automatic rank selection via cumulative shrinkage. The framework jointly models baseline and time-varying effects with subject-specific random-effects cores, yielding a cohesive mean-covariance tensor representation that remains computationally scalable through parallel core updates. The method demonstrates strong performance on simulated data and yields novel insights into local brain changes associated with Alzheimer's disease progression in ADNI-3 diffusion MRI data, highlighting its potential for broad applicability in multi-dimensional longitudinal studies across domains.
Abstract
We introduce a novel longitudinal mixed model for analyzing complex multidimensional functional data, addressing challenges such as high-resolution, structural complexities, and computational demands. Our approach integrates dimension reduction techniques, including basis function representation and Tucker tensor decomposition, to model complex functional (e.g., spatial and temporal) variations, group differences, and individual heterogeneity while drastically reducing model dimensions. The model accommodates multiplicative random effects whose marginalization yields a novel Tucker-decomposed covariance-tensor framework. To ensure scalability, we employ semi-orthogonal mode matrices implemented via a novel graph-Laplacian-based smoothness prior with low-rank approximation, leading to an efficient posterior sampling method. A cumulative shrinkage strategy promotes sparsity and enables semiautomated rank selection. We establish theoretical guarantees for posterior convergence and demonstrate the method's effectiveness through simulations, showing significant improvements over existing techniques. Applying the method to Alzheimer's Disease Neuroimaging Initiative (ADNI) neuroimaging data reveals novel insights into local brain changes associated with disease progression, highlighting the method's practical utility for studying cognitive decline and neurodegenerative conditions.