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

DOI10.1214/11-AOAS483zbMATH Open1223.62160arXiv1107.4192WikidataQ35540778 ScholiaQ35540778MaRDI QIDQ641071FDOQ641071

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)

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.

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

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

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)

Cites Work

• ScaLAPACK Users' Guide
• Title not available (Why is that?)
• Smoothness priors analysis of time series

Cited In (6)

• Parallel Minimum p-Norm Solution of the Neuromagnetic Inverse Problem for Realistic Signals Using Exact Hessian-Vector Products
• Evaluating stationarity via change-point alternatives with applications to fMRI data
• Bayesian spatiotemporal modeling for inverse problems
• Dynamic filtering of static dipoles in magnetoencephalography
• Quantifying time-varying sources in magnetoencephalography -- a discrete approach
• Hybrid projection methods for large-scale inverse problems with mixed Gaussian priors

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)