Miguel Angel Mañanas

A comparative study of automatic techniques for ocular artifact reduction in spontaneous EEG signals based on clinical target variables: A simulation case

Eye movement artifacts represent a critical issue for quantitative electroencephalography (EEG) analysis and a number of mathematical approaches have been proposed to reduce their contribution in EEG recordings. The aim of this paper was to objectively and quantitatively evaluate the performance of ocular filtering methods with respect to spectral target variables widely used in clinical and functional EEG studies. In particular the following methods were applied: regression analysis and some blind source separation (BSS) techniques based on second-order statistics (PCA, AMUSE and SOBI) and on higher-order statistics (JADE, INFOMAX and FASTICA). Considering blind source decomposition methods, a completely automatic procedure of BSS based on logical rules related to spectral and topographical information was proposed in order to identify the components related to ocular interference. The automatic procedure was applied in different montages of simulated EEG and electrooculography (EOG) recordings: a full montage with 19 EEG and 2 EOG channels, a reduced one with only 6 EEG leads and a third one where EOG channels were not available. Time and frequency results in all of them indicated that AMUSE and SOBI algorithms preserved and recovered more brain activity than the other methods mainly at anterior regions. In the case of full montage: (i) errors were lower than 5% for all spectral variables at anterior sites; and (ii) the highest improvement in the signal-to-artifact (SAR) ratio was obtained up to 40dB at these anterior sites. Finally, we concluded that second-order BSS-based algorithms (AMUSE and SOBI) provided an effective technique for eye movement removal even when EOG recordings were not available or when data length was short.

High-density surface EMG maps from upper-arm and forearm muscles

BackgroundsEMG signal has been widely used in different applications in kinesiology and rehabilitation as well as in the control of human-machine interfaces. In general, the signals are recorded with bipolar electrodes located in different muscles. However, such configuration may disregard some aspects of the spatial distribution of the potentials like location of innervation zones and the manifestation of inhomogineties in the control of the muscular fibers. On the other hand, the spatial distribution of motor unit action potentials has recently been assessed with activation maps obtained from High Density EMG signals (HD-EMG), these lasts recorded with arrays of closely spaced electrodes. The main objective of this work is to analyze patterns in the activation maps, associating them with four movement directions at the elbow joint and with different strengths of those tasks. Although the activation pattern can be assessed with bipolar electrodes, HD-EMG maps could enable the extraction of features that depend on the spatial distribution of the potentials and on the load-sharing between muscles, in order to have a better differentiation between tasks and effort levels.MethodsAn experimental protocol consisting of isometric contractions at three levels of effort during flexion, extension, supination and pronation at the elbow joint was designed and HD-EMG signals were recorded with 2D electrode arrays on different upper-limb muscles. Techniques for the identification and interpolation of artifacts are explained, as well as a method for the segmentation of the activation areas. In addition, variables related to the intensity and spatial distribution of the maps were obtained, as well as variables associated to signal power of traditional single bipolar recordings. Finally, statistical tests were applied in order to assess differences between information extracted from single bipolar signals or from HD-EMG maps and to analyze differences due to type of task and effort level.ResultsSignificant differences were observed between EMG signal power obtained from single bipolar configuration and HD-EMG and better results regarding the identification of tasks and effort levels were obtained with the latter. Additionally, average maps for a population of 12 subjects were obtained and differences in the co-activation pattern of muscles were found not only from variables related to the intensity of the maps but also to their spatial distribution.ConclusionsIntensity and spatial distribution of HD-EMG maps could be useful in applications where the identification of movement intention and its strength is needed, for example in robotic-aided therapies or for devices like powered- prostheses or orthoses. Finally, additional data transformations or other features are necessary in order to improve the performance of tasks identification.

Ocular Reduction in EEG Signals Based on Adaptive Filtering, Regression and Blind Source Separation

Quantitative electroencephalographic (EEG) analysis is very useful for diagnosing dysfunctional neural states and for evaluating drug effects on the brain, among others. However, the bidirectional contamination between electrooculographic (EOG) and cerebral activities can mislead and induce wrong conclusions from EEG recordings. Different methods for ocular reduction have been developed but only few studies have shown an objective evaluation of their performance. For this purpose, the following approaches were evaluated with simulated data: regression analysis, adaptive filtering, and blind source separation (BSS). In the first two, filtered versions were also taken into account by filtering EOG references in order to reduce the cancellation of cerebral high frequency components in EEG data. Performance of these methods was quantitatively evaluated by level of similarity, agreement and errors in spectral variables both between sources and corrected EEG recordings. Topographic distributions showed that errors were located at anterior sites and especially in frontopolar and lateral–frontal regions. In addition, these errors were higher in theta and especially delta band. In general, filtered versions of time-domain regression and of adaptive filtering with RLS algorithm provided a very effective ocular reduction. However, BSS based on second order statistics showed the highest similarity indexes and the lowest errors in spectral variables.

