Multimodal integration of EEG and MEG data: a simulation study with variable signal-to-noise ratio and number of sensors

Exploring the impact of manual and automatic EEG pre-processing methods on interpersonal neural synchrony measures in parent-infant hyperscanning studies

BACKGROUND

Electroencephalograph (EEG) hyperscanning allows studying Interpersonal Neural Synchrony (INS) between two or more individuals across different social conditions, including parent-infant interactions. Signal pre-processing is crucial to optimize computation of INS estimates; however, few attempts have been made at comparing the impact of different dyadic EEG data pre-processing methods on INS estimates.

NEW METHODS

EEG data collected on 31 mother-infant dyads (8-10 months) engaged in a Face-to-Face Still-Face Procedure were pre-processed with two versions of the same pipeline, the "automated" and the "manual". Cross-frequency PLV in the theta (3-5 Hz, 4-7 Hz) and alpha (6-9 Hz, 8-12 Hz) frequency bands were computed after automated and manual pre-processing and compared through Pearson's correlations and Repeated Measures ANOVAs.

RESULTS

PLVs computed in the theta, but not alpha, frequency band were significantly higher after automated pre-processing than after manual pre-processing. Moreover, the automated pipeline rejected a significantly lower percentage of ICs and epochs compared to the manual pipeline.

COMPARISON WITH EXISTING METHODS

While no direct comparison with existing dyadic EEG data pre-processing pipelines was made, this is the first study assessing the impact of different methodological decisions, particularly of the degree of pre-processing automatization, on cross-frequency PLV computed on a dataset of parent-infant dyads.

CONCLUSIONS

Non-directional phase-based INS indexes such as the PLV seem to be affected by the degree of automatization of the pre-processing pipeline. Future research should strive for standardization of dyadic EEG pre-processing methods.

Volumetric mesoscopic electrophysiology: a new imaging modality for the nonhuman primate

The primate brain is a densely interconnected organ whose function is best understood by recording from the entire structure in parallel, rather than parts of it in sequence. However, available methods either have limited temporal resolution (functional magnetic resonance imaging; fMRI), limited spatial resolution (macroscopic electroencephalography), or a limited field of view (microscopic electrophysiology). To address this need, we developed a volumetric, mesoscopic recording approach (MePhys) by tessellating the volume of a monkey hemisphere with 992 electrode contacts that were distributed across 62 chronically implanted multielectrode shafts. We showcase the scientific promise of MePhys by describing the functional interactions of local field potentials between the more than 300,000 simultaneously recorded pairs of electrodes. We find that a subanesthetic dose of ketamine-believed to mimic certain aspects of psychosis-can create a pronounced state of functional disconnection and prevent the formation of stable large-scale intrinsic states. We conclude that MePhys provides a new and fundamentally distinct window into brain function whose unique profile of strengths and weaknesses complements existing approaches in synergistic ways.NEW & NOTEWORTHY We created a new imaging modality for the nonhuman primate, mesoscopic electrophysiology, or MePhys by sampling local field potentials (LFPs) in a dense three-dimensional grid from across the volume of one entire hemisphere. MePhys combines the millisecond temporal resolution of electrophysiology with the large field of view and millimeter spatial resolution of functional magnetic resonance imaging (fMRI). MePhys' unique profile of strengths and limitations makes it an ideal imaging method for the nonhuman primate brain observatories of the future.

Volumetric mesoscopic electrophysiology: a new imaging modality for the non-human primate

The primate brain is a densely interconnected organ whose function is best understood by recording from the entire structure in parallel, rather than parts of it in sequence. However, available methods either have limited temporal resolution (functional magnetic resonance imaging), limited spatial resolution (macroscopic electroencephalography), or a limited field of view (microscopic electrophysiology). To address this need, we developed a volumetric, mesoscopic recording approach ( MePhys ) by tessellating the volume of a monkey hemisphere with 992 electrode contacts that were distributed across 62 chronically implanted multi-electrode shafts. We showcase the scientific promise of MePhys by describing the functional interactions of local field potentials between the more than 300,000 simultaneously recorded pairs of electrodes. We find that a subanesthetic dose of ketamine -believed to mimic certain aspects of psychosis- can create a pronounced state of functional disconnection and prevent the formation of stable large-scale intrinsic states. We conclude that MePhys provides a new and fundamentally distinct window into brain function whose unique profile of strengths and weaknesses complements existing approaches in synergistic ways.

Spatial fidelity of MEG/EEG source estimates: A general evaluation approach

Low spatial resolution is often cited as the most critical limitation of magneto- and electroencephalography (MEG and EEG), but a unifying framework for quantifying the spatial fidelity of M/EEG source estimates has yet to be established; previous studies have focused on linear estimation methods under ideal scenarios without noise. Here we present an approach that quantifies the spatial fidelity of M/EEG estimates from simulated patch activations over the entire neocortex superposed on measured resting-state data. This approach grants more generalizability in the evaluation process that allows for, e.g., comparing linear and non-linear estimates in the whole brain for different signal-to-noise ratios (SNR), number of active sources and activation waveforms. Using this framework, we evaluated the MNE, dSPM, sLORETA, eLORETA, and MxNE methods and found that the spatial fidelity varies significantly with SNR, following a largely sigmoidal curve whose shape varies depending on which aspect of spatial fidelity that is being quantified and the source estimation method. We believe that these methods and results will be useful when interpreting M/EEG source estimates as well as in methods development.

Mechanisms of face specificity - Differentiating speed and accuracy in face cognition by event-related potentials of central processing

Given the crucial role of face recognition in social life, it is hardly surprising that cognitive processes specific for faces have been identified. In previous individual differences studies, the speed (measured in easy tasks) and accuracy (difficult tasks) of face cognition (FC, involving perception and recognition of faces) have been shown to form distinct abilities, going along with divergent factorial structures. This result has been replicated, but remained unexplained. To fill this gap, we first parameterized the sub-processes underlying speed vs. accuracy in easy and difficult memory tasks for faces and houses in a large sample. Then, we analyzed event-related potentials (ERPs) extracted from the EEG by using residue iteration decomposition (RIDE), yielding a central (C) component that is comparable to a purified P300. Structural equation modeling (SEM) was applied to estimate face specificity of C component latencies and amplitudes. If performance in easy tasks relies on purely general processes that are insensitive to stimulus content, there should be no specificity of individual differences in the latency recorded in easy tasks. However, in difficult tasks specificity was expected. Results indicated that, contrary to our predictions, specificity occurred in the C component latency of both speed-based and accuracy-based measures, but was stronger in accuracy. Further analyses suggested specific relationships between the face-related C latency and FC ability. Finally, we detected specificity in RTs of easy tasks when single tasks were modeled, but not when multiple tasks were jointly modeled. This suggests that the mechanisms leading to face specificity in performance speed are distinct across tasks.

