Evaluation of different measures of functional connectivity using a neural mass model

We use a neural mass model to address some important issues in characterising functional integration among remote cortical areas using magnetoencephalography or electroencephalography (MEG or EEG). In a previous paper [Neuroimage (in press)], we showed how the coupling among cortical areas can modulate the MEG or EEG spectrum and synchronise oscillatory dynamics. In this work, we exploit the model further by evaluating different measures of statistical dependencies (i.e., functional connectivity) among MEG or EEG signals that are mediated by neuronal coupling. We have examined linear and nonlinear methods, including phase synchronisation. Our results show that each method can detect coupling but with different sensitivity profiles that depended on (i) the frequency specificity of the interaction (broad vs. narrow band) and (ii) the nature of the coupling (linear vs. nonlinear). Our analyses suggest that methods based on the concept of generalised synchronisation are the most sensitive when interactions encompass different frequencies (broadband analyses). In the context of narrow-band analyses, mutual information was found to be the most sensitive way to disclose frequency-specific couplings. Measures based on generalised synchronisation and phase synchronisation are the most sensitive to nonlinear coupling. These different sensitivity profiles mean that the choice of coupling measures can have dramatic effects on the cortical networks identified. We illustrate this using a single-subject MEG study of binocular rivalry and highlight the greater recovery of statistical dependencies among cortical areas in the beta band when mutual information is used.