Development of a variational scheme for model inversion of multi-area model of brain. Part I: simulation evaluation

We previously developed an integrated model of the brain within a single cortical area for functional Magnetic Resonance Imaging (fMRI), electroencephalography (EEG), and magnetoencephalography (MEG) using an extended neural mass model (ENMM). We then extended ENMM from a single-area to a multi-area model to develop a neural mass model of the entire brain. To this end, we derived a nonlinear state-space representation of the multi-area model. In Parts I and II of these two companion papers (henceforth called Part I and Part II), we develop and evaluate a variational Bayesian expectation maximization (VBEM) method to estimate parameters of multi-area ENMM (MEN) using E/MEG data. In Part I, we derive a state-space representation of MEN and use VBEM method for model inversion (parameter estimation). We evaluate and validate performance of VBEM method for model inversion of MEN using simulation studies in various signal-to-noise ratios. Details of VBEM method are presented in Part II. The proposed approach provides a useful technique for analyzing effective connectivity using non-invasive EEG and MEG methods.