Empirical Bayes evaluation of fused EEG-MEG source reconstruction: Application to auditory mismatch evoked responses

We here turn the general and theoretical question of the complementarity of EEG and MEG for source reconstruction, into a practical empirical one. Precisely, we address the challenge of evaluating multimodal data fusion on real data. For this purpose, we build on the flexibility of Parametric Empirical Bayes, namely for EEG-MEG data fusion, group level inference and formal hypothesis testing. The proposed approach follows a two-step procedure by first using unimodal or multimodal inference to derive a cortical solution at the group level; and second by using this solution as a prior model for single subject level inference based on either unimodal or multimodal data. Interestingly, for inference based on the same data (EEG, MEG or both), one can then formally compare, as alternative hypotheses, the relative plausibility of the two unimodal and the multimodal group priors. Using auditory data, we show that this approach enables to draw important conclusions, namely on (i) the superiority of multimodal inference, (ii) the greater spatial sensitivity of MEG compared to EEG, (iii) the ability of EEG data alone to source reconstruct temporal lobe activity, (iv) the usefulness of EEG to improve MEG based source reconstruction. Importantly, we largely reproduce those findings over two different experimental conditions. We here focused on Mismatch Negativity (MMN) responses for which generators have been extensively investigated with little homogeneity in the reported results. Our multimodal inference at the group level revealed spatio-temporal activity within the supratemporal plane with a precision which, to our knowledge, has never been achieved before with non-invasive recordings.

Wearable neuroimaging: Combining and contrasting magnetoencephalography and electroencephalography

One of the most severe limitations of functional neuroimaging techniques, such as magnetoencephalography (MEG), is that participants must maintain a fixed head position during data acquisition. This imposes restrictions on the characteristics of the experimental cohorts that can be scanned and the experimental questions that can be addressed. For these reasons, the use of 'wearable' neuroimaging, in which participants can move freely during scanning, is attractive. The most successful example of wearable neuroimaging is electroencephalography (EEG), which employs lightweight and flexible instrumentation that makes it useable in almost any experimental setting. However, EEG has major technical limitations compared to MEG, and therefore the development of wearable MEG, or hybrid MEG/EEG systems, is a compelling prospect. In this paper, we combine and compare EEG and MEG measurements, the latter made using a new generation of optically-pumped magnetometers (OPMs). We show that these new second generation commercial OPMs, can be mounted on the scalp in an 'EEG-like' cap, enabling the acquisition of high fidelity electrophysiological measurements. We show that these sensors can be used in conjunction with conventional EEG electrodes, offering the potential for the development of hybrid MEG/EEG systems. We compare concurrently measured signals, showing that, whilst both modalities offer high quality data in stationary subjects, OPM-MEG measurements are less sensitive to artefacts produced when subjects move. Finally, we show using simulations that OPM-MEG offers a fundamentally better spatial specificity than EEG. The demonstrated technology holds the potential to revolutionise the utility of functional brain imaging, exploiting the flexibility of wearable systems to facilitate hitherto impractical experimental paradigms.

A Time-Varying Connectivity Analysis from Distributed EEG Sources: A Simulation Study

Time-varying connectivity analysis based on sources reconstructed using inverse modeling of electroencephalographic (EEG) data is important to understand the dynamic behaviour of the brain. We simulated cortical data from a visual spatial attention network with a time-varying connectivity structure, and then simulated the propagation to the scalp to obtain EEG data. Distributed EEG source modeling using sLORETA was applied. We compared different dipole (representing a source) selection strategies based on their time series in a region of interest. Next, we estimated multivariate autoregressive (MVAR) parameters using classical Kalman filter and general linear Kalman filter approaches followed by the calculation of partial directed coherence (PDC). MVAR parameters and PDC values for the selected sources were compared with the ground-truth. We found that the best strategy to extract the time series of a region of interest was to select a dipole with time series showing the highest correlation with the average time series in the region of interest. Dipole selection based on power or based on the largest singular value offer comparable alternatives. Among the different Kalman filter approaches, the use of a general linear Kalman filter was preferred to estimate PDC based connectivity except when only a small number of trials are available. In the latter case, the classical Kalman filter can be an alternative.

Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy

Fusion of electroencephalography (EEG) and magnetoencephalography (MEG) data using maximum entropy on the mean method (MEM-fusion) takes advantage of the complementarities between EEG and MEG to improve localization accuracy. Simulation studies demonstrated MEM-fusion to be robust especially in noisy conditions such as single spike source localizations (SSSL). Our objective was to assess the reliability of SSSL using MEM-fusion on clinical data. We proposed to cluster SSSL results to find the most reliable and consistent source map from the reconstructed sources, the so-called consensus map. Thirty-four types of interictal epileptic discharges (IEDs) were analyzed from 26 patients with well-defined epileptogenic focus. SSSLs were performed on EEG, MEG, and fusion data and consensus maps were estimated using hierarchical clustering. Qualitative (spike-to-spike reproducibility rate, SSR) and quantitative (localization error and spatial dispersion) assessments were performed using the epileptogenic focus as clinical reference. Fusion SSSL provided significantly better results than EEG or MEG alone. Fusion found at least one cluster concordant with the clinical reference in all cases. This concordant cluster was always the one involving the highest number of spikes. Fusion yielded highest reproducibility (SSR EEG = 55%, MEG = 71%, fusion = 90%) and lowest localization error. Also, using only few channels from either modality (21EEG + 272MEG or 54EEG + 25MEG) was sufficient to reach accurate fusion. MEM-fusion with consensus map approach provides an objective way of finding the most reliable and concordant generators of IEDs. We, therefore, suggest the pertinence of SSSL using MEM-fusion as a valuable clinical tool for presurgical evaluation of epilepsy.

