Bayesian modeling and analysis of the magnetoencephalography and electroencephalography modalities provide a flexible framework for introducing prior information complementary to the measured data. This prior information is often qualitative in nature, making the translation of the available information into a computational model a challenging task. We propose a generalized gamma family of hyperpriors which allows the impressed currents to be focal and we advocate a fast and efficient iterative algorithm, the iterative alternating sequential algorithm for computing maximum a posteriori (MAP) estimates. Furthermore, we show that for particular choices of the scalar parameters specifying the hyperprior, the algorithm effectively approximates popular regularization strategies such as the minimum current estimate and the minimum support estimate. The connection between priorconditioning and adaptive regularization methods is also pointed out. The posterior densities are explored by means of a Markov chain Monte Carlo strategy suitable for this family of hypermodels. The computed experiments suggest that the known preference of regularization methods for superficial sources over deep sources is a property of the MAP estimators only, and that estimation of the posterior mean in the hierarchical model is better adapted for localizing deep sources.

No abstract

Electroencephalography (EEG) is a non-invasive imaging modality in which a primary current density generated by the neural activity in the brain is to be reconstructed based on external electric potential measurements. This paper focuses on the finite element method (FEM) from both forward and inverse aspects. The goal is to establish a clear correspondence between the lowest order Raviart-Thomas basis functions and dipole sources as well as to show that the adopted FEM approach is computationally effective. Each basis function is associated with a dipole moment and a location. Four candidate locations are tested. Numerical experiments cover two different spherical multilayer head models, four mesh resolutions and two different forward simulation approaches, one based on FEM and one based on the boundary element method (BEM) with standard dipoles as sources. The forward simulation accuracy is examined through column-and matrix-wise relative errors as well as through performance in inverse dipole localization. A closed-form approximation of dipole potential was used as the reference forward simulation. The results suggest that the present approach is comparable or superior to BEM and to the recent FEM based subtraction approach regarding both accuracy, computation time and accessibility of implementation.