Model Estimation of Cerebral Hemodynamics Between Blood Flow and Volume Changes: A Data-Based Modeling Approach

It is well known that there is a dynamic relationship between cerebral blood flow (CBF) and cerebral blood volume (CBV). With increasing applications of functional MRI, where the blood oxygen-level-dependent signals are recorded, the understanding and accurate modeling of the hemodynamic relationship between CBF and CBV becomes increasingly important. This study presents an empirical and data-based modeling framework for model identification from CBF and CBV experimental data. It is shown that the relationship between the changes in CBF and CBV can be described using a parsimonious autoregressive with exogenous input model structure. It is observed that neither the ordinary least-squares (LS) method nor the classical total least-squares (TLS) method can produce accurate estimates from the original noisy CBF and CBV data. A regularized total least-squares (RTLS) method is thus introduced and extended to solve such an error-in-the-variables problem. Quantitative results show that the RTLS method works very well on the noisy CBF and CBV data. Finally, a combination of RTLS with a filtering method can lead to a parsimonious but very effective model that can characterize the relationship between the changes in CBF and CBV.

[1]  Steve A. Billings,et al.  Term and variable selection for non-linear system identification , 2004 .

[2]  Sabine Van Huffel,et al.  Total least squares and errors-in-variables modeling , 2007, Signal Process..

[3]  Stephen A. Billings,et al.  A three-compartment model of the hemodynamic response and oxygen delivery to brain , 2005, NeuroImage.

[4]  Sheng Chen,et al.  Orthogonal least squares methods and their application to non-linear system identification , 1989 .

[5]  Gene H. Golub,et al.  Tikhonov Regularization and Total Least Squares , 1999, SIAM J. Matrix Anal. Appl..

[6]  S. Billings,et al.  Piecewise linear identification of non-linear systems , 1987 .

[7]  Karl J. Friston,et al.  Analysis of fMRI Time-Series Revisited—Again , 1995, NeuroImage.

[8]  Naoki Saito,et al.  Simultaneous noise suppression and signal compression using a library of orthonormal bases and the minimum-description-length criterion , 1994, Defense, Security, and Sensing.

[9]  Karl J. Friston,et al.  Analysis of fMRI Time-Series Revisited , 1995, NeuroImage.

[10]  Sabine Van Huffel,et al.  Total least squares problem - computational aspects and analysis , 1991, Frontiers in applied mathematics.

[11]  S. Billings,et al.  Correlation based model validity tests for non-linear models , 1986 .

[12]  G. Golub,et al.  Regularized Total Least Squares Based on Quadratic Eigenvalue Problem Solvers , 2004 .

[13]  Gene H. Golub,et al.  An analysis of the total least squares problem , 1980, Milestones in Matrix Computation.

[14]  Stephen A. Billings,et al.  Improved model identification for non-linear systems using a random subsampling and multifold modelling (RSMM) approach , 2009, Int. J. Control.

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

[16]  L. A. Aguirre,et al.  Dynamical effects of overparametrization in nonlinear models , 1995 .

[17]  Stephen A. Billings,et al.  Nonlinear model validation using correlation tests , 1994 .

[18]  Sheng Chen,et al.  Representations of non-linear systems: the NARMAX model , 1989 .

[19]  M. Peruggia.  Total Least Squares and Errors-in-Variables Modeling: Analysis, Algorithms and Applications , 2003 .

[20]  S. Billings,et al.  Interpretation of non-linear frequency response functions , 1990 .

[21]  S. A. Billings,et al.  Structure Detection and Model Validity Tests in the Identification of Nonlinear Systems , 1983 .

[22]  Stephen A. Billings,et al.  Model structure selection using an integrated forward orthogonal search algorithm assisted by squared correlation and mutual information , 2008, Int. J. Model. Identif. Control..

[23]  S. Billings,et al.  Recursive algorithm for computing the frequency response of a class of non-linear difference equation models , 1989 .

[24]  Stephen A. Billings,et al.  Model validation tests for multivariable nonlinear models including neural networks , 1995 .

[25]  V. Torczon,et al.  Direct search methods: then and now , 2000 .

[26]  Stephen A. Billings,et al.  Recurrent radial basis function networks for adaptive noise cancellation , 1995, Neural Networks.

[27]  J. Mayhew,et al.  Concurrent Optical Imaging Spectroscopy and Laser-Doppler Flowmetry: The Relationship between Blood Flow, Oxygenation, and Volume in Rodent Barrel Cortex , 2001, NeuroImage.

[28]  S. Billings,et al.  Output Frequency Characteristics of Nonlinear Systems , 1996 .

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

[30]  Nikolas P. Galatsanos,et al.  Regularized constrained total least-squares image restoration , 1994, Other Conferences.