Noninvasive Electromagnetic Source Imaging and Granger Causality Analysis: An Electrophysiological Connectome (eConnectome) Approach

OBJECTIVE

Combined source-imaging techniques and directional connectivity analysis can provide useful information about the underlying brain networks in a noninvasive fashion. Source-imaging techniques have been used successfully to either determine the source of activity or to extract source time-courses for Granger causality analysis, previously. In this work, we utilize source-imaging algorithms to both find the network nodes [regions of interest (ROI)] and then extract the activation time series for further Granger causality analysis. The aim of this work is to find network nodes objectively from noninvasive electromagnetic signals, extract activation time-courses, and apply Granger analysis on the extracted series to study brain networks under realistic conditions.

METHODS

Source-imaging methods are used to identify network nodes and extract time-courses and then Granger causality analysis is applied to delineate the directional functional connectivity of underlying brain networks. Computer simulations studies where the underlying network (nodes and connectivity pattern) is known were performed; additionally, this approach has been evaluated in partial epilepsy patients to study epilepsy networks from interictal and ictal signals recorded by EEG and/or Magnetoencephalography (MEG).

RESULTS

Localization errors of network nodes are less than 5 mm and normalized connectivity errors of ∼20% in estimating underlying brain networks in simulation studies. Additionally, two focal epilepsy patients were studied and the identified nodes driving the epileptic network were concordant with clinical findings from intracranial recordings or surgical resection.

CONCLUSION

Our study indicates that combined source-imaging algorithms with Granger causality analysis can identify underlying networks precisely (both in terms of network nodes location and internodal connectivity).

SIGNIFICANCE

The combined source imaging and Granger analysis technique is an effective tool for studying normal or pathological brain conditions.

EEG-MEG Integration Enhances the Characterization of Functional and Effective Connectivity in the Resting State Network

At the sensor level many aspects, such as spectral power, functional and effective connectivity as well as relative-power-ratio ratio (RPR) and spatial resolution have been comprehensively investigated through both electroencephalography (EEG) and magnetoencephalography (MEG). Despite this, differences between both modalities have not yet been systematically studied by direct comparison. It remains an open question as to whether the integration of EEG and MEG data would improve the information obtained from the above mentioned parameters. Here, EEG (64-channel system) and MEG (275 sensor system) were recorded simultaneously in conditions with eyes open (EO) and eyes closed (EC) in 29 healthy adults. Spectral power, functional and effective connectivity, RPR, and spatial resolution were analyzed at five different frequency bands (delta, theta, alpha, beta and gamma). Networks of functional and effective connectivity were described using a spatial filter approach called the dynamic imaging of coherent sources (DICS) followed by the renormalized partial directed coherence (RPDC). Absolute mean power at the sensor level was significantly higher in EEG than in MEG data in both EO and EC conditions. At the source level, there was a trend towards a better performance of the combined EEG+MEG analysis compared with separate EEG or MEG analyses for the source mean power, functional correlation, effective connectivity for both EO and EC. The network of coherent sources and the spatial resolution were similar for both the EEG and MEG data if they were analyzed separately. Results indicate that the combined approach has several advantages over the separate analyses of both EEG and MEG. Moreover, by a direct comparison of EEG and MEG, EEG was characterized by significantly higher values in all measured parameters in both sensor and source level. All the above conclusions are specific to the resting state task and the specific analysis used in this study to have general conclusion multi-center studies would be helpful.

Effects of reconstructed magnetic field from sparse noisy boundary measurements on localization of active neural source

Localization of active neural source (ANS) from measurements on head surface is vital in magnetoencephalography. As neuron-generated magnetic fields are extremely weak, significant uncertainties caused by stochastic measurement interference complicate its localization. This paper presents a novel computational method based on reconstructed magnetic field from sparse noisy measurements for enhanced ANS localization by suppressing effects of unrelated noise. In this approach, the magnetic flux density (MFD) in the nearby current-free space outside the head is reconstructed from measurements through formulating the infinite series solution of the Laplace's equation, where boundary condition (BC) integrals over the entire measurements provide "smooth" reconstructed MFD with the decrease in unrelated noise. Using a gradient-based method, reconstructed MFDs with good fidelity are selected for enhanced ANS localization. The reconstruction model, spatial interpolation of BC, parametric equivalent current dipole-based inverse estimation algorithm using reconstruction, and gradient-based selection are detailed and validated. The influences of various source depths and measurement signal-to-noise ratio levels on the estimated ANS location are analyzed numerically and compared with a traditional method (where measurements are directly used), and it was demonstrated that gradient-selected high-fidelity reconstructed data can effectively improve the accuracy of ANS localization.

Electroencephalography-based real-time cortical monitoring system that uses hierarchical Bayesian estimations for the brain-machine interface

In this study, a real-time cortical activity monitoring system was constructed, which could estimate cortical activities every 125 milliseconds over 2,240 vertexes from 64 channel electroencephalography signals through the Hierarchical Bayesian estimation that uses functional magnetic resonance imaging data as its prior information. Recently, functional magnetic resonance imaging has mostly been used in the neurofeedback field because it allows for high spatial resolution. However, in functional magnetic resonance imaging, the time for the neurofeedback information to reach the patient is delayed several seconds because of its poor temporal resolution. Therefore, a number of problems need to be solved to effectively implement feedback training paradigms in patients. To address this issue, this study used a new cortical activity monitoring system that improved both spatial and temporal resolution by using both functional magnetic resonance imaging data and electroencephalography signals in conjunction with one another. This system is advantageous as it can improve applications in the fields of real-time diagnosis, neurofeedback, and the brain-machine interface.

A multi-subject, multi-modal human neuroimaging dataset

We describe data acquired with multiple functional and structural neuroimaging modalities on the same nineteen healthy volunteers. The functional data include Electroencephalography (EEG), Magnetoencephalography (MEG) and functional Magnetic Resonance Imaging (fMRI) data, recorded while the volunteers performed multiple runs of hundreds of trials of a simple perceptual task on pictures of familiar, unfamiliar and scrambled faces during two visits to the laboratory. The structural data include T1-weighted MPRAGE, Multi-Echo FLASH and Diffusion-weighted MR sequences. Though only from a small sample of volunteers, these data can be used to develop methods for integrating multiple modalities from multiple runs on multiple participants, with the aim of increasing the spatial and temporal resolution above that of any one modality alone. They can also be used to integrate measures of functional and structural connectivity, and as a benchmark dataset to compare results across the many neuroimaging analysis packages. The data are freely available from https://openfmri.org/.

