Some code for fMRI