State-space solutions to the dynamic magnetoencephalography inverse problem using high performance computing

DOI10.1214/11-AOAS483MaRDI QIDQ641071zbMATHWikidataFDO

Authors Christopher J. Long, Patrick L. Purdon, Simona Temereanca, Neil U. Desai, Matti S. Hämäläinen, Emery N. Brown

Publication date 21 October 2011

Published in The Annals of Applied Statistics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1107.4192

zbMATH Keywords

Kalman filtersource localizationfixed interval smoother

Mathematics Subject Classification ID

Bayesian inference (62F15) Applications of statistics to biology and medical sciences; meta analysis (62P10) Inference from stochastic processes and prediction (62M20) Sums of independent random variables; random walks (60G50) Neural biology (92C20) Biomedical imaging and signal processing (92C55) Numerical algorithms for specific classes of architectures (65Y10)

Abstract: Determining the magnitude and location of neural sources within the brain that are responsible for generating magnetoencephalography (MEG) signals measured on the surface of the head is a challenging problem in functional neuroimaging. The number of potential sources within the brain exceeds by an order of magnitude the number of recording sites. As a consequence, the estimates for the magnitude and location of the neural sources will be ill-conditioned because of the underdetermined nature of the problem. One well-known technique designed to address this imbalance is the minimum norm estimator (MNE). This approach imposes an $L^2$ regularization constraint that serves to stabilize and condition the source parameter estimates. However, these classes of regularizer are static in time and do not consider the temporal constraints inherent to the biophysics of the MEG experiment. In this paper we propose a dynamic state-space model that accounts for both spatial and temporal correlations within and across candidate intracortical sources. In our model, the observation model is derived from the steady-state solution to Maxwell's equations while the latent model representing neural dynamics is given by a random walk process.

Recommendations

• A statistical approach to the inverse problem in magnetoencephalography
• A spatio-temporal solution for the EEG/MEG inverse problem using group penalization methods
• Particle filters for magnetoencephalography
• Bayesian preconditioned CGLS for source separation in MEG time series
• A two-way regularization method for MEG source reconstruction

Cites work

• scientific article; zbMATH DE number 432498 (Why is no real title available?)
• ScaLAPACK Users' Guide
• Smoothness priors analysis of time series

Cited in

(9)

• Estimating the number of sources in magnetoencephalography using spiked population eigenvalues
• Quantifying time-varying sources in magnetoencephalography -- a discrete approach
• Recent trends and advances in solving the inverse problem for EEG source localization
• Hybrid projection methods for large-scale inverse problems with mixed Gaussian priors
• Parallel Minimum p-Norm Solution of the Neuromagnetic Inverse Problem for Realistic Signals Using Exact Hessian-Vector Products
• Bayesian spatiotemporal modeling for inverse problems
• Evaluating stationarity via change-point alternatives with applications to fMRI data
• Dynamic filtering of static dipoles in magnetoencephalography
• A statistical approach to the inverse problem in magnetoencephalography

Describes a project that uses

Uses Software

• ScaLAPACK

This page was built for publication: State-space solutions to the dynamic magnetoencephalography inverse problem using high performance computing

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q641071)