How to use fMRI functional localizers to improve EEG/MEG source estimation

EEG and MEG have excellent temporal resolution, but the estimation of the neural sources that generate the signals recorded by the sensors is a difficult, ill-posed problem. The high spatial resolution of functional MRI makes it an ideal tool to improve the localization of the EEG/MEG sources using data fusion. However, the combination of the two techniques remains challenging, as the neural generators of the EEG/MEG and BOLD signals might in some cases be very different. Here we describe a data fusion approach that was developed by our team over the last decade in which fMRI is used to provide source constraints that are based on functional areas defined individually for each subject. This mini-review describes the different steps that are necessary to perform source estimation using this approach. It also provides a list of pitfalls that should be avoided when doing fMRI-informed EEG/MEG source imaging. Finally, it describes the advantages of using a ROI-based approach for group-level analysis and for the study of sensory systems.

Beamformer source analysis and connectivity on concurrent EEG and MEG data during voluntary movements

Electroencephalography (EEG) and magnetoencephalography (MEG) are the two modalities for measuring neuronal dynamics at a millisecond temporal resolution. Different source analysis methods, to locate the dipoles in the brain from which these dynamics originate, have been readily applied to both modalities alone. However, direct comparisons and possible advantages of combining both modalities have rarely been assessed during voluntary movements using coherent source analysis. In the present study, the cortical and sub-cortical network of coherent sources at the finger tapping task frequency (2–4 Hz) and the modes of interaction within this network were analysed in 15 healthy subjects using a beamformer approach called the dynamic imaging of coherent sources (DICS) with subsequent source signal reconstruction and renormalized partial directed coherence analysis (RPDC). MEG and EEG data were recorded simultaneously allowing the comparison of each of the modalities separately to that of the combined approach. We found the identified network of coherent sources for the finger tapping task as described in earlier studies when using only the MEG or combined MEG+EEG whereas the EEG data alone failed to detect single sub-cortical sources. The signal-to-noise ratio (SNR) level of the coherent rhythmic activity at the tapping frequency in MEG and combined MEG+EEG data was significantly higher than EEG alone. The functional connectivity analysis revealed that the combined approach had more active connections compared to either of the modalities during the finger tapping (FT) task. These results indicate that MEG is superior in the detection of deep coherent sources and that the SNR seems to be more vital than the sensitivity to theoretical dipole orientation and the volume conduction effect in the case of EEG.

The evolution of a disparity decision in human visual cortex

We used fMRI-informed EEG source-imaging in humans to characterize the dynamics of cortical responses during a disparity-discrimination task. After the onset of a disparity-defined target, decision-related activity was found within an extended cortical network that included several occipital regions of interest (ROIs): V4, V3A, hMT+ and the Lateral Occipital Complex (LOC). By using a response-locked analysis, we were able to determine the timing relationships in this network of ROIs relative to the subject's behavioral response. Choice-related activity appeared first in the V4 ROI almost 200 ms before the button press and then subsequently in the V3A ROI. Modeling of the responses in the V4 ROI suggests that this area provides an early contribution to disparity discrimination. Choice-related responses were also found after the button-press in ROIs V4, V3A, LOC and hMT+. Outside the visual cortex, choice-related activity was found in the frontal and temporal poles before the button-press. By combining the spatial resolution of fMRI-informed EEG source imaging with the ability to sort out neural activity occurring before, during and after the behavioral manifestation of the decision, our study is the first to assign distinct functional roles to the extra-striate ROIs involved in perceptual decisions based on disparity, the primary cue for depth.

Adaptive cluster analysis approach for functional localization using magnetoencephalography

In this paper we propose an agglomerative hierarchical clustering Ward's algorithm in tandem with the Affinity Propagation algorithm to reliably localize active brain regions from magnetoencephalography (MEG) brain signals. Reliable localization of brain areas with MEG has been difficult due to variations in signal strength, and the spatial extent of the reconstructed activity. The proposed approach to resolve this difficulty is based on adaptive clustering on reconstructed beamformer images to find locations that are consistently active across different participants and experimental conditions with high spatial resolution. Using data from a human reaching task, we show that the method allows more accurate and reliable localization from MEG data alone without using functional magnetic resonance imaging (fMRI) or any other imaging techniques.

Increasing the accuracy of electromagnetic inverses using functional area source correlation constraints

Estimating cortical current distributions from electroencephalographic (EEG) or magnetoencephalographic data is a difficult inverse problem whose solution can be improved by the addition of priors on the associated neural responses. In the context of visual activation studies, we propose a new approach that uses a functional area constrained estimator (FACE) to increase the accuracy of the reconstructions. It derives the source correlation matrix from a segmentation of the cortex into areas defined by retinotopic maps of the visual field or by functional localizers obtained independently by fMRI. These areas are computed once for each individual subject and the associated estimators can therefore be reused for any new study on the same participant. The resulting FACE reconstructions emphasize the activity of sources within these areas or enforce their intercorrelations. We used realistic Monte-Carlo simulations to demonstrate that this approach improved our estimates of a diverse set of source configurations. Reconstructions obtained from a real EEG dataset demonstrate that our priors improve the localization of the cortical areas involved in horizontal disparity processing.

The primary somatosensory cortex largely contributes to the early part of the cortical response elicited by nociceptive stimuli

Research on the cortical sources of nociceptive laser-evoked brain potentials (LEPs) began almost two decades ago (Tarkka and Treede, 1993). Whereas there is a large consensus on the sources of the late part of the LEP waveform (N2 and P2 waves), the relative contribution of the primary somatosensory cortex (S1) to the early part of the LEP waveform (N1 wave) is still debated. To address this issue we recorded LEPs elicited by the stimulation of four limbs in a large population (n=35). Early LEP generators were estimated both at single-subject and group level, using three different approaches: distributed source analysis, dipolar source modeling, and probabilistic independent component analysis (ICA). We show that the scalp distribution of the earliest LEP response to hand stimulation was maximal over the central-parietal electrodes contralateral to the stimulated side, while that of the earliest LEP response to foot stimulation was maximal over the central-parietal midline electrodes. Crucially, all three approaches indicated hand and foot S1 areas as generators of the earliest LEP response. Altogether, these findings indicate that the earliest part of the scalp response elicited by a selective nociceptive stimulus is largely explained by activity in the contralateral S1, with negligible contribution from the secondary somatosensory cortex (S2).

