Characterization of early partial seizure onset: frequency, complexity and entropy

Characterization of antiseizure medications effects on the EEG neurodynamic by fractal dimension

Objectives

An important challenge in epilepsy is to define biomarkers of response to treatment. Many electroencephalography (EEG) methods and indices have been developed mainly using linear methods, e.g., spectral power and individual alpha frequency peak (IAF). However, brain activity is complex and non-linear, hence there is a need to explore EEG neurodynamics using nonlinear approaches. Here, we use the Fractal Dimension (FD), a measure of whole brain signal complexity, to measure the response to anti-seizure therapy in patients with Focal Epilepsy (FE) and compare it with linear methods.

Materials

Twenty-five drug-responder (DR) patients with focal epilepsy were studied before (t1, named DR-t1) and after (t2, named DR-t2) the introduction of the anti-seizure medications (ASMs). DR-t1 and DR-t2 EEG results were compared against 40 age-matched healthy controls (HC).

Methods

EEG data were investigated from two different angles: frequency domain-spectral properties in δ, θ, α, β, and γ bands and the IAF peak, and time-domain-FD as a signature of the nonlinear complexity of the EEG signals. Those features were compared among the three groups.

Results

The δ power differed between DR patients pre and post-ASM and HC (DR-t1 vs. HC, p < 0.01 and DR-t2 vs. HC, p < 0.01). The θ power differed between DR-t1 and DR-t2 (p = 0.015) and between DR-t1 and HC (p = 0.01). The α power, similar to the δ, differed between DR patients pre and post-ASM and HC (DR-t1 vs. HC, p < 0.01 and DR-t2 vs. HC, p < 0.01). The IAF value was lower for DR-t1 than DR-t2 (p = 0.048) and HC (p = 0.042). The FD value was lower in DR-t1 than in DR-t2 (p = 0.015) and HC (p = 0.011). Finally, Bayes Factor analysis showed that FD was 195 times more likely to separate DR-t1 from DR-t2 than IAF and 231 times than θ.

Discussion

FD measured in baseline EEG signals is a non-linear brain measure of complexity more sensitive than EEG power or IAF in detecting a response to ASMs. This likely reflects the non-oscillatory nature of neural activity, which FD better describes.

Conclusion

Our work suggests that FD is a promising measure to monitor the response to ASMs in FE.

Clinical Sensitivity of Fractal Neurodynamics

Among the significant advances in the understanding of the organization of the neuronal networks that coordinate the body and brain, their complex nature is increasingly important, resulting from the interaction between the very large number of constituents strongly organized hierarchically and at the same time with "self-emerging." This awareness drives us to identify the measures that best quantify the "complexity" that accompanies the continuous evolutionary dynamics of the brain. In this chapter, after an introductory section (Sect. 15.1), we examine how the Higuchi fractal dimension is able to perceive physiological processes (15.2), neurological (15.3) and psychiatric (15.4) disorders, and neuromodulation effects (15.5), giving a mention of other methods of measuring neuronal electrical activity in addition to electroencephalography, such as magnetoencephalography and functional magnetic resonance. Conscious that further progress will support a deeper understanding of the temporal course of neuronal activity because of continuous interaction with the environment, we conclude confident that the fractal dimension has begun to uncover important features of the physiology of brain activity and its alterations.

Technical feasibility of multimodal imaging in neonatal hypoxic-ischemic encephalopathy from an ovine model to a human case series

