Most existing analytical techniques for EEG-fMRI data need specific assumptions about the hemodynamic response function (HRF). These assumptions may not be appropriate when the HRF varies from subject to subject or from region to region. In this article, we introduce a deconvolution method for EEG-fMRI activation detection, which can be implemented with voxel-specific HRFs. A comparison of performance is made between three fixed HRFs and the deconvolution method under the framework of the general linear model. The main results are as follows: (1) the volume of detected regions from the deconvolved HRFs is larger. (2) In some subjects, the deconvolution technique can find areas of activation that have not been detected with the three fixed HRFs at our threshold of significance. (3) Deconvolution obtained higher adjusted coefficients of multiple determination compared to those obtained with the three fixed HRFs. The results suggest that the fixed HRF methods may not be the most appropriate for the analysis of epileptic activity with EEG-fMRI, and the deconvolution method may be a better choice.