Mapping the signal-to-noise-ratios of cortical sources in magnetoencephalography and electroencephalography

Although magnetoencephalography (MEG) and electroencephalography (EEG) have been available for decades, their relative merits are still debated. We examined regional differences in signal-to-noise-ratios (SNRs) of cortical sources in MEG and EEG. Data from four subjects were used to simulate focal and extended sources located on the cortical surface reconstructed from high-resolution magnetic resonance images. The SNR maps for MEG and EEG were found to be complementary. The SNR of deep sources was larger in EEG than in MEG, whereas the opposite was typically the case for superficial sources. Overall, the SNR maps were more uniform for EEG than for MEG. When using a noise model based on uniformly distributed random sources on the cortex, the SNR in MEG was found to be underestimated, compared with the maps obtained with noise estimated from actual recorded MEG and EEG data. With extended sources, the total area of cortex in which the SNR was higher in EEG than in MEG was larger than with focal sources. Clinically, SNR maps in a patient explained differential sensitivity of MEG and EEG in detecting epileptic activity. Our results emphasize the benefits of recording MEG and EEG simultaneously.

MEG and EEG data fusion: simultaneous localisation of face-evoked responses

We present an empirical Bayesian scheme for distributed multimodal inversion of electromagnetic forward models of EEG and MEG signals. We used a generative model with common source activity and separate error components for each modality. Under this scheme, the weightings of error for each modality, relative to source components, are estimated automatically from the data, by optimising the model-evidence. This obviates the need for arbitrary user-defined weightings. To evaluate the scheme, we acquired three types of data simultaneously from twelve participants: total magnetic flux (as recorded by 102 magnetometers), orthogonal in-plane gradients of the magnetic field (as recorded by 204 planar gradiometers) and voltage differences in the electrical field (recorded by 70 electrodes). We assessed the relative precision of each sensor-type in terms of signal-to-noise ratio (SNR); using empirical sample variances and optimised estimators from the generative model. We then compared the localisation of face-evoked responses, using each modality separately, with that obtained by their "fusion" under the common generative model. Finally, we quantified the conditional precisions of the source estimates using their posterior covariance, confirming that EEG can improve MEG-based source reconstructions.

Estimation of effective and functional cortical connectivity from neuroelectric and hemodynamic recordings

In this paper, different linear and nonlinear methodologies for the estimation of cortical connectivity from neuroelectric and hemodynamic measurements are reviewed and applied on common data set in order to highlight similarities and differences in the results. Different effective and functional connectivity methods were applied to motor and cognitive data sets, including structural equation modeling (SEM), directed transfer function (DTF), partial directed coherence (PDC), and direct directed transfer function (dDTF). Comparisons were made between the results in order to understand if, for a same dataset, effective and functional connectivity estimators can return the same cortical connectivity patterns. An application of a nonlinear method [phase synchronization index (PSI)] to similar executed and imagined movements was also reviewed. Connectivity patterns estimated with the use of the neuroelectric information and of the information from the multimodal integration of neuroelectric and hemodynamic data were also compared. Results suggests that the estimation of the cortical connectivity patterns performed with the linear methods (SEM, DTF, PDC, dDTF) or with the nonlinear method (PSI) on movement related potentials returned similar cortical networks. Differences in cortical connectivity were noted between the patterns estimated with the use of multimodal integration and those estimated by using only the neuroelectric data.

Motor cortical responsiveness to attempted movements in tetraplegia: evidence from neuroelectrical imaging

OBJECTIVE

The maintenance of a motor cortical program in the temporal domain is relevant to current neuroinformatic efforts to use non-invasive EEG signals to control neuroprosthetic devices designed to restore natural movements of paralyzed body parts. Here we use an advance neuroelectrical imaging approach to examine the motor cortical responsiveness in human tetraplegia.

METHODS

High resolution-electroencephalographic (EEG) recordings were performed in five subjects with tetraplegia due to chronic, complete spinal cord injuries (SCIs) while they attempted self-generated movements of a plegic body part (foot), and in five healthy subjects executing simple foot movements.

RESULTS

Self-generated movement attempts induced significant EEG sources of activity in a set of motor-related areas (including the primary motor area, MI) similar to what observed during the preparatory stages of movement execution (control subjects). Functional connectivity showed a preferential interaction between the "non-primary" motor areas and the putative MI foot site, as estimated for both motor execution and attempt. Under this latter condition however, it could be observed an "enlargement" of the functional network by including the left superior parietal cortex.

CONCLUSIONS

Our findings indicate the existence of a functional circuit subserving the attempted motion in SCI subjects that encompasses a set of areas known to play a role in motor execution, yet reveals differences in the functional interaction between these areas.

SIGNIFICANCE

The understanding of changes in the motor circuitry is relevant to current neuroinformatic efforts to use non-invasive EEG signals to control neuroprosthetic devices designed to benefit paralyzed persons.

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.

Estimate of causality between independent cortical spatial patterns during movement volition in spinal cord injured patients

