Recursive penalized least squares solution for dynamical inverse problems of EEG generation

In the dynamical inverse problem of electroencephalogram (EEG) generation where a specific dynamics for the electrical current distribution is assumed, we can impose general spatiotemporal constraints onto the solution by casting the problem into a state space representation and assuming a specific class of parametric models for the dynamics. The Akaike Bayesian Information Criterion (ABIC), which is based on the Type II likelihood, was used to estimate the parameters and evaluate the model. In addition, dynamic low‐resolution brain electromagnetic tomography (LORETA), a new approach for estimating the current distribution is introduced. A recursive penalized least squares (RPLS) step forms the main element of our implementation. To obtain improved inverse solutions, dynamic LORETA exploits both spatial and temporal information, whereas LORETA uses only spatial information. A considerable improvement in performance compared to LORETA was found when dynamic LORETA was applied to simulated EEG data, and the new method was applied also to clinical EEG data. Hum. Brain Mapp. 21:221–235, 2004. © 2004 Wiley‐Liss, Inc.

[1]  A K Louis,et al.  Efficient algorithms for the regularization of dynamic inverse problems: I. Theory , 2002 .

[2]  R. Pascual-Marqui Review of methods for solving the EEG inverse problem , 1999 .

[3]  Genshiro Kitagawa,et al.  State Space Modeling of Time Series , 1994 .

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

[5]  Rolando Grave de Peralta Menendez,et al.  Backus and gilbert method for vector fields , 1999 .

[6]  G. Backus,et al.  The Resolving Power of Gross Earth Data , 1968 .

[7]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

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

[9]  N. L. Johnson,et al.  Multivariate Analysis , 1958, Nature.

[10]  Rob S. MacLeod,et al.  Inverse electrocardiography by simultaneous imposition of multiple constraints , 1999, IEEE Transactions on Biomedical Engineering.

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

[12]  Rolando Grave de Peralta,et al.  Comparison of Algorithms for the Localization of Focal Sources: Evaluation with simulated data and analysis of experimental data. , 2002 .

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

[14]  Karl J. Friston,et al.  Systematic Regularization of Linear Inverse Solutions of the EEG Source Localization Problem , 2002, NeuroImage.

[15]  K. Matsuura,et al.  Selective minimum-norm solution of the biomagnetic inverse problem , 1995, IEEE Transactions on Biomedical Engineering.

[16]  C. Michel,et al.  Linear inverse solutions with optimal resolution kernels applied to electromagnetic tomography , 1997, Human brain mapping.

[17]  Grace Wahba,et al.  Spline Models for Observational Data , 1990 .

[18]  J. Riera,et al.  Electric lead field for a piecewise homogeneous volume conductor model of the head , 1998, IEEE Transactions on Biomedical Engineering.

[19]  E. Rodin,et al.  Dipole sources of the human alpha rhythm , 2005, Brain Topography.

[20]  Erkki Somersalo,et al.  State Estimation in Time‐Varying Electrical Impedance Tomography , 1999 .

[21]  P. Nunez,et al.  Electric fields of the brain , 1981 .

[22]  Jari P. Kaipio,et al.  Dynamic reconstruction in SPET , 1997, Proceedings of the 19th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. 'Magnificent Milestones and Emerging Opportunities in Medical Engineering' (Cat. No.97CH36136).

[23]  H. Akaike SEASONAL ADJUSTMENT BY A BAYESIAN MODELING , 1980 .

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

[25]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[26]  Manfred Fuchs,et al.  Numerical aspects of spatio-temporal current density reconstruction from EEG-/MEG-data , 2001, IEEE Transactions on Medical Imaging.

[27]  David S. Stoffer Time series analysis by state space models: J. Durbin and S. J. Koopman; Oxford University Press, Oxford, 2001, pp 253 + xvii, ISBN: 0 19 852354 8 , 2003, Autom..

[28]  Carsten H. Wolters,et al.  Efficient algorithms for the regularization of dynamic inverse problems: II. Applications , 2002 .

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

[30]  Rolando J. Biscay,et al.  Frequency domain models of the EEG , 2005, Brain Topography.

[31]  Jari P. Kaipio,et al.  Fixed-lag smoothing and state estimation in dynamic electrical impedance tomography , 2001 .

[32]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[33]  Morris Grenfell Davies Frequency‐Domain Models for Room Response , 2006 .