EEG Oscillatory Power and Complexity for Epileptic Seizure Detection

: Monitoring patients at risk of epileptic seizure is critical for optimal treatment and ensuing the reduction of seizure risk and complications. In general, seizure detection is done manually in hospitals and involves time-consuming visual inspection and interpretation by experts of electroencephalography (EEG) recordings. The purpose of this study is to investigate the pertinence of band-limited spectral power and signal complexity in order to discriminate between seizure and seizure-free EEG brain activity. The signal complexity and spectral power are evaluated in ﬁve frequency intervals, namely, the delta, theta, alpha, beta, and gamma bands, to be used as EEG signal feature representation. Classiﬁcation of seizure and seizure-free data was performed by prevalent potent classiﬁers. Substantial comparative performance evaluation experiments were performed on a large EEG data record of 341 patients in the Temple University Hospital EEG seizure database. Based on statistically validated criteria, results show the efﬁciency of band-limited spectral power and signal complexity when using random forest and gradient-boosting decision tree classiﬁers (95% of the area under the curve (AUC) and 91% for both F-measure and accuracy). These results support the use of these automatic classiﬁcation schemes to assist the practicing neurologist interpret EEG records more accurately and without tedious visual inspection. and an F-measure of 90.76%, with a decrease of 1%, compared to the experiments done without the feature selection step. There was no signiﬁcant change in the performance of the three classiﬁers RF, GBDT, and k-NN between the experiments done with and without univariate feature selection. SVM accuracy increased from 78.72% to 81.92%.


Introduction
Epilepsy is a neurological disorder characterized by recurrent epileptic seizures. It affects as many as 50 million people of all ages worldwide. Uncontrolled seizures can lead to a disruption of the nervous system and physical risks such as injury and even death. It is generally accepted that better seizure control is essential for better brain health in the long term. Brain activity can be recorded by one of several means, such as electroencephalography (EEG), magnetoencephalography (MEG), and functional magnetic resonance imaging (fMRI). Of these, EEG is the most prevalent because it is simple to effect, low cost, and has high temporal resolution. Its interpretation is commonly used in the study of various brain disorders such as epilepsy, autism, and attention deficit hyperactivity. It is also the most common method in seizure detection and prediction [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18].
Visual inspection and interpretation of EEG signals are tedious and complex tasks done by experts, which justify research and development of computer-processing methods of EEG that use EEG data representation and machine-learning technics.
Investigations of EEG patterns of seizure have shown that, in general, seizure events have been observed within a wide range of EEG bandwidths. Likewise, their amplitude and using of invasive intracranial EEG records. Furthermore, previous research has used a patient-specific problem formulation and experimental test in which the model is trained and tested on a single patient. Despite attempts to create algorithms to obtain values as accurate as possible, no single one has yet gained widespread acceptance in clinical practice, and the results of studies are far from substituting manual interpretation of the EEG [35,36]. Indeed, despite the numerous commercially available seizure detection devices, there are several issues that must be taken into consideration before such devices can be used in everyday clinical practice. According to a study published in 2021 [37], a comparative study was conducted between 23 devices on the market for seizure detection/alerting. This study showed that these devices rely largely on movement detectors, autonomic change detectors, heart rates, or eye tracking to detect seizures. According to the study's authors, almost no commercially available seizure detection methods use EEG except for Epihunter, a wearable headband device that is connected to a smartphone and is specifically designed for automatic absence seizure detection, which is a specific type of seizure. Thus, the authors concluded that commercially available seizure detection devices do not have the capability of detecting multiple seizure types. Another recent study [38] compared the sensitivity of three well-known commercial seizure detection softwares (Besa, Encevis, and Persyst) to determine whether they correctly detected seizures over long-term videoelectroencephalography monitoring (VEM). Based on results from 81 unseen patients, the researchers found a sensitivity of 67.6% for Besa, 77.8% for Encevis, and 81.6% for Persyst on a patient-by-patient basis. They suggested that the false alarm rate needs to be improved. Furthermore, in another study published in 2020 [39], the authors concluded that most commercially available methods focus on using non-electroencephalography EEG signals, which is perceived as a major limitation given the inability of these devices to detect all seizures types. Therefore, it is crucial to develop an automated system for the detection of epileptic seizures based on EEG signals. To develop such an automated system, significant further research and experimentation are necessary so as to develop seizure detection methods that are clinically relevant, high-performing, and statistically validated. Therefore, testing the clinical relevance of state-of-the-art findings through the analysis of a large number of EEG data collated from an accurate representation of clinical situations and collected from multiple sites will be necessary.
The purpose of this study is to investigate a new data representation model that simultaneously exploits signal oscillatory power in frequency bands and signal complexity as EEG. This representation serves prevailing potent ML algorithms for seizure vs. seizurefree classification. We will rigorously evaluate the proposed EEG data representation features on the publicly available large dataset of Temple University Hospital (TUSZ) collected in clinical settings [40].
The main contribution of this paper can be summarized as follows: (1) An analytical framework for seizure and seizure-free EEG classification is proposed and validated on a new large EEG seizure corpus. (2) The feature significance level is investigated by a univariate data analysis. (3) The performance of complexity measures and oscillatory power is analyzed individually and then combined, as EEG signal representation features, to classify normal and seizure data. (4) The results of four supervised methods are compared: RF, GBDT, SVM, and k-NN. (5) The impact of adding a feature selection step on the performance of classifiers is investigated. (6) The contribution of each channel to the performance of classifiers is examined.