Drug effect on EEG connectivity assessed by linear and nonlinear couplings

Quantitative analysis of human electroencephalogram (EEG) is a valuable method for evaluating psychopharmacological agents. Although the effects of different drug classes on EEG spectra are already known, interactions between brain locations remain unclear. In this work, cross mutual information function and appropriate surrogate data were applied to assess linear and nonlinear couplings between EEG signals. The main goal was to evaluate the pharmacological effects of alprazolam on brain connectivity during wakefulness in healthy volunteers using a cross‐over, placebo‐controlled design. Eighty‐five pairs of EEG leads were selected for the analysis, and connectivity was evaluated inside anterior, central, and posterior zones of the scalp. Connectivity between these zones and interhemispheric connectivity were also measured. Results showed that alprazolam induced significant changes in EEG connectivity in terms of information transfer in comparison with placebo. Trends were opposite depending on the statistical characteristics: decreases in linear connectivity and increases in nonlinear couplings. These effects were generally spread over the entire scalp. Linear changes were negatively correlated, and nonlinear changes were positively correlated with drug plasma concentrations; the latter showed higher correlation coefficients. The use of both linear and nonlinear approaches revealed the importance of assessing changes in EEG connectivity as this can provide interesting information about psychopharmacological effects. Hum Brain Mapp, 2010.

Evaluation of Respiratory Muscles Activity by Means of Cross Mutual Information Function at Different Levels of Ventilatory Effort

Analysis of respiratory muscles activity is an effective technique for the study of pulmonary diseases such as obstructive sleep apnea syndrome (OSAS). Respiratory diseases, especially those associated with changes in the mechanical properties of the respiratory apparatus, are often associated with disruptions of the normally highly coordinated contractions of respiratory muscles. Due to the complexity of the respiratory control, the assessment of OSAS related dysfunctions by linear methods are not sufficient. Therefore, the objective of this study was the detection of diagnostically relevant nonlinear complex respiratory mechanisms. Two aims of this work were: 1) to assess coordination of respiratory muscles contractions through evaluation of interactions between respiratory signals and myographic signals through nonlinear analysis by means of cross mutual information function (CMIF); 2) to differentiate between functioning of respiratory muscles in patients with OSAS and in normal subjects. Electromyographic (EMG) and mechanomyographic (MMG) signals were recorded from three respiratory muscles: genioglossus, sternomastoid and diaphragm. Inspiratory pressure and flow were also acquired. All signals were measured in eight patients with OSAS and eight healthy subjects during an increased respiratory effort while awake. Several variables were defined and calculated from CMIF in order to describe correlation between signals. The results indicate different nonlinear couplings of respiratory muscles in both populations. This effect is progressively more evident at higher levels of respiratory effort.

Study of myographic signals from sternomastoid muscle in patients with chronic obstructive pulmonary disease

Analysis of the respiratory muscle activity is a promising technique for diagnosis of respiratory diseases, such as chronic obstructive pulmonary disease (COPD). The sternomastoid muscle (SMM) was selected to study the activity of respiratory muscles due to its accessibility in order to define a noninvasive analysis. The aims of this work are two: analyze the relationship between the SMM function and pulmonary obstruction, and study the influence of spectral estimator on frequency parameters related with the muscle activity. For the first goal, we propose the analysis of vibromyographic and electromyographic signals from the SMM to study the muscle function during two ventilatory tests. Activity of SMM was found by means of several indexes: root-mean-square (rms) values, mean and median frequencies, and ratio between high and low-frequency components. For the second goal, spectral analysis was performed by means of nonparametric methods: Correlogram and Welch periodogram, and parametric methods: autoregressive (AR), moving average (MA), and ARMA models. It is deduced that these indexes show muscle activity and certain fatigue of the SMM, whose muscle function depends on the level of pulmonary obstruction, and they depend a lot on spectral estimator being the more suitable an AR model with high order.

Stress assessment based on EEG univariate features and functional connectivity measures

