Abstract:Task functional magnetic resonance imaging (fMRI) is a type of neuroimaging data used to identify areas of the brain that activate during specific tasks or stimuli. These data are conventionally modeled using a massive univariate approach across all data locations, which ignores spatial dependence at the cost of model power. We previously developed and validated a spatial Bayesian model leveraging dependencies along the cortical surface of the brain in order to improve accuracy and power. This model utilizes stochastic partial differential equation spatial priors with sparse precision matrices to allow for appropriate modeling of spatially-dependent activations seen in the neuroimaging literature, resulting in substantial increases in model power. Our original implementation relies on the computational efficiencies of the integrated nested Laplace approximation (INLA) to overcome the computational challenges of analyzing high-dimensional fMRI data while avoiding issues associated with variational Bayes implementations. However, this requires significant memory resources, extra software, and software licenses to run. In this article, we develop an exact Bayesian analysis method for the general linear model, employing an efficient expectation-maximization algorithm to find maximum a posteriori estimates of task-based regressors on cortical surface fMRI data. Through an extensive simulation study of cortical surface-based fMRI data, we compare our proposed method to the existing INLA implementation, as well as a conventional massive univariate approach employing ad-hoc spatial smoothing. We also apply the method to task fMRI data from the Human Connectome Project and show that our proposed implementation produces similar results to the validated INLA implementation. Both the INLA and EM-based implementations are available through our open-source BayesfMRI R package.

An efficient approach to characterize spatio-temporal dependence in cortical surface fMRI data

A Bayesian General Linear Modeling Approach to Cortical Surface fMRI Data Analysis

A robust model for spatiotemporal dependencies

Fast Bayesian estimation of brain activation with cortical surface fMRI data using EM

Spurious correlations in surface-based functional brain imaging

Analysis of fMRI Time-Series Revisited

A novel method for sparse dynamic functional connectivity analysis from resting-state fMRI

Towards a predictive spatio-temporal representation of brain data

Fast Bayesian estimation of brain activation with cortical surface and subcortical fMRI data using EM

Self-similar correlation function in brain resting-state fMRI

A spatial mixture approach to inferring sub-ROI spatio-temporal patterns from rapid event-related fMRI data.

fMRI Signal Restoration Using a Spatio-Temporal Markov Random Field Preserving Transitions

A framework for registration, statistical characterization and classification of cortically constrained functional imaging data

Spatio-temporally separable non-linear latent factor learning: an application to somatomotor cortex fMRI data

A nonstationary nonparametric Bayesian approach to dynamically modeling effective connectivity in functional magnetic resonance imaging experiments

Time Series Analysis of fMRI Data: Spatial Modelling and Bayesian Computation

Recent Developments in Complex and Spatially Correlated Functional Data

fMRI-Kernel Regression: A Kernel-based Method for Pointwise Statistical Analysis of rs-fMRI for Population Studies

On inferring structural connectivity from brain functional-MRI data

Bayesian Analysis of fMRI data with Spatially-Varying Autoregressive Orders

Dynamic Adaptive Spatio-temporal Graph Convolution for fMRI Modelling