Hypoxic-ischemic encephalopathy (HIE) secondary to perinatal asphyxia occurs when the brain does not receive enough oxygen and blood. A surrogate marker for "intact survival" is necessary for the successful management of HIE. The severity of HIE can be classified based on clinical presentation, including the presence of seizures, using a clinical classification scale called Sarnat staging; however, Sarnat staging is subjective, and the score changes over time. Furthermore, seizures are difficult to detect clinically and are associated with a poor prognosis. Therefore, a tool for continuous monitoring on the cot side is necessary, for example, an electroencephalogram (EEG) that noninvasively measures the electrical activity of the brain from the scalp. Then, multimodal brain imaging, when combined with functional near-infrared spectroscopy (fNIRS), can capture the neurovascular coupling (NVC) status. In this study, we first tested the feasibility of a low-cost EEG-fNIRS imaging system to differentiate between normal, hypoxic, and ictal states in a perinatal ovine hypoxia model. Here, the objective was to evaluate a portable cot-side device and perform autoregressive with extra input (ARX) modeling to capture the perinatal ovine brain states during a simulated HIE injury. So, ARX parameters were tested with a linear classifier using a single differential channel EEG, with varying states of tissue oxygenation detected using fNIRS, to label simulated HIE states in the ovine model. Then, we showed the technical feasibility of the low-cost EEG-fNIRS device and ARX modeling with support vector machine classification for a human HIE case series with and without sepsis. The classifier trained with the ovine hypoxia data labeled ten severe HIE human cases (with and without sepsis) as the "hypoxia" group and the four moderate HIE human cases as the "control" group. Furthermore, we showed the feasibility of experimental modal analysis (EMA) based on the ARX model to investigate the NVC dynamics using EEG-fNIRS joint-imaging data that differentiated six severe HIE human cases without sepsis from four severe HIE human cases with sepsis. In conclusion, our study showed the technical feasibility of EEG-fNIRS imaging, ARX modeling of NVC for HIE classification, and EMA that may provide a biomarker of sepsis effects on the NVC in HIE.

A complex systems perspective on psychedelic brain action

Recent findings suggesting the potential transdiagnostic efficacy of psychedelic-assisted therapy have fostered the need to deepen our understanding of psychedelic brain action. Functional neuroimaging investigations have found that psychedelics reduce the functional segregation of large-scale brain networks. However, beyond this general trend, findings have been largely inconsistent. We argue here that a perspective based on complexity science that foregrounds the distributed, interactional, and dynamic nature of brain function may render these inconsistencies intelligible. We propose that psychedelics induce a mode of brain function that is more dynamically flexible, diverse, integrated, and tuned for information sharing, consistent with greater criticality. This 'meta' perspective has the potential to unify past findings and guide intuitions toward compelling mechanistic models.

Dynamic training of a novelty classifier algorithm for real-time detection of early seizure onset

OBJECTIVE

To develop an adaptive framework for seizure detection in real-time that is practical to use in the Epilepsy Monitoring Unit (EMU) as a warning signal, and whose output helps characterize epileptiform activity.

METHODS

Our algorithm was tested on intracranial EEG from epilepsy patients admitted to the EMU for presurgical evaluation. Our framework uses a one-class Support Vector Machine (SVM) that is being trained dynamically according to past activity in all available channels to classify the novelty of the current activity. In this study we compared multiple configurations using a one-class SVM to assess if there is significance over specific neural features or electrode locations.

RESULTS

Our results show that the algorithm reaches a sensitivity of 87% for early-onset seizure detection and of 97.7% as a generic seizure detection.

CONCLUSIONS

Our algorithm is capable of running in real-time and achieving a high performance for early seizure-onset detection with a low false positive rate and robustness in detection of different type of seizure-onset patterns.

SIGNIFICANCE

This algorithm offers a solution to warning systems in the EMU as well as a tool for seizure characterization during post-hoc analysis of intracranial EEG data for surgical resection of the epileptogenic network.

Performance Analysis of Entropy Methods in Detecting Epileptic Seizure from Surface Electroencephalograms

Physiological signals like Electrocardiography (ECG) and Electroencephalography (EEG) are complex and nonlinear in nature. To retrieve diagnostic information from these, we need the help of nonlinear methods of analysis. Entropy estimation is a very popular approach in the nonlinear category, where entropy estimates are used as features for signal classification and analysis. In this study, we analyze and compare the performances of four entropy methods; namely Distribution entropy (DistEn), Shannon entropy (ShanEn), Renyi entropy (RenEn) and LempelZiv complexity (LempelZiv) as classification features to detect epileptic seizure (ES) from surface Electroencephalography (sEEG) signal. Experiments were conducted on sEEG data from 23 subjects, obtained from the CHB-MIT database of PhysioNet. ShanEn, RenEn and LempelZiv entropy are found to be potential features for accurate and consistent detection of ES from sEEG, across multiple channels and subjects.

CEPS: An Open Access MATLAB Graphical User Interface (GUI) for the Analysis of Complexity and Entropy in Physiological Signals

