Bayesian hierarchical multi-subject multiscale analysis of functional MRI data

We develop a methodology for Bayesian hierarchical multi-subject multiscale analysis of functional Magnetic Resonance Imaging (fMRI) data. We begin by modeling the brain images temporally with a standard general linear model. After that, we transform the resulting estimated standardized regression coefficient maps through a discrete wavelet transformation to obtain a sparse representation in the wavelet space. Subsequently, we assign to the wavelet coefficients a prior that is a mixture of a point mass at zero and a Gaussian white noise. In this mixture prior for the wavelet coefficients, the mixture probabilities are related to the pattern of brain activity across different resolutions. To incorporate this information, we assume that the mixture probabilities for wavelet coefficients at the same location and level are common across subjects. Furthermore, we assign for the mixture probabilities a prior that depends on a few hyperparameters. We develop an empirical Bayes methodology to estimate the hyperparameters and, as these hyperparameters are shared by all subjects, we obtain precise estimated values. Then we carry out inference in the wavelet space and obtain smoothed images of the regression coefficients by applying the inverse wavelet transform to the posterior means of the wavelet coefficients. An application to computer simulated synthetic data has shown that, when compared to single-subject analysis, our multi-subject methodology performs better in terms of mean squared error. Finally, we illustrate the utility and flexibility of our multi-subject methodology with an application to an event-related fMRI dataset generated by Postle (2005) through a multi-subject fMRI study of working memory related brain activation.

[1]  Bradley R. Postle,et al.  Delay-period Activity in the Prefrontal Cortex: One Function Is Sensory Gating , 2005, Journal of Cognitive Neuroscience.

[2]  Herbert K. H. Lee,et al.  Multiscale Modeling: A Bayesian Perspective , 2007 .

[3]  Sadasivan Puthusserypady,et al.  Analysis of fMRI Data With Drift: Modified General Linear Model and Bayesian Estimator , 2008, IEEE Transactions on Biomedical Engineering.

[4]  Emery N. Brown,et al.  Spatiotemporal wavelet analysis for functional MRI , 2004, NeuroImage.

[5]  François G. Meyer Wavelet-based estimation of a semiparametric generalized linear model of fMRI time-series , 2003, IEEE Transactions on Medical Imaging.

[6]  Marco A. R. Ferreira,et al.  Dynamic multiscale spatiotemporal models for Gaussian areal data , 2011 .

[7]  Jos B. T. M. Roerdink,et al.  Denoising functional MR images: a comparison of wavelet denoising and Gaussian smoothing , 2004, IEEE Transactions on Medical Imaging.

[8]  Lee M. Harrison,et al.  A Bayesian spatiotemporal model for very large data sets , 2010, NeuroImage.

[9]  Bradley P. Carlin,et al.  Bayesian Methods for Data Analysis , 2008 .

[10]  L. Fahrmeir,et al.  Spatial Bayesian Variable Selection With Application to Functional Magnetic Resonance Imaging , 2007 .

[11]  Dani Gamerman,et al.  Bayesian spatiotemporal model of fMRI data , 2010, NeuroImage.

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

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

[14]  William D. Penny,et al.  Bayesian fMRI data analysis with sparse spatial basis function priors , 2007, NeuroImage.

[15]  Laurent Risser,et al.  Spatially adaptive mixture modeling for analysis of fMRI time series , 2009, NeuroImage.

[16]  Thierry Blu,et al.  WSPM or How to Obtain Statistical Parametric Maps using Shift-Invariant Wavelet Processing , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[17]  Brian Caffo,et al.  A Bayesian hierarchical framework for spatial modeling of fMRI data , 2008, NeuroImage.

[18]  S. Zaroubi,et al.  Complex denoising of MR data via wavelet analysis: application for functional MRI. , 2000, Magnetic resonance imaging.

[19]  Michael S. Beauchamp,et al.  FMRI group analysis combining effect estimates and their variances , 2012, NeuroImage.

[20]  Cyril Riddell,et al.  A linear wavelet filter for parametric imaging with dynamic PET , 2003, IEEE Transactions on Medical Imaging.

[21]  R W Cox,et al.  AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. , 1996, Computers and biomedical research, an international journal.

[22]  M. Vannucci,et al.  Covariance structure of wavelet coefficients: theory and models in a Bayesian perspective , 1999 .

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