Static hemodynamic or neuroelectric images of brain regions activated during particular tasks do not convey the information of how these regions communicate to each other. Cortical connectivity estimation aims at describing these interactions as connectivity patterns which hold the direction and strength of the information flow between cortical areas. In this study, we attempted to estimate the causality between distributed cortical systems during a movement volition task in preparation for execution of simple movements by a group of normal healthy subjects and by a group of Spinal Cord Injured (SCI) patients. To estimate the causality between the spatial distributed patterns of cortical activity in the frequency domain, we applied a series of processing steps on the recorded EEG data. From the high-resolution EEG recordings we estimated the cortical waveforms for the regions of interest (ROIs), each representing a selected sensor group population. The solutions of the linear inverse problem returned a series of cortical waveforms for each ROI considered and for each trial analyzed. For each subject, the cortical waveforms were then subjected to Independent Component Analysis (ICA) pre-processing. The independent components obtained by the application of the ThinICA algorithm were further processed by a Partial Directed Coherence algorithm, in order to extract the causality between spatial cortical patterns of the estimated data. The source-target cortical dependencies found in the group of normal subjects were relatively similar in all frequency bands analyzed. For the normal subjects we observed a common source pattern in an ensemble of cortical areas including the right parietal and right lip primary motor areas and bilaterally the primary foot and posterior SMA areas. The target of this cortical network, in the Granger-sense of causality, was shown to be a smaller network composed mostly by the primary foot motor areas and the posterior SMA bilaterally. In the case of the SCI population, both the source and the target cortical patterns had larger sizes than in the normal population. The source cortical areas included always the primary foot and lip motor areas, often bilaterally. In addition, the right parietal area and the bilateral premotor area 6 were also involved. Again, the patterns remained substantially stable across the different frequency bands analyzed. The target cortical patterns observed in the SCI population had larger extensions when compared to the normal ones, since in most cases they involved the bilateral activation of the primary foot movement areas as well as the SMA, the primary lip areas and the parietal cortical areas.

A novel integrated MEG and EEG analysis method for dipolar sources

The ability of magnetoencephalography (MEG) to accurately localize neuronal currents and obtain tangential components of the source is largely due to MEG's insensitivity to the conductivity profile of the head tissues. However, MEG cannot reliably detect the radial component of the neuronal current. In contrast, the localization accuracy of electroencephalography (EEG) is not as good as MEG, but EEG can detect both the tangential and radial components of the source. In the present study, we investigated the conductivity dependence in a new approach that combines MEG and EEG to accurately obtain, not only the location and tangential components, but also the radial component of the source. In this approach, the source location and tangential components are obtained from MEG alone, and optimal conductivity values of the EEG model are estimated by best-fitting EEG signal, while precisely matching the tangential components of the source in EEG and MEG. Then, the radial components are obtained from EEG using the previously estimated optimal conductivity values. Computer simulations testing this integrated approach demonstrated two main findings. First, there are well-organized optimal combinations of the conductivity values that provide an accurate fit to the combined MEG and EEG data. Second, the radial component, in addition to the location and tangential components, can be obtained with high accuracy without needing to know the precise conductivity profile of the head. We then demonstrated that this new approach performed reliably in an analysis of the 20-ms component from human somatosensory responses elicited by electric median-nerve stimulation.

Comparison of different cortical connectivity estimators for high-resolution EEG recordings

The aim of this work is to characterize quantitatively the performance of a body of techniques in the frequency domain for the estimation of cortical connectivity from high-resolution EEG recordings in different operative conditions commonly encountered in practice. Connectivity pattern estimators investigated are the Directed Transfer Function (DTF), its modification known as direct DTF (dDTF) and the Partial Directed Coherence (PDC). Predefined patterns of cortical connectivity were simulated and then retrieved by the application of the DTF, dDTF, and PDC methods. Signal-to-noise ratio (SNR) and length (LENGTH) of EEG epochs were studied as factors affecting the reconstruction of the imposed connectivity patterns. Reconstruction quality and error rate in estimated connectivity patterns were evaluated by means of some indexes of quality for the reconstructed connectivity pattern. The error functions were statistically analyzed with analysis of variance (ANOVA). The whole methodology was then applied to high-resolution EEG data recorded during the well-known Stroop paradigm. Simulations indicated that all three methods correctly estimated the simulated connectivity patterns under reasonable conditions. However, performance of the methods differed somewhat as a function of SNR and LENGTH factors. The methods were generally equivalent when applied to the Stroop data. In general, the amount of available EEG affected the accuracy of connectivity pattern estimations. Analysis of 27 s of nonconsecutive recordings with an SNR of 3 or more ensured that the connectivity pattern could be accurately recovered with an error below 7% for the PDC and 5% for the DTF. In conclusion, functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in most EEG recordings by combining high-resolution EEG techniques, linear inverse estimation of the cortical activity, and frequency domain multivariate methods such as PDC, DTF, and dDTF.

Assessing cortical functional connectivity by partial directed coherence: simulations and application to real data

The aim of this paper is to test a technique called partial directed coherence (PDC) and its modification (squared PDC; sPDC) for the estimation of human cortical connectivity by means of simulation study, in which both PDC and sPDC were studied by analysis of variance. The statistical analysis performed returned that both PDC and sPDC are able to estimate correctly the imposed connectivity patterns when data exhibit a signal-to-noise ratio of at least 3 and a length of at least 27 s of nonconsecutive recordings at 250 Hz of sampling rate, equivalent, more generally, to 6750 data samples.

Motor-related cortical dynamics to intact movements in tetraplegics as revealed by high-resolution EEG

We explored the cortical dynamics during movements of an unaffected body part in tetraplegic subjects with chronic spinal cord injury (SCI). The aims were to find out whether the intact movements were associated with a physiological time-varying pattern of activity in the motor-related cortical areas and whether the primary motor area (MI) activation followed a somatotopic distribution. Event-related potentials to self-initiated lip movements were analyzed by means of cortical source imaging of EEG recorded from seven tetraplegic subjects and seven control subjects. Regions of interest (ROIs) were selected on individual MRI and the time-varying electrophysiologic activity (cortical current density, CCD) was estimated on these ROIs and subjected to across-subject analysis. A significant, bilateral movement-related pattern of MI activation was detected during motor task execution in SCI patients as well as in controls. The site of local maxima activation displayed a symmetrical discrete distribution within MI, consistently with a putative somatotopic lip representation, in all the subjects. The supplementary motor area proper (SMAp) was always coactivated with MI and coactivation was characterized by a time course with typical premotion and motion phases over both motor areas. A clear-cut temporal delay between the SMAp and MI activation did not occur either in SCI patients or in controls. These findings obtained with noninvasive neuroelectrical source imaging document that in chronic SCI subjects "executive" motor areas are engaged with a preserved temporal and spatial pattern during preparation and execution of intact movements.

MEG source localization under multiple constraints: An extended Bayesian framework

