Using voxel-specific hemodynamic response function in EEG-fMRI data analysis: An estimation and detection model

Research groups who study epileptic spikes with simultaneous EEG-fMRI have used mostly the general linear model (GLM). A shortcoming of the GLM is that the specification of a simple hemodynamic response function (HRF) may lead to biased results. Other methods, which predict the hemodynamic response from the measured data, have been termed "recognition models". The merit of recognition models lies in the power of estimating the region-specific or voxel-specific HRF. We propose an approach that merges these two models in a general framework: estimate the HRF on the training data sets, and applying the estimated HRF on the other part of the data sets. The merit of this framework is that it can utilize the advantages of both models. A comparison of performance is made between the GLM with three fixed HRFs and the new model with voxel-specific HRFs. The main results are as follows: (1) in 18 of the 21 patients, the new model has a higher adjusted coefficient of multiple determination than the GLM with fixed HRF; (2) in some subjects, with the new model, we found areas of activation that had not been detected with the three fixed HRFs at our threshold of significance. The results suggest that the new model can do better than the fixed HRF GLM for the analysis of epileptic activity with EEG-fMRI.