Cookies on this website
We use cookies to ensure that we give you the best experience on our website. If you click 'Continue' we'll assume that you are happy to receive all cookies and you won't see this message again. Click 'Find out more' for information on how to change your cookie settings.

FMRI modelling requires flexible haemodynamic response function (HRF) modelling, with the HRF being allowed to vary spatially and between subjects. To achieve this flexibility, voxelwise parameterised HRFs have been proposed; however, inference on such models is very slow. An alternative approach is to use basis functions allowing inference to proceed in the more manageable General Linear Model (GLM) framework. However, a large amount of the subspace spanned by the basis functions produces nonsensical HRF shapes. In this work we propose a technique for choosing a basis set, and then the means to constrain the subspace spanned by the basis set to only include sensible HRF shapes. Penny et al. showed how Variational Bayes can be used to infer on the GLM for FMRI. Here we extend the work of Penny et al. to give inference on the GLM with constrained HRF basis functions and with spatial Markov Random Fields on the autoregressive noise parameters. Constraining the subspace spanned by the basis set allows for far superior separation of activating voxels from nonactivating voxels in FMRI data. We use spatial mixture modelling to produce final probabilities of activation and demonstrate increased sensitivity on an FMRI dataset.

Original publication




Journal article



Publication Date





1748 - 1761


Oxford Centre for Functional Magnetic Resonance Imaging of the Brain, University of Oxford, John Radcliffe Hospital, Headington, Oxford OX3 9DU, UK.


Brain, Humans, Image Interpretation, Computer-Assisted, Magnetic Resonance Imaging, Brain Mapping, Models, Statistical, Linear Models, Bayes Theorem, Markov Chains, Sensitivity and Specificity, Pattern Recognition, Visual, Pain Threshold, Evoked Potentials, Mathematical Computing, Computer Simulation, Thermosensing