Modeling the electroencephalogram by means of spatial spline smoothing and temporal autoregression

A Dynamic Window Recognition Algorithm for SSVEP-Based Brain–Computer Interfaces Using a Spatio-Temporal Equalizer

The past decade has witnessed rapid development in the field of brain–computer interfaces (BCIs). While the performance is no longer the biggest bottleneck in the BCI application, the tedious training process and the poor ease-of-use have become the most significant challenges. In this study, a spatio-temporal equalization dynamic window (STE-DW) recognition algorithm is proposed for steady-state visual evoked potential (SSVEP)-based BCIs. The algorithm can adaptively control the stimulus time while maintaining the recognition accuracy, which significantly improves the information transfer rate (ITR) and enhances the adaptability of the system to different subjects. Specifically, a spatio-temporal equalization algorithm is used to reduce the adverse effects of spatial and temporal correlation of background noise. Based on the theory of multiple hypotheses testing, a stimulus termination criterion is used to adaptively control the dynamic window. The offline analysis which used a benchmark dataset and an offline dataset collected from 16 subjects demonstrated that the STE-DW algorithm is superior to the filter bank canonical correlation analysis (FBCCA), canonical variates with autoregressive spectral analysis (CVARS), canonical correlation analysis (CCA) and CCA reducing variation (CCA-RV) algorithms in terms of accuracy and ITR. The results show that in the benchmark dataset, the STE-DW algorithm achieved an average ITR of 134 bits/min, which exceeds the FBCCA, CVARS, CCA and CCA-RV. In off-line experiments, the STE-DW algorithm also achieved an average ITR of 116 bits/min. In addition, the online experiment also showed that the STE-DW algorithm can effectively expand the number of applicable users of the SSVEP-based BCI system. We suggest that the STE-DW algorithm can be used as a reliable identification algorithm for training-free SSVEP-based BCIs, because of the good balance between ease of use, recognition accuracy, ITR and user applicability.

Protocol Design Challenges in the Detection of Awareness in Aware Subjects Using EEG Signals

Recent studies have evidenced serious difficulties in detecting covert awareness with electroencephalography-based techniques both in unresponsive patients and in healthy control subjects. This work reproduces the protocol design in two recent mental imagery studies with a larger group comprising 20 healthy volunteers. The main goal is assessing if modifications in the signal extraction techniques, training-testing/cross-validation routines, and hypotheses evoked in the statistical analysis, can provide solutions to the serious difficulties documented in the literature. The lack of robustness in the results advises for further search of alternative protocols more suitable for machine learning classification and of better performing signal treatment techniques. Specific recommendations are made using the findings in this work.

Assessing direct paths of intracortical causal information flow of oscillatory activity with the isolated effective coherence (iCoh)

Functional connectivity is of central importance in understanding brain function. For this purpose, multiple time series of electric cortical activity can be used for assessing the properties of a network: the strength, directionality, and spectral characteristics (i.e., which oscillations are preferentially transmitted) of the connections. The partial directed coherence (PDC) of Baccala and Sameshima (2001) is a widely used method for this problem. The three aims of this study are: (1) To show that the PDC can misrepresent the frequency response under plausible realistic conditions, thus defeating the main purpose for which the measure was developed; (2) To provide a solution to this problem, namely the "isolated effective coherence" (iCoh), which consists of estimating the partial coherence under a multivariate autoregressive model, followed by setting all irrelevant associations to zero, other than the particular directional association of interest; and (3) To show that adequate iCoh estimators can be obtained from non-invasively computed cortical signals based on exact low resolution electromagnetic tomography (eLORETA) applied to scalp EEG recordings. To illustrate the severity of the problem with the PDC, and the solution achieved by the iCoh, three examples are given, based on: (1) Simulated time series with known dynamics; (2) Simulated cortical sources with known dynamics, used for generating EEG recordings, which are then used for estimating (with eLORETA) the source signals for the final connectivity assessment; and (3) EEG recordings in rats. Lastly, real human recordings are analyzed, where the iCoh between six cortical regions of interest are calculated and compared under eyes open and closed conditions, using 61-channel EEG recordings from 109 subjects. During eyes closed, the posterior cingulate sends alpha activity to all other regions. During eyes open, the anterior cingulate sends theta-alpha activity to other frontal regions.

From EEG dependency multichannel matching pursuit to sparse topographic EEG decomposition

In this work, we present a multichannel EEG decomposition model based on an adaptive topographic time-frequency approximation technique. It is an extension of the Matching Pursuit algorithm and called dependency multichannel matching pursuit (DMMP). It takes the physiologically explainable and statistically observable topographic dependencies between the channels into account, namely the spatial smoothness of neighboring electrodes that is implied by the electric leadfield. DMMP decomposes a multichannel signal as a weighted sum of atoms from a given dictionary where the single channels are represented from exactly the same subset of a complete dictionary. The decomposition is illustrated on topographical EEG data during different physiological conditions using a complete Gabor dictionary. Further the extension of the single-channel time-frequency distribution to a multichannel time-frequency distribution is given. This can be used for the visualization of the decomposition structure of multichannel EEG. A clustering procedure applied to the topographies, the vectors of the corresponding contribution of an atom to the signal in each channel produced by DMMP, leads to an extremely sparse topographic decomposition of the EEG.

