Investigation of BOLD responses to epileptic spikes: a combined modeling and data analysis approach