A Bayesian hierarchical framework for spatial modeling of fMRI data

Applications of functional magnetic resonance imaging (fMRI) have provided novel insights into the neuropathophysiology of major psychiatric, neurological, and substance abuse disorders and their treatments. Modern activation studies often compare localized task-induced changes in brain activity between experimental groups. Complementary approaches consider the ensemble of voxels constituting an anatomically defined region of interest (ROI) or summary statistics, such as means or quantiles, of the ROI. In this work, we present a Bayesian extension of voxel-level analyses that offers several notable benefits. Among these, it combines whole-brain voxel-by-voxel modeling and ROI analyses within a unified framework. Secondly, an unstructured variance/covariance matrix for regional mean parameters allows for the study of inter-regional (long-range) correlations, and the model employs an exchangeable correlation structure to capture intra-regional (short-range) correlations. Estimation is performed using Markov Chain Monte Carlo (MCMC) techniques implemented via Gibbs sampling. We apply our Bayesian hierarchical model to two novel fMRI data sets: one considering inhibitory control in cocaine-dependent men and the second considering verbal memory in subjects at high risk for Alzheimer's disease.

[1]  N. L. Johnson,et al.  Multivariate Analysis , 1958, Nature.

[2]  Mark W. Woolrich,et al.  Constrained linear basis sets for HRF modelling using Variational Bayes , 2004, NeuroImage.

[3]  N. Volkow,et al.  The neural basis of addiction: a pathology of motivation and choice. , 2005, The American journal of psychiatry.

[4]  Karl J. Friston,et al.  Human Brain Function , 1997 .

[5]  R. Royall.  On the Probability of Observing Misleading Statistical Evidence , 2000 .

[6]  Nikos Makris,et al.  Automatically parcellating the human cerebral cortex. , 2004, Cerebral cortex.

[7]  H. Garavan,et al.  Dissociable Executive Functions in the Dynamic Control of Behavior: Inhibition, Error Detection, and Correction , 2002, NeuroImage.

[8]  R. Tibshirani,et al.  Empirical bayes methods and false discovery rates for microarrays , 2002, Genetic epidemiology.

[9]  T. Goldberg,et al.  Brain regions underlying response inhibition and interference monitoring and suppression , 2006, The European journal of neuroscience.

[10]  Mark S. Cohen,et al.  Patterns of brain activation in people at risk for Alzheimer's disease. , 2000, The New England journal of medicine.

[11]  Karl J. Friston,et al.  Bayesian fMRI time series analysis with spatial priors , 2005, NeuroImage.

[12]  Jeff Harrison,et al.  Applied Bayesian Forecasting and Time Series Analysis , 1994 .

[13]  Alan C. Evans,et al.  A Linear Spatial Correlation Model, with Applications to Positron Emission Tomography , 1991 .

[14]  James G. Scott,et al.  An exploration of aspects of Bayesian multiple testing , 2006 .

[15]  Marina Vannucci,et al.  Wavelet Thresholding with Bayesian False Discovery Rate Control , 2005, Biometrics.

[16]  R. Poldrack,et al.  Cortical and Subcortical Contributions to Stop Signal Response Inhibition: Role of the Subthalamic Nucleus , 2006, The Journal of Neuroscience.

[17]  Mark W. Woolrich,et al.  Fully Bayesian spatio-temporal modeling of FMRI data , 2004, IEEE Transactions on Medical Imaging.

[18]  Dietmar Cordes,et al.  Hierarchical clustering to measure connectivity in fMRI resting-state data. , 2002, Magnetic resonance imaging.

[19]  Korbinian Brodmann,et al.  Brodmann's localization in the cerebral cortex : the principles of comparative localisation in the cerebral cortex based on cytoarchitectonics , 2006 .

[20]  Rajan S. Patel,et al.  A Bayesian approach to determining connectivity of the human brain , 2006, Human brain mapping.