The biological response to stress originates in the brain but involves different biochemical and physiological effects. Many common clinical methods to assess stress are based on the presence of specific hormones and on features extracted from different signals, including electrocardiogram, blood pressure, skin temperature, or galvanic skin response. The aim of this paper was to assess stress using EEG-based variables obtained from univariate analysis and functional connectivity evaluation. Two different stressors, the Stroop test and sleep deprivation, were applied to 30 volunteers to find common EEG patterns related to stress effects. Results showed a decrease of the high alpha power (11 to 12 Hz), an increase in the high beta band (23 to 36 Hz, considered a busy brain indicator), and a decrease in the approximate entropy. Moreover, connectivity showed that the high beta coherence and the interhemispheric nonlinear couplings, measured by the cross mutual information function, increased significantly for both stressors, suggesting that useful stress indexes may be obtained from EEG-based features.

Analyzing dynamic EMG and VMG signals of respiratory muscles

A promising technique is described for evaluating ventilatory disease by studying activity and fatigue in the sternocleidomastoid muscle. We analyze dynamic muscular function in time and frequency domains during two respiratory load tests at different levels of ventilation.

A Real-Time Method for Decoding the Neural Drive to Muscles Using Single-Channel Intra-Muscular EMG Recordings

The neural command from motor neurons to muscles - sometimes referred to as the neural drive to muscle - can be identified by decomposition of electromyographic (EMG) signals. This approach can be used for inferring the voluntary commands in neural interfaces in patients with limb amputations. This paper proposes for the first time an innovative method for fully automatic and real-time intramuscular EMG (iEMG) decomposition. The method is based on online single-pass density-based clustering and adaptive classification of bivariate features, using the concept of potential measure. No attempt was made to resolve superimposed motor unit action potentials. The proposed algorithm was validated on sets of simulated and experimental iEMG signals. Signals were recorded from the biceps femoris long-head, vastus medialis and lateralis and tibialis anterior muscles during low-to-moderate isometric constant-force and linearly-varying force contractions. The average number of missed, duplicated and erroneous clusters for the examined signals was [Formula: see text], [Formula: see text], and [Formula: see text], respectively. The average decomposition accuracy (defined similar to signal detection theory but without using True Negatives in the denominator) and coefficient of determination (variance accounted for) for the cumulative discharge rate estimation were [Formula: see text], and [Formula: see text], respectively. The time cost for processing each 200[Formula: see text]ms iEMG interval was [Formula: see text] (21-97)[Formula: see text]ms. However, computational time generally increases over time as a function of frames/signal epochs. Meanwhile, the incremental accuracy defined as the accuracy of real-time analysis of each signal epoch, was [Formula: see text]% for epochs recorded after initial one second. The proposed algorithm is thus a promising new tool for neural decoding in the next-generation of prosthetic control.

A Hybrid Computer-aided-diagnosis System for Prediction of Breast Cancer Recurrence (HPBCR) Using Optimized Ensemble Learning.

Cancer is a collection of diseases that involves growing abnormal cells with the potential to invade or spread to the body. Breast cancer is the second leading cause of cancer death among women. A method for 5-year breast cancer recurrence prediction is presented in this manuscript. Clinicopathologic characteristics of 579 breast cancer patients (recurrence prevalence of 19.3%) were analyzed and discriminative features were selected using statistical feature selection methods. They were further refined by Particle Swarm Optimization (PSO) as the inputs of the classification system with ensemble learning (Bagged Decision Tree: BDT). The proper combination of selected categorical features and also the weight (importance) of the selected interval-measurement-scale features were identified by the PSO algorithm. The performance of HPBCR (hybrid predictor of breast cancer recurrence) was assessed using the holdout and 4-fold cross-validation. Three other classifiers namely as supported vector machines, DT, and multilayer perceptron neural network were used for comparison. The selected features were diagnosis age, tumor size, lymph node involvement ratio, number of involved axillary lymph nodes, progesterone receptor expression, having hormone therapy and type of surgery. The minimum sensitivity, specificity, precision and accuracy of HPBCR were 77%, 93%, 95% and 85%, respectively in the entire cross-validation folds and the hold-out test fold. HPBCR outperformed the other tested classifiers. It showed excellent agreement with the gold standard (i.e. the oncologist opinion after blood tumor marker and imaging tests, and tissue biopsy). This algorithm is thus a promising online tool for the prediction of breast cancer recurrence.