Next Article in Journal
Air Quality and Cancer Prevalence Trends across the Sub-Saharan African Regions during 2005–2020
Previous Article in Journal
Students’ Burnout at University: The Role of Gender and Worker Status
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epilepsy Seizures Prediction Based on Nonlinear Features of EEG Signal and Gradient Boosting Decision Tree

School of Communications and Information Engineering, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(18), 11326; https://doi.org/10.3390/ijerph191811326
Submission received: 9 August 2022 / Revised: 5 September 2022 / Accepted: 6 September 2022 / Published: 9 September 2022

Abstract

:
Epilepsy is a common neurological disorder with sudden and recurrent seizures. Early prediction of seizures and effective intervention can significantly reduce the harm suffered by patients. In this paper, a method based on nonlinear features of EEG signal and gradient boosting decision tree (GBDT) is proposed for early prediction of epilepsy seizures. First, the EEG signals were divided into two categories: those that had seizures onset over a period of time (represented by InT) and those that did not. Second, the noise in the EEG was removed using complementary ensemble empirical mode decomposition (CEEMD) and wavelet threshold denoising. Third, the nonlinear features of the two categories of EEG were extracted, including approximate entropy, sample entropy, permutation entropy, spectral entropy and wavelet entropy. Fourth, a GBDT classifier with random forest as the initial result was designed to distinguish the two categories of EEG. Fifth, a two-step “k of n” method was used to reduce the number of false alarms. The proposed method was evaluated on 13 patients’ EEG data from the CHB-MIT Scalp EEG Database. Based on ten-fold cross validation, the average accuracy was 91.76% when the InT was taken at 30 min, and 38 out of 39 seizures were successfully predicted. When the InT was taken for 40 min, the average accuracy was 92.50% and all 42 seizures selected were successfully predicted. The results indicate the effectiveness of the proposed method for predicting epilepsy seizures.

1. Introduction

Epilepsy is one of the most common neurological disorders in humans, affecting approximately 60 million people worldwide who are unable to live normally [1]. Approximately 70–80% of epilepsy patients are well controlled with medication or surgery, but about 20–30% of patients are still not effectively treated [2]. During seizures, patients experience loss of consciousness, involuntary convulsions and mental abnormalities. If a patient has a seizure while exercising or driving, it is likely to cause further harm and place a heavy burden on the patient and their family [3]. Therefore, if seizures can be predicted and interfered with in a timely manner, it can reduce the harm caused by seizures and improve the quality of life for patients.
Electroencephalographic (EEG) signals directly record the electrical activity in the cerebral cortex and contain a wealth of information about brain activity, which are widely used in the fields of brain disease diagnosis and medical rehabilitation. Depending on the method of acquisition, there are two types of EEG: intracranial EEG (iEEG) and scalp EEG (sEEG). Although sEEG is susceptible to interference from external factors, its non-invasive and convenient acquisition properties make it one of the most important tools in epilepsy diagnosis and treatment. With the continuous development of signal processing and artificial intelligence (AI), an increasing number of researchers have been working on the early prediction of epilepsy seizures using EEG, with some success.
At present, methods for epilepsy seizures prediction using EEG can be divided into two categories: threshold-based prediction methods and classification-based methods [3]. The former involves identifying patterns of changes in certain features of the EEG prior to seizures and using this pattern to determine whether a seizure is likely to occur. Salvatierra et al. [4] calculated the spectral power ratio (SPR) of the EEG in the hour before a seizure and plotted the SPR curve to predict seizures using the pattern of change in SPR. The method was evaluated on the CHB-MIT Scalp EEG Database with an optimal prediction time of 60 min and a prediction success rate of 100%. The threshold-based prediction method is intuitive and easy to understand, but the predictive effect of the same threshold varies greatly from patient to patient, making it difficult to establish stable models. In contrast, classification-based methods are more stable and effective, and are the main methods used for seizures prediction. This type of method extracts several features from the interictal and pre-ictal EEG signals, and uses machine learning or deep learning to identify the pre-ictal EEG to achieve seizures prediction. Zhang et al. [5] decomposed the patient’s EEG into eight frequency bands, extracted features for each band using a common spatial pattern (CSP), and then used a convolutional neural network (CNN) to classify the pre-ictal and interictal EEGs. This method obtained an average sensitivity of 0.92, an average accuracy of 0.90, and an average false alarm rate of 0.12/h. Usman et al. [6] combined handcrafted EEG features (mean, kurtosis, approximate entropy, etc.) with features automatically extracted using a three-layer CNN, and designed an ensemble classifier with SVM, CNN, and LSTM as base classifiers for EEG classification, which achieved a sensitivity of 96.28% and an average prediction time of 33 min. As patients with epilepsy spend most of their time in an interictal state, the amount of interictal data is much larger than the amount of pre-ictal data. Therefore, classification-based methods often suffer from class imbalance problems, which affect the performance of the model. Slicing and re-splicing EEG data or using generative adversarial networks (GAN) can increase the amount of pre-ictal data and, thus, mitigate the impact of the imbalanced class problem [6,7].
In this work, the epileptic EEG signal was divided into two categories; category A has seizures after a time interval, while category B has no seizures after a time interval. This interval is referred to as the interval time (InT) and can be either 30 or 40 min. The data in category B were taken from an interval of 4 consecutive hours without a seizure and had the same amount of data as the data in category A. This approach addressed the class imbalance and made full use of the patient’s EEG data. The two categories of EEG signals were first pre-processed by complementary ensemble empirical mode decomposition (CEEMD) and wavelet threshold denoising. Then, the approximate entropy, sample entropy, permutation entropy, spectral entropy and wavelet entropy were extracted as features and classified using a gradient boosting decision tree (GBDT) with random forest as the initial result.
The rest of this paper is organized as follows: Section 2 describes the dataset, pre-processing, feature extraction, classification and post-processing methods, Section 3 presents the results, and the discussion and conclusions are in Section 4 and Section 5 respectively.