Dataset
EEG data were drawn from the TUH EEG seizure database (TUSZ), which is a part of Temple University Hospital EEG corpus, the largest open-source EEG corpus. The latter comprises more than 16,986 sessions of EEG recordings collected from 10,874 unique subjects [40,41]. The version of the database used in this work was v1.5.1, released in March 2020. In this study, the dataset consists of EEG signals collected from 341 patients, providing Appl. Sci. 2022, 12, 4181 4 of 14 a substantial amount of data for seizure detection investigation. The EEG segments were labelled by a team of neuroscience and bioengineering students who underwent several months of intensive training to gain annotation skills. They have worked closely with a team of neurologists at Temple Hospital to understand their workflow and their clinical needs while labeling EEG segments. There have been numerous revisions of the corpus during the annotation process (two or three annotators per file system), each time applying more refinement and criteria to enhance the clarity and accuracy of the data. Annotations are created by viewing files with open-source software known as an EDF viewer or by using a customized annotation tool developed by the team, which provides time and frequency domain visualizations. The annotation was made available publicly, except for an evaluation corpus not released to the public, to be allocated for research competition purposes. The dataset used contains 886 sessions that were broken down into 7634 files, of which 1780 are seizure activities of different lengths, in seconds, for a total of 40.40 h. The data were collected in real-time clinical environments, including an intensive care unit, an epilepsy monitoring unit, emergency rooms, cardiac intensive care, and surgical and respiratory intensive care units [41]. Most of the EEG recordings have at least 19 electrodes corresponding to the international standard 10/20 system and range from one second to one hour in duration. The sampling rate varied between 250 Hz and 500 Hz. We resampled all data at 256 Hz. In this study, we discarded seizure events that lasted less than our sliding window of 20 s. Eight types of seizures were present in the TUSZ database: focal nonspecific seizure, generalized non-specific seizure, complex partial seizure, simple partial seizure, tonic-clonic seizure, absence seizure, tonic seizure, and myoclonic seizure. Since our goal was not to detect the type of seizure but to determine the presence or absence of seizures, we combined all these types under one label, "seizure." The dataset description is summarized in Table 1, and Figure 1 shows an example of seizure and seizure-free EEG epochs. Detailed information about the database annotation can be found in the following reference [42].

Methodology
An overview of this study's generalized ML-processing strategy is shown in Figure  2. It gives the main steps for data pre-processing, feature extraction, and, finally, classification, which includes training and evaluation. The TUSZ EEG input data consist of

Methodology
An overview of this study's generalized ML-processing strategy is shown in Figure 2. It gives the main steps for data pre-processing, feature extraction, and, finally, classification, which includes training and evaluation. The TUSZ EEG input data consist of pruned EEG recordings, and all uninteresting portions of the data (muscle artifacts, noise, electrode movements, and eye blinks) were purged from the cortical EEG signals. The steps involved in our proposed approach for the classification of epileptic EEG seizures are as follows:

