  Functional Magnetic Resonance Imaging (fMRI) maps cerebral activation in
response to stimuli but this activation is often difficult to detect,
especially in low-signal contexts and single-subject studies. Accurate
activation detection can be guided by the fact that very few voxels are, in
reality, truly activated and that these voxels are spatially localized, but it
is challenging to incorporate both these facts. We address these twin
challenges to single-subject and low-signal fMRI by developing a
computationally feasible and methodologically sound model-based approach,
implemented in the R package MixfMRI, that bounds the a priori expected
proportion of activated voxels while also incorporating spatial context. An
added benefit of our methodology is the ability to distinguish voxels and
regions having different intensities of activation. Our suggested approach is
evaluated in realistic two- and three-dimensional simulation experiments as
well as on multiple real-world datasets. Finally, the value of our suggested
approach in low-signal and single-subject fMRI studies is illustrated on a
sports imagination experiment that is often used to detect awareness and
improve treatment in patients in persistent vegetative state (PVS). Our ability
to reliably distinguish activation in this experiment potentially opens the
door to the adoption of fMRI as a clinical tool for the improved treatment and
therapy of PVS survivors and other patients.