[31]  Sabine Van Huffel,et al.  Recent advances in total least squares techniques and errors-in-variables modeling , 1997 .

[32]  Karl J. Friston,et al.  Characterizing the Response of PET and fMRI Data Using Multivariate Linear Models , 1997, NeuroImage.

[33]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[34]  S. A. Billings,et al.  Experimental design and identifiability for non-linear systems , 1987 .

[35]  J. Mayhew,et al.  Changes in Blood Flow, Oxygenation, and Volume Following Extended Stimulation of Rodent Barrel Cortex , 2002, NeuroImage.

[36]  Anestis Antoniadis,et al.  Model selection using wavelet decomposition and applications , 1997 .

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

[38]  I. J. Leontaritis,et al.  Model selection and validation methods for non-linear systems , 1987 .

[39]  Ling Xia,et al.  Truncated Total Least Squares: A New Regularization Method for the Solution of ECG Inverse Problems , 2008, IEEE Transactions on Biomedical Engineering.

[40]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[41]  Stephen A. Billings,et al.  Generalized multiscale radial basis function networks , 2007, Neural Networks.

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

[43]  S. Huffel,et al.  Total Least Squares and Errors-in-Variables Modeling : Analysis, Algorithms and Applications , 2002 .

[44]  Ying Zheng,et al.  The Hemodynamic Impulse Response to a Single Neural Event , 2003, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[45]  Gene H. Golub,et al.  Some modified matrix eigenvalue problems , 1973, Milestones in Matrix Computation.

[46]  Stephen A. Billings,et al.  An adaptive orthogonal search algorithm for model subset selection and non-linear system identification , 2008, Int. J. Control.

[47]  Tohru Ozaki,et al.  fMRI activation maps based on the NN-ARx model , 2004, NeuroImage.

[48]  J. Mayhew,et al.  A Model of the Dynamic Relationship between Blood Flow and Volume Changes during Brain Activation , 2004, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[49]  David M. Simpson,et al.  Multivariate dynamic analysis of cerebral blood flow regulation in humans , 2000, IEEE Transactions on Biomedical Engineering.

[50]  Lennart Ljung,et al.  System Identification: Theory for the User , 1987 .

[51]  Szi-Wen Chen,et al.  A two-stage discrimination of cardiac arrhythmias using a total least squares-based Prony modeling algorithm , 2000, IEEE Trans. Biomed. Eng..

[52]  Ying Zheng,et al.  Spectroscopic Analysis of Changes in Remitted Illumination: The Response to Increased Neural Activity in Brain , 1999, NeuroImage.

[53]  Stephen M. Smith,et al.  Temporal Autocorrelation in Univariate Linear Modeling of FMRI Data , 2001, NeuroImage.

[54]  Stephen A. Billings,et al.  Constructing an overall dynamical model for a system with changing design parameter properties , 2008, Int. J. Model. Identif. Control..

[55]  Nikolas P. Galatsanos,et al.  Regularized constrained total least-squares image restoration , 1994, Other Conferences.

[56]  Stergios B. Fotopoulos,et al.  Discrete-Time Dynamic Models , 2001, Technometrics.

[57]  I. J. Leontaritis,et al.  Input-output parametric models for non-linear systems Part II: stochastic non-linear systems , 1985 .

[58]  K. Saito,et al.  Tooth shape reconstruction from ct images using spline Curves , 2007, 2007 International Conference on Wavelet Analysis and Pattern Recognition.

[59]  Stephen A. Billings,et al.  Sparse Model Identification Using a Forward Orthogonal Regression Algorithm Aided by Mutual Information , 2007, IEEE Transactions on Neural Networks.

[60]  Karl J. Friston,et al.  Bayesian Estimation of Dynamical Systems: An Application to fMRI , 2002, NeuroImage.

[61]  Georgios D. Mitsis,et al.  Nonlinear modeling of the dynamic effects of arterial pressure and CO2 variations on cerebral blood flow in healthy humans , 2004, IEEE Trans. Biomed. Eng..

[62]  Jeffrey C. Lagarias,et al.  Convergence Properties of the Nelder-Mead Simplex Method in Low Dimensions , 1998, SIAM J. Optim..

[63]  Sabine Van Huffel,et al.  Total Least Squares and Errors-in-variables Modeling , 2007, Comput. Stat. Data Anal..

[64]  S. Billings,et al.  Identification of Nonlinear Parameter-Dependent Common-Structured models to accommodate varying experimental conditions and design parameter properties , 2006 .

[65]  Sabine Van Huffel,et al.  Overview of total least-squares methods , 2007, Signal Process..

[66]  Patrizia Baraldi,et al.  An ARX model-based approach to trial by trial identification of fMRI-BOLD responses , 2007, NeuroImage.