TY - JOUR
T1 - A spatio-temporal nonparametric Bayesian variable selection model of fMRI data for clustering correlated time courses
AU - Zhang, Linlin
AU - Guindani, Michele
AU - Versace, Francesco
AU - Vannucci, Marina
N1 - Funding Information:
MG and FV were partially supported by the University of Texas MD Anderson Cancer Center's Cancer Center Support Grant ( NIH CA016672 ).
PY - 2014/7/15
Y1 - 2014/7/15
N2 - In this paper we present a novel wavelet-based Bayesian nonparametric regression model for the analysis of functional magnetic resonance imaging (fMRI) data. Our goal is to provide a joint analytical framework that allows to detect regions of the brain which exhibit neuronal activity in response to a stimulus and, simultaneously, infer the association, or clustering, of spatially remote voxels that exhibit fMRI time series with similar characteristics. We start by modeling the data with a hemodynamic response function (HRF) with a voxel-dependent shape parameter. We detect regions of the brain activated in response to a given stimulus by using mixture priors with a spike at zero on the coefficients of the regression model. We account for the complex spatial correlation structure of the brain by using a Markov random field (MRF) prior on the parameters guiding the selection of the activated voxels, therefore capturing correlation among nearby voxels. In order to infer association of the voxel time courses, we assume correlated errors, in particular long memory, and exploit the whitening properties of discrete wavelet transforms. Furthermore, we achieve clustering of the voxels by imposing a Dirichlet process (DP) prior on the parameters of the long memory process. For inference, we use Markov Chain Monte Carlo (MCMC) sampling techniques that combine Metropolis-Hastings schemes employed in Bayesian variable selection with sampling algorithms for nonparametric DP models. We explore the performance of the proposed model on simulated data, with both block- and event-related design, and on real fMRI data.
AB - In this paper we present a novel wavelet-based Bayesian nonparametric regression model for the analysis of functional magnetic resonance imaging (fMRI) data. Our goal is to provide a joint analytical framework that allows to detect regions of the brain which exhibit neuronal activity in response to a stimulus and, simultaneously, infer the association, or clustering, of spatially remote voxels that exhibit fMRI time series with similar characteristics. We start by modeling the data with a hemodynamic response function (HRF) with a voxel-dependent shape parameter. We detect regions of the brain activated in response to a given stimulus by using mixture priors with a spike at zero on the coefficients of the regression model. We account for the complex spatial correlation structure of the brain by using a Markov random field (MRF) prior on the parameters guiding the selection of the activated voxels, therefore capturing correlation among nearby voxels. In order to infer association of the voxel time courses, we assume correlated errors, in particular long memory, and exploit the whitening properties of discrete wavelet transforms. Furthermore, we achieve clustering of the voxels by imposing a Dirichlet process (DP) prior on the parameters of the long memory process. For inference, we use Markov Chain Monte Carlo (MCMC) sampling techniques that combine Metropolis-Hastings schemes employed in Bayesian variable selection with sampling algorithms for nonparametric DP models. We explore the performance of the proposed model on simulated data, with both block- and event-related design, and on real fMRI data.
KW - Bayesian nonparametric
KW - Dirichlet process prior
KW - Discrete wavelet transform
KW - FMRI
KW - Long memory errors
KW - Markov random field prior
UR - http://www.scopus.com/inward/record.url?scp=84899692109&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84899692109&partnerID=8YFLogxK
U2 - 10.1016/j.neuroimage.2014.03.024
DO - 10.1016/j.neuroimage.2014.03.024
M3 - Article
C2 - 24650600
AN - SCOPUS:84899692109
SN - 1053-8119
VL - 95
SP - 162
EP - 175
JO - NeuroImage
JF - NeuroImage
ER -