Abstract:All neuroimaging packages can handle group analysis with t-tests or general linear modeling (GLM). However, they are quite hamstrung when there are multiple within-subject factors or when quantitative covariates are involved in the presence of a within-subject factor. In addition, sphericity is typically assumed for the variance-covariance structure when there are more than two levels in a within-subject factor. To overcome such limitations in the traditional AN(C)OVA and GLM, we adopt a multivariate modeling (MVM) approach to analyzing neuroimaging data at the group level with the following advantages: a) there is no limit on the number of factors as long as sample sizes are deemed appropriate; b) quantitative covariates can be analyzed together with within-subject factors; c) when a within-subject factor is involved, three testing methodologies are provided: traditional univariate testing (UVT) with sphericity assumption (UVT-UC) and with correction when the assumption is violated (UVT-SC), and within-subject multivariate testing (MVT-WS); d) to correct for sphericity violation at the voxel level, we propose a hybrid testing (HT) approach that achieves equal or higher power via combining traditional sphericity correction methods (Greenhouse-Geisser and Huynh-Feldt) with MVT-WS. To validate the MVM methodology, we performed simulations to assess the controllability for false positives and power achievement. A real FMRI dataset was analyzed to demonstrate the capability of the MVM approach. The methodology has been implemented into an open source program 3dMVM in AFNI, and all the statistical tests can be performed through symbolic coding with variable names instead of the tedious process of dummy coding. Our data indicates that the severity of sphericity violation varies substantially across brain regions. The differences among various modeling methodologies were addressed through direct comparisons between the MVM approach and some of the GLM implementations in the field, and the following two issues were raised: a) the improper formulation of test statistics in some univariate GLM implementations when a within-subject factor is involved in a data structure with two or more factors, and b) the unjustified presumption of uniform sphericity violation and the practice of estimating the variance-covariance structure through pooling across brain regions.

Joint Generalized Models for Multidimensional Outcomes: A Case Study of Neuroscience Data from Multimodalities

A Joint Multimodal Group Analysis Framework for Modeling Corticomuscular Activity.

Joint model for survival and multivariate sparse functional data with application to a study of Alzheimer's Disease

Three-part Joint Modeling Methods for Complex Functional Data Mixed with Zero-and-one-inflated Proportions and Zero-Inflated Continuous Outcomes with Skewness

MULTIVARIATE VARYING COEFFICIENT MODEL FOR FUNCTIONAL RESPONSES

Exploration and Application of a Muscle Fatigue Assessment Model Based on NMF for Multi-Muscle Synergistic Movements

Modeling Heterogeneity and Dependence for Analysis of Neuronal Data.

Modeling Multivariate Mixed-Response Functional Data

A Tridirectional Method for Corticomuscular Coupling Analysis in Parkinson'S Disease

Applications of multivariate modeling to neuroimaging group analysis: a comprehensive alternative to univariate general linear model

Bayesian Joint Modeling of Multiple Brain Functional Networks

Latent Multimodal Functional Graphical Model Estimation

A three-step multimodal analysis framework for modeling corticomuscular activity with application to Parkinson's disease

Extracting Stable Control Information from EMG Signals to Drive a Musculoskeletal Model - A Preliminary Study

Multimodal Predictive Modeling: Scalable Imaging Informed Approaches to Predict Future Brain Health

A Generalized Regression Model for Region of Interest Analysis of Fmri Data.

Joint analysis of multivariate longitudinal, imaging, and time-to-event data

Assessing Time-Dependent Association Between Scalp EEG and Muscle Activation: A Functional Random-Effects Model Approach

Joint Modeling of An Outcome Variable and Integrated Omic Datasets Using GLM-PO2PLS

Flexible joint models for multivariate longitudinal and time-to-event data using multivariate functional principal components

An Approach to Extract Nonlinear Muscle Synergies from sEMG through Multi-Model Learning