Simultaneously acquired functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) data hold great promise for localizing the spatial source of epileptiform events detected in the EEG trace. Despite a number of studies applying this method, there has been no independent and systematic validation of the approach. The present study uses a nonparametric method to show that interictal discharges lead to a blood oxygen level dependent (BOLD) response that is significantly different to that obtained by examining random 'events'. We also use this approach to examine the optimization of analysis strategy for detecting these BOLD responses. Two patients with frequent epileptiform events and a healthy control were studied. The fMRI data for each patient were analyzed using a model derived from the timings of the epileptiform events detected on EEG during fMRI scanning. Twenty sets of random pseudoevents were used to generate a null distribution representing the level of chance correlation between the EEG events and fMRI data. The same pseudoevents were applied to control data. We demonstrate that it is possible to detect blood oxygen level-dependent (BOLD) changes related to interictal discharges with specific and independent knowledge about the reliability of this activation. Biologically generated events complicate the fMRI-EEG experiment. Our proposed validation examines whether identified events have an associated BOLD response beyond chance and allows optimization of analysis strategies. This is an important step beyond standard analysis. It informs clinical interpretation because it permits assessment of the reliability of the connection between interictal EEG events and the BOLD response to those events.