We studied single-event and average BOLD responses to EEG interictal epileptic discharges (IEDs) in four patients with focal epilepsy, using continuous EEG-fMRI during 80-min sessions. The detection of activated areas was performed by comparing the BOLD signal at each voxel to a model of the expected signal. Since little is known about the BOLD response to IEDs, we modeled it with the response to brief auditory events (G. H., NeuroImage 9, 416-429). For each activated area, we then obtained the time course of the BOLD signal for the complete session and computed the actual average hemodynamic response function (HRF) to IEDs. In two of four patients, we observed clear BOLD responses to single IEDs. The average response was composed of a positive lobe peaking between 6 and 7 s in all patients and a negative undershoot in three patients. There were important variations in amplitude and shape between average HRFs across patients. The average HRF presented a wider positive lobe than the Glover model in three patients and a longer undershoot in two. There was a remarkable similarity in the shape of the HRF across areas for patients presenting multiple activation sites. There was no clear correlation between the amplitude of individual BOLD responses and the amplitude of the corresponding EEG spike. The possibility of a longer HRF could be used to improve statistical detection of activation in simultaneous EEG-fMRI. The variability in average HRFs across patients could reflect in part different pathophysiological mechanisms.