BACKGROUND

We developed CEPS as an open access MATLAB^® GUI (graphical user interface) for the analysis of Complexity and Entropy in Physiological Signals (CEPS), and demonstrate its use with an example data set that shows the effects of paced breathing (PB) on variability of heart, pulse and respiration rates. CEPS is also sufficiently adaptable to be used for other time series physiological data such as EEG (electroencephalography), postural sway or temperature measurements.

METHODS

Data were collected from a convenience sample of nine healthy adults in a pilot for a larger study investigating the effects on vagal tone of breathing paced at various different rates, part of a development programme for a home training stress reduction system.

RESULTS

The current version of CEPS focuses on those complexity and entropy measures that appear most frequently in the literature, together with some recently introduced entropy measures which may have advantages over those that are more established. Ten methods of estimating data complexity are currently included, and some 28 entropy measures. The GUI also includes a section for data pre-processing and standard ancillary methods to enable parameter estimation of embedding dimension m and time delay τ ('tau') where required. The software is freely available under version 3 of the GNU Lesser General Public License (LGPLv3) for non-commercial users. CEPS can be downloaded from Bitbucket. In our illustration on PB, most complexity and entropy measures decreased significantly in response to breathing at 7 breaths per minute, differentiating more clearly than conventional linear, time- and frequency-domain measures between breathing states. In contrast, Higuchi fractal dimension increased during paced breathing.

CONCLUSIONS

We have developed CEPS software as a physiological data visualiser able to integrate state of the art techniques. The interface is designed for clinical research and has a structure designed for integrating new tools. The aim is to strengthen collaboration between clinicians and the biomedical community, as demonstrated here by using CEPS to analyse various physiological responses to paced breathing.

Localization of the epileptogenic zone based on ictal stereo-electroencephalogram: Brain network and single-channel signal feature analysis

Accurate localization of the epileptogenic zone (EZ) is crucial for refractory focal epilepsy patients to achieve freedom from seizures following epilepsy surgery. In this study, ictal stereo-electroencephalography data from 35 patients with refractory focal epilepsy were analyzed. Effective networks based on partial directed coherence were analyzed, and a gray level co-occurrence matrix was applied to extract the time-varying features of the in-degree. These features, combined with the single-channel signal time-frequency features, including approximate entropy and line length, were used to localize the EZ based on a cluster algorithm. For all seizure-free patients (n = 23), the proposed method was effective in identifying the clinical-EZ-contacts and clinical-EZ-blocks, with an F1-score of 62.47 % and 72.18 %, respectively. The sensitivity was 96.00 % for the clinical-EZ-block identification, which provided the information for the decision-making of clinicians, prompting clinicians to focus on the identified EZ-blocks and their nearby contacts. The agreement between the EZ identified by the proposed method and the clinical-EZ was worse for non-seizure-free patients (n = 12) than for seizure-free patients. Furthermore, our method provided better results than using only brain network or single-channel signal features. This suggests that combining these complementary features can facilitate more accurate localization of the EZ.

The complexity of clinically-normal sinus-rhythm ECGs is decreased in equine athletes with a diagnosis of paroxysmal atrial fibrillation

Equine athletes have a pattern of exercise which is analogous to human athletes and the cardiovascular risks in both species are similar. Both species have a propensity for atrial fibrillation (AF), which is challenging to detect by ECG analysis when in paroxysmal form. We hypothesised that the proarrhythmic background present between fibrillation episodes in paroxysmal AF (PAF) might be detectable by complexity analysis of apparently normal sinus-rhythm ECGs. In this retrospective study ECG recordings were obtained during routine clinical work from 82 healthy horses and from 10 horses with a diagnosis of PAF. Artefact-free 60-second strips of normal sinus-rhythm ECGs were converted to binary strings using threshold crossing, beat detection and a novel feature detection parsing algorithm. Complexity of the resulting binary strings was calculated using Lempel-Ziv (‘76 & ‘78) and Titchener complexity estimators. Dependence of Lempel-Ziv ‘76 and Titchener T-complexity on the heart rate in ECG strips obtained at low heart rates (25–60 bpm) and processed by the feature detection method was found to be significantly different in control animals and those diagnosed with PAF. This allows identification of horses with PAF from sinus-rhythm ECGs with high accuracy.

