Abstract:The aim of this work is to develop a spatial model for multi-subject fMRI data. While there has been much work on univariate modeling of each voxel for singleand multi-subject data, and some work on spatial modeling for single-subject data, there has been no work on spatial models that explicitly account for intersubject variability in activation location. We use a Bayesian hierarchical spatial model to fit the data. At the first level we model “population centers” that mark the centers of regions of activation. For a given population center each subject may have zero or more associated ”individual components”. While most previous work uses Gaussian mixtures for the activation shape, we instead use Gaussian mixtures for the probability that a voxel belongs to an activated region, assuming homogeneous mean intensity within a region. Our approach incorporates the unknown number of mixture components into the model as a parameter whose posterior distribution is estimated by reversible jump Markov Chain Monte Carlo. We demonstrate our method with a fMRI study of visual working memory and show dramatically better precision of localization with our method relative to the standard massunivariate method. Although we are motivated by fMRI data, this model could easily be modified to handle other types of imaging data. Bayesian Spatial Modeling of fMRI data: A Multiple-Subject Analysis Lei Xu1∗ Timothy D. Johnson1 Thomas E. Nichols2,3 1Department of Biostatistics School of Public Health University of Michigan Ann Arbor, MI 48109 2GlaxoSmithKline Clinical Imaging Centre Imperial College London Hammersmith Hospital, W12 0NN United Kingdom 3FMRIB Centre Oxford University Oxford, OX3 9DU United Kingdom Summary. The aim of this work is to develop a spatial model for multi-subject fMRI data. While there has been much work on univariate modeling of each voxel for singleand multi-subject data, and some work on spatial modeling for single-subject data, there has been no work on spatial models that explicitly account for intersubject variability in activation location. We use a Bayesian hierarchical spatial model to fit the data. At the first level we model “population centers” that mark the centers of regions of activation. For a given population center each subject may have zero or more associated “individual components”. While most previous work uses Gaussian mixtures for the activation shape, we instead use Gaussian mixtures for the probability that a voxel belongs to an activated region, assuming homogeneous mean intensity within a region. Our approach incorporates the unknown number of mixture components into the model as a parameter whose posterior distribution is estimated by reversible jump Markov Chain Monte Carlo. We demonstrate our method with a fMRI study of visual working memory and show dramatically better precision of localization with our method relative to the standard mass-univariate method. Although we are motivated by fMRI data, this model could easily be modified to handle other types of imaging data.

Novel Models for Multiple Dependent Heteroskedastic Time Series

Bayesian Analysis of fMRI data with Spatially-Varying Autoregressive Orders

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

Covariate-guided Bayesian mixture model for multivariate time series

Modelling inter-shot variability for robust temporal sub-sampling of dynamic, GABA-edited MR spectroscopy data

Prototype-Based Spatio-Temporal Probabilistic Modelling of fMRI Data.

Analysis of fMRI Time-Series Revisited

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

Evaluating stationarity via change-point alternatives with applications to fMRI data

Bayesian models for functional magnetic resonance imaging data analysis

Modelling variability in dynamic functional brain networks using embeddings

Regression Analysis of Asynchronous Longitudinal Functional and Scalar Data

A Bayesian non-stationary heteroskedastic time series model for multivariate critical care data

Exponential Family Graphical Models: Correlated Replicates and Unmeasured Confounders, with Applications to fMRI Data

Structured Estimation of Heterogeneous Time Series

Multiscale Adaptive Marginal Analysis of Longitudinal Neuroimaging Data with Time-Varying Covariates

Semi-parametric inference for large-scale data with temporally dependent noise

Bayesian Spatial Modeling of Fmri Data: A Multiple Subject Analysis.

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

Single-trial Variable Model for Event-Related Fmri Data Analysis