Modeling of region-specific fMRI BOLD neurovascular response functions in rat brain reveals residual differences that correlate with the differences in regional evoked potentials.