Automated classification of five seizure onset patterns from intracranial electroencephalogram signals

OBJECTIVE

The electroencephalographic (EEG) signals contain information about seizures and their onset location. There are several seizure onset patterns reported in the literature, and these patterns have clinical significance. In this work, we propose a system to automatically classify five seizure onset patterns from intracerebral EEG signals.

METHODS

The EEG was segmented by clinicians indicating the start and end time of each seizure onset pattern, the channels involved at onset and the seizure onset pattern. Twelve features that represent the time domain characteristics and signal complexity were extracted from 663 seizures channels of 24 patients. The features were used for classification of the patterns with support vector machine - Error-Correcting Output Codes (SVM-ECOC). Three patient groups with a similar number of seizure segments were created, and one group was used for testing and the rest for training. This test was repeated by rotating the testing and training data.

RESULTS

The feature space formed by both time domain and multiscale sample entropy features perform well in classification of the data. An overall accuracy of 80.7% was obtained with these features and a linear kernel of SVM-ECOC.

CONCLUSIONS

The seizure onset patterns consist of varied time and complexity characteristics. It is possible to automatically classify various seizure onset patterns very similarly to visual classification.

SIGNIFICANCE

The proposed system could aid the medical team in assessing intracerebral EEG by providing an objective classification of seizure onset patterns.

Resting-state brain entropy in right temporal lobe epilepsy and its relationship with alertness

BACKGROUND

To date, no functional MRI (fMRI) studies have focused on brain entropy in right temporal lobe epilepsy (rTLE) patients. Here, we characterized brain entropy (BEN) alterations in patients with rTLE using resting-state functional MRI(rs-fMRI) and explored the relationship between BEN and alertness.

METHOD

Thirty-one rTLE patients and 33 controls underwent MRI scanning to investigate differences in BEN and resting-state functional connectivity (rs-FC) in regions of interest (ROIs) between patients and controls. Correlation analyses were performed to examine relationships between the BEN of each ROI and alertness reaction times (RTs) in rTLE patients.

RESULTS

Compared with controls, the BEN of rTLE patients was significantly increased in the right middle temporal gyrus, inferior temporal gyrus, and other regions of the left hemisphere and significantly decreased in the right middle frontal gyrus and left supplementary motor area (p < .05). The rs-FCs between the ROIs (at p < .01, with the left superior parietal lobule and right precentral gyrus defined as ROI1 and ROI2, respectively) and the whole brain showed an increasing trend in rTLE patients. In addition, the BEN of ROI2 was associated with the intrinsic alertness and phasic alertness RTs of patients with rTLE.

CONCLUSIONS

Our findings suggest that BEN is altered in patients with rTLE and that decreased BEN in the right precentral gyrus is positively related to intrinsic and phasic alertness; the abnormal FC in the brain regions with altered entropy suggests a reconstruction of brain functional connectivity. These findings suggest that BEN mapping may provide a useful tool for probing brain mechanisms related to TLE.

Ultra Broad Band Neural Activity Portends Seizure Onset in a Rat Model of Epilepsy

Epilepsy affects over 70 million people worldwide and 30% of patients' seizures cannot be controlled with medications, motivating the development of alternative therapies such as electrical stimulation. Current stimulation strategies attempt to stop seizures after they start, but none aim to prevent seizures altogether. Preventing seizures requires knowing when the brain is entering a preictal state (i.e., approaching seizure onset). Here we show that such preictal activity can be detected by an informative neural signal that progressively and monotonically changes as the brain approaches a seizure event. Specifically, we use local field potentials (LFP) from a rat model of epilepsy to develop an innovative measure of signal novelty relative to nonseizure activity, that shows the presence of progressive neural dynamics in an ultra broad band (4 Hz - 5 kHz). The measure is extracted from functional connectivity features computed from the LFPs which are used as an input to a one-class Support Vector Machine (SVM). The SVM outputs a scalar signal which quantifies how novel the current activity looks relative to baseline (non-seizure) activity and shows a progression towards seizure onset minutes ahead of time. The use of ultra broad band multivariate features into the SVM results in a novelty signal that has a significantly higher slope in the progression to seizure onset when compared to using power in conventional frequency bands (4 - 500 Hz) on individual channels as input features to the SVM. Functional connectivity in conjunction with the SVM is a strategy that generates a new measurement of novelty that can be used by closed-loop systems for seizure forecasting and prevention.