Methodology
An overview of this study's generalized ML-processing strategy is shown in Figure  2. It gives the main steps for data pre-processing, feature extraction, and, finally, classification, which includes training and evaluation. The TUSZ EEG input data consist of pruned EEG recordings, and all uninteresting portions of the data (muscle artifacts, noise, electrode movements, and eye blinks) were purged from the cortical EEG signals. The steps involved in our proposed approach for the classification of epileptic EEG seizures are as follows: Figure 2. A block diagram of this study's automatic seizure detection method. A total of 19 pruned EEG channels from the TUSZ corpus are resampled, filtered, re-referenced, and segmented. Features are subsequently extracted and used as input for modelling and testing four classifiers separately (RF, GB, k-NN, and SVM).

Pre-Processing
In the staging phase, a 60 Hz infinite impulse response notch filter was applied to attenuate the power line, followed by a bandpass filter with lower and higher cut-off frequencies of 0.5 Hz and 75 Hz, respectively. This was followed by resampling to 256 Hz Figure 2. A block diagram of this study's automatic seizure detection method. A total of 19 pruned EEG channels from the TUSZ corpus are resampled, filtered, re-referenced, and segmented. Features are subsequently extracted and used as input for modelling and testing four classifiers separately (RF, GB, k-NN, and SVM).

Pre-Processing
In the staging phase, a 60 Hz infinite impulse response notch filter was applied to attenuate the power line, followed by a bandpass filter with lower and higher cut-off frequencies of 0.5 Hz and 75 Hz, respectively. This was followed by resampling to 256 Hz and re-referencing to the average of all electrodes. From the remaining EEG signal, 20 s fixed-length epochs were segmented.

Feature Extraction
Nine signal complexity measures were derived from the pre-processed EEG: multipermutation entropy (4 levels), sample entropy, wavelet entropies, logarithmic entropy, Shannon entropy, and spectral entropy. The permutation entropy is a complexity measure with a low computational cost for time series based on comparing neighboring values using the distribution of order patterns [16]. The sample entropy quantifies the regularity in the EEG signal regardless of its length [14]. The discrete wavelet transform is used to compute the wavelet entropies [5], whereas the wavelet packet decomposition is used to compute log energy entropy and Shannon entropy to measure the degree of uncertainty in the signal and to evaluate the dynamical order of the signal [4]. The spectral entropy is calculated using the normalized power spectral distribution of the EEG signal [14]. For the frequency domain analysis, we performed the Welch method, which takes an average of the periodograms obtained using fast Fourier transform (FFT). We calculated the absolute power density and relative power density within each of the five frequency sub-bands-delta (1-4 Hz), theta (4-8 Hz), alpha (8)(9)(10)(11)(12)(13)(14), beta (13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30), and gamma (30-70 Hz)-to increase the feature vector size to 19 per channel. The absolute power of a band can be calculated as the sum of all of its power values in that frequency range, whereas the relative power (RP) index for each band was determined based on the absolute power in each frequency band, expressed as a percentage of the absolute power (AP) totaled across all frequency bands. Both absolute and relative PSD analysis are essential for achieving accurate brain analysis. Feature values were averaged using a scalar method before serving as input to the classifier. To investigate whether a subset of the extracted features had a greater significance compared to the entire set, we subdivided the features based on their category: frequency domain and information theory (see Table 2). A univariate feature selection method based on the one-way ANOVA test was also applied for comparison to investigate the effect of reducing the vector size. We evaluated the statistical difference of each feature between the two classes: seizure and seizure-free. In order to identify highly significant features, we sorted features based on their p-values, which correspond to the probability of the results being observed if the null hypothesis H0 is true. In this work, the null hypothesis H0 implies that there is no difference between the means of two groups (seizure and seizure-free records). In general, the lower the p-value, the higher the reliability of results. In medical applications, 0.05, 0.01, or 0.001 are often the recommended threshold values. The features' level of significance was implemented using Python (3.8.0), specifically the Sklearn library.

