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. (C) 2004 Wiley-Liss, Inc.