The application of Lempel-Ziv and Titchener complexity analysis for equine telemetric electrocardiographic recordings

The analysis of equine electrocardiographic (ECG) recordings is complicated by the absence of agreed abnormality classification criteria. We explore the applicability of several complexity analysis methods for characterization of non-linear aspects of electrocardiographic recordings. We here show that complexity estimates provided by Lempel-Ziv ’76, Titchener’s T-complexity and Lempel-Ziv ’78 analysis of ECG recordings of healthy Thoroughbred horses are highly dependent on the duration of analysed ECG fragments and the heart rate. The results provide a methodological basis and a feasible reference point for the complexity analysis of equine telemetric ECG recordings that might be applied to automate detection of equine arrhythmias in equine clinical practice.

Multivariate Matching Pursuit Decomposition and Normalized Gabor Entropy for Quantification of Preictal Trends in Epilepsy

Quantification of the complexity of signals recorded concurrently from multivariate systems, such as the brain, plays an important role in the study and characterization of their state and state transitions. Multivariate analysis of the electroencephalographic signals (EEG) over time is conceptually most promising in unveiling the global dynamics of dynamical brain disorders such as epilepsy. We employed a novel methodology to study the global complexity of the epileptic brain en route to seizures. The developed measures of complexity were based on Multivariate Matching Pursuit (MMP) decomposition of signals in terms of time–frequency Gabor functions (atoms) and Shannon entropy. The measures were first validated on simulation data (Lorenz system) and then applied to EEGs from preictal (before seizure onsets) periods, recorded by intracranial electrodes from eight patients with temporal lobe epilepsy and a total of 42 seizures, in search of global trends of complexity before seizures onset. Out of five Gabor measures of complexity we tested, we found that our newly defined measure, the normalized Gabor entropy (NGE), was able to detect statistically significant (p < 0.05) nonlinear trends of the mean global complexity across all patients over 1 h periods prior to seizures’ onset. These trends pointed to a slow decrease of the epileptic brain’s global complexity over time accompanied by an increase of the variance of complexity closer to seizure onsets. These results show that the global complexity of the epileptic brain decreases at least 1 h prior to seizures and imply that the employed methodology and measures could be useful in identifying different brain states, monitoring of seizure susceptibility over time, and potentially in seizure prediction.

Bio-Signal Complexity Analysis in Epileptic Seizure Monitoring: A Topic Review

Complexity science has provided new perspectives and opportunities for understanding a variety of complex natural or social phenomena, including brain dysfunctions like epilepsy. By delving into the complexity in electrophysiological signals and neuroimaging, new insights have emerged. These discoveries have revealed that complexity is a fundamental aspect of physiological processes. The inherent nonlinearity and non-stationarity of physiological processes limits the methods based on simpler underlying assumptions to point out the pathway to a more comprehensive understanding of their behavior and relation with certain diseases. The perspective of complexity may benefit both the research and clinical practice through providing novel data analytics tools devoted for the understanding of and the intervention about epilepsies. This review aims to provide a sketchy overview of the methods derived from different disciplines lucubrating to the complexity of bio-signals in the field of epilepsy monitoring. Although the complexity of bio-signals is still not fully understood, bundles of new insights have been already obtained. Despite the promising results about epileptic seizure detection and prediction through offline analysis, we are still lacking robust, tried-and-true real-time applications. Multidisciplinary collaborations and more high-quality data accessible to the whole community are needed for reproducible research and the development of such applications.

Application of Higuchi's fractal dimension from basic to clinical neurophysiology: A review

BACKGROUND AND OBJECTIVE

For more than 20 years, Higuchi's fractal dimension (HFD), as a nonlinear method, has occupied an important place in the analysis of biological signals. The use of HFD has evolved from EEG and single neuron activity analysis to the most recent application in automated assessments of different clinical conditions. Our objective is to provide an updated review of the HFD method applied in basic and clinical neurophysiological research.

METHODS