2. Materials and Methods

2.1. Datasets

The EEG data used in this paper were obtained from the CHB-MIT Scalp EEG Database [7], which contains sEEG data from 22 epilepsy patients aged 3–22. Each patient containsbetween 9 and 42 EEG data files, recorded at a sampling rate of 256 Hz and 16 bit resolution, with electrodes named and placed according to the international 10–20 system. An annotation file in the dataset described the start and end time of each seizure. Hardware limitations resulted in gaps between consecutively-numbered EEG files, during which the signals were not recorded; in most cases, the gaps were 10 s or less.

2.2. Data Segmentation

For each patient, EEG data were divided into two categories, A and B. For a seizure, the period before the onset of the seizure was defined as the interval time (InT), and the EEG data from 10 min before the InT was included as part of the category A data. It should be noted that if the time interval between two seizures was short, then the category A data for the second seizure might contain the ictal or pre-ictal EEG of the first seizure, which would make the category A data ‘impure’. Therefore, seizures that were the subjects of this study were not permitted have any other seizures in the two hours prior, and seizures that did not meet this condition were excluded from this study. Additionally, if there was a gap in category A data for which no data were recorded, this A data was discarded and the corresponding seizure was not studied.
Suppose a patient has acquired n segments of category A data as described above, and n segments of category B data were obtained for that patient in order to keep the amount of data consistent between the two categories of EEG. Among all the patient’s EEGs, n sets of 4 consecutive hours without seizures were sought, and 10 min of the first 3 h of each set were randomly selected as a segment of category B data, which ensured that no seizures occurred in each segment of category B data after IT. No gap of no data recording was allowed in each category B data. A 10 min segment of data was then divided into 60 segments of 10 s in length, and the subsequent analysis was carried out in 10 s segments. The data segmentation process is shown in Figure 1.
In some previous studies, researchers have only used the EEG signal from the period before each seizure to analyze and process the signal, for example, defining the first 30 min before a seizure as the pre-ictal period, and the 30–90 min as the interictal period, and identifying the pre-ictal period to predict seizures. In this paper, category A data were obtained in the period before each seizure, while category B data were obtained in the period away from the seizures, which has rarely been used in previous studies, making full use of the patient’s EEG in the dataset. In addition, because the duration of the interictal is longer than the that of the pre-ictal, in some previous studies, the amount of data for the interictal period was higher than the pre-ictal period, which led to class imbalance problems. In this paper, the two categories of data, A and B, which need to be classified and identified, had the same amount of data, avoiding the class imbalance problem.
If the prediction model triggers an alarm, patients can intervene (e.g., take medication) for the upcoming seizure within the InT. Too short an InT may leave patients ill-prepared, while too long may cause stress as patients wait for the upcoming seizure. Therefore, the InT can be set for 30 or 40 min, which gives patients enough time to be well protected and at the same time reduces their anxiety and stress.
After screening, EEG data from 13 patients in the dataset were selected for the experiment. The selected EEG data of all patients contained 23 channels. A total of 39 seizures were selected for the experiment with an InT of 30 min, while 42 seizures were selected with an InT of 40 min. When different values of InT were taken, the amount of category A data for the same patient meeting the above criteria might be different, so that the number of selected seizures appears to be different. Details of the selected patients are shown in Table 1. In the case of patient chb01, with an InT of 30 min, four seizures were selected in the manner described above, so four segments of category A, as well as B EEG data, were divided, each with a recording time of 10 min.

2.3. Pre-Processing

