Quantification of the benefit from integrating MEG and EEG data in minimum l2-norm estimation

Source current estimation from electromagnetic (MEG and EEG) signals is an ill-posed problem that often produces blurry or inaccurately positioned estimates. The two modalities have distinct factors limiting the resolution, e.g., MEG cannot detect radially oriented sources, while EEG is sensitive to accuracy of the head model. This makes combined EEG+MEG estimation techniques desirable, but different acquisition noise statistics, complexity of the head models, and lack of pertinent metrics all complicate the assessment of the resulting improvements. We investigated analytically the effect of including EEG recordings in MEG studies versus the addition of new MEG channels when computing noise-normalized minimum l(2)-norm estimates. Three-compartment boundary-element forward models were constructed using structural MRI scans for four subjects. Singular value analysis of the resulting forward models predicted better performance of the EEG+MEG case in the form of higher matrix rank. MNE inverse operators for EEG, MEG and EEG+MEG were constructed using the sensor noise covariance estimated from data. Metrics derived from the resolution matrices predicted higher spatial resolution in EEG+MEG as compared to MEG due to decreased spread (lower spatial dispersion, higher resolution index) with no reduction in dipole localization error. The effect was apparent in all source locations, with increased magnitude for deep areas such as the cingulate cortex. We were also able to corroborate the results for the somatosensory cortex using median nerve responses.