[24]  Simon P. Wilson,et al.  Bayesian spatiotemporal model of fMRI data using transfer functions , 2010, NeuroImage.

[25]  Dimitri Van De Ville,et al.  Integrated wavelet processing and spatial statistical testing of fMRI data , 2004, NeuroImage.

[26]  B. Whitcher,et al.  Multiple hypothesis mapping of functional MRI data in orthogonal and complex wavelet domains , 2005, IEEE Transactions on Signal Processing.

[27]  Yodchanan Wongsawat Trends on Wavelet-Based Functional MRI for Activation Detection , 2009 .

[28]  Mark W. Woolrich,et al.  Mixture models with adaptive spatial regularization for segmentation with an application to FMRI data , 2005, IEEE Transactions on Medical Imaging.

[29]  Nouna Kettaneh,et al.  Statistical Modeling by Wavelets , 1999, Technometrics.

[30]  Ludwig Fahrmeir,et al.  Assessing brain activity through spatial bayesian variable selection , 2003, NeuroImage.

[31]  Gholam-Ali Hossein-Zadeh,et al.  Fixed and random effect analysis of multi-subject fMRI data using wavelet transform , 2009, Journal of Neuroscience Methods.

[32]  Karl J. Friston,et al.  Variational Bayesian inference for fMRI time series , 2003, NeuroImage.

[33]  Gareth J. Barker,et al.  Bayesian wavelet-based analysis of functional magnetic resonance time series. , 2009, Magnetic resonance imaging.

[34]  R Baumgartner,et al.  Wavelet domain de-noising of time-courses in MR image sequences. , 2000, Magnetic resonance imaging.

[35]  Bertrand Thirion,et al.  A fully Bayesian approach to the parcel-based detection-estimation of brain activity in fMRI , 2008, NeuroImage.

[36]  Ghassan Hamarneh,et al.  Detecting Brain Activation in fMRI Using Group Random Walker , 2010, MICCAI.

[37]  Michael Unser,et al.  Statistical analysis of functional MRI data in the wavelet domain , 1998, IEEE Transactions on Medical Imaging.

[38]  Robert Tibshirani,et al.  The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd Edition , 2001, Springer Series in Statistics.

[39]  Rainer Hinz,et al.  Multi-resolution Bayesian regression in PET dynamic studies using wavelets , 2006, NeuroImage.

[40]  Mark W. Woolrich,et al.  Combined spatial and non-spatial prior for inference on MRI time-series , 2009, NeuroImage.

[41]  E. Bullmore,et al.  Wavelet-Generalized Least Squares: A New BLU Estimator of Linear Regression Models with 1/f Errors , 2002, NeuroImage.

[42]  Juan Zhou,et al.  Learning effective brain connectivity with dynamic Bayesian networks , 2007, NeuroImage.

[43]  M. Unser,et al.  Surfing the brain , 2006, IEEE Engineering in Medicine and Biology Magazine.

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

[45]  M. Clyde,et al.  Multiple shrinkage and subset selection in wavelets , 1998 .

[46]  Karl J. Friston,et al.  Diffusion-based spatial priors for functional magnetic resonance images , 2008, NeuroImage.

[47]  Jeffrey S. Morris,et al.  AUTOMATED ANALYSIS OF QUANTITATIVE IMAGE DATA USING ISOMORPHIC FUNCTIONAL MIXED MODELS, WITH APPLICATION TO PROTEOMICS DATA. , 2011, The annals of applied statistics.

[48]  Christopher R. Genovese,et al.  A Bayesian Time-Course Model for Functional Magnetic Resonance Imaging Data , 2000 .

[49]  William H. Press,et al.  Numerical recipes in C , 2002 .

[50]  B. Silverman,et al.  Wavelet thresholding via a Bayesian approach , 1998 .

[51]  E. Bullmore,et al.  Wavelets and functional magnetic resonance imaging of the human brain , 2004, NeuroImage.

[52]  Thomas Vincent,et al.  Impact of the joint detection-estimation approach on random effects group studies in FMRI , 2011, 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[53]  Ghassan Hamarneh,et al.  Group MRF for fMRI activation detection , 2010, 2010 IEEE Computer Society Conference on Computer Vision and Pattern Recognition.