We took extra care to investigate the possibility of pulsation artifacts confounding our results, by constructing fMRI regressors from BCG-artifact pulse timing. This timing was determined using a peak detection algorithm on temporal lobe EEG channels, which are the channels most strongly contaminated with BCG artifact. We fit this regressor to the functional data, and thresholded the resulting maps at an uncorrected p = 0.005. We do not report any results that overlap with these pulse correlates.