Training and Evaluation
Comparative supervised classification was carried out with four classifiers: random forest (RF), gradient-boosting decision tree (GBDT), support vector machine (SVM), and k-nearest neighbors (k-NN). The four classifiers were chosen based on findings from previous studies, which showed that they had superior classification compared to others [15,26,27,43,44] The random forest algorithm, an ensemble learning method, involves multiple un-pruned decision tree classifiers. Every tree is an individual classifier, built at the training stage using randomly selected attributes from the original data and at each node to determine the best split. Classification is determined based on majority voting [45]. The gradient-boosting decision tree classifier runs an iterative algorithm to build multiple trees sequentially, where each tree learns and updates its model from the errors of all preceding ones [46]. The support vector machine, an extensively adopted supervised learning classifier, has been widely and successfully used in EEG binary classification problems and specifically in the medical field (automatic seizure detection, mental task classification, epileptic EEG classification, emotion recognition, etc.). The k-nearest neighbors algorithm assigns an unlabeled sample to the most frequently occurring class among the k-nearest labeled training samples. The Euclidean distance is used as the distance metric for k-NN [47].
These proposed ML techniques are based on the most widely used algorithms in previous studies and in various classification problems, as these algorithms have better predictive power and are designed to perform better than linear algorithms, especially for complex non-linearly separable EEG data. Additionally, these techniques have the great advantage of having fewer parameters and hyper-parameters, which makes optimization easier. The model was developed with a view toward being sufficiently transferable to enable the integration in the connected objects for real-time detection.
To avoid overfitting, for each experiment, the dataset was partitioned randomly into training and testing portions, using 10-fold cross-validation repeated 100 times so that onefold served for testing and the remaining k-1 folds were used for training. The parameter settings of the classifiers in all experiments were as follows: for RF, the number of trees was 100, for GBDT it was 100 as well, for SVM the kernel type was linear, and for k-NN the number of nearest neighbors was 5.
To handle imbalanced data, which may occur due to uneven data representation of classes in the training dataset, with a ratio of seizure to seizure-free of approximately 1:4, we corrected the weights of minority and majority classes according to the distribution of the classes in the entire training set. Furthermore, the weights were inversely proportional to the frequency of classes in the dataset. In a final step, the accuracy and the area under the curve (AUC), an effective metric that combines sensitivity and specificity for classification performance evaluation, were computed for each subset of features and also for the entire set. An F-measure representing the harmonic mean of precision and recall was used as well to study performance. For clarity, results are given per 20 s window size input vector and with respect to the three set of features: complexity, power, and a combination of the two. We compared the results with and without the feature selection step.

Univariate Data Analysis
Each feature was subjected to a univariate data analysis to identify its level of significance with the studied classes, i.e., normal or seizure class, and to investigate its discrimination capability in characterizing seizure records. Tables 3 and 4 show the normalized values of the extracted measures (mean ± standard deviation) for seizure and seizure-free segments. Significant p-values (obtained using the ANOVA test) indicate a high discrimination capability in characterizing seizure segments. The values displayed in Tables 3 and 4 correspond to the highest p-value per feature. It is evident from Table 3 that seizure-free records had lower complexity values compared to epileptic seizure records. Results in Table 4 show that the relative power in the alpha band had the lowest p-value. The results also show that oscillatory power change features had higher levels of importance than complexity measures (Tables 3 and 4), which indicates that oscillatory power change in all channels can contribute to better classification. In addition, Shannon entropy and relative power in the alpha band had lower p-values, which indicates that they had higher levels of importance than other features of this category.

Performance of Group of Features Extracted from All Channels without Feature Selection
We evaluated classification performance corresponding to the three sets of features: complexity measures (9 features per channel), oscillatory power (10 features per channel), and a combination of the two (19 features per channel). Table 5 displays the performance of the three sets of features. Classification of seizure versus seizure-free reached an accuracy and F-measure of 90.68% and 91.05%, respectively, using the complexity measures, and 90.95% and 91.33%, respectively, using relative and absolute power in the alpha, beta, gamma, theta, and delta bands. A high AUC value of 95%, an F-measure of 91.41%, and an accuracy of 91.07% were obtained using the RF classifier when the entire set of features was given as classification input. RF and GBDT outperformed SVM and k-NN in all experiments. There was only a slight difference in performance between RF and GBDT.