[21]  Rajan S. Patel,et al.  Determining hierarchical functional networks from auditory stimuli fMRI , 2006, Human brain mapping.

[22]  C Gössl,et al.  Bayesian Spatiotemporal Inference in Functional Magnetic Resonance Imaging , 2001, Biometrics.

[23]  B. Caffo,et al.  Familial risk for Alzheimer's disease alters fMRI activation patterns. , 2006, Brain : a journal of neurology.

[24]  Howard J Federoff,et al.  Current Treatment for Alzheimer Disease and Future Prospects , 2003, Alzheimer disease and associated disorders.

[25]  Jeffrey D Blume,et al.  Likelihood methods for measuring statistical evidence , 2002, Statistics in medicine.

[26]  F. Bowman,et al.  Spatiotemporal Models for Region of Interest Analyses of Functional Neuroimaging Data , 2007 .

[27]  P. Skudlarski,et al.  Detection of functional connectivity using temporal correlations in MR images , 2002, Human brain mapping.

[28]  R Brookmeyer,et al.  Projections of Alzheimer's disease in the United States and the public health impact of delaying disease onset. , 1998, American journal of public health.

[29]  Karl J. Friston,et al.  Statistical parametric maps in functional imaging: A general linear approach , 1994 .

[30]  P S Albert,et al.  A generalized estimating equations approach for spatially correlated binary data: applications to the analysis of neuroimaging data. , 1995, Biometrics.

[31]  Alan C. Evans,et al.  A general statistical analysis for fMRI data , 2000, NeuroImage.

[32]  D. Berry,et al.  Bayesian perspectives on multiple comparisons , 1999 .

[33]  N. Tzourio-Mazoyer,et al.  Automated Anatomical Labeling of Activations in SPM Using a Macroscopic Anatomical Parcellation of the MNI MRI Single-Subject Brain , 2002, NeuroImage.

[34]  H. Lüders,et al.  Functional connectivity in the human language system: a cortico-cortical evoked potential study. , 2004, Brain : a journal of neurology.

[35]  F. DuBois Bowman,et al.  Identifying spatial relationships in neural processing using a multiple classification approach , 2004, NeuroImage.

[36]  J D Watson,et al.  Nonparametric Analysis of Statistic Images from Functional Mapping Experiments , 1996, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[37]  R. Royall.  Statistical Evidence: A Likelihood Paradigm , 1997 .

[38]  F Dubois Bowman,et al.  Spatio-temporal modeling of localized brain activity. , 2005, Biostatistics.

[39]  Galin L. Jones,et al.  Fixed-Width Output Analysis for Markov Chain Monte Carlo , 2006, math/0601446.

[40]  Y. Miyashita,et al.  Common inhibitory mechanism in human inferior prefrontal cortex revealed by event-related functional MRI. , 1999, Brain : a journal of neurology.

[41]  Karl J. Friston,et al.  Functional Connectivity: The Principal-Component Analysis of Large (PET) Data Sets , 1993, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[42]  E. Taylor,et al.  Neuropsychological analyses of impulsiveness in childhood hyperactivity , 2001, British Journal of Psychiatry.

[43]  Yasumasa Matsuda,et al.  A Spatio-temporal Regression Model for the Analysis of Functional MRI Data , 2002, NeuroImage.

[44]  Karl J. Friston,et al.  Classical and Bayesian Inference in Neuroimaging: Theory , 2002, NeuroImage.

[45]  Alan C. Evans,et al.  A Three-Dimensional Statistical Analysis for CBF Activation Studies in Human Brain , 1992, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[46]  K. Worsley,et al.  Local Maxima and the Expected Euler Characteristic of Excursion Sets of χ 2, F and t Fields , 1994, Advances in Applied Probability.

[47]  A. Dale,et al.  Whole Brain Segmentation Automated Labeling of Neuroanatomical Structures in the Human Brain , 2002, Neuron.