José David López

Algorithmic procedures for Bayesian MEG/EEG source reconstruction in SPM

The MEG/EEG inverse problem is ill-posed, giving different source reconstructions depending on the initial assumption sets. Parametric Empirical Bayes allows one to implement most popular MEG/EEG inversion schemes (Minimum Norm, LORETA, etc.) within the same generic Bayesian framework. It also provides a cost-function in terms of the variational Free energy—an approximation to the marginal likelihood or evidence of the solution. In this manuscript, we revisit the algorithm for MEG/EEG source reconstruction with a view to providing a didactic and practical guide. The aim is to promote and help standardise the development and consolidation of other schemes within the same framework. We describe the implementation in the Statistical Parametric Mapping (SPM) software package, carefully explaining each of its stages with the help of a simple simulated data example. We focus on the Multiple Sparse Priors (MSP) model, which we compare with the well-known Minimum Norm and LORETA models, using the negative variational Free energy for model comparison. The manuscript is accompanied by Matlab scripts to allow the reader to test and explore the underlying algorithm.

High precision anatomy for MEG

Precise MEG estimates of neuronal current flow are undermined by uncertain knowledge of the head location with respect to the MEG sensors. This is either due to head movements within the scanning session or systematic errors in co-registration to anatomy. Here we show how such errors can be minimized using subject-specific head-casts produced using 3D printing technology. The casts fit the scalp of the subject internally and the inside of the MEG dewar externally, reducing within session and between session head movements. Systematic errors in matching to MRI coordinate system are also reduced through the use of MRI-visible fiducial markers placed on the same cast. Bootstrap estimates of absolute co-registration error were of the order of 1 mm. Estimates of relative co-registration error were < 1.5 mm between sessions. We corroborated these scalp based estimates by looking at the MEG data recorded over a 6 month period. We found that the between session sensor variability of the subjects evoked response was of the order of the within session noise, showing no appreciable noise due to between-session movement. Simulations suggest that the between-session sensor level amplitude SNR improved by a factor of 5 over conventional strategies. We show that at this level of coregistration accuracy there is strong evidence for anatomical models based on the individual rather than canonical anatomy; but that this advantage disappears for errors of greater than 5 mm. This work paves the way for source reconstruction methods which can exploit very high SNR signals and accurate anatomical models; and also significantly increases the sensitivity of longitudinal studies with MEG.

Cortical dynamics and subcortical signatures of motor-language coupling in Parkinson’s disease

Impairments of action language have been documented in early stage Parkinson’s disease (EPD). The action-sentence compatibility effect (ACE) paradigm has revealed that EPD involves deficits to integrate action-verb processing and ongoing motor actions. Recent studies suggest that an abolished ACE in EPD reflects a cortico-subcortical disruption, and recent neurocognitive models highlight the role of the basal ganglia (BG) in motor-language coupling. Building on such breakthroughs, we report the first exploration of convergent cortical and subcortical signatures of ACE in EPD patients and matched controls. Specifically, we combined cortical recordings of the motor potential, functional connectivity measures, and structural analysis of the BG through voxel-based morphometry. Relative to controls, EPD patients exhibited an impaired ACE, a reduced motor potential, and aberrant frontotemporal connectivity. Furthermore, motor potential abnormalities during the ACE task were predicted by overall BG volume and atrophy. These results corroborate that motor-language coupling is mainly subserved by a cortico-subcortical network including the BG as a key hub. They also evince that action-verb processing may constitute a neurocognitive marker of EPD. Our findings suggest that research on the relationship between language and motor domains is crucial to develop models of motor cognition as well as diagnostic and intervention strategies.

Discrimination of cortical laminae using MEG

Typically MEG source reconstruction is used to estimate the distribution of current flow on a single anatomically derived cortical surface model. In this study we use two such models representing superficial and deep cortical laminae. We establish how well we can discriminate between these two different cortical layer models based on the same MEG data in the presence of different levels of co-registration noise, Signal-to-Noise Ratio (SNR) and cortical patch size. We demonstrate that it is possible to make a distinction between superficial and deep cortical laminae for levels of co-registration noise of less than 2 mm translation and 2° rotation at SNR > 11 dB. We also show that an incorrect estimate of cortical patch size will tend to bias layer estimates. We then use a 3D printed head-cast (Troebinger et al., 2014) to achieve comparable levels of co-registration noise, in an auditory evoked response paradigm, and show that it is possible to discriminate between these cortical layer models in real data.

Automatic recognition of anuran species based on syllable identification

Abstract Monitoring of biological populations is well known for being a complex task that involves high operational costs, unknown reproductive intervals of the studied species, and difficult visualization of isolated individuals (due to their mimetic and cryptic capabilities). Therefore, the development of new methodologies able to measure quantities of individuals in specific biological populations without direct contact is desired. Species and individual recognition, based on acoustic analysis of their calls (Bioacoustics), is possible for many animals and has proven to be a useful tool in the study and monitoring of animal species. In this paper, an unsupervised methodology for anuran automatic identification is proposed; it is based on the use of a fuzzy classifier and Mel Frequency Cepstral Coefficients. This methodology is able to detect species not presented in the training stage, although they belong to different populations. Additionally, correlations among species of the same genus can be determined through the similarities of their calls. For testing the proposed method, two different datasets with species from the northeastern Colombia (Choco and Antioquia departments with 103 and 813 mating calls respectively) were used. In validation tests performed, accuracies between 99.38% and 100% were achieved in all species by applying the proposed methodology to both datasets. Thirteen different species of anurans in both datasets were correctly identified.

