Back

A Test for Confounding in Coupling of Multimodal Neuroimaging Data

Hao, Y.; Vandekar, S.; Alexander-Bloch, A.; Satterthwaite, T.; White, B.; Shinohara, R.

2026-09-21 neuroscience
10.64898/2026.09.14.750774 bioRxiv
Show abstract

Multimodal neuroimaging studies often involve comparisons between brain maps. Recently, statistical methods have been proposed to quantify and assess spatial correspondence between two modalities. The simple permutation-based inter-modal correspondence (SPICE) test evaluates whether within-subject correspondence exceeds chance, where the null distribution is constructed by permuting subject labels for one modality. Despite its easy implementation and minimal spatial assumptions, the critical assumption underlying permutation analysis, that subjects are exchangeable under the null, may be violated when covariates such as age, sex, or disease status systematically alter brain map distributions. This violation can produce results analogous to Simpson's paradox, where apparent population-level correspondence reflects between-group differences rather than genuine within-subject coupling. We propose a formal U-statistic-based test for such covariate effects, enabling both diagnostic evaluation of assumption violations and scientific discovery. Using synthetic and semi-synthetic neuroimaging data, we demonstrate well-controlled Type I error and high statistical power. We apply our method to test for confounding effects due to age and sex using real data from two pairs of imaging modalities in the Philadelphia Neurodevelopmental Cohort. Our framework increases the rigor and interpretability of intermodal coupling analyses, with broad implications for neuroimaging studies in heterogeneous populations, especially in developmental, aging, and disease-focused research.

Matching journals

The top 2 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.