Bayesian Estimation of Dynamical Systems: An Application to fMRI

This paper presents a method for estimating the conditional or posterior distribution of the parameters of deterministic dynamical systems. The procedure conforms to an EM implementation of a Gauss-Newton search for the maximum of the conditional or posterior density. The inclusion of priors in the estimation procedure ensures robust and rapid convergence and the resulting conditional densities enable Bayesian inference about the model parameters. The method is demonstrated using an input-state-output model of the hemodynamic coupling between experimentally designed causes or factors in fMRI studies and the ensuing BOLD response. This example represents a generalization of current fMRI analysis models that accommodates nonlinearities and in which the parameters have an explicit physical interpretation. Second, the approach extends classical inference, based on the likelihood of the data given a null hypothesis about the parameters, to more plausible inferences about the parameters of the model given the data. This inference provides for confidence intervals based on the conditional density.

[1]  Geoffrey E. Hinton,et al.  A View of the Em Algorithm that Justifies Incremental, Sparse, and other Variants , 1998, Learning in Graphical Models.

[2]  Karl J. Friston,et al.  Classical and Bayesian Inference in fMRI , 2004 .

[3]  M. Moskowitz,et al.  Importance of Nitric Oxide Synthase Inhibition to the Attenuated Vascular Responses Induced by Topical L-Nitroarginine during Vibrissal Stimulation , 1994, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[4]  M. Raichle,et al.  The Effects of Changes in PaCO2 Cerebral Blood Volume, Blood Flow, and Vascular Mean Transit Time , 1974, Stroke.

[5]  R. Buxton,et al.  Dynamics of blood flow and oxygenation changes during brain activation: The balloon model , 1998, Magnetic resonance in medicine.

[7]  B. Rosen,et al.  Evidence of a Cerebrovascular Postarteriole Windkessel with Delayed Compliance , 1999, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[8]  Karl J. Friston,et al.  Nonlinear event‐related responses in fMRI , 1998, Magnetic resonance in medicine.

[9]  L. Fahrmeir,et al.  Multivariate statistical modelling based on generalized linear models , 1994 .

[10]  Karl J. Friston,et al.  Characterizing the Relationship between BOLD Contrast and Regional Cerebral Blood Flow Measurements by Varying the Stimulus Presentation Rate , 1997, NeuroImage.

[11]  Winson Taam,et al.  Nonlinear System Analysis and Identification From Random Data , 1991 .

[12]  Françoise Lamnabhi-Lagarrigue,et al.  An algebraic approach to nonlinear functional expansions , 1983 .

[13]  G. Crelier,et al.  Linear coupling between cerebral blood flow and oxygen consumption in activated human cortex. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[14]  Tom Leonard Bayesian methods for binomial data , 1972 .

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

[16]  Ying Zheng,et al.  A Model of the Hemodynamic Response and Oxygen Delivery to Brain , 2002, NeuroImage.

[17]  Karl J. Friston Regulation of rCBF by diffusible signals: An analysis of constraints on diffusion and elimination , 1995 .

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

[19]  C. Büchel,et al.  Modulation of connectivity in visual pathways by attention: cortical interactions evaluated with structural equation modelling and fMRI. , 1997, Cerebral cortex.

[20]  Karl J. Friston,et al.  Nonlinear Responses in fMRI: The Balloon Model, Volterra Kernels, and Other Hemodynamics , 2000, NeuroImage.

[21]  Dewen Hu,et al.  An Evaluation of Linear Model Analysis Techniques for Processing Images of Microcirculation Activity , 1998, NeuroImage.

[22]  D. Harville Maximum Likelihood Approaches to Variance Component Estimation and to Related Problems , 1977 .

[23]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[24]  D. G. Simpson,et al.  The Statistical Analysis of Discrete Data , 1989 .

[25]  Rick Chappell,et al.  Multivariate Statistical Modelling Based on Generalized Linear Models , 2003 .

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