A data-driven model of the generation of human EEG based on a spatially distributed stochastic wave equation

We discuss a model for the dynamics of the primary current density vector field within the grey matter of human brain. The model is based on a linear damped wave equation, driven by a stochastic term. By employing a realistically shaped average brain model and an estimate of the matrix which maps the primary currents distributed over grey matter to the electric potentials at the surface of the head, the model can be put into relation with recordings of the electroencephalogram (EEG). Through this step it becomes possible to employ EEG recordings for the purpose of estimating the primary current density vector field, i.e. finding a solution of the inverse problem of EEG generation. As a technique for inferring the unobserved high-dimensional primary current density field from EEG data of much lower dimension, a linear state space modelling approach is suggested, based on a generalisation of Kalman filtering, in combination with maximum-likelihood parameter estimation. The resulting algorithm for estimating dynamical solutions of the EEG inverse problem is applied to the task of localising the source of an epileptic spike from a clinical EEG data set; for comparison, we apply to the same task also a non-dynamical standard algorithm.

[1]  赤池 弘次,et al.  Statistical analysis and control of dynamic systems , 1988 .

[2]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[3]  Arnold Neumaier,et al.  Solving Ill-Conditioned and Singular Linear Systems: A Tutorial on Regularization , 1998, SIAM Rev..

[4]  Tohru Ozaki,et al.  Recursive penalized least squares solution for dynamical inverse problems of EEG generation , 2004, Human brain mapping.

[5]  H. Akaike.  A new look at the statistical model identification , 1974 .

[6]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[7]  R Grave de Peralta Menendez,et al.  Imaging the electrical activity of the brain: ELECTRA , 2000, Human brain mapping.

[8]  R. Engle.  Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation , 1982 .

[9]  Tohru Ozaki,et al.  GARCH modelling of covariance in dynamical estimation of inverse solutions , 2004 .

[10]  Jens Timmer,et al.  Handbook of Time Series Analysis , 2006 .

[11]  Jorge J. Riera,et al.  EEG-distributed inverse solutions for a spherical head model , 1998 .

[12]  G. Tunnicliffe Wilson,et al.  Statistical analysis and control of dynamic systems , 1988 .

[13]  Richard M. Leahy,et al.  Electromagnetic brain mapping , 2001, IEEE Signal Process. Mag..

[14]  Michael A. Arbib,et al.  Topics in Mathematical System Theory , 1969 .

[15]  J. P. Ary,et al.  Location of Sources of Evoked Scalp Potentials: Corrections for Skull and Scalp Thicknesses , 1981, IEEE Transactions on Biomedical Engineering.

[16]  H. Akaike.  Markovian representation of stochastic processes and its application to the analysis of autoregressive moving average processes , 1974 .

[17]  F. H. Lopes da Silva,et al.  Model of brain rhythmic activity , 1974, Kybernetik.

[18]  P A Robinson,et al.  Estimation of multiscale neurophysiologic parameters by electroencephalographic means , 2004, Human brain mapping.

[19]  Melvin J. Hinich,et al.  Time Series Analysis by State Space Methods , 2001 .

[20]  R D Pascual-Marqui,et al.  Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. , 2002, Methods and findings in experimental and clinical pharmacology.

[21]  G. Kitagawa.  Smoothness priors analysis of time series , 1996 .

[22]  Hirotugu Akaike,et al.  Likelihood and the Bayes procedure , 1980 .

[23]  D. Lehmann,et al.  Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. , 1994, International journal of psychophysiology : official journal of the International Organization of Psychophysiology.

[24]  P. Hansen,et al.  Methods and Applications of Inversion , 2000 .

[25]  Nelson J. Trujillo-Barreto,et al.  Realistically Coupled Neural Mass Models Can Generate EEG Rhythms , 2007, Neural Computation.

[26]  A. Jazwinski.  Stochastic Processes and Filtering Theory , 1970 .

[27]  Tohru Ozaki,et al.  A solution to the dynamical inverse problem of EEG generation using spatiotemporal Kalman filtering , 2004, NeuroImage.

[28]  Walter J. Freeman,et al.  Neurodynamics: An Exploration in Mesoscopic Brain Dynamics , 2000, Perspectives in Neural Computing.

[29]  Manfred Deistler,et al.  Linear Models for Multivariate Time Series , 2006 .

[30]  Arthur W. Toga,et al.  A Probabilistic Atlas of the Human Brain: Theory and Rationale for Its Development The International Consortium for Brain Mapping (ICBM) , 1995, NeuroImage.

[31]  T. Bollerslev,et al.  Generalized autoregressive conditional heteroskedasticity , 1986 .

[32]  Tohru Ozaki,et al.  Modelling non-stationary variance in EEG time series by state space GARCH model , 2006, Comput. Biol. Medicine.