Implications of the Rician distribution for fMRI generalized likelihood ratio tests

Nonlinear Bayesian estimation of BOLD signal under non-Gaussian noise

Modeling the blood oxygenation level dependent (BOLD) signal has been a subject of study for over a decade in the neuroimaging community. Inspired from fluid dynamics, the hemodynamic model provides a plausible yet convincing interpretation of the BOLD signal by amalgamating effects of dynamic physiological changes in blood oxygenation, cerebral blood flow and volume. The nonautonomous, nonlinear set of differential equations of the hemodynamic model constitutes the process model while the weighted nonlinear sum of the physiological variables forms the measurement model. Plagued by various noise sources, the time series fMRI measurement data is mostly assumed to be affected by additive Gaussian noise. Though more feasible, the assumption may cause the designed filter to perform poorly if made to work under non-Gaussian environment. In this paper, we present a data assimilation scheme that assumes additive non-Gaussian noise, namely, the e-mixture noise, affecting the measurements. The proposed filter MAGSF and the celebrated EKF are put to test by performing joint optimal Bayesian filtering to estimate both the states and parameters governing the hemodynamic model under non-Gaussian environment. Analyses using both the synthetic and real data reveal superior performance of the MAGSF as compared to EKF.

Ricean over Gaussian modelling in magnitude fMRI Analysis-Added Complexity with Negligible Practical Benefits

It is well-known that Gaussian modeling 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 ROC curve methodology show that tests derived using Ricean modeling are substantially superior over the Gaussian-based activation tests only for SNRs below 0.6, i.e SNR values far lower than those encountered in fMRI as currently practiced.

Bayesian estimation of changes in transverse relaxation rates

Although the biasing of R(2)* estimates by assuming magnitude MR data to be normally distributed has been described, the effect on changes in R(2)* (DeltaR(2)*), such as induced by a paramagnetic contrast agent, has not been reported. In this study, two versions of a novel Bayesian maximum a posteriori approach for estimating DeltaR(2)* are described and evaluated: one that assumes normally distributed data and the other, Rice-distributed data. The approach enables the robust, voxelwise determination of the uncertainty in DeltaR(2)* estimates and provides a useful statistical framework for quantifying the probability that a pixel has been significantly enhanced. This technique was evaluated in vivo, using ultrasmall superparamagnetic iron oxide particles in orthotopic murine prostate tumors. It is shown that assuming magnitude data to be normally distributed causes DeltaR(2)* to be underestimated when signal-to-noise ratio is modest. However, the biasing effect is less than is found in R(2)* estimates, implying that the simplifying assumption of normally distributed noise is more justifiable when evaluating DeltaR(2)* compared with when evaluating precontrast R(2)* values.

A short note on the analysis of distance measurements by electron paramagnetic resonance

In electron paramagnetic resonance (EPR) distance distributions between site-directedly attached spin labels in soft matter are obtained by measuring their dipole-dipole interaction. The analysis of these distance distributions can be misleading particularly for broad distributions, because the most probable distance deviates from the distance between the most probable label positions. The current manuscript studies this effect using numerically generated spin label positions, molecular dynamics simulations, and experimental data of a model system. An approach involving Rice distributions is proposed to overcome this problem.

Regression Models for Identifying Noise Sources in Magnetic Resonance Images

Stochastic noise, susceptibility artifacts, magnetic field and radiofrequency inhomogeneities, and other noise components in magnetic resonance images (MRIs) can introduce serious bias into any measurements made with those images. We formally introduce three regression models including a Rician regression model and two associated normal models to characterize stochastic noise in various magnetic resonance imaging modalities, including diffusion-weighted imaging (DWI) and functional MRI (fMRI). Estimation algorithms are introduced to maximize the likelihood function of the three regression models. We also develop a diagnostic procedure for systematically exploring MR images to identify noise components other than simple stochastic noise, and to detect discrepancies between the fitted regression models and MRI data. The diagnostic procedure includes goodness-of-fit statistics, measures of influence, and tools for graphical display. The goodness-of-fit statistics can assess the key assumptions of the three regression models, whereas measures of influence can isolate outliers caused by certain noise components, including motion artifacts. The tools for graphical display permit graphical visualization of the values for the goodness-of-fit statistic and influence measures. Finally, we conduct simulation studies to evaluate performance of these methods, and we analyze a real dataset to illustrate how our diagnostic procedure localizes subtle image artifacts by detecting intravoxel variability that is not captured by the regression models.

BOLD Noise Assumptions in fMRI

This paper discusses the assumption of Gaussian noise in the blood-oxygenation-dependent (BOLD) contrast for functional MRI (fMRI). In principle, magnitudes in MRI images follow a Rice distribution. We start by reviewing differences between Rician and Gaussian noise. An analytic expression is derived for the null (resting-state) distribution of the difference between two Rician distributed images. This distribution is shown to be symmetric, and an exact expression for its standard deviation is derived. This distribution can be well approximated by a Gaussian, with very high precision for high SNR, and high precision for lower SNR. Tests on simulated and real MR images show that subtracting the time-series mean in fMRI yields asymmetrically distributed temporal noise. Subtracting a resting-state time series from the first results in symmetric and nearly Gaussian noise. This has important consequences for fMRI analyses using standard statistical tests.