EEG signals contains artifacts such as EOG and EMG signals, reducing the signal-to-noise ratio of EEG and, thus, adversely affecting the classification of the two categories of signals [8]. Therefore, a pre-processing method based on CEEMD and wavelet denoising was used for artefact removal.
CEEMD is an adaptive time-frequency analysis method for nonlinear and non-stationary signals that decomposes EEG signals into a number of intrinsic mode functions (IMFs) and a residual signal. The method adds a pair of white noises of opposite amplitude to the source signal as auxiliary noise when decomposing, in order to eliminate the excess auxiliary white noise remaining in the reconstructed signal and to reduce the number of iterations required for the decomposition. The wavelet transform of the EEG signal can be used to obtain a set of wavelet coefficients, where the larger coefficients correspond to the EEG signal and the smaller coefficients correspond to the noise.
The EEG signal was decomposed using CEEMD to obtain several IMFs, and the spectrum corresponding to each IMF was calculated, as shown in Figure 2, from which we can see that IMF1 and IMF2 contained frequency components greater than 30 Hz. Since the frequency of EEG usually lies in the range of 0.5–30 Hz, IMF1 and IMF2 were considered to contain noise. The “db4” wavelet and 4-layers decomposition were used for IMF1 and IMF2 for wavelet denoise with “sqtwolog” threshold method and soft threshold function. IMF1 and IMF2 after denoising were reconstructed with the remaining IMFs (IMF3-IMF10) to obtain the denoised EEG signals. The IMFs containing noise were not discarded directly but were used again for signal reconstruction after wavelet denoising, which avoided excessive loss of useful signal while removing noise. Figure 3 shows the EEG signal before and after pre-processing, with a signal duration of 10 s. The noise of high frequency in the EEG was suppressed without excessive loss of information.

2.4. Feature Extraction

The term “entropy” was originally derived from thermodynamics and was first proposed by Clausius in 1865. To date, entropy theory has been developed for more than 100 years, and several entropy features have become common in the field of EEG signal analysis. Researchers have applied various entropy parameters to the prediction of epilepsy seizures and achieved some results [9,10,11].

2.4.1. Entropy Based on Time Domain: Approximate Entropy, Sample Entropy, Permutation Entropy

Approximate entropy (ApEn) is used to describe the unpredictability or randomness of a finite-length signal and is an important nonlinear signal feature. The computation involves embedding a signal into phase space, and estimates the rate of increment in the number of phase space patterns within a predefined value, r , when the embedding dimension of phase space increases from m to m + 1 [12]. For a time series x ( i )   ( 1     i     N ) of finite length, N , to reconstitute the N     m + 1 vectors X m ( i ) , follow the form:
X m i = x i ,   x i + 1 ,     x i + m     1 ,   i = 1 , 2 ,     N   m + 1
where m is the embedding dimension. Let C i m r be the probability that any vector X m ( j ) is within r of X m ( i ) . It can be defined as:
C i m r = 1 N   m + 1 j = 0 N   m + 1 Θ d ij m r ,   i ,   j = 1 , 2 ,   N   m + 1
where d is the distance between vector X m ( j )   and X m ( i ) , defined as:
d ij m = d X i m , X j m = max x i + k x j + k ,   k = 0 ,   1 ,   ,   m  
where Θ is the Heaviside function. Define parameter Φ m ( r ) as:
Φ m ( r ) = ( N   m + 1 ) 1 i = 1 N   m + 1 lnC i m r
Change the embedding dimension to m + 1 , repeat the above steps to calculate   Φ m + 1 ( r ) :
Φ m + 1 ( r ) = ( N   m ) 1 i = 1 N   m lnC i m + 1 ( r )  
Then the approximate entropy of the signal is defined by:
ApEn ( m ,   r ,   N ) = Φ m ( r )   Φ m + 1 ( r )
Approximate entropy is susceptible to data length ( N ), similar tolerance ( r ) and the embedding dimension ( m ). In this paper, we chose m as 2, and r as 0.2 times the standard deviation of the signal.
Sample entropy (SampEn) has a similar meaning to approximate entropy in that it measures the complexity of a time series by measuring the probability of generating a new pattern in the signal. Similar to approximate entropy, the computation of sample entropy first requires calculating   C i m r . When the embedding dimension is m , define a parameter B m ( r ) as:
B m ( r ) = ( N   m ) 1 i = 1 N   m C i m ( r )  
Increase the embedding dimension to m + 1 , define A m ( r ) as:
A m ( r ) = ( N   m ) 1 i = 1 N   m C i m ( r )  
Finally, the sample entropy can be estimated by:
SampEn ( r ,   m ,   N ) = ln A m ( r ) B m ( r )  
Permutation entropy (PEn) is an entropy parameter that measures the complexity of a time series and is characterized by simplicity of calculation and resistance to noise interference. Given a time series x ( i )   ( 1     i     N ) . Firstly, the reconstruction time series:
X i = x i , x i + τ , , x i + m   1 τ ,   i = 1 ,   2 ,     N   ( m 1 ) τ  
where τ is time delay, and m is the embedding dimension.
The elements in each of the reconstructed components X i are rearranged in ascending order:
x i + j 1 1 τ x i + j 2 1 τ x i + j m 1 τ  
where j 1 ,   j 2 ,   j m denotes the index of the column in which each element of the reconstructed component is located. As a result, each reconstructed component of a time series yields a sequence of symbols j 1 ,   j 2 ,   j m . There are m! different sequences in the m -dimensional phase space; let the number of occurrences of the j th sequence be n j , then its probability of occurrence p j can be calculated by:
p j = n j j = 1 m ! n j
According to the formula of Shannon entropy, there are:
H x ( m ) = j = 1 m ! p j lnp j  
Finally, the permutation entropy can be estimated by:
PEn = H x ( m ) ln ( m ! )