To use Electroencephalography (EEG) and Magnetoencephalography (MEG) as functional brain 3D imaging techniques, identifiable distributed source models are required. The reconstruction of EEG/MEG sources rests on inverting these models and is ill-posed because the solution does not depend continuously on the data and there is no unique solution in the absence of prior information or constraints. We have described a general framework that can account for several priors in a common inverse solution. An empirical Bayesian framework based on hierarchical linear models was proposed for the analysis of functional neuroimaging data [Friston, K., Penny, W., Phillips, C., Kiebel, S., Hinton, G., Ashburner, J., 2002. Classical and Bayesian inference in neuroimaging: theory. NeuroImage 16, 465-483] and was evaluated recently in the context of EEG [Phillips, C., Mattout, J., Rugg, M.D., Maquet, P., Friston, K., 2005. An empirical Bayesian solution to the source reconstruction problem in EEG. NeuroImage 24, 997-1011]. The approach consists of estimating the expected source distribution and its conditional variance that is constrained by an empirically determined mixture of prior variance components. Estimation uses Expectation-Maximization (EM) to give the Restricted Maximum Likelihood (ReML) estimate of the variance components (in terms of hyperparameters) and the Maximum A Posteriori (MAP) estimate of the source parameters. In this paper, we extend the framework to compare different combinations of priors, using a second level of inference based on Bayesian model selection. Using Monte-Carlo simulations, ReML is first compared to a classic Weighted Minimum Norm (WMN) solution under a single constraint. Then, the ReML estimates are evaluated using various combinations of priors. Both standard criterion and ROC-based measures were used to assess localization and detection performance. The empirical Bayes approach proved useful as: (1) ReML was significantly better than WMN for single priors; (2) valid location priors improved ReML source localization; (3) invalid location priors did not significantly impair performance. Finally, we show how model selection, using the log-evidence, can be used to select the best combination of priors. This enables a global strategy for multiple prior-based regularization of the MEG/EEG source reconstruction.

Principal components analysis of Laplacian waveforms as a generic method for identifying ERP generator patterns: II. Adequacy of low-density estimates

OBJECTIVE

To evaluate the comparability of high- and low-density surface Laplacian estimates for determining ERP generator patterns of group data derived from a typical ERP sample size and paradigm.

METHODS

High-density ERP data (129 sites) recorded from 17 adults during tonal and phonetic oddball tasks were converted to a 10-20-system EEG montage (31 sites) using spherical spline interpolations. Current source density (CSD) waveforms were computed from the high- and low-density, but otherwise identical, ERPs, and correlated at corresponding locations. CSD data were submitted to separate covariance-based, unrestricted temporal PCAs (Varimax of covariance loadings) to identify and effectively summarize temporally and spatially overlapping CSD components. Solutions were compared by correlating factor loadings and scores, and by plotting ANOVA F statistics derived from corresponding high- and low-resolution factor scores using representative sites.

RESULTS

High- and low-density CSD waveforms, PCA solutions, and F statistics were remarkably similar, yielding correlations of .9 < or = r < or = .999 between waveforms, loadings, and scores for almost all comparisons at low-density locations except for low-signal CSD waveforms at occipital sites. Each of the first 10 high-density factors corresponded precisely to one factor of the first 10 low-density factors, with each 10-factor set accounting for the meaningful CSD variance (> 91.6%).

CONCLUSIONS

Low-density surface Laplacian estimates were shown to be accurate approximations of high-density CSDs at these locations, which adequately and quite sufficiently summarized group data. Moreover, reasonable approximations of many high-density scalp locations were obtained for group data from interpolations of low-density data. If group findings are the primary objective, as typical for cognitive ERP research, low-resolution CSD topographies may be as efficient, given the effective spatial smoothing when averaging across subjects and/or conditions.

SIGNIFICANCE

Conservative recommendations for restricting surface Laplacians to high-density recordings may not be appropriate for all ERP research applications, and should be re-evaluated considering objective, costs and benefits.

Whole-head MEG analysis of cortical spatial organization from unilateral stimulation of median nerve in both hands: No complete hemispheric homology

We examined the contralateral hemispheric cortical activity in MEG (151 ch) after unilateral median nerve stimulation of the right and left hand in twenty healthy right-handed subjects. The goal was to establish parameters to describe cortical activity of the hemispheric responses and to study the potential ability to assess differences in volunteers and patients. We focused on the within-subject similarity and differences between evoked fields in both hands. Cortical activity was characterized by the overlay display of waveforms (CWP), number of peak stages, loci of focal maxima and minima in each stage, 3D topographic maps and exemplified equivalent current dipole characteristics. The paired-wise test was used to analyze the hemispheric differences. The waveform morphology was unique across the subjects, similar CWPs were noted in both hemispheres of the individual. The contralateral hemispheric responses showed a well defined temporal-spatial activation of six to seven stages in the 500 ms window. Consistently (in over 80% of subjects), the six stages across the subjects were 20M, 30M, 50M, 70M, 90M, and 150M. A 240M was present in the left hemisphere (LH) in 15/20 subjects and in the right hemisphere (RH) in 10/20. Statistics of the latencies and amplitudes of these seven stages were calculated. Our results indicated that the latency was highly consistent and exhibited no statistical mean difference for all stages. Furthermore, no mean amplitude differences between both hemispheres at each stage were found. The patterns of magnetic fields in both hemispheres were consistent in 70% of the subjects. A laterality index (L.I.) was used for defining the magnetic field amplitude differences between two hemispheres for each individual. Overall, the absolute amplitude of the brain responses was larger in the left than in the right hemisphere in the majority of subjects (16/20), yet a significant portion (4/20) exhibited right dominance of the N20m activity. Each individual exhibited a unique CWP, there was reliable consistency of peak latencies and mean amplitudes in median nerve MEG. Nevertheless, this study indicates the limitations of using the intact hemisphere responses to compare with those from the affected (brain) side and suggests caution in assuming full homology in the cortical organization of both hemispheres. This study provides some results to address clinical issues like which parameter describes individual differences best. Whether a genuine difference is found or whether any difference may simply represent the variability encountered in a normal population.

Estimation of the effective and functional human cortical connectivity with structural equation modeling and directed transfer function applied to high-resolution EEG

