Bayesian inverse analysis of neuromagnetic data using cortically constrained multiple dipoles.
A recently introduced Bayesian model for magnetoencephalographic (MEG) data consistently localized multiple simulated dipoles with the help of marginalization of spatiotemporal background noise covariance structure in the analysis [Jun et al., (2005): Neuroimage 28:84-98]. Here, we elaborated this model to include subject's individual brain surface reconstructions with cortical location and orientation constraints. To enable efficient Markov chain Monte Carlo sampling of the dipole locations, we adopted a parametrization of the source space surfaces with two continuous variables (i.e., spherical angle coordinates). Prior to analysis, we simplified the likelihood by exploiting only a small set of independent measurement combinations obtained by singular value decomposition of the gain matrix, which also makes the sampler significantly faster. We analyzed both realistically simulated and empirical MEG data recorded during simple auditory and visual stimulation. The results show that our model produces reasonable solutions and adequate data fits without much manual interaction. However, the rigid cortical constraints seemed to make the utilized scheme challenging as the sampler did not switch modes of the dipoles efficiently. This is problematic in the presence of evidently highly multimodal posterior distribution, and especially in the relative quantitative comparison of the different modes. To overcome the difficulties with the present model, we propose the use of loose orientation constraints and combined model of prelocalization utilizing the hierarchical minimum-norm estimate and multiple dipole sampling scheme.
Auranen
Toni
T
Laboratory of Computational Engineering, Helsinki University of Technology, Espoo, Finland. Toni.Auranen@hut.fi
Nummenmaa
Aapo
A
Hämäläinen
Matti S
MS
Jääskeläinen
Iiro P
IP
Lampinen
Jouko
J
Vehtari
Aki
A
Sams
Mikko
M
D000161:Acoustic Stimulation; D000465:Algorithms; D001499:Bayes Theorem; D001931:Brain Mapping / Q000379:methods*; D002540:Cerebral Cortex / Q000033:anatomy & histology* / Q000502:physiology*; D003198:Computer Simulation; D005071:Evoked Potentials / Q000502:physiology; D006801:Humans; D007091:Image Processing, Computer-Assisted / Q000379:methods*; D015225:Magnetoencephalography / Q000379:methods*; D008390:Markov Chains; D008959:Models, Neurological; D009010:Monte Carlo Method; D016011:Normal Distribution; D010775:Photic Stimulation; D011336:Probability; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D012815:Signal Processing, Computer-Assisted