2.4.2. Entropy Based on Frequency Domain: Spectral Entropy

Spectral entropy (SpEn) is the Shannon entropy of the power spectral density of the signal [13]. It is calculated as follows:
SpEn = j = 1 p f logp f  
where p f is the relative power of the component with frequency f .

2.4.3. Entropy Based on Time-Frequency Domain: Wavelet Entropy

Wavelet decomposition is performed on a segment of the signal with a scale of j and a window length of   N j , then the wavelet energy E j N j at each scale j is defined as:
E j N j = k d j ( k ) 2  
where d j ( k ) are the wavelet coefficients. The total energy of all scales is:
E tot = j E j N j  
The percentage of E j N j in E tot is:
p j N j = E j N j E tot  
Finally, the wavelet entropy of the signal is calculated as follows:
WEn = j p j N j logp j N j  
The five entropy parameters described above were calculated for each of the 23 channels for each 10 s segment of data for all patients. Assuming a patient had a total of n selected seizures, then they had 2 n segments of category A and B data. The feature extraction process generated a feature matrix of 2 n   ×   60 rows and 23   ×   ( 5 + 1 ) columns, where the last column was labeled as the target value of 1 for category A data, and 0 for category B. For example, the feature extraction for patient chb01 resulted in a feature matrix of size 480   ×   116 .

2.5. Classification

Ensemble learning is the process of constructing and combining multiple learners to accomplish specific learning tasks. The general structure is to generate a set of “base learners”, which are then combined with a certain strategy to produce the final result. Based on the way base learners are generated, there are two broad categories of ensemble learning methods: bagging, in which the base learners are built in parallel and do not affect each other, and boosting, in which the base learners are built one by one and the previous base learner has an impact on the creation of the next. Ensemble learning methods based on boosting (e.g., AdaBoost, GBDT, XGBoost, etc.) have been widely used in many fields.
Figure 4 illustrates the general modeling process of boosting. For each sample in the training set   ( x i , y i ) (where x i is the feature vector and   y i is the target), the base learner gives its evaluation result f x i and calculates a loss function L ( y i , f x i ) , then, adaptively influences the construction of the next base learner by a certain rule. Finally, the evaluation results of the ensemble learner are derived from all base learners.
Assuming that the evaluation result of sample x i under the t th base learner is f t x i and a total of T base learners are constructed, the output result H x i of the boosting ensemble learner is the weighted sum of all base learners, which can be expressed as:
H x i = t = 1 T ϕ t f t x i
where ϕ t denotes the weight of the t th base learner. It is worth noting that the final evaluation result is not calculated after all the base learners are constructed, but is computed over iterations as the base learners are constructed. That is, after the first t base learner is constructed, there are:
H t x i = H t 1 x i + ϕ t f t x i
GBDT is one of the most representative algorithms in boosting and one of the most consistently performing machine learning algorithms in practical application scenarios. GBDT uses CART decision trees as base learners and features negative gradients of the loss function of each base learner as residuals, and builds the base learners by continuously fitting the residuals. For the sample ( x i , y i ) , the t th base learner evaluates the result of f t x i , the loss function is L ( y i , f t x i ) , then the residual is the negative gradient of the loss function, i.e.,
r it = [ L ( y i ,   f t ( x i ) ) f t ( x i ) ] f t ( x i ) = f t 1 ( x i )  
Use the residuals to replace the labels of the samples to obtain new samples ( x i , r it ) , then use the new samples to construct the next base learner. This step is repeated until the number of base learners reaches the preset value.
L ( y i ,   f t ( x i ) ) = ( y i logp ( x i ) + ( 1 y i ) log ( 1 p ( x i ) ) )  
where p ( x i ) is the result of the evaluation based on the first t evaluation result of the base learner H t x i and the probability value calculated by the Sigmoid function, i.e.,
p x i = Sigmoid ( H t x i )  
According to Equation (21), the first base classifier is constructed with:
H 1 x i = H 0 x i + ϕ 1 f 1 x i  
where H 0 x i is called the “initial result” of GBDT and is often taken as zero. In order to improve the classification effect of GBDT, a random forest (RF) classifier was trained first, and each sample was given a prediction result and its corresponding probability by the RF classifier. This probability was taken as the initial result of GBDT, and then the “RF + GBDT” classification model was constructed.
The feature matrices calculated for all patients in the feature extraction session were fed separately into a GBDT classifier with random forest as the initial result to construct a patient-specific EEG signal classification model. Each 10 s segment of the EEG was a sample, and the output of the classifier resulted in a target value of 0 or 1. Ten-fold cross validation procedures were followed to evaluate the generalization performance of the proposed approach.

2.6. Post-Processing

