Bayesian deconvolution fMRI data using bilinear dynamical systems

In Penny et al. [Penny, W., Ghahramani, Z., Friston, K.J. 2005. Bilinear dynamical systems. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360(1457) 983-993], a particular case of the Linear Dynamical Systems (LDSs) was used to model the dynamic behavior of the BOLD response in functional MRI. This state-space model, called bilinear dynamical system (BDS), is used to deconvolve the fMRI time series in order to estimate the neuronal response induced by the different stimuli of the experimental paradigm. The BDS model parameters are estimated using an expectation-maximization (EM) algorithm proposed by Ghahramani and Hinton [Ghahramani, Z., Hinton, G.E. 1996. Parameter Estimation for Linear Dynamical Systems. Technical Report, Department of Computer Science, University of Toronto]. In this paper we introduce modifications to the BDS model in order to explicitly model the spatial variations of the haemodynamic response function (HRF) in the brain using a non-parametric approach. While in Penny et al. [Penny, W., Ghahramani, Z., Friston, K.J. 2005. Bilinear dynamical systems. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360(1457) 983-993] the relationship between neuronal activation and fMRI signals is formulated as a first-order convolution with a kernel expansion using basis functions (typically two or three), in this paper, we argue in favor of a spatially adaptive GLM in which a local non-parametric estimation of the HRF is performed. Furthermore, in order to overcome the overfitting problem typically associated with simple EM estimates, we propose a full Variational Bayes (VB) solution to infer the BDS model parameters. We demonstrate the usefulness of our model which is able to estimate both the neuronal activity and the haemodynamic response function in every voxel of the brain. We first examine the behavior of this approach when applied to simulated data with different temporal and noise features. As an example we will show how this method can be used to improve interpretability of estimates from an independent component analysis (ICA) analysis of fMRI data. We finally demonstrate its use on real fMRI data in one slice of the brain.

[1]  Matthew J. Beal,et al.  The variational Bayesian EM algorithm for incomplete data: with application to scoring graphical model structures , 2003 .

[2]  P. Jezzard,et al.  Sources of distortion in functional MRI data , 1999, Human brain mapping.

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

[4]  M. D’Esposito,et al.  The variability of human BOLD hemodynamic responses , 1998, NeuroImage.

[5]  Karl J. Friston,et al.  Multivariate Autoregressive Modelling of fMRI time series , 2003 .

[6]  J. Q. Smith,et al.  1. Bayesian Statistics 4 , 1993 .

[7]  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.

[8]  Michael I. Jordan.  Learning in Graphical Models , 1999, NATO ASI Series.

[9]  S. Petersen,et al.  Characterizing the Hemodynamic Response: Effects of Presentation Rate, Sampling Procedure, and the Possibility of Ordering Brain Activity Based on Relative Timing , 2000, NeuroImage.

[10]  O. Witte,et al.  Functional Mapping of the Human Brain , 2000 .

[11]  R. Turner,et al.  Detecting Latency Differences in Event-Related BOLD Responses: Application to Words versus Nonwords and Initial versus Repeated Face Presentations , 2002, NeuroImage.

[12]  Karl J. Friston,et al.  Assessing the significance of focal activations using their spatial extent , 1994, Human brain mapping.

[13]  Karl J. Friston,et al.  Correcting for non-sphericity in imaging data using classical and Bayesian approaches , 2001, NeuroImage.

[14]  David Barber,et al.  Unified Inference for Variational Bayesian Linear Gaussian State-Space Models , 2006, NIPS.

[15]  S C Strother,et al.  Commentary and Opinion: I. Principal Component Analysis, Variance Partitioning, and “Functional Connectivity” , 1995, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[16]  Matthew J. Beal.  Variational algorithms for approximate Bayesian inference , 2003 .

[17]  Stephen M. Smith,et al.  Probabilistic independent component analysis for functional magnetic resonance imaging , 2004, IEEE Transactions on Medical Imaging.

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

[19]  Zoubin Ghahramani,et al.  The variational Kalman smoother , 2001 .

[20]  J B Poline,et al.  Joint detection-estimation of brain activity in functional MRI: a Multichannel Deconvolution solution , 2005, IEEE Transactions on Signal Processing.

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

[22]  S Makeig,et al.  Analysis of fMRI data by blind separation into independent spatial components , 1998, Human brain mapping.

[23]  Stephen M. Smith,et al.  Improved Optimization for the Robust and Accurate Linear Registration and Motion Correction of Brain Images , 2002, NeuroImage.

[24]  Karl J. Friston,et al.  The choice of basis functions in event-related fMRI , 2001, NeuroImage.

[25]  Michael Brady,et al.  Improved Optimization for the Robust and Accurate Linear Registration and Motion Correction of Brain Images , 2002, NeuroImage.

[26]  Mari Ostendorf,et al.  Cross-Validation EM Training for Robust Parameter Estimation , 2007, 2007 IEEE International Conference on Acoustics, Speech and Signal Processing - ICASSP '07.

[27]  Shinichi Nakajima,et al.  Variational Bayes Solution of Linear Neural Networks and Its Generalization Performance , 2007, Neural Computation.

[28]  Mark D'Esposito,et al.  Variation of BOLD hemodynamic responses across subjects and brain regions and their effects on statistical analyses , 2004, NeuroImage.

[29]  Karl J. Friston.  Imaging neuroscience: principles or maps? , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Gabriele Lohmann,et al.  Bayesian second-level analysis of functional magnetic resonance images , 2003, NeuroImage.

[31]  Karl J. Friston,et al.  Bilinear dynamical systems , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.