Performance of Group of Features Extracted from All Channels with Feature Selection
For this study, we executed three sets of experiments to investigate the effect of the recommended thresholds. Consistent with previous findings [48,49], we found that reducing the features by the criteria that their p-values be less than 0.001 led to improved performance. All features with a p-value higher than 0.001 were discarded from the feature vector as non-descriptive. Classification performance was evaluated using AUC, F-measure, and accuracy. Table 6 displays the performance of the three sets of features. RF was able to achieve an accuracy of 90.30% and an F-measure of 90.76%, with a decrease of 1%, compared to the experiments done without the feature selection step. There was no significant change in the performance of the three classifiers RF, GBDT, and k-NN between the experiments done with and without univariate feature selection. SVM accuracy increased from 78.72% to 81.92%.  Table 7 displays the results of RF-and GBDT-based classification using the combined complexity and power features for the 19 EEG channels. A total of 19 features was considered per run. The best AUC, F-measure, and accuracy of 94%, 89.96%, and 89.51% respectively, were obtained with the RF classifier in the Pz channel.

ROC Analysis of RF and GBDT Classifiers
The performance of this study's proposed approach was measured using the receiver operating curve (ROC) analysis and AUC metric. A higher AUC designates better performance. Figure 3 displays a comparison between ROC curves when the entire set of features is used. The AUC for each fold is shown in the caption of each figure, and the mean of the AUC was computed. Figure 3 illustrates the ranking performance of 100 iterations of 10fold cross-validation. The ROC curve summarizes the results and shows that performance was good for all folds, giving an average AUC equal to 95% for RF and 94% for GBDT. The RF classifier yielded better performance for each of the 10 folds compared to GBDT. It is remarkable for both classifiers that the fluctuations of the AUC were smaller among the different folds, indicating the stability of the feature set.

ROC Analysis of RF and GBDT Classifiers
The performance of this study's proposed approach was measured using the receiver operating curve (ROC) analysis and AUC metric. A higher AUC designates better performance. Figure 3 displays a comparison between ROC curves when the entire set of features is used. The AUC for each fold is shown in the caption of each figure, and the mean of the AUC was computed. Figure 3 illustrates the ranking performance of 100 iterations of 10-fold cross-validation. The ROC curve summarizes the results and shows that performance was good for all folds, giving an average AUC equal to 95% for RF and 94% for GBDT. The RF classifier yielded better performance for each of the 10 folds compared to GBDT. It is remarkable for both classifiers that the fluctuations of the AUC were smaller among the different folds, indicating the stability of the feature set.

Comparison of ROC Curves and Accuracies of Classifiers
A comparison between the ROC curves of the four classifiers is displayed in Figure  4. All four classifiers achieved a good AUC. RF produced the highest AUC of 96% and

Comparison of ROC Curves and Accuracies of Classifiers
A comparison between the ROC curves of the four classifiers is displayed in Figure 4. All four classifiers achieved a good AUC. RF produced the highest AUC of 96% and differed slightly from GDBT at 95%. RF outperformed k-NN and SVM in all experiments. The 10-fold cross-validation was performed and the mean and standard deviation of the cross-validation scores over 100 iterations were computed for each classifier. A comparison between accuracies achieved by the classifiers is shown in Figure 5. It is remarkable that the performance of RF classification was superior in both ROC and accuracy. An average accuracy of 90.90% was reached using RF, followed by 89.87%, 89.16%, and 78.72% with GBDT, k-NN, and SVM, respectively.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 11 of 15 differed slightly from GDBT at 95%. RF outperformed k-NN and SVM in all experiments. The 10-fold cross-validation was performed and the mean and standard deviation of the cross-validation scores over 100 iterations were computed for each classifier. A comparison between accuracies achieved by the classifiers is shown in Figure 5. It is remarkable that the performance of RF classification was superior in both ROC and accuracy. An average accuracy of 90.90% was reached using RF, followed by 89.87%, 89.16%, and 78.72% with GBDT, k-NN, and SVM, respectively.

Discussion
This study addressed the problem of seizure vs. seizure-free EEG record classification. The classifier design and validation used data from 341 subjects in the TUSZ database, a larger amount of data than previously used.

