A novel smoothing-based goodness-of-fit test of covariance for multivariate sparse functional data

Abstract: Accurately specifying covariance structures is critical for valid inference in longitudinal and functional data analysis, particularly when data are sparsely observed. In this study, we develop a global goodness-of-fit test to assess parametric covariance structures in multivariate sparse functional data. Our contribution is twofold. First, we extend the univariate goodness-of-fit test proposed by Chen et al. (2019) to better accommodate sparse data by improving error variance estimation and applying positive semi-definite smoothing to covariance estimation. These corrections ensure appropriate Type I error control under sparse designs. Second, we introduce a multivariate extension of the improved test that jointly evaluates covariance structures across multiple outcomes, employing novel test statistics based on the maximum and $\ell_2$ norms to account for inter-outcome dependencies and enhance statistical power. Through extensive simulation studies, we demonstrate that the proposed methods maintain proper Type I error rates and achieve greater power than univariate tests with multiple testing adjustments. Applications to longitudinal neuroimaging and clinical data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) and the Parkinson's Progression Marker Initiative (PPMI) illustrate the practical utility of the proposed methods for evaluating covariance structures in sparse multivariate longitudinal data.