1,230 research outputs found
Neural Connectivity with Hidden Gaussian Graphical State-Model
The noninvasive procedures for neural connectivity are under questioning.
Theoretical models sustain that the electromagnetic field registered at
external sensors is elicited by currents at neural space. Nevertheless, what we
observe at the sensor space is a superposition of projected fields, from the
whole gray-matter. This is the reason for a major pitfall of noninvasive
Electrophysiology methods: distorted reconstruction of neural activity and its
connectivity or leakage. It has been proven that current methods produce
incorrect connectomes. Somewhat related to the incorrect connectivity
modelling, they disregard either Systems Theory and Bayesian Information
Theory. We introduce a new formalism that attains for it, Hidden Gaussian
Graphical State-Model (HIGGS). A neural Gaussian Graphical Model (GGM) hidden
by the observation equation of Magneto-encephalographic (MEEG) signals. HIGGS
is equivalent to a frequency domain Linear State Space Model (LSSM) but with
sparse connectivity prior. The mathematical contribution here is the theory for
high-dimensional and frequency-domain HIGGS solvers. We demonstrate that HIGGS
can attenuate the leakage effect in the most critical case: the distortion EEG
signal due to head volume conduction heterogeneities. Its application in EEG is
illustrated with retrieved connectivity patterns from human Steady State Visual
Evoked Potentials (SSVEP). We provide for the first time confirmatory evidence
for noninvasive procedures of neural connectivity: concurrent EEG and
Electrocorticography (ECoG) recordings on monkey. Open source packages are
freely available online, to reproduce the results presented in this paper and
to analyze external MEEG databases
A two-way regularization method for MEG source reconstruction
The MEG inverse problem refers to the reconstruction of the neural activity
of the brain from magnetoencephalography (MEG) measurements. We propose a
two-way regularization (TWR) method to solve the MEG inverse problem under the
assumptions that only a small number of locations in space are responsible for
the measured signals (focality), and each source time course is smooth in time
(smoothness). The focality and smoothness of the reconstructed signals are
ensured respectively by imposing a sparsity-inducing penalty and a roughness
penalty in the data fitting criterion. A two-stage algorithm is developed for
fast computation, where a raw estimate of the source time course is obtained in
the first stage and then refined in the second stage by the two-way
regularization. The proposed method is shown to be effective on both synthetic
and real-world examples.Comment: Published in at http://dx.doi.org/10.1214/11-AOAS531 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Dynamic filtering of static dipoles in magnetoencephalography
We consider the problem of estimating neural activity from measurements
of the magnetic fields recorded by magnetoencephalography. We exploit
the temporal structure of the problem and model the neural current as a
collection of evolving current dipoles, which appear and disappear, but whose
locations are constant throughout their lifetime. This fully reflects the physiological
interpretation of the model.
In order to conduct inference under this proposed model, it was necessary
to develop an algorithm based around state-of-the-art sequential Monte
Carlo methods employing carefully designed importance distributions. Previous
work employed a bootstrap filter and an artificial dynamic structure
where dipoles performed a random walk in space, yielding nonphysical artefacts
in the reconstructions; such artefacts are not observed when using the
proposed model. The algorithm is validated with simulated data, in which
it provided an average localisation error which is approximately half that of
the bootstrap filter. An application to complex real data derived from a somatosensory
experiment is presented. Assessment of model fit via marginal
likelihood showed a clear preference for the proposed model and the associated
reconstructions show better localisation
- âŠ