Discussion
This study addressed the problem of seizure vs. seizure-free EEG record classification. The classifier design and validation used data from 341 subjects in the TUSZ database, a larger amount of data than previously used.
Several studies have highlighted the advantage of complexity measures as EEG signal description in seizure detection [2,4,5,[9][10][11]13,15,19,20], and others showed the importance, as well, of signal power in different frequency bands [3,7,8,19]. We inquired into the performance of these features with the larger TUSZ database to gain some understanding of the effect of the experimentation database size on seizure detection potency. We followed this with the study of a new method that exploits both types of features. We investigated different feature-grouping options, followed by machine-learning classification of selected features to identify seizure segments.
Two feature subsets of prevailing features, as well as their combination, were investigated to evaluate their seizure detection potency. In contrast to earlier studies that focused on patient-specific seizure detection, this study addressed patient-independent detection, using data from 341 subjects. Signal complexity and power both yielded good classification, with signal power providing slightly better results. When used jointly, signal complexity and power gave superior performance: 95% AUC with RF classification, 94% with GBDT, 86% with SVM, and 88% with k-NN. RF outperformed other classifiers systematically in all experiments. This may be explained by the fact that decisions made by the RF classifier correspond to unpruned and diverse trees that lead to high resolution in feature space. Random operations in the training and voting procedures of RF also contribute to better classification by addressing the issue of overfitting [50].
This study showed the classification significance of the features at the channel level. Unlike most others, which limited features to one channel (e.g., BONN dataset) or six channels (e.g., Freiburg dataset), the TUSZ database used in validating the method in this study has more than 19 channels corresponding to the 10-20 standard system. This allowed us to investigate the performance of the method using only one channel at each run. Results showed that a maximum distinction between seizure and seizure-free records was obtained from the Pz channel, with a corresponding AUC of 94%. In addition, the method results show that the parietal and frontal-parietal regions were effective in extracting features that can discriminate seizure from seizure-free records, wherein the performance for each channel obtained was above 90%. Additionally, this study addressed the effect of adding a feature selection step prior to classification. The analysis removed the less informative features using the p-value test: If the p-value was greater than 0.001, the feature extracted from a specified channel was discarded from the group considered. Only a slight decrease in performance occurred for RF and GBDT with feature selection, keeping in mind that the method is computationally more efficient when running on a smaller set of features.
The results are consistent with findings in a previous study [51] that compared random forest variable selection methods for classification. Indeed, the latter concluded that there was practically no difference in prediction error rates between all methods. However, it was obvious that SVM was more sensitive to the number of features, with a decrease in accuracy rate as the number of features increased.
It may be important to note that the good performance of all classifiers is an indirect indication that the feature set is stable in addition to being representative of EEG data for the problem of seizure detection. An analysis of the significance level of a set of entropy properties determined that seizure records have invariably lower signal complexity than seizure-free records. This finding agrees with previous studies that concluded that the EEG signal during seizure is less complex than when seizure free, and therefore, a reduction in signal information content and complexity can be inferred [11,12,14,15]. Additionally, an analysis of the significance level of relative and absolute signal power in frequency bands, typically, delta (1-4 Hz), theta (4-8 Hz), alpha (8)(9)(10)(11)(12)(13)(14), beta (13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30), and gamma , indicated that these power features can contribute to higher classification, explaining in part their successful use in [3,7,8,19]. A higher level of importance for oscillatory power-based feature vectors in comparison to complexity measures was apparent and could justify their outperformance in classification.
It is noteworthy to mention that one of the limitations of this study is related to the medical frequency sub-bands that were not considered in this work. This is the case of beta frequency Beta-1, 13-15 Hz; Beta-2, 15-18 Hz; Beta-3, 18-25 Hz; and Hi-Beta, 25-30 Hz. Recent reports [52,53] indicate that additional frequency band analysis is beneficial for epilepsy detection and shows the impact of the frequency sub-bands to the epileptic EEG classification accuracy, and the obtained results revealed several frequency sub-band combinations that achieved high classification accuracy, including the medical frequency sub-bands.

Conclusions
This study investigated the use of various ML methods for seizure vs. seizure-free EEG record classification. The classification performance of three sets of extracted featurescomplexity, oscillatory power, and a combination of the two-was evaluated. Our results showed that abnormalities in a seizure EEG record could be distinguished from a seizurefree EEG record by using its complexity measures and oscillatory power, collected as one feature vector fitted to a random forest classifier. This study's proposed approach could provide a prominent contribution to the development of a fully automated seizure detection system. In a future work, we will expand our study to include pre-ictal and ictal. As there are eight types of seizures present in the TUSZ database, we plan to classify seizures by type, incorporating gender as a predictor, since there are certain differences associated with epilepsy patterns between genders. In addition, we will include a cross-database evaluation to validate the effectiveness of our proposed method and to confirm its generalization ability in seizure classification.