The inverse problem of MEG aims at estimating electromagnetic cerebral activity from measurements of the magnetic fields outside the head. After formulating the problem within the Bayesian framework, a hierarchical conditionally Gaussian prior model is introduced, including a physiologically inspired prior model that takes into account the preferred directions of the source currents. The hyperparameter vector consists of prior variances of the dipole moments, assumed to follow a non-conjugate gamma distribution with variable scaling and shape parameters. A point estimate of both dipole moments and their variances can be computed using an iterative alternating sequential updating algorithm, which is shown to be globally convergent. The numerical solution is based on computing an approximation of the dipole moments using a Krylov subspace iterative linear solver equipped with statistically inspired preconditioning and a suitable termination rule. The shape parameters of the model are shown to control the focality, and furthermore, using an empirical Bayes argument, it is shown that the scaling parameters can be naturally adjusted to provide a statistically well justified depth sensitivity scaling. The validity of this interpretation is verified through computed numerical examples. Also, a computed example showing the applicability of the algorithm to analyze realistic time series data is presented.

A hierarchical Krylov–Bayes iterative inverse solver for MEG with physiological preconditioning / Calvetti, D; Pascarella, A; Pitolli, Francesca; Somersalo, E; Vantaggi, Barbara. - In: INVERSE PROBLEMS. - ISSN 0266-5611. - ELETTRONICO. - 31:12(2015). [10.1088/0266-5611/31/12/125005]

A hierarchical Krylov–Bayes iterative inverse solver for MEG with physiological preconditioning

PITOLLI, Francesca;VANTAGGI, Barbara
2015

Abstract

The inverse problem of MEG aims at estimating electromagnetic cerebral activity from measurements of the magnetic fields outside the head. After formulating the problem within the Bayesian framework, a hierarchical conditionally Gaussian prior model is introduced, including a physiologically inspired prior model that takes into account the preferred directions of the source currents. The hyperparameter vector consists of prior variances of the dipole moments, assumed to follow a non-conjugate gamma distribution with variable scaling and shape parameters. A point estimate of both dipole moments and their variances can be computed using an iterative alternating sequential updating algorithm, which is shown to be globally convergent. The numerical solution is based on computing an approximation of the dipole moments using a Krylov subspace iterative linear solver equipped with statistically inspired preconditioning and a suitable termination rule. The shape parameters of the model are shown to control the focality, and furthermore, using an empirical Bayes argument, it is shown that the scaling parameters can be naturally adjusted to provide a statistically well justified depth sensitivity scaling. The validity of this interpretation is verified through computed numerical examples. Also, a computed example showing the applicability of the algorithm to analyze realistic time series data is presented.
2015
Bayesian hierarchical model; brain activity; magnetoencephalography (MEG); prior information; sparsity
01 Pubblicazione su rivista::01a Articolo in rivista
A hierarchical Krylov–Bayes iterative inverse solver for MEG with physiological preconditioning / Calvetti, D; Pascarella, A; Pitolli, Francesca; Somersalo, E; Vantaggi, Barbara. - In: INVERSE PROBLEMS. - ISSN 0266-5611. - ELETTRONICO. - 31:12(2015). [10.1088/0266-5611/31/12/125005]
File allegati a questo prodotto
File Dimensione Formato  
Vantaggi_hierarchical-Krylov–Bayes_2015.pdf

solo gestori archivio

Tipologia: Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 1.72 MB
Formato Adobe PDF
1.72 MB Adobe PDF   Contatta l'autore

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11573/845910
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 28
  • ???jsp.display-item.citation.isi??? 26
social impact