Quantitative multichannel EEG measure predicting the optimal weaning from ventilator in ICU patients with acute respiratory failure

OBJECTIVE

The objective of this study was to develop a novel quantitative multichannel EEG (qEEG) based analysis method, called Global Field Damping Time (GFDT), in order to detect potential EEG changes of patients admitted to the ICU with acute respiratory failure, and correlate them to the patients' recovery outcome predicting the optimal time-point to disconnect the patient from mechanical ventilation.

METHODS

Twenty-nine adult patients with acute respiratory failure out of 98 admitted to the Intensive Care Unit of Saint Paul General Hospital were enrolled, and among them only 15 completed the study. The patients were classified in 3 groups according to their outcome after 3 months follow-up. The patients were intubated with fraction of inspired oxygen (FiO2) of 100%. Neurological Deficit Scores (NDS) were measured 24 h after intubation to assess patients' neurological condition. Twenty-four hours after patient's intubation, FiO2 was decreased to 40% (weaning session), followed by a 5 min early recovery session, a 5 min recovery 1 session and a 5 min recovery 2 session. EEG recordings were performed during this experimental procedure. Multichannel EEG segments were processed and fitted into a multivariate autoregressive (mAR) model, and single channel EEG segments into a scalar autoregressive (sAR) model. The mAR and the sAR models of arbitrary order p were decomposed into mp and p oscillators and relaxators, respectively. Damping time of each oscillator and each relaxator, and the Global Field Damping Time (GFDT) as a weighted damping time were estimated for both mAR and sAR models.

RESULTS

A statistically significant increase of mAR model's GFDT during the weaning session was observed in the subjects of all groups. Comparing the 3 patients' groups, statistically significant differences for mAR model's GFDT were observed for the weaning and early recovery session. Linear regression analysis between NDS and mean mAR model's GFDT showed statistical significance during weaning session, early recovery session, and recovery 1 session. There was no statistical significance for SaO2 in the regression analysis with NDS. The sAR model's GFDT presented worst results in comparison with the mAR modelling GFDT in the identification of hypoxic conditions during weaning session and in the discrimination of patients with acute respiratory failure according to their neurological outcome.

CONCLUSIONS

Global Field Damping Time as correlated to the patients' neurological outcome appears to be a simple, compact, and substantial novel indicator of cerebral hypoxia and a potential predictor of the optimal time-point to disconnect the patient from the ventilator.

SIGNIFICANCE

Quantitative EEG seems to be an important tool for ICU clinicians assisting them to decide for the patients' optimal time-point to disconnect the patient from the ventilator.

Random field-union intersection tests for EEG/MEG imaging

Electrophysiological (EEG/MEG) imaging challenges statistics by providing two views of the same spatiotemporal data: topographic and tomographic. Until now, statistical tests for these two situations have developed separately. This work introduces statistical tests for assessing simultaneously the significance of spatiotemporal event-related potential/event-related field (ERP/ERF) components and that of their sources. The test for detecting a component at a given time instant is provided by a Hotelling's T(2) statistic. This statistic is constructed in such a manner to be invariant to any choice of reference and is based upon a generalized version of the average reference transform of the data. As a consequence, the proposed test is a generalization of the well-known Global Field Power statistic. Consideration of tests at all time instants leads to a multiple comparison problem addressed by the use of Random Field Theory (RFT). The Union-Intersection (UI) principle is the basis for testing hypotheses about the topographic and tomographic distributions of such ERP/ERF components. The performance of the method is illustrated with actual EEG recordings obtained from a visual experiment of pattern reversal stimuli.

Topographic time-frequency decomposition of the EEG

Frequency-transformed EEG resting data has been widely used to describe normal and abnormal brain functional states as function of the spectral power in different frequency bands. This has yielded a series of clinically relevant findings. However, by transforming the EEG into the frequency domain, the initially excellent time resolution of time-domain EEG is lost. The topographic time-frequency decomposition is a novel computerized EEG analysis method that combines previously available techniques from time-domain spatial EEG analysis and time-frequency decomposition of single-channel time series. It yields a new, physiologically and statistically plausible topographic time-frequency representation of human multichannel EEG. The original EEG is accounted by the coefficients of a large set of user defined EEG like time-series, which are optimized for maximal spatial smoothness and minimal norm. These coefficients are then reduced to a small number of model scalp field configurations, which vary in intensity as a function of time and frequency. The result is thus a small number of EEG field configurations, each with a corresponding time-frequency (Wigner) plot. The method has several advantages: It does not assume that the data is composed of orthogonal elements, it does not assume stationarity, it produces topographical maps and it allows to include user-defined, specific EEG elements, such as spike and wave patterns. After a formal introduction of the method, several examples are given, which include artificial data and multichannel EEG during different physiological and pathological conditions.