Different brain imaging devices are presently available to provide images of the human functional cortical activity, based on hemodynamic, metabolic or electromagnetic measurements. However, static images of brain regions activated during particular tasks do not convey the information of how these regions are interconnected. The concept of brain connectivity plays a central role in the neuroscience, and different definitions of connectivity, functional and effective, have been adopted in literature. While the functional connectivity is defined as the temporal coherence among the activities of different brain areas, the effective connectivity is defined as the simplest brain circuit that would produce the same temporal relationship as observed experimentally among cortical sites. The structural equation modeling (SEM) is the most used method to estimate effective connectivity in neuroscience, and its typical application is on data related to brain hemodynamic behavior tested by functional magnetic resonance imaging (fMRI), whereas the directed transfer function (DTF) method is a frequency-domain approach based on both a multivariate autoregressive (MVAR) modeling of time series and on the concept of Granger causality. This study presents advanced methods for the estimation of cortical connectivity by applying SEM and DTF on the cortical signals estimated from high-resolution electroencephalography (EEG) recordings, since these signals exhibit a higher spatial resolution than conventional cerebral electromagnetic measures. To estimate correctly the cortical signals, we used a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual MRI, a distributed source model and a regularized linear inverse source estimates of cortical current density. Before the application of SEM and DTF methodology to the cortical waveforms estimated from high-resolution EEG data, we performed a simulation study, in which different main factors (signal-to-noise ratio, SNR, and simulated cortical activity duration, LENGTH) were systematically manipulated in the generation of test signals, and the errors in the estimated connectivity were evaluated by the analysis of variance (ANOVA). The statistical analysis returned that during simulations, both SEM and DTF estimators were able to correctly estimate the imposed connectivity patterns under reasonable operative conditions, that is, when data exhibit an SNR of at least 3 and a LENGTH of at least 75 s of nonconsecutive EEG recordings at 64 Hz of sampling rate. Hence, effective and functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in any practical EEG recordings, by combining high-resolution EEG techniques and linear inverse estimation with SEM or DTF methods. We conclude that the estimation of cortical connectivity can be performed not only with hemodynamic measurements, but also with EEG signals treated with advanced computational techniques.

Estimation of the Cortical Connectivity by High-Resolution EEG and Structural Equation Modeling: Simulations and Application to Finger Tapping Data

Today, the concept of brain connectivity plays a central role in the neuroscience. While functional connectivity is defined as the temporal coherence between the activities of different brain areas, the effective connectivity is defined as the simplest brain circuit that would produce the same temporal relationship as observed experimentally between cortical sites. The most used method to estimate effective connectivity in neuroscience is the structural equation modeling (SEM), typically used on data related to the brain hemodynamic behavior. However, the use of hemodynamic measures limits the temporal resolution on which the brain process can be followed. The present research proposes the use of the SEM approach on the cortical waveforms estimated from the high-resolution EEG data, which exhibits a good spatial resolution and a higher temporal resolution than hemodynamic measures. We performed a simulation study, in which different main factors were systematically manipulated in the generation of test signals, and the errors in the estimated connectivity were evaluated by the analysis of variance (ANOVA). Such factors were the signal-to-noise ratio and the duration of the simulated cortical activity. Since SEM technique is based on the use of a model formulated on the basis of anatomical and physiological constraints, different experimental conditions were analyzed, in order to evaluate the effect of errors made in the a priori model formulation on its performances. The feasibility of the proposed approach has been shown in a human study using high-resolution EEG recordings related to finger tapping movements.

Evaluation of different cortical source localization methods using simulated and experimental EEG data

Different cortical source localization methods have been developed to directly link the scalp potentials with the cortical activities. Up to now, these methods are the only possible solution to noninvasively investigate cortical activities with both high spatial and time resolutions. However, the application of these methods is hindered by the fact that they have not been rigorously evaluated nor compared. In this paper, the performances of several source localization methods (moving dipoles, minimum Lp norm, and low resolution tomography (LRT) with Lp norm, p equal to 1, 1.5, and 2) were evaluated by using simulated scalp EEG data, scalp somatosensory evoked potentials (SEPs), and upper limb motor-related potentials (MRPs) obtained on human subjects (all with 163 scalp electrodes). By using simulated EEG data, we first evaluated the source localization ability of the above methods quantitatively. Subsequently, the performance of the various methods was evaluated qualitatively by using experimental SEPs and MRPs. Our results show that the overall LRT Lp norm method with p equal to 1 has a better source localization ability than any of the other investigated methods and provides physiologically meaningful reconstruction results. Our evaluation results provide useful information for choosing cortical source localization approaches for future EEG/MEG studies.

Assessing cortical functional connectivity by linear inverse estimation and directed transfer function: simulations and application to real data

OBJECTIVE

To test a technique called Directed Transfer Function (DTF) for the estimation of human cortical connectivity, by means of simulation study and human study, using high resolution EEG recordings related to finger movements.

METHODS

The method of the Directed Transfer Function (DTF) is a frequency-domain approach, based on a multivariate autoregressive modeling of time series and on the concept of Granger causality. Since the spreading of the potential from the cortex to the sensors makes it difficult to infer the relation between the spatial patterns on the sensor space and those on the cortical sites, we propose the use of the DTF method on cortical signals estimated from high resolution EEG recordings, which exhibit a higher spatial resolution than conventional cerebral electromagnetic measures. The simulation study was followed by an analysis of variance (ANOVA) of the results obtained for different levels of Signal to Noise Ratio (SNR) and temporal length, as they have been systematically imposed on simulated signals. The whole methodology was then applied to high resolution EEG data recorded during a visually paced finger movement.

RESULTS

The statistical analysis performed returns that during simulations, DTF is able to estimate correctly the imposed connectivity patterns under reasonable operative conditions, i.e. when data exhibit a SNR of at least 3 and a length of at least 75 s of non-consecutive recordings at 64 Hz of sampling rate, equivalent, more generally, to 4800 data samples.

CONCLUSIONS

Functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in any practical EEG recordings, by combining high resolution EEG techniques, linear inverse estimation and the DTF method.

SIGNIFICANCE

The estimation of cortical connectivity can be performed not only with hemodynamic measurements, by using functional MRI recordings, but also with modern EEG recordings treated with advanced computational techniques.