Detecting Brain Activations in Functional Magnetic Resonance Imaging (fMRI) Experiments with a Maximum Cross-Correlation Statistic

Abstract

Various statistical models have been proposed to analyze fMRI data. The usual goal is to make inferences about the eects that are related to an external stimulus. The primary focus of this paper is on those statisti- cal methods that enable one to detect 'signicantly activated' regions of the brain due to event-related stimuli. Most of these methods share a common property, requiring estimation of the hemodynamic response function (HRF) as part of the deterministic component of the statistical model. We propose and investigate a new approach that does not require HRF ts to detect 'activated' voxels. We argue that the method not only avoids tting a specic HRF, but still takes into account that the unknown response is delayed and smeared in time. This method also adapts to dierential re- sponses of the BOLD response across dierent brain regions and experimen- tal sessions. The maximum cross-correlation between the kernel-smoothed stimulus sequence and shifted (lagged) values of the observed response is the proposed test statistic. Using our recommended approach we show through realistic simulations and with real data that we obtain better sensitivity than simple correlation methods using default values of SPM2. The simulation experiment incor- porates dierent HRFs empirically determined from real data. The noise models are also dierent AR(3) ts and fractional Gaussians estimated from real data. We conclude that our proposed method is more powerful than simple correlation procedures, because of its robustness to variation in the HRF.