This article summarizes and critically reviews a broad literature and major findings concerning the applications of HFD for measuring the complexity of neuronal activity during different neurophysiological conditions. The source of information used in this review comes from the PubMed, Scopus, Google Scholar and IEEE Xplore Digital Library databases.

RESULTS

The review process substantiated the significance, advantages and shortcomings of HFD application within all key areas of basic and clinical neurophysiology. Therefore, the paper discusses HFD application alone, combined with other linear or nonlinear measures, or as a part of automated methods for analyzing neurophysiological signals.

CONCLUSIONS

The speed, accuracy and cost of applying the HFD method for research and medical diagnosis make it stand out from the widely used linear methods. However, only a combination of HFD with other nonlinear methods ensures reliable and accurate analysis of a wide range of neurophysiological signals.

EEG analysis via multiscale Lempel-Ziv complexity for seizure detection

Robust seizure detection and seizure prediction continues to be a challenge. Lempel-Ziv Complexity (LZC) is one of the features that has shown to be relevant in seizure detection. Recent work has shown that augmenting LZC can be beneficial to emphasize variations in amplitude or frequency when analyzing biomedical signals. In this paper, we present a first look into evaluating the feasibility of using a recently proposed feature stemmed from LZC, namely the Multiscale Lempel-Ziv Complexity (MLZC) for seizure detection. MLZC does not allow the high-frequency signal components to be overwhelmed by the low frequency signal components when calculating complexity values. We compare MLZC and LZC for identifying seizures for three cases and show MLZC can provide a clear separation between non-ictal and ictal periods for all three cases using a single threshold over 7 recordings and 7 seizures per patient, whereas LZC provided such a clear separation for only one of the patients.

Multifractal Analysis and Relevance Vector Machine-Based Automatic Seizure Detection in Intracranial EEG

Automatic seizure detection technology is of great significance for long-term electroencephalogram (EEG) monitoring of epilepsy patients. The aim of this work is to develop a seizure detection system with high accuracy. The proposed system was mainly based on multifractal analysis, which describes the local singular behavior of fractal objects and characterizes the multifractal structure using a continuous spectrum. Compared with computing the single fractal dimension, multifractal analysis can provide a better description on the transient behavior of EEG fractal time series during the evolvement from interictal stage to seizures. Thus both interictal EEG and ictal EEG were analyzed by multifractal formalism and their differences in the multifractal features were used to distinguish the two class of EEG and detect seizures. In the proposed detection system, eight features (α0, αmin, αmax, Δα, f(α min ), f(α max ), Δf and R) were extracted from the multifractal spectrums of the preprocessed EEG to construct feature vectors. Subsequently, relevance vector machine (RVM) was applied for EEG patterns classification, and a series of post-processing operations were used to increase the accuracy and reduce false detections. Both epoch-based and event-based evaluation methods were performed to appraise the system's performance on the EEG recordings of 21 patients in the Freiburg database. The epoch-based sensitivity of 92.94% and specificity of 97.47% were achieved, and the proposed system obtained a sensitivity of 92.06% with a false detection rate of 0.34/h in event-based performance assessment.

Ordinal symbolic analysis and its application to biomedical recordings

Ordinal symbolic analysis opens an interesting and powerful perspective on time-series analysis. Here, we review this relatively new approach and highlight its relation to symbolic dynamics and representations. Our exposition reaches from the general ideas up to recent developments, with special emphasis on its applications to biomedical recordings. The latter will be illustrated with epilepsy data.

Multiscale Lempel-Ziv complexity for EEG measures

OBJECTIVE

To demonstrate that the classical calculation of Lempel-Ziv complexity (LZC) has an important limitation when applied to EEGs with rapid rhythms, and to propose a multiscale approach that overcomes this limitation.

METHODS

We have evaluated, both with simulated and real EEGs, whether LZC calculation neglects functional characteristics of rapid EEG rhythms. In addition, we have proposed a procedure to obtain multiple binarization sequences that yield a spectrum of LZC, and we have explored whether complexity would be better captured using this computation.

RESULTS

In our simulated signals, classical LZC did not capture modulations of a rapid component when a slower component of more amplitude was included in the signal. In real EEGs from healthy participants with eyes closed and eyes open, classical LZC calculation failed to show any difference between these two conditions. However, a multiscale LZC showed that complexity was lower for eyes closed than for eyes open conditions.

CONCLUSIONS

As hypothesized, our new approximation captures the complexity of series with fast components masked by slower rhythms.

SIGNIFICANCE

The method we introduce significantly improves LZC calculation, and it allows a better characterization of complexity of EEG signals.

Robust deep network with maximum correntropy criterion for seizure detection

Effective seizure detection from long-term EEG is highly important for seizure diagnosis. Existing methods usually design the feature and classifier individually, while little work has been done for the simultaneous optimization of the two parts. This work proposes a deep network to jointly learn a feature and a classifier so that they could help each other to make the whole system optimal. To deal with the challenge of the impulsive noises and outliers caused by EMG artifacts in EEG signals, we formulate a robust stacked autoencoder (R-SAE) as a part of the network to learn an effective feature. In R-SAE, the maximum correntropy criterion (MCC) is proposed to reduce the effect of noise/outliers. Unlike the mean square error (MSE), the output of the new kernel MCC increases more slowly than that of MSE when the input goes away from the center. Thus, the effect of those noises/outliers positioned far away from the center can be suppressed. The proposed method is evaluated on six patients of 33.6 hours of scalp EEG data. Our method achieves a sensitivity of 100% and a specificity of 99%, which is promising for clinical applications.

A boosted cascade for efficient epileptic seizure detection

Seizure detection from electroencephalogram (EEG) plays an important role for epilepsy therapy. Due to the diversity of seizure EEG patterns between different individuals, multiple features are necessary for high accuracy since a single feature could hardly encode all types of epileptiform discharges. However, a large feature set inevitably causes the increase of the computational cost. This paper proposes a boosted cascade chain to obtain both high detection performance and high computational efficiency. Sixteen features that are widely used in seizure detection are implemented. Considering the sequential characteristics of EEG signals, the features are extracted on each 1-second segment and its former three segments. Thus, a total of 64 features are used to construct a feature pool. Based on the feature pool, Real AdaBoost is used to select a group of effective features, on which weak classifiers are learned to assemble a strong classifier. The strong classifier is transformed to a cascade classifier by reordering the weak classifiers and learning a threshold for each weak classifier. The cascade classifier still has the similar classification strength to the original strong classifier. More importantly, it is able to reject easy non-seizure samples by the first a few weak classifiers in the cascade, thus high computational efficiency can be obtained. To evaluate our method, 90.6-hour EEG signals from four patients are tested. The experimental results show that our method can achieve an average accuracy of 95.31% and an average detection rate of 91.29% with the false positive rate of 4.68%. On average, only about 4 features are used. Compared with support vector machine (SVM), our method is much more efficient with the similar detection performance.

Seizure localization using three-dimensional surface projections of intracranial EEG power

Intracranial EEG (icEEG) provides a critical road map for epilepsy surgery but it has become increasingly difficult to interpret as technology has allowed the number of icEEG channels to grow. Borrowing methods from neuroimaging, we aimed to simplify data analysis and increase consistency between reviewers by using 3D surface projections of intracranial EEG poweR (3D-SPIER). We analyzed 139 seizures from 48 intractable epilepsy patients (28 temporal and 20 extratemporal) who had icEEG recordings, epilepsy surgery, and at least one year of post-surgical follow-up. We coregistered and plotted icEEG β frequency band signal power over time onto MRI-based surface renderings for each patient, to create color 3D-SPIER movies. Two independent reviewers interpreted the icEEG data using visual analysis vs. 3D-SPIER, blinded to any clinical information. Overall agreement rates between 3D-SPIER and icEEG visual analysis or surgery were about 90% for side of seizure onset, 80% for lobe, and just under 80% for sublobar localization. These agreement rates were improved when flexible thresholds or frequency ranges were allowed for 3D-SPIER, especially for sublobar localization. Interestingly, agreement was better for patients with good surgical outcome than for patients with poor outcome. Localization using 3D-SPIER was measurably faster and considered qualitatively easier to interpret than visual analysis. These findings suggest that 3D-SPIER could be an improved diagnostic method for presurgical seizure localization in patients with intractable epilepsy and may also be useful for mapping normal brain function.