Script demonstrating differences in GLM modelling approaches between SPM and FSL:
import numpy as np
from nipy.modalities.fmri.glm import GeneralLinearModel
# Simulated data (20 'rest' volumes followed by 20 'task' volumes)
Y = np.hstack((np.random.normal(4200, 50, size=20),
np.random.normal(4300, 50, size=20)))