Ricean over Gaussian modelling in magnitude fMRI analysis—added complexity with negligible practical benefits

Abstract

<p>It is well known that Gaussian modelling of functional magnetic resonance imaging (fMRI) magnitude time-course data, which are truly Rice distributed, constitutes an approximation, especially at low signal-to-noise ratios (SNRs). Based on this fact, previous work has argued that Rice-based activation tests show superior performance over their Gaussian-based counterparts at low SNRs and should be preferred in spite of the attendant additional computational and estimation burden. Here, we revisit these past studies and, after identifying and removing their underlying limiting assumptions and approximations, provide a more comprehensive comparison. Our experimental evaluations using Receiver Operating Characteristic (ROC) curve methodology show that tests derived using Ricean modelling are substantially superior over the Gaussian-based activation tests only for SNRs below 0.6, that is, SNR values far lower than those encountered in fMRI as currently practiced</p>