Space-time statistical model for functional MRI image sequences. In: Duncan J, Gindi G, editors. Information Processing in Medical Imaging. Berlin: Springer Berlin Heidelberg; 1997. pp. 285-98. [DOI: 10.1007/3-540-63046-5_22]

A Spatio-temporal Model for Detecting the Effect of Cocaine Use Disorder on Functional Connectivity

Drug addiction can lead to many health-related problems and social concerns. Researchers are interested in the association between long-term drug usage and abnormal functional connectivity. Functional connectivity obtained from functional magnetic resonance imaging data promotes a variety of fundamental understandings in such association. Due to the complex correlation structure and large dimensionality, the modeling and analysis of the functional connectivity from neuroimage are challenging. By proposing a spatio-temporal model for multi-subject neuroimage data, we incorporate voxel-level spatio-temporal dependencies of whole-brain measurements to improve the accuracy of statistical inference. To tackle large-scale spatio-temporal neuroimage data, we develop a computational efficient algorithm to estimate the parameters. Our method is used to first identify functional connectivity, and then detect the effect of cocaine use disorder (CUD) on functional connectivity between different brain regions. The functional connectivity identified by our spatio-temporal model matches existing studies on brain networks, and further indicates that CUD may alter the functional connectivity in the medial orbitofrontal cortex subregions and the supplementary motor areas.

Modeling the spatial and temporal dependence in FMRI data

Functional magnetic resonance imaging (fMRI) data sets are large and characterized by complex dependence structures driven by highly sophisticated neurophysiology and aspects of the experimental designs. Typical analyses investigating task-related changes in measured brain activity use a two-stage procedure in which the first stage involves subject-specific models and the second-stage specifies group (or population) level parameters. Customarily, the first-level accounts for temporal correlations between the serial scans acquired during one scanning session. Despite accounting for these correlations, fMRI studies often include multiple sessions and temporal dependencies may persist between the corresponding estimates of mean neural activity. Further, spatial correlations between brain activity measurements in different locations are often unaccounted for in statistical modeling and estimation. We propose a two-stage, spatio-temporal, autoregressive model that simultaneously accounts for spatial dependencies between voxels within the same anatomical region and for temporal dependencies between a subject's estimates from multiple sessions. We develop an algorithm that leverages the special structure of our covariance model, enabling relatively fast and efficient estimation. Using our proposed method, we analyze fMRI data from a study of inhibitory control in cocaine addicts.

Unified framework for robust estimation of brain networks from FMRI using temporal and spatial correlation analyses

There is a rapidly growing interest in the neuroimaging field to use functional magnetic resonance imaging (fMRI) to explore brain networks, i.e., how regions of the brain communicate with one another. This paper presents a general and novel statistical framework for robust and more complete estimation of brain functional connectivity from fMRI based on correlation analyses and hypothesis testing. In addition to the ability of examining the correlations with each individual seed as in the standard and existing methods, the proposed framework can detect functional interactions by simultaneously examining multiseed correlations via multiple correlation coefficients. Spatially structured noise in fMRI is also taken into account during the identification of functional interconnection networks through noncentral F hypothesis tests. The associated issues for the multiple testing and the effective degrees-of-freedom are considered as well. Furthermore, partial multiple correlations are introduced and formulated to measure any additional task-induced but not stimulus-locked relation over brain regions so that we can take the analysis of functional connectivity closer to the characterization of direct functional interactions of the brain. Evaluation for accuracy and advantages, and comparisons of the new approaches in the presented general framework are performed using both realistic synthetic data and in vivo fMRI data.

Fully Bayesian spatio-temporal modeling of FMRI data

We present a fully Bayesian approach to modeling in functional magnetic resonance imaging (FMRI), incorporating spatio-temporal noise modeling and haemodynamic response function (HRF) modeling. A fully Bayesian approach allows for the uncertainties in the noise and signal modeling to be incorporated together to provide full posterior distributions of the HRF parameters. The noise modeling is achieved via a nonseparable space-time vector autoregressive process. Previous FMRI noise models have either been purely temporal, separable or modeling deterministic trends. The specific form of the noise process is determined using model selection techniques. Notably, this results in the need for a spatially nonstationary and temporally stationary spatial component. Within the same full model, we also investigate the variation of the HRF in different areas of the activation, and for different experimental stimuli. We propose a novel HRF model made up of half-cosines, which allows distinct combinations of parameters to represent characteristics of interest. In addition, to adaptively avoid over-fitting we propose the use of automatic relevance determination priors to force certain parameters in the model to zero with high precision if there is no evidence to support them in the data. We apply the model to three datasets and observe matter-type dependence of the spatial and temporal noise, and a negative correlation between activation height and HRF time to main peak (although we suggest that this apparent correlation may be due to a number of different effects).

Comparison of filtering methods for fMRI datasets

When studying complex cognitive tasks using functional magnetic resonance imaging (fMRI) one often encounters weak signal responses. These weak responses are corrupted by noise and artifacts of various sources. Preprocessing of the raw data before the application of test statistics helps to extract the signal and can vastly improve signal detection. Artifact sources and algorithms to handle them are discussed. In an empirical approach targeted to yield an optimal recovery of the hemodynamic response, we implemented a test bed for baseline correction and noise-filtering methods. A known signal is modulated onto foreground patches obtained from event-related fMRI experiments. Quantitative performance measures are defined to optimize the characteristics of a given filter and to compare their results. Marked improvements in the sensitivity and selectivity are achieved by optimized filtering. Examples using real data underline the usefulness of this preprocessing sequence.