A general Bayesian treatment for MEG source reconstruction incorporating lead field uncertainty

There is uncertainty introduced when a cortical surface based model derived from an anatomical MRI is used to reconstruct neural activity with MEG data. This is a specific case of a problem with uncertainty in parameters on which M/EEG lead fields depend non-linearly. Here we present a general mathematical treatment of any such problem with a particular focus on co-registration. We use a Metropolis search followed by Bayesian Model Averaging over multiple sparse prior source inversions with different headlocation/orientation parameters. Based on MEG data alone we can locate the cortex to within 4 mm at empirically realistic signal to noise ratios. We also show that this process gives improved posterior distributions on the estimated current distributions, and can be extended to make inference on the locations of local maxima by providing confidence intervals for each source.

SisFall: A Fall and Movement Dataset

Research on fall and movement detection with wearable devices has witnessed promising growth. However, there are few publicly available datasets, all recorded with smartphones, which are insufficient for testing new proposals due to their absence of objective population, lack of performed activities, and limited information. Here, we present a dataset of falls and activities of daily living (ADLs) acquired with a self-developed device composed of two types of accelerometer and one gyroscope. It consists of 19 ADLs and 15 fall types performed by 23 young adults, 15 ADL types performed by 14 healthy and independent participants over 62 years old, and data from one participant of 60 years old that performed all ADLs and falls. These activities were selected based on a survey and a literature analysis. We test the dataset with widely used feature extraction and a simple to implement threshold based classification, achieving up to 96% of accuracy in fall detection. An individual activity analysis demonstrates that most errors coincide in a few number of activities where new approaches could be focused. Finally, validation tests with elderly people significantly reduced the fall detection performance of the tested features. This validates findings of other authors and encourages developing new strategies with this new dataset as the benchmark.

Bayesian model selection of template forward models for EEG source reconstruction.

Several EEG source reconstruction techniques have been proposed to identify the generating neuronal sources of electrical activity measured on the scalp. The solution of these techniques depends directly on the accuracy of the forward model that is inverted. Recently, a parametric empirical Bayesian (PEB) framework for distributed source reconstruction in EEG/MEG was introduced and implemented in the Statistical Parametric Mapping (SPM) software. The framework allows us to compare different forward modeling approaches, using real data, instead of using more traditional simulated data from an assumed true forward model. In the absence of a subject specific MR image, a 3-layered boundary element method (BEM) template head model is currently used including a scalp, skull and brain compartment. In this study, we introduced volumetric template head models based on the finite difference method (FDM). We constructed a FDM head model equivalent to the BEM model and an extended FDM model including CSF. These models were compared within the context of three different types of source priors related to the type of inversion used in the PEB framework: independent and identically distributed (IID) sources, equivalent to classical minimum norm approaches, coherence (COH) priors similar to methods such as LORETA, and multiple sparse priors (MSP). The resulting models were compared based on ERP data of 20 subjects using Bayesian model selection for group studies. The reconstructed activity was also compared with the findings of previous studies using functional magnetic resonance imaging. We found very strong evidence in favor of the extended FDM head model with CSF and assuming MSP. These results suggest that the use of realistic volumetric forward models can improve PEB EEG source reconstruction.

A multiobjective-based switching topology for hierarchical model predictive control applied to a hydro-power valley

In a Hierarchical Model Predictive Control (H-MPC) framework, this paper explores suitable time-variant structures for the hierarchies of different local MPC controllers. The idea is to adapt to different operational conditions by changing the importance of the local controllers. This is done by defining the level of the hierarchy they belong to, and solving within each level the local MPC problem using the information provided by the higher levels at the current time step and the predicted information from the lower levels obtained in the previous time step. As selecting a hierarchy results in a combinatorial optimization problem, it is explicitly solved for a limited number of relevant topologies only and then the switching between topologies is defined with a multiobjective optimizer, so as to decide the best H-MPC scheme according to the expected performance. A comparison with fixed-topology H-MPC controllers is done, showing the effectiveness of the proposed approach for the power control of a hydro-power valley.

Artificial Neural Networks as an alternative to traditional fall detection methods

Falls are common events among older adults and may have serious consequences. Automatic fall detection systems are becoming a popular tool to rapidly detect such events, helping family or health personal to rapidly help the person that falls. This paper presents the results obtained in the process of testing a new fall detection method, based on Artificial Neural Networks (ANN). This method intends to improve fall detection accuracy, by avoiding the traditional threshold - based fall detection methods, and introducing ANN as a suitable option on this application.Also ANN have low computational cost, this characteristic makes them easy to implement on a portable device, comfortable to be wear by the patient.