Merlin Keller, Alexis Roche, Alan Tucholka and Bertrand Thirion (2008). Dealing with spatial normalization errors in fMRI group inference using hierarchical modeling. Vol. 18, No. 4, 1357-1374.

Statistica Sinica 18(2008), 1357-1374

DEALING WITH SPATIAL NORMALIZATION ERRORS

IN fMRI GROUP INFERENCE

USING HIERARCHICAL MODELING

Merlin Keller $^{1, 2}$ , Alexis Roche

, Alan Tucholka $^{1, 2}$

and Bertrand Thirion $^{1, 2}$

CEA, Neurospin, Gif-sur-Yvette, France and INRIA, Saclay, France

Abstract: An important challenge in neuroimaging multi-subject studies is to take into account that different brains cannot be aligned perfectly. To this end, we extend the classical mass univariate model for group analysis to incorporate uncertainty on localization by introducing, for each subject, a spatial ``jitter'' variable to be marginalized out. We derive a Bayes factor to test for the mean population effect's sign in each voxel of a search volume, and discuss a Gibbs sampler to compute it. This Bayes factor, which generalizes the classical -statistic, may be combined with a permutation test in order to control the frequentist false positive rate. Results on both simulated and experimental data suggest that this test may outperform conventional mass univariate tests in terms of detection power, while limiting the problem of overestimating the size of activity clusters.

Key words and phrases: Group analysis, hierarchical modeling, mixed effects, spatial uncertainty, Bayes factor, Metropolis within Gibbs, permutation test.