Likelihood-based hypothesis tests for brain activation detection from MRI data disturbed by colored noise: a simulation study

Adaptive Brain Activity Detection in Structured Interference and Partially Homogeneous Locally Correlated Disturbance

OBJECTIVE

In this paper, we aim to address the problem of subspace detection in the presence of locally-correlated complex Gaussian noise and interference. For applications like brain activity detection using functional magnetic resonance imaging (fMRI) data where the noise is possibly locally correlated, using the sample covariance estimator is not a suitable choice due to significant dependency of its accuracy on the number of observations.

METHODS

In this study, we take advantage of an assumed banded structure in the covariance matrix to model the local dependence in the noise and propose a new covariance estimation approach. In particular, we use the idea of fac-torizing the joint likelihood function into a few conditional likelihood terms and maximizing each term independently of the others. This process leads to an explicit estimator for banded covariance matrices which requires fewer observations to achieve the same accuracy as the sample covari-ance. This estimate is then fed into an adaptive matched filter, two-step Rao and two-step Wald tests for detection.

RESULTS

Simulation results reveal the superiority of the proposed methods over well known classical detectors. Finally, the proposed methods are applied to functional magnetic resonance imaging (fMRI) data to localize neural activities in the brain.

CONCLUSION

The proposed method can offer better activation maps in terms of accuracy and spatial smoothness.

SIGNIFICANCE

The proposed methods can be seen as alternatives for standard detection approaches which are not perfectly aligned with the properties of fMRI data.

A decomposition model and voxel selection framework for fMRI analysis to predict neural response of visual stimuli

This paper presents an fMRI signal analysis methodology using geometric mean curve decomposition (GMCD) and mutual information-based voxel selection framework. Previously, the fMRI signal analysis has been conducted using empirical mean curve decomposition (EMCD) model and voxel selection on raw fMRI signal. The erstwhile methodology loses frequency component, while the latter methodology suffers from signal redundancy. Both challenges are addressed by our methodology in which the frequency component is considered by decomposing the raw fMRI signal using geometric mean rather than arithmetic mean and the voxels are selected from EMCD signal using GMCD components, rather than raw fMRI signal. The proposed methodologies are adopted for predicting the neural response. Experimentations are conducted in the openly available fMRI data of six subjects, and comparisons are made with existing decomposition models and voxel selection frameworks. Subsequently, the effect of degree of selected voxels and the selection constraints are analyzed. The comparative results and the analysis demonstrate the superiority and the reliability of the proposed methodology.

On the definition of signal-to-noise ratio and contrast-to-noise ratio for FMRI data

Signal-to-noise ratio, the ratio between signal and noise, is a quantity that has been well established for MRI data but is still subject of ongoing debate and confusion when it comes to fMRI data. fMRI data are characterised by small activation fluctuations in a background of noise. Depending on how the signal of interest and the noise are identified, signal-to-noise ratio for fMRI data is reported by using many different definitions. Since each definition comes with a different scale, interpreting and comparing signal-to-noise ratio values for fMRI data can be a very challenging job. In this paper, we provide an overview of existing definitions. Further, the relationship with activation detection power is investigated. Reference tables and conversion formulae are provided to facilitate comparability between fMRI studies.

Fractional-order time series models for extracting the haemodynamic response from functional magnetic resonance imaging data

The postprocessing of functional magnetic resonance imaging (fMRI) data to study the brain functions deals mainly with two objectives: signal detection and extraction of the haemodynamic response. Signal detection consists of exploring and detecting those areas of the brain that are triggered due to an external stimulus. Extraction of the haemodynamic response deals with describing and measuring the physiological process of activated regions in the brain due to stimulus. The haemodynamic response represents the change in oxygen levels since the brain functions require more glucose and oxygen upon stimulus that implies a change in blood flow. In the literature, different approaches to estimate and model the haemodynamic response have been proposed. These approaches can be discriminated in model structures that either provide a proper representation of the obtained measurements but provide no or a limited amount of physiological information, or provide physiological insight but lacks a proper fit to the data. In this paper, a novel model structure is studied for describing the haemodynamics in fMRI measurements: fractional models. We show that these models are flexible enough to describe the gathered data with the additional merit of providing physiological information.

A generalized form of context-dependent psychophysiological interactions (gPPI): a comparison to standard approaches

Functional MRI (fMRI) allows one to study task-related regional responses and task-dependent connectivity analysis using psychophysiological interaction (PPI) methods. The latter affords the additional opportunity to understand how brain regions interact in a task-dependent manner. The current implementation of PPI in Statistical Parametric Mapping (SPM8) is configured primarily to assess connectivity differences between two task conditions, when in practice fMRI tasks frequently employ more than two conditions. Here we evaluate how a generalized form of context-dependent PPI (gPPI; http://www.nitrc.org/projects/gppi), which is configured to automatically accommodate more than two task conditions in the same PPI model by spanning the entire experimental space, compares to the standard implementation in SPM8. These comparisons are made using both simulations and an empirical dataset. In the simulated dataset, we compare the interaction beta estimates to their expected values and model fit using the Akaike information criterion (AIC). We found that interaction beta estimates in gPPI were robust to different simulated data models, were not different from the expected beta value, and had better model fits than when using standard PPI (sPPI) methods. In the empirical dataset, we compare the model fit of the gPPI approach to sPPI. We found that the gPPI approach improved model fit compared to sPPI. There were several regions that became non-significant with gPPI. These regions all showed significantly better model fits with gPPI. Also, there were several regions where task-dependent connectivity was only detected using gPPI methods, also with improved model fit. Regions that were detected with all methods had more similar model fits. These results suggest that gPPI may have greater sensitivity and specificity than standard implementation in SPM. This notion is tempered slightly as there is no gold standard; however, data simulations with a known outcome support our conclusions about gPPI. In sum, the generalized form of context-dependent PPI approach has increased flexibility of statistical modeling, and potentially improves model fit, specificity to true negative findings, and sensitivity to true positive findings.