A two-step ‘k of n’ method was used to post-process the classified data to produce the final seizure prediction results while reducing the false alarm rate. For every 10 min segment of EEG data (containing 60 segments of 10 s in length), the classifier output a sequence of 60 numbers with a predictive target of 0 or 1. First, if three or more of the six consecutive numbers in the sequence were 1, then those six numbers were recorded as 1, otherwise they were recorded as 0. Once this step was completed, a sequence of 10 numbers was obtained. Secondly, if 5 or more of the numbers in a sequence of 10 were 1, the 10 numbers was recorded as 1, otherwise 0 was recorded. At this point, a prediction value of 0 or 1 was generated for every 10 min of data. A prediction value of 1 predicted the segment of EEG as category A, where a seizure was considered to occur after a period of time and an alarm was given; conversely, it was predicted as category B, where no seizure was considered to occur after a period of time. Figure 5 illustrates the process of post-processing.

3. Results

3.1. Evaluation Metrics

Based on the output of the GBDT classifier, the confusion matrix was drawn and used as the basis for calculating accuracy, sensitivity, false positive rate (FPR) and F1-score as the evaluation metrics of the classifier. Equations (26)–(29) are the formulas for the four metrics.
Accuracy = TP + TN TP + TN + FP + FN
Sensitivity = TP TP + FN
FPR = FP FP + TN
F 1 - score = 2 TP 2 TP + FN + FP
The post-processing session gave the prediction results for each 10 min data segment. If a segment of category A was correctly predicted, this was considered a successful prediction of a seizure, and conversely, this seizure was considered missed; if a segment of category B was incorrectly predicted as category A, this was considered a false alarm. The number of successful and missed predictions and false alarms were calculated for each patient based on the post-processing results.

3.2. Result

3.2.1. Classification Result

The classification evaluation metrics at InT of 30 and 40 min are shown in Table 2 and Table 3, respectively. The average accuracy for both categories of EEG at InT of 30 min was 91.76%, the sensitivity was 91.87%, the FPR was 0.083 and the F1-score was 91.78%. The average accuracy at InT of 40 min was 92.50%, the sensitivity was 91.90%, the FPR was 0.069 and F1-score was 92.37%. Using support vector machines (SVM), random forest (RF) and GBDT with zero initial results as classifier separately, the results (shown in Table 4) indicated that the “RF+GBDT” classification method used in this paper outperformed the above three classifiers.

3.2.2. Prediction Result

The seizure prediction results are shown in Table 5 and Table 6. A total of 39 seizures were selected at an InT of 30 min, and 38 were predicted, with one missed seizure and one false alarm. A total of 42 seizures selected were successfully predicted at an InT of 40 min, with one false alarm. The results show that the proposed method can provide accurate early prediction of seizures with a low number of false alarms.

4. Discussion

4.1. Comparison with Other Approaches

Table 7 shows the comparison between the proposed seizure prediction method and some of the existing methods, all of which use data from the CHB-MIT Scalp EEG Database. In comparison, the method proposed in this paper outperformed some of the existing methods and was slightly inferior to some methods that use deep learning.

4.2. Practical Application

In this paper, an epilepsy prediction method based on the nonlinear features of EEG signal and gradient boosting decision tree was proposed. The experimental results illustrate that this method can effectively classify EEG signals and achieve early prediction of seizures with a low false alarm rate. This method can provide early warning of seizures and reduce the harm suffered by epilepsy patients, which has some practical application prospects.
Patients can have sufficient EEG signals collected in hospitals and other institutions, and their EEG are subjected to a series of analyses and processing by the relevant researchers, including data segmentation, pre-processing, feature extraction, classification and post-processing, to form a patient-specific seizure prediction model. Once an early alarm is issued by the model, patients, their families or healthcare professionals can intervene in a timely manner (with medication or electrical stimulation, etc.) to prevent the seizure or reduce its severity, thus, greatly reducing the physical and psychological harm suffered by the patient.
In order to achieve the practical application of prediction models, the algorithms used in the models must have the ability of strong real-time analysis, which requires the design of algorithms with low complexity, short computing time and high accuracy. Making more accurate predictions with fewer features, speeding up the computing time of classifiers, and promoting the practical application of this approach will be future research directions.

4.3. Limitation of the Current Study

In this paper, a series of analyses of EEG signals from epileptic patients in the CHB-MIT Scalp EEG Dataset, including pre-processing, feature extraction and classification, were carried out to predict epilepsy seizures, early. However, there are some limitations to our current study.
First, the seizure prediction approach proposed in this paper relied on the annotation of seizures in the dataset. As described in Section 2.2, the 10 min prior to the InT before a certain seizure, was divided into category A data. In other words, the acquisition of category A data relied on the annotations about seizures, which also makes this approach still somewhat far from practical application. If the annotations of seizures could be ignored, and a more comprehensive analysis of the patient’s EEG could be carried out to identify the pattern of changes in EEG prior to seizures in each patient, it would be possible to predict whether a seizure will occur for any segment of the EEG without annotation. The ultimate goal of this and related studies is to analyze the EEG of diagnosed patients in order to predict seizures early and help them improve their quality of life. To achieve this, it is necessary to address the limitation of relying on annotations.
Second, it is still unknown what impact the choice of InT will have on the performance of the model. In this paper, the InT was only set to 30 and 40 min and no further analysis of the impact of the InT was undertaken. As described in Section 2.2, the choice of InT in this paper was made in order to give patients sufficient time to take the intervention without increasing their anxiety. From a practical application point of view, our current study prefers to improve the model performance as much as possible with a fixed InT. However, it is undeniably meaningful to analyze the impact of InT. For example, will a further increase in InT invalidate the model? Is there a threshold value of InT at which the model performs best? Addressing these issues will further refine the prediction model.

5. Conclusions

In this paper, a method based on nonlinear features of EEG signals and GBDT was proposed for early prediction of epilepsy seizures. To avoid the class imbalance problem, the EEG signals were divided into two categories: with and without seizures after a period of time. Based on CEEMD and wavelet denoising, the two categories of EEG signals were pre-processed to remove noise. Then, approximate entropy, sample entropy, permutation entropy, spectral entropy and wavelet entropy features were calculated; and a GBDT classifier with random forest as the initial result was used to classify the two categories of EEG. Finally, a two-step “k of n” method was used to obtain the prediction results. The proposed method was evaluated on the CHB-MIT Scalp EEG Database, and the results showed that it could accurately classify the two categories of EEG signals with a low false alarm rate. In addition, our “RF + GBDT” classification method outperformed other traditional machine learning classifiers. The promising results suggest the effectiveness of the proposed approach for early prediction of seizures, which can help patients to intervene in advance, to reduce the harm they suffer. Meanwhile, this approach has some limitations in that it relies on annotations about seizures and the impact of InT on the model is not yet clear. Further study and optimization is warranted.

Author Contributions

Conceptualization, X.X.; methodology, M.L.; validation, M.L; writing—original draft preparation, M.L.; writing—review and editing, X.X. and T.X; supervision, X.X. and T.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Excellent Youth Foundation of Jiangsu Scientific Committee (grant number: BK20211538), the National Science Foundations of China (grant number: No. 61991431), the National Basic Research Program of China (grant number: 2018YFA0209101) and Postgraduate Research & Practice Innovation Program of Jiangsu Province (grant number: KYCX22_0966).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of The Affiliated Brain Hospital of Nanjing Medical University on 20 March 2020 (protocol code: 2020-KY052-01).

Informed Consent Statement

Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are openly available in CHB-MIT Scalp EEG Database at (https://physionet.org/content/chbmit/1.0.0/, accessed on 10 March 2022), reference number [7].

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Rahman, R.; Varnosfaderani, S.M.; Makke, O.; Sarhan, N.J.; Asano, E.; Luat, A.; Alhawari, M.; IEEE. Comprehensive Analysis of EEG Datasets for Epileptic Seizure Prediction. In Proceedings of the IEEE International Symposium on Circuits and Systems (IEEE ISCAS), Daegu, Korea, 22–28 May 2021. [Google Scholar]
  2. Maimaiti, B.; Meng, H.M.; Lv, Y.D.; Qiu, J.Q.; Zhu, Z.P.; Xie, Y.Y.; Li, Y.; Cheng, Y.; Zhao, W.X.; Liu, J.Y.; et al. An Overview of EEG-based Machine Learning Methods in Seizure Prediction and Opportunities for Neurologists in this Field. Neuroscience 2022, 481, 197–218. [Google Scholar] [CrossRef] [PubMed]
  3. Han, C.; Peng, F.; Chen, C.; Li, W.; Zhang, X.; Wang, X.; Zhou, W. Research progress of epileptic seizure predictions based on electroencephalogram signals. Sheng Wu Yi Xue Gong Cheng Xue Za Zhi = Chin. J. Biomed. Eng. 2021, 38, 1193–1202. [Google Scholar]
  4. Salvatierra, N.; Sakanishi, R.; Flores, C. Epileptic Seizure Prediction from Scalp EEG Using Ratios of Spectral Power. In Proceedings of the 2020 IEEE Engineering International Research Conference (EIRCON), Lima, Peru, 21–23 October 2020. [Google Scholar]
  5. Zhang, Y.; Guo, Y.; Yang, P.; Chen, W.; Lo, B. Epilepsy Seizure Prediction on EEG Using Common Spatial Pattern and Convolutional Neural Network. IEEE J. Biomed. Health Inform. 2020, 24, 465–474. [Google Scholar] [CrossRef] [PubMed]
  6. Usman, S.M.; Khalid, S.; Bashir, S. A deep learning based ensemble learning method for epileptic seizure prediction. Comput. Biol. Med. 2021, 136, 12. [Google Scholar]
  7. CHB-MIT Scalp EEG Database. Available online: https://physionet.org/content/chbmit/1.0.0/ (accessed on 10 March 2022).
  8. Suhail, T.A.; Indiradevi, K.P.; Suhara, E.M.; Poovathinal, S.A.; Anitha, A. Performance Analysis of Mother Wavelet Functions and Thresholding Methods for Denoising EEG Signals during Cognitive Tasks. In Proceedings of the 2020 International Conference on Power, Instrumentation, Control and Computing (PICC), Thrissur, India, 17–19 December 2020. [Google Scholar]
  9. Niu, Y.; Cao, R.; Wang, H.Y.; Li, C.G.; Zhou, M.N.; Guo, Y.X.; Wang, B.; Yan, P.F.; Xiang, J. Permutation Fuzzy Entropy-An Index for the Analysis of Epileptic Electroencephalogram. J. Med. Imaging Health Inform. 2019, 9, 637–645. [Google Scholar] [CrossRef]
  10. Sukriti; Chakraborty, M.; Mitra, D.; IEEE. Epilepsy Seizure Detection using Non-linear and DWT-based Features. In Proceedings of the 4th IEEE International Conference on Wireless Communications Signal Processing and Networking (WiSPNET), Chennai, India, 21–23 March 2019. [Google Scholar]
  11. Zhang, J.; Wei, Z.C.; Zou, J.Z.; Fu, H. Automatic epileptic EEG classification based on differential entropy and attention model. Eng. Appl. Artif. Intell. 2020, 96, 10. [Google Scholar] [CrossRef]
  12. Bai, Y.; Li, X. Nonlinear Neural Dynamics. In EEG Signal Processing and Feature Extraction; Hu, L., Zhang, Z., Eds.; Springer Nature: Singapore, 2019. [Google Scholar]
  13. Yan, J.; Li, J.; Xu, H.; Yu, Y.; Pan, L.; Cheng, X.; Tan, S. EEG Seizure Prediction Based on Empirical Mode Decomposition and Convolutional Neural Network. In Proceedings of the 14th International Conference on Brain Informatics (BI), Virtual Event, 17–19 September 2021. [Google Scholar]
  14. Alotaiby, T.N.; Alshebeili, S.A.; Alotaibi, F.M.; Alrshoud, S.R. Epileptic Seizure Prediction Using CSP and LDA for Scalp EEG Signals. Comput. Intell. Neurosci. 2017, 2017, 1240323. [Google Scholar] [CrossRef] [PubMed]
  15. Agboola, H.; Solebo, C.; Aribike, D.; Lesi, A.; Susu, A. Seizure Prediction with Adaptive Feature Representation Learning. J. Neurol. Neurosci. 2019, 10, 294. [Google Scholar] [CrossRef]
  16. Zhang, Q.Z.; Ding, J.; Kong, W.Z.; Liu, Y.; Wang, Q.; Jiang, T.J. Epilepsy prediction through optimized multidimensional sample entropy and Bi-LSTM. Biomed. Signal Process. Control 2021, 64, 9. [Google Scholar] [CrossRef]
  17. Rusnac, A.L.; Grigore, O.; IEEE. Intelligent Seizure Prediction System Based on Spectral Entropy. In Proceedings of the 14th International Symposium on Signals, Circuits and Systems (ISSCS), Iasi, Romania, 11–12 July 2019. [Google Scholar]
  18. Daoud, H.; Bayoumi, M.A. Efficient Epileptic Seizure Prediction Based on Deep Learning. IEEE Trans. Biomed. Circuits Syst. 2019, 13, 804–813. [Google Scholar] [CrossRef] [PubMed]
  19. Jana, R.; Mukherjee, I. Deep learning based efficient epileptic seizure prediction with EEG channel optimization. Biomed. Signal Process. Control 2021, 68, 8. [Google Scholar] [CrossRef]
  20. Yan, J.Z.; Li, J.N.; Xu, H.X.; Yu, Y.C.; Xu, T.Y. Seizure Prediction Based on Transformer Using Scalp Electroencephalogram. Appl. Sci. 2022, 12, 4158. [Google Scholar] [CrossRef]
Figure 1. Segmentation of EEG data: (a) Category A data, (b) Category B data.
Figure 1. Segmentation of EEG data: (a) Category A data, (b) Category B data.
Ijerph 19 11326 g001
Figure 2. The IMFs and their corresponding spectrum.
Figure 2. The IMFs and their corresponding spectrum.
Ijerph 19 11326 g002
Figure 3. Comparison of EEG signal before and after pre-processing.
Figure 3. Comparison of EEG signal before and after pre-processing.
Ijerph 19 11326 g003
Figure 4. Boosting algorithm modeling process.
Figure 4. Boosting algorithm modeling process.
Ijerph 19 11326 g004
Figure 5. Two-step “k of n” post-processing.
Figure 5. Two-step “k of n” post-processing.
Ijerph 19 11326 g005
Table 1. Details of selected patients.
Table 1. Details of selected patients.
Patient IDGender *AgeNumber of
Total Seizures
Number of Seizures Selected
InT = 30 minInT = 40 min
chb01F11744
chb02M11322
chb04M22422
chb05F7524
chb06F1.51055
chb07F14.5333
chb09F10433
chb10M3745
chb12F24034
chb14F9822
chb15M162054
chb22F9322
chb23F6722
Total----1213942
* Gender: F (female), M (male).
Table 2. Classification evaluation metrics (InT = 30 min).
Table 2. Classification evaluation metrics (InT = 30 min).
Patient IDAcc *Sen *FPR *F1-Score
chb0199.37%99.58%0.00899.38%
chb0297.92%99.17%0.03397.94%
chb0495.83%95.00%0.03395.80%
chb0597.08%97.50%0.03397.10%
chb0681.00%78.67%0.16380.68%
chb0794.72%96.11%0.06794.79%
chb0995.56%95.00%0.03995.53%
chb1092.71%92.92%0.07592.72%
chb1292.78%93.89%0.08392.86%
chb1476.67%79.17%0.25877.24%
chb1581.33%80.67%0.18081.21%
chb2292.50%92.50%0.07592.50%
chb2395.42%94.17%0.03395.36%
Average91.76%91.87%0.08391.78%
* Acc: accuracy; Sen: sensitivity; FPR: false positive rate.
Table 3. Classification evaluation metrics (InT = 40 min).
Table 3. Classification evaluation metrics (InT = 40 min).
Patient IDAcc *Sen *FPR *F1-Score
chb0199.79%100.00%0.00499.79%
chb0295.83%96.67%0.05095.87%
chb0497.08%95.83%0.01797.05%
chb0595.83%97.08%0.05495.88%
chb0681.33%76.33%0.13780.35%
chb0794.72%97.22%0.07894.85%
chb0995.56%96.67%0.05695.60%
chb1093.00%90.67%0.04792.83%
chb1282.29%82.92%0.18382.40%
chb1481.67%77.50%0.14280.87%
chb1597.08%97.08%0.02997.08%
chb2292.50%90.83%0.05892.37%
chb2395.83%95.83%0.04295.83%
Average92.50%91.90%0.06992.37%
* Acc: accuracy; Sen: sensitivity; FPR: false positive rate.
Table 4. Comparison of evaluation metrics of different classifiers.
Table 4. Comparison of evaluation metrics of different classifiers.
Classification MethodInT = 30 minInT = 40 min
Acc *Sen *FPR *Acc *Sen *FPR *
SVM89.41%88.33%0.09589.66%85.62%0.063
RF90.55%90.55%0.09491.45%90.77%0.079
GBDT   ( H 0 x i = 0 )90.33%90.18%0.09591.15%90.68%0.084
RF + GBDT92.00%91.87%0.08392.50%91.90%0.069
* Acc: accuracy; Sen: sensitivity; FPR: false positive rate.
Table 5. Prediction results (InT = 30 min).
Table 5. Prediction results (InT = 30 min).
Patient IDNumber of SeizuresNumber of PredictionsNumber of Missed PredictionsNumber of False Alarms
chb014400
chb022200
chb042200
chb052200
chb065500
chb073300
chb093300
chb104400
chb123300
chb142200
chb155411
chb222200
chb232200
Total393811
Table 6. Prediction results (InT = 40 min).
Table 6. Prediction results (InT = 40 min).
Patient IDNumber of SeizuresNumber of PredictionsNumber of Missed PredictionsNumber of False Alarms
chb014400
chb022200
chb042200
chb054400
chb065500
chb073300
chb093300
chb105500
chb124401
chb142200
chb154400
chb222200
chb23 2200
Total424201
Table 7. Comparison of existing methods.
Table 7. Comparison of existing methods.
AuthorsYearClassifierAcc *Sen *FPR *
Alotaiby et al. [14]2017LDA-89%0.390
Agboola et al. [15]2019SVM-87.26%0.080
Zhang et al. [16]2021Bi-LSTM80.09%-0.260
Rusnac et al. [17]2019MLP91.14%91.37%0.090
Daoud et al. [18]2019DCAE +
Bi-LSTM + CS
99.66%99.72%0.004
Jana et al. [19]2021CNN99.47%97.83%0.0764
Yan et al. [20]2022Transformer-96.01%0.047
This work
(InT = 30 min)
-RF + GBDT91.76%91.87%0.083
This work
(InT = 40 min)
-RF + GBDT92.50%91.90%0.069
* Acc: accuracy; Sen: sensitivity; FPR: false positive rate.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, X.; Lin, M.; Xu, T. Epilepsy Seizures Prediction Based on Nonlinear Features of EEG Signal and Gradient Boosting Decision Tree. Int. J. Environ. Res. Public Health 2022, 19, 11326. https://doi.org/10.3390/ijerph191811326

AMA Style

Xu X, Lin M, Xu T. Epilepsy Seizures Prediction Based on Nonlinear Features of EEG Signal and Gradient Boosting Decision Tree. International Journal of Environmental Research and Public Health. 2022; 19(18):11326. https://doi.org/10.3390/ijerph191811326

Chicago/Turabian Style

Xu, Xin, Maokun Lin, and Tingting Xu. 2022. "Epilepsy Seizures Prediction Based on Nonlinear Features of EEG Signal and Gradient Boosting Decision Tree" International Journal of Environmental Research and Public Health 19, no. 18: 11326. https://doi.org/10.3390/ijerph191811326

APA Style

Xu, X., Lin, M., & Xu, T. (2022). Epilepsy Seizures Prediction Based on Nonlinear Features of EEG Signal and Gradient Boosting Decision Tree. International Journal of Environmental Research and Public Health, 19(18), 11326. https://doi.org/10.3390/ijerph191811326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop