Generalisability of epileptiform patterns across time and patients
==================================================================

* Hamid Karimi-Rouzbahani
* Aileen McGonigal

## Abstract

Complexity of epileptogenic zone (EZ) localisation contributes to failure of surgical resection to produce seizure freedom. This is to some extent a result of distinct patterns of epileptiform activity between (i.e., interictal) and during seizures (i.e., ictal) and their diversity across patients. This often leads to suboptimal localisation based on inspection of electroencephalography (EEG) features. We asked two open questions. First, whether neural signal reflecting epileptogenicity would be generalisable from interictal to ictal time window within each patient. This would be critical for patients who are monitored in hospital without having a seizure to help with EZ localisation, and more generally for understanding the predictive power of resting state (interictal) EEG data in determining EZ. Second, whether epileptiform patterns would generalise across patients, and if so, which aspects of those patterns are the most generalisable.

We used an intracranial EEG dataset that included fifty-five patients with lesional and non-lesional pathology, who had subsequently undergone cortical resection in frontal or temporal lobe with different levels of seizure freedom. We extracted a large set of simple to complex features from stereo-EEG (SEEG) and electrocorticographic (ECoG) neural signals recorded during interictal and ictal time windows. We fed those features to decision tree classifiers for EZ localisation and to quantify the diversity of ictal and interictal epileptiform patterns through a cross-time and cross-patient generalisation procedure.

We observed significant evidence (*Bayes factor* >> 10) for generalisability of patterns from interictal to ictal time windows across patients, which were dominantly reflected in signal power and high-frequency network-based features. Majority of patients showed consistent patterns of epileptogenicity across interictal and ictal time windows, reflected in above-chance area-under-curve (mean *AUC =* 0.6). We observed significant evidence (*Bayes factor* >> 10) that signal features of epileptogenic regions could generalise across patients in both interictal and ictal time windows with significant evidence for higher generalisability in ictal than interictal time window (mean *AUC* 0.75 vs. 0.59; *Bayes factor* >> 10). While signal power and moment features were the most contributory to the cross-patient generalisation in the interictal window, signal complexity features were the most contributory in the ictal window.

These results provide new insights about features of epileptic neural activity that generalise across interictal-ictal time windows and patients, which can have implications for both qualitative and quantitative EZ localisation. The explainable machine-learning pipeline developed here can guide future developments in epilepsy investigations.

Keywords
*   quantitative analysis
*   epilepsy surgery
*   intracranial EEG
*   epileptiform activity
*   feature extraction

## Introduction

There are over 50 million people with epilepsy worldwide1. Anti-seizure medications cannot adequately control the disorder in about 30% of cases2. If the epilepsy is considered focal (i.e., seizures arising from part of one hemisphere3), those with drug-resistant focal epilepsy may undergo presurgical evaluation to detect areas involved in the generation of seizure activity, which may require intracranial electroencephalography (EEG) in some. These areas can be collectively referred to as the epileptogenic zone (EZ), a term that was conceptually developed from stereo-electroencephalography (SEEG)4, a method of intracerebral recording based on multiple depth electrodes. The EZ is considered as the region of primary seizure organization5. After localisation, if the clinical risk-benefit ratio is deemed favourable for a specific patient, the EZ can be removed and/or disconnected through surgical resection or laser-based ablation. Despite great progress in use of multimodal approaches (e.g., magnetic resonance imaging (MRI), electroencephalography (EEG), positron emission tomography (PET) scans, etc.) and wealth of clinical expertise, precise localisation of the EZ often remains difficult and may lead to failure to achieve seizure freedom6,7.

Quantification methods have shown great promise in localising the EZ through quantification of intracranial EEG signals8–10(see11 for review). These methods generally investigate either the interictal or the ictal time window (see Supplementary Table 1 for an overview). In the ictal window, low voltage fast activity (LVFA), baseline shift, rhythmic spikes/spike-waves and preictal low frequency spiking, were found to be the most prevalent epileptiform activities12. These characteristics were successfully extracted from signals and used for EZ localisation in several studies13–19. In the interictal window, the traditional epileptiform characteristics include interictal spikes/discharges20 and high-frequency oscillations (HFOs21) with a debate on more efficacy of one over the other, and ultimately possible increased predictive EZ by measuring their co-occurrence22. Modelling approaches based on patient-specific imaging data can predict spatial extent of epileptogenicity23.

The relation between electrical seizure onset and electrical disturbances detectable between seizures is of great clinical and neuroscientific importance and yet remains incompletely known20. Interictal spiking is a heterogeneous phenomenon that reflects the involvement of different neuronal networks and mechanisms (e.g., synaptic conductance) in different regions of an epileptic brain24 and has shown predictive value in differentiating underlying neuropathological substrates25. In fact, the original concept of “epileptic focus” arose not from seizure recordings, but from observations of interictal spiking during ECoG5. Basing surgical excision extent on the zone of interictal spiking measured peri-operatively by ECoG was thought to help improve likelihood of surgical outcome, albeit with variable accuracy26. In pioneering SEEG work27, the regions involved by interictal spiking typically showed (partial) overlap with the zone of primary seizure organization. Observations from clinical data of this type as well as animal models suggest that interictal and ictal signal features may reflect some similar neuronal mechanisms24. Apart from spikes and HFOs, many other linear and nonlinear signal features including randomness, power and entropy have shown success in localising EZ in the interictal window28–35. While many of the classical methods for EZ localisation relied on univariate/single-channel signal activity, there has been a shift to multivariate/network-based localisation10,35–42, which aligns with the conceptualisation of epilepsy as a network disorder43–45 and has shown better performance than univariate methods10,38,42,46–48.

Despite the large literature on EZ localisation, using various signal analysis approaches applied to both interictal and ictal windows (Supplementary Table 1), the correspondence between the two windows has remained unclear. This might be because of the distinct, pre-defined sets of features which were looked for in the two windows and which appear visually different (e.g., interictal spikes and ictal LVFA). Nonetheless, if there are signal features which are shared between the two windows, interictal activity, which represents most of the patient’s brain state and is generally easier to record, could potentially be sufficient to predict EZ. As a first step to explore this, we looked at a heterogeneous group of epilepsies studied with intracranial EEG, available in an open-access dataset46,48. We compared a large battery of explainable signal characteristics, ranging from simple single-channel to computationally complex network-based features, from both interictal and ictal time windows, to see which features generalise across the two time-windows using the data from each individual patient (i.e., within-individual across-time generalisation1). Moreover, to see if there are features which are shared between individuals, we also evaluated the generalisability of features across patients within each individual time window (i.e., across-individual within-time generalisation). Finally, we evaluated the effect of surgical outcome (seizure-free/not seizure-free), EZ, pathology of epilepsy (lesional/non-lesional) and type of recording (SEEG/ECoG) on the EZ localisation performances.

## Materials and methods

### Dataset

This study uses a well-structured open-access intracranial dataset which brings together data from multiple centres46,48. The dataset includes 57 patients who had been implanted with either subdural grid/strip (termed “electrocorticography” (ECoG)47) or SEEG as their presurgical workup, and subsequently treated with surgical resection or laser ablation. Clinically determined seizure onset channels were provided, as well as marking of channels which overlap with the resection/ablation zone, which was rigorously determined by segmenting the resection cavity. Two patients’ data were excluded from our analyses as one had no interictal and the other no ictal recordings. Each patient had 2 interictal recordings and between 1 to 5 (*mean =* 3.7) ictal recordings/seizures (110 interictal and 204 ictal recordings over all patients). The interictal data was selected from awake brain activities determined both by the selection of day-time epochs (8 am – 8 pm) and the use of a custom non-REM sleep detector (explained in detail in Bernabei et al.,48). The interictal data were at least 2 hours before the beginning of a seizure and at least 2 hours after a subclinical seizure, 6 hours after a focal seizure and 12 hours after a generalised seizure. The sampling frequencies of the signals varied across patients and ranged from 256 to 1024 Hz. We adjusted the sampling rate of all datasets to 256 Hz across patients. The details of the patients included in the analyses are provided in Supplementary Table 2. Epileptogenic zones/resected areas ranged from frontal, frontoparietal, mesiofrontal, temporal, mesiotemporal, parietal and insular areas.

### Pre-processing

We used a 5-minute signal from each interictal recording (10 minutes per patient) and a 58-second signal from each ictal recording (−30 to +28 seconds around the time of seizure onset). Bad channels, as marked in the dataset, were excluded from analyses. An average of 105.6 contacts (*std* = 38.04) per patient remained for analysis after bad channels were removed from the dataset. There was an average of 114.2 (*std* = 41.2) and 88.8 (*std* = 25.3) channels recorded in patients implanted with SEEG and ECoG, respectively. Among these, an average of 12.87% (*std* = 11.1%) of channels were in the EZ/resected area in each patient. We applied no filtering or artefact removal on the dataset. As the low- and high-frequency noise is shared across both groups of contacts, and as classifiers rely on the differences between classes rather than similarities49, we did not apply filters. Moreover, by not applying any filters, we allowed easier replication of results in future studies as any choice of filters can potentially affect the results in some way50.

### Feature extraction

We quantified the signal patterns by extracting 34 mathematically distinct features. Features were extracted in 2s non-overlapping sliding windows along the interictal and ictal signals as in previous studies51–54. This led to 14 pre- and 14 post-seizure onset time windows in the ictal period excluding the last window. To quantify changes to neural activities upon seizure onset, we normalised the extracted post-seizure onset data by the pre-seizure onset data using equation (1): ![Formula][1]</img>  where *Pre* and *Post* refer to the arrays of 14 feature values extracted from neural data. This led to 14 normalised ictal feature values which were used for analysis. In interictal data, we down sampled the number of extracted feature samples (n ∼ 150) to 28 samples using the Matlab “resample” function. This led to approximately equal number of data points in interictal and ictal windows, from 110 interictal and 204 ictal recordings. A range of simple to complex signal features were extracted. All these features have been previously used to quantify EEG patterns55–57 and the reader is referred to Supplementary Text 1 and the mentioned publications for details. Briefly, we extracted four categories of signal features to obtain a relatively comprehensive view of signal characteristics. These include the signal *moment* features, nonlinear *complexity* features, *frequency*-domain features and *network*-based features.

### Multivariate pattern classification

We used a standard multivariate pattern classification procedure to localise EZ (i.e., to discriminate epileptogenic/resected and non-epileptogenic/non-resected contacts). We use the term “epileptiform” patterns/activities in a general sense to refer to any patterns which discriminated epileptogenic and non-epileptogenic contacts. Accordingly, the classification performance indicates how discriminable were the signal patterns across these two sets of contacts. We quantified the classification performance by area-under-the-curve (*AUC*) to provide a comprehensive, threshold-free classification performance49. As in recent studies48,58, we used decision tree (DT) classifiers, and each contact was treated like an observation in classification. Our DT classifiers used a random forest algorithm with 50 bags of feature combinations. DT classifiers are well suited for nonlinear feature classifications and provide insights into feature contributions. This method also provides a “feature contribution” metric by permuting the observation/contact labels in each feature separately and quantifying its effect on performance-contribution is in inverse proportion to performance drop. We performed three distinct types of classifications - one within patient and time (non-generalisation) and two which involved generalisation either across time or patients. In all three analyses, we classified epileptogenic and non-epileptogenic contacts (i.e., EZ localisation). In the non-generalisation classification, we performed the classification within the interictal and ictal time window separately for each patient using a 10-fold cross-validation procedure. In the cross-time generalisation, in each patient, we trained the classifier using the data from the two interictal recordings and tested the classifier using all the ictal recordings (*mean =* 3.7). In the cross-patient generalisation, in each time window (interictal/ictal), we trained the classifier using the data from all patients minus one and tested the classifier using the data from the left-out patient and repeated this procedure until every patient was used once for testing the classifier.

To equalise the number of epileptogenic to non-epileptogenic contacts (12.87% vs. 77.13% on average, respectively) which is essential for avoiding bias toward one class in classification49, we used a down-sampling procedure in our analyses and repeated every classification of data 1000 times before averaging the results.

To generate chance-level performances, against which we could evaluate the validity of our true classification performances, we shuffled (epileptogenic/non-epileptogenic) contact labels 1000 times and recalculated the classification performance leading to 1000 chance-level classification results.

### Statistical analysis

We used a Bayes Factor analysis for statistical inference59. We compared the levels of *AUCs* against chance-level *AUCs* as well as evaluated main effects on classifications. We used standard rules of thumb for interpreting levels of evidence60,61: Bayes factors between 3 and 10 and between 1/10 and 1/3 were interpreted as evidence for the alternative and null hypotheses, respectively. Bayes factors > 10 and < 1/10 were interpreted as significant evidence for the alternative and null hypotheses, respectively. We considered the Bayes factors which fell between 1/3 and 3 as insufficient evidence either way. Insufficient evidence means that no conclusions can be made about difference between a pair of variables.

To evaluate the evidence for the null and alternative hypotheses of at-chance and above-chance classification, respectively, we compared the classification rates in each analysis and the classification rates obtained from the chance-level classification results (e.g., panel A in all figures). For that, we performed an unpaired Bayes factor *t-test* for alternative (i.e., difference from chance; H1) and the null (i.e., no difference from chance; H0) hypotheses. To evaluate the evidence for the null and alternative hypotheses of difference between classification levels across analyses (e.g., Interictal vs. Ictal), we compared the classification rates obtained from each of those analyses using paired Bayes factor *t-test*. To evaluate the main effects of resection outcome, EZ, pathology (lesional/non-lesional) and type of recording (SEEG/ECoG), we used a Bayes factor *ANOVA*, with these four factors as independent variables and classification/generalisation *AUC* as the dependent variable. For statistical power in *ANOVA*, we excluded patients with insular, frontoparietal, parietal and mesiofrontal resection which were under-sampled (n < 3). The priors for all Bayes factor analyses were determined based on Jeffrey-Zellner-Siow priors62,63 which are from the Cauchy distribution based on the effect size that is initially calculated in the algorithm using *t-test*59.

### Data and code availability

The dataset used in this study was from previous studies and is available at [https://openneuro.org/datasets/ds004100/versions/1.1.3](https://openneuro.org/datasets/ds004100/versions/1.1.3). The code developed for this project is available at [https://github.com/HamidKarimi-Rouzbahani/Intracranial\_EEG_generalisation](https://github.com/HamidKarimi-Rouzbahani/Intracranial_EEG_generalisation).

## Results

We used a multivariate pattern analysis approach on features extracted from intracranial SEEG/ECoG data in patients with epilepsy to address two main questions. First, we wondered if there were similarities between the epileptiform patterns which discriminated epileptogenic from non-epileptogenic areas in interictal and ictal time windows. Second, we wondered how generalisable epileptiform patterns were across patients.

### How discriminable are epileptogenic and non-epileptogenic contacts?

As an initial step in our analyses, we quantified the discriminability of epileptogenic and non-epileptogenic contacts. This was done for each patient and time window (interictal and ictal) separately. There was significant evidence (*BF* >> 10) for above-chance *AUC* which averaged to 0.97 (*std* = 0.03) in the interictal and 0.91 (*std* = 0.07) in ictal time windows, respectively (Fig. 1A). These showed that our multi-feature classification pipeline could robustly differentiate epileptogenic from non-epileptogenic contacts.

![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/12/2023.08.29.23294708/F1.medium.gif)

[Figure 1](http://medrxiv.org/content/early/2023/12/12/2023.08.29.23294708/F1)

Figure 1 Classification of contacts with and without epileptogenic activities in each patient.
**(A)** Area Under Curve (*AUC*) of classification performance for interictal and ictal classifications. Box plots show the distribution of data, its quartiles and median and whiskers indicate the maximum and minimum of the data over patients. Each dot indicates the data from one patient. Numbers below the bars indicate Bayesian evidence for the difference between true and null classification performances. Horizontal dashed line refers to theoretical chance-level classification (0.5). Bayes factor reflecting evidence for the difference between interictal and ictal classifications are also shown. *BF* >> 10 indicates *BF* > 100. **(B)** Contribution of each feature to the classification performance shown in A, calculated using random permutation. Yellow, pink, green and purple dots indicate moment, complexity, frequency and network-based features. **(C)** *Pearson* correlation between interictal and ictal classification performances across patients with each dot showing one patient. **(D)** *Pearson* correlation between interictal and ictal feature contributions across features with each dot showing one feature. Correlation and the corresponding *p* values are shown on top of panel **C** and **D** with the dashed line showing the theoretical perfect positive correlation.

Interestingly, there was significant evidence (*BF* >> 10) for higher classification in the interictal than ictal time window (Fig. 1A). This points to the importance of interictal neural activities in localisation of EZ. While the classification performances were high across all patients (*AUC* > 0.8), there was no correlation (*r =* −0. 17, *p* = 0.24; *Pearson*) between the level of *AUCs* in interictal and ictal time windows across patients (Fig. 1C). This suggests that patients with the clearest separation between epileptogenic and non-epileptogenic contacts in interictal window did not necessarily show the clearest separation between those contacts in the ictal window and vice versa.

We then evaluated the contribution of each feature to the performance (Fig. 1B). In interictal data, *variance* from the moment features, *Hjorth mobility* from the complexity features, *beta-band powe*r from the frequency features, and *gamma-band coherence* from the network features were among the most contributory features. In ictal data, *variance* from the moment features, *approximate entropy* from the complexity features, *gamma-band power* from the frequency features, and *beta-band coherence* from the network features were among the most contributory features. There was significant correlation between the features’ contributions across the two windows (*r =* 0.73, *p*<0.0001; *Pearson*; Fig. 1D) suggesting that similar sets of features dominantly contributed to the EZ localisation across interictal and ictal time windows.

Next, we evaluated the effect of outcome, region of resection (EZ), pathology and recording type on the classification results in each window (Supplementary Fig. 1). In the interictal data, there was evidence (1/10 < *BF* < 1/3) against any effect of outcome, region of resection, pathology and recording type on classification results. In the ictal data, there was significant evidence (*BF* < 1/10) against any effect of outcome, evidence (3 < *BF* < 10) for an effect of region of resection, insufficient evidence (1/3 < *BF* < 3) for an effect of pathology and evidence (1/10 < *BF* < 1/3) against any effect of recording type on classification results (Supplementary Fig. 1). To check the direction of region of resection effect, we used Bayes-factor *t-test* which showed insufficient evidence (1 < *BF* < 3) for higher classification in patients where the epileptogenic zone/resection was in temporal than frontal and mesiotemporal area (Supplementary Fig. 1).

Our classifications used all signal features simultaneously. To check if any individual feature could predict the resection outcome, we performed a direct comparison (unpaired Bayes factor *t-test*) between feature contributions in patients who became seizure-free (Engel I) vs. not seizure-free (Engel II-IV) outcomes (Supplementary Fig. 2). In interictal data, there was evidence (3 < *BF* < 10) that signal *median* contributed to better EZ localisation in patients who became seizure-free vs. those who did not. However, as *median* was among the least contributory features overall (c.f., Fig. 1B), we prefer not to put too much weight on this result. In ictal data, there was insufficient evidence (1/3 < *BF* < 3) for any feature to predict resection outcome.

As our features relied on signal patterns which were relatively sustained, compared to transient patterns such as interictal spikes or HFOs, we wondered whether accurate classification was possible using even shorter time windows. To test this, we repeated the classifications using the earliest, the middle and the latest 2-second time window of data in interictal and ictal data separately. Interestingly, we found significant evidence (*BF* >> 10) for above-chance *AUC* in both interictal and ictal time, with significant evidence (*BF* >> 10) for higher classification in ictal than interictal data (Supplementary Fig. 3). This repeated the pattern observed when using all windows of data in interictal and ictal periods (c.f., Fig. 1A).

### Do epileptiform patterns generalise across time windows?

We showed that a correlated set of features contributed to EZ localisation in both interictal and ictal windows (Fig. 1C), which might point to shared neural mechanisms underlying signal patterns in both time windows. We wondered if we could localise the EZ in the ictal window based on patterns of interictal activities. To test this, we trained our classifiers on interictal data and tested them on ictal data for each patient separately. We observed that, while the performance was lower (mean *AUC* = 0.60, *std* = 0.1; Fig. 2A) than those obtained by training and testing the classifiers within each time window separately (c.f., Fig. 1A), there was still significant evidence for above-chance cross-time generalisation (interictal to ictal; *BF* >> 10).

![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/12/2023.08.29.23294708/F2.medium.gif)

[Figure 2](http://medrxiv.org/content/early/2023/12/12/2023.08.29.23294708/F2)

Figure 2 Classification of contacts with and without epileptogenic activities for each patient across time windows.
**(A)** *AUC* of cross-time generalisation performance for interictal-to-ictal and ictal-to-interictal generalisations. Box plots show the distribution of data, its quartiles and median and whiskers indicate the maximum and minimum of the data over patients. Each dot indicates the data from one patient. Numbers below the bars indicate Bayesian evidence for the difference between true and null generalisation performances. Horizontal dashed line refers to theoretical chance-level generalisation (0.5). Bayes factor reflecting evidence for the difference between interictal-to-ictal and ictal-to-interictal generalisations are shown. *BF* >> 10 indicates *BF* > 100. **(B)** Contribution of each feature to the generalisation performance shown in A, calculated using random permutation. Yellow, pink, green and purple dots indicate moment, complexity, frequency and network-based features. **(C)** *Pearson* linear correlation between interictal-to-ictal and ictal-to-interictal generalisation performances across patients with each dot showing one patient. **(D)** *Pearson* linear correlation between interictal-to-ictal and ictal-to-interictal feature contributions across features with each dot showing one feature. Correlation and the corresponding *p* values are shown on top of panel **C** and **D** with the dashed line showing the theoretical perfect positive correlation.

Feature contribution results showed an advantage for *variance* from the moment features, *Higuchi fractal dimension* from the complexity features, *delta-band power* from the frequency features, and *gamma-band coherence* from network features (Fig. 2B). There was evidence (1/10 < *BF* < 1/3) against any effect of outcome, region of resection, pathology and insufficient evidence (1/3 < *BF* < 3) for an effect of recording type on the cross-time generalisation results (Supplementary Fig. 4).

Note that we only used two 5-minute windows of interictal recording to train the classifiers, which is relatively short given usual large clinical datasets of interictal activities recorded over several days. The generalisation performance is predicted to improve with higher volumes and more varied sets of training data. To check if increasing the diversity of the training data could improve the generalisation performance, we trained the classifiers using the ictal data and tested them using the interictal data. This would provide the classifiers with a more diverse training set as ictal data were obtained from more recordings than interictal recordings (3.7 vs. 2). Numbers of observations were equalised between interictal and ictal time windows.

While there was significant evidence (*BF* >> 10) for above-chance generalisation performance when training on the ictal data, there was insufficient evidence (*BF* = 0.91) for higher cross-time generalisation when the training data was from the ictal than interictal time windows. Therefore, while a more diverse dataset seems to have improved the classification (shifted the mean *AUC* from 0.6 up to 0.63), more data is needed to establish an improvement effect. There was significant correlation (*r =* 0.60, *p* < 0.001, *Pearson*; Fig. 2C) between the level of performance in interictal-to-ictal and ictal-to-interictal generalisations across patients. This shows that patients who showed the best generalisation from interictal to ictal windows also showed the highest generalisation in the opposite direction. This suggests that each patient has a certain level of similarity between interictal and ictal epileptiform patterns. There was significant correlation between the features’ contributions across the interictal-to-ictal and ictal-to-interictal generalisations (*r =* 0.56, *p* < 0.001; *Pearson*; Fig. 2D) suggesting that generalisable epileptiform patterns were reflected in similar sets of features no matter if generalising from interictal to ictal or vice versa.

In interictal-to-ictal generalisation data, there was evidence or significant evidence (*BF* > 3) that features of *Katz fractal dimension*, *energy ratio*, *theta-band power* led to poorer EZ localisation in patients with seizure-free vs. not seizure-free outcome (Supplementary Fig. 5). In ictal-to-interictal generalisation, this pattern was repeated for features of *energy ratio* and *delta-band power*. These suggest that specific features such as *energy ratio* might be more informative for EZ localisation when they show differences between their interictal and ictal patterns (i.e., as reflected in lower cross-time generalisability; Fig. 2B). This might mean that, patients whose signals’ *energy ratio* changes from interictal to ictal windows (e.g., through a significant increase) have a higher chance for their EZ to be localised; energy ratio change was indeed the basis for the original Epileptogenicity Index method13.

### Do epileptiform patterns generalise across patients?

So far, our analyses focused on within-time classification of contacts and cross-time generalisation of classifications both done within each patient. A clinically important aspect is to ascertain the generalisability of epileptiform patterns across patients, and test the feasibility of using the data from previous patients to help localise the EZ in new out-of-sample patients. To test this, we trained classifiers on the data from all patients minus one and tested the classifiers on the data from the left-out patient. This was done separately for interictal and ictal windows.

In interictal data, there was significant evidence (*BF* >> 10) for above chance cross-patient generalisation (Fig. 3A), which suggests that there were interictal epileptiform patterns which had similarities across patients. We evaluated the features’ contribution to the generalisation (Fig. 3B). Results showed an advantage for *kurtosis* from the moment features, *Higuchi fractal dimension* from the complexity features, *beta-band power* from the frequency features, and *gamma-band coherence* from the network features. There was significant evidence (*BF* < 1/10) against any effect of outcome, evidence (1/10 < *BF* < 1/3) against any effect of region of resection, insufficient evidence (1/3 < *BF* < 3) for an effect of pathology and significant evidence (*BF* > 10) for an effect of recording type on the cross-patient generalisation results (Supplementary Fig. 6). There was significant evidence (*BF* = 12) for higher generalisation to test patients with resection in mesiotemporal than temporal and significant evidence (*BF* >> 10) for higher generalisation to test patients with SEEG than ECoG recording (Supplementary Fig. 6).

![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/12/2023.08.29.23294708/F3.medium.gif)

[Figure 3](http://medrxiv.org/content/early/2023/12/12/2023.08.29.23294708/F3)

Figure 3 Classification of contacts with and without epileptogenic activities across patients.
**(A)** *AUC* of cross-patient generalisation performance for interictal and ictal generalisations. Box plots show the distribution of data, its quartiles and median and whiskers indicate the maximum and minimum of the data over patients. Each dot indicates the data from one patient. Numbers below the bars indicate Bayesian evidence for the difference between true and null generalisation performances. Horizontal dashed line refers to theoretical chance-level generalisation (0.5). Bayes factor reflecting evidence for the difference between interictal and ictal generalisations are shown. *BF* >> 10 indicates *BF* > 100. **(B)** Contribution of each feature to the generalisation performance shown in A, calculated using random permutation. Yellow, pink, green and purple dots indicate moment, complexity, frequency and network-based features. **(C)** *Pearson* linear correlation between interictal and ictal generalisation performances across patients with each dot showing one patient. **(D)** *Pearson* linear correlation between interictal-to-ictal and ictal-to-interictal feature contributions across features with each dot showing one feature. Correlation and the corresponding *p* values are shown on top of panel **C** and **D** with the dashed line showing the theoretical perfect positive correlation.

In ictal data, there was significant evidence (*BF* >> 10) for above-chance cross-patient generalisation (Fig. 3A) which suggests that there were ictal epileptiform patterns which had similarities across patients. Results showed the highest contribution to generalisation by *skewness* from the moment features (however very low compared to other features), *Hjorth mobility* from the complexity features, *mean frequency* from the frequency features, and *gamma-band coherence* from the network features. There was evidence (1/10 < *BF* < 1/3) against any effect of outcome, insufficient evidence (1/3 < *BF* < 3) for an effect of region of resection and pathology and significant evidence (*BF* > 10) for an effect of recording type on the cross-patient generalisation. There was significant evidence (*BF* >> 10) for higher generalisation to test patients with resection in mesiotemporal than temporal and significant evidence (*BF* >> 10) for higher generalisation to test patients with SEEG than ECoG recording (Supplementary Fig. 6).

There was significant correlation (*r =* 0.73, *p* < 0.001; *Pearson*) between the level of generalisation in interictal and ictal time windows across patients (Fig. 3C). This suggests that if a patient’s data (i.e., testing set) has similarities to the pool of other patients’ data (i.e., training set), this will reflect in both interictal and ictal generalisations. On the other hand, patients with very distinct epileptiform activity patterns show this distinction across both interictal and ictal time windows. There was no correlation (*r =* −0.23, *p* = 0.19; *Pearson*) between informative features in the interictal and ictal time windows. This suggests that distinct sets of features contributed to cross-patient generalisation in interictal and ictal time windows.

In interictal data, there was significant evidence (*BF* >> 10; Bayes factor *t-test*) that *approximate entropy* led to better EZ localisation in patients who became seizure-free vs. those who did not (Supplementary Fig. 7). In ictal data, there was evidence or significant evidence (*BF* > 3) that *mean* and *median* led to better EZ localisation in patients who became seizure-free vs. those who did not and *median frequency* led to the opposite pattern. As these were all among the least contributory features to the overall cross-patient generalisation performance we avoid over-interpreting them.

## Discussion

Precise localisation of epileptogenic zone has remained a challenging problem. It has recently been highlighted that data and code sharing are fundamental to moving computational epilepsy studies towards clinical translation11. The current work uses one of the few open-access SEEG/ECoG datasets which provides epochs of interictal and ictal activities with meticulous labelling of electrode contacts, resection volume and clinical information including surgical outcome46,48. Using a rigorous ML-based pipeline, our study shows the feasibility of establishing generalisability of patterns within individuals from interictal to ictal periods, and across individuals during both interictal and ictal time windows. These results make several contributions to EZ localisation, explained below.

First, to establish that our ML-based method could discriminate areas with and without epileptiform activities, we classified epileptogenic and non-epileptogenic areas (EZ localisation, based on clinician labelling of the dataset) within each patient and found significant differences between the two classes of areas. At the individual patient level, we observed a remarkable EZ localisation performance in the interictal time window, which counter-intuitively surpassed that obtained in the ictal time window. Therefore, while majority of available studies have been developed to localise EZ during the ictal time window13–19, this study finds significant information within the interictal signal that can be utilised for EZ localisation. It is important to note that our approach did not pre-select any specific feature (such as spikes, HFO) as biomarkers *a priori*, but rather analysed the ensemble of the neural signal over time, which appears to contain predictive information beyond these well-known features, even when analysing time-windows as short as 2 seconds (c.f., Supplementary Fig. 3). This agrees with reports showing that interictal patterns are relatively stable over time36,40,64 (also see65,66).

While systematic comparison between interictal and ictal epileptiform patterns are rare in the literature, our result aligns with a surface EEG study in children with MRI-visible lesions which found better predictive value of interictal rather than ictal data67. We observed that the most informative features in interictal data included *variance*, *beta-band power*, *correlation* and *gamma-band coherence*, which in order support previous studies finding information in multiscale entropy in the gamma band28,30 (a complexity measure) and relative entropy34(a network measure). Interictal network studies have shown a gradient of within-area connectivity decreasing progressively from epileptogenic regions to propagation regions to non-involved zones39,41,54, providing a proxy for EZ localisation. The prominence of network measures in all of our analyses also aligns with a recent study showing greater information in network measures compared to univariate power-based measures in EZ localisation48. While we cannot perfectly equalise the interictal and ictal data for fair comparison as there are systematic differences in their collection time, number of epochs and potential artefacts in ictal signals, these results provide valuable insights into the richness of interictal activity patterns for EZ localisation.

The remarkable classification performance observed especially in the interictal window (> 0.90; Fig. 1A) supports the value of ML-based interictal localisation methods, which here used multiple features. While these features could have overlapped in their selectivity, they worked in synergy to detect as much non-overlapping information as possible. DT classifiers are good at combining distinct combinations of features to generate representational spaces in which classes can be separated. We also tested Support Vector Machine and Linear Discriminant Analysis classifiers, but both provided poorer classification. Our approach of using a range of simple-to-complex features is different from older localisation methods which used one main feature such as high-to-low frequency energy ratio13,16, and aligns with later studies which have combined several features for EZ localisation17,31,68,69, temporal detection of seizures70 and quantification of seizure severity71.

Second, we found prediction power in activity patterns of interictal data to localise EZ in the ictal period. While previous studies have localised EZ in interictal and ictal windows, the interictal-ictal correspondence has not been systematically investigated. In eleven patients with epilepsy who had been implanted with ECoG, one study found that the template of connectivity-based ictal epileptogenic areas could be helpful in informing the localisation of EZ interictally72. One recent study investigated the fine-grained timing and direction of interictal and ictal discharges using microelectrode grids, and suggested that interictal discharges are traveling waves that traverse the same path as seizure discharges73. After confirming a consistent temporal ordering of discharges in interictal and ictal windows, another study developed a novel source localisation method based on wave propagation, which successfully localised the EZ74. Our work evaluates the generalisability of a large battery of epileptiform features from interictal to ictal time windows and vice versa. The drop in cross-time generalisation performances (Fig. 2A) compared to the non-generalisation analysis (c.f., Fig. 1A) is supported by known differences in visualisable epileptogenic patterns from the interictal to the ictal time window75,76. In our study, the observed above-chance cross-time generalisation had not been necessarily predictable. It could be the case that epileptogenic areas would show higher value of a particular feature (e.g., power) than non-epileptogenic areas in the interictal window with this pattern flipping in the ictal window. This would have been detectable by our machine learning pipeline and would have reflected in below-chance (AUC < 0.5) generalisation performance. The potential of using interictal recording to predict EZ localisation is significant because many patients have insufficient or sometimes no seizures during their one/two-week hospitalisation for EZ localisation. This is an important limitation for visual localisation of EZ based on electrical patterns during seizures and for training ML algorithms, which, like humans, need enough samples to learn and localise epileptiform patterns from the data.

The third contribution of this work is showing that, despite clear inter-subject differences, there were patterns of epileptiform activities which were shared across patients. Machine learning allowed us to train the classifiers using data from one set of patients and test the generalisability of patterns to the data from a new out-of-sample patient. While these performances were expectedly lower in the cross-patient generalisation than within-patient classification (Fig. 3A vs. Fig. 1A), this result is promising and informative. The decrease in performance can be explained by large differences across patients’ data including epilepsy characteristics as well as distinct sampling of the brain, recording type, etc. Moreover, significant inter-subject differences may be present in terms of patient-specific epileptogenic “signatures”, the features of which are detectable across both interictal and ictal time windows for that individual77. On the other hand, a few studies showed that specific patterns can be generalised across patients, but only evaluated it in either interictal or ictal time window (Supplementary Table 1). A universal repertoire of seizure patterns across species has previously been observed, which suggests that some invariant properties characterise seizures under different physiological and pathological conditions78. Here, we showed that *ictal* epileptogenic patterns, especially those captured by complexity features, were more generalisable across subjects than *interictal* patterns. The present results also showed that SEEG recordings provided advantageous generalisability compared to ECoG (Supplementary Fig. 6). This is likely explained by SEEG’s greater sampling of a wider range of brain structures, which contributes to a more consistent sampling of the brain across patients79. These cross-patient generalisable patterns make it possible and desirable to use them on new out-of-sample data, to potentially build on these results by testing larger datasets. To facilitate future testing in novel datasets, we have shared our Matlab scripts.

One of the main concerns when using artificial intelligence in applications such as EZ localisation is the explainability of algorithms. Lack of knowledge about how a ML algorithm decides why a contact is classified as “epileptogenic” makes the algorithm less trustworthy for clinicians80, who may not be able to validate if a specific feature of a signal is indicative of epileptogenicity or whether the algorithm is simply wrong. Methodologies incorporating explainable features can mitigate the explainability issue and provide complementary insights into the growing body of work in EZ localisation and seizure prediction, which tend to adopt unexplainable ML algorithms such as deep neural networks81,82. Every individual feature used in our work has clear mathematical definition and has been validated in previous quantification analysis of neural data55,56,83,84. We also quantified the contribution of each feature in our analyses, thus avoiding the “black box” effect encountered when using algorithms such as deep neural networks. Accordingly, our proposed pipeline can be added as a primary feature extractor to prediction pipelines to make them more explainable to humans. It is of note that, while the mathematical definition of each of our features are clear, the neurophysiological correlates of these features needs to be sought for in the future.

There are several future directions which can facilitate the translation of this work to clinical practice. One can come through the improvement of the classification and generalisation performance. We used relatively short time windows of interictal (5-minute windows) and ictal (1-minute windows) data, both of which can be lengthened to potentially improve the classification performance. We did not apply any filtering or artefact removal, as ML classification algorithms are mainly sensitive to distinct patterns between classes (i.e., contacts with and without epileptiform patterns) rather than patterns which are common between classes (e.g., line noise). Nonetheless, one future direction would be to test if application of filters or artefact removal algorithms in the pre-processing stage can improve classification performance. Rather than an optimised work, the current study was only a feasibility effort to establish the generalisability of patterns of activities across time within patients and within time across patients. Another future direction would be to test the generalisability of the classification pipeline to datasets from other centres, which have undergone meticulous clinical evaluation and labelling. Finally, it would be interesting to evaluate the generalisability of the proposed pipeline to non-invasive modalities such as scalp EEG and magnetoencephalography (MEG). As the features extracted here are not modality-specific and rely on characteristics of time series, these methods can be applied to sensor- or source-space E/MEG data. The generalisability of these methods to non-invasive modalities is significant because current gold-standard invasive methods of SEEG and ECoG suffer from incomplete or sub-optimal sampling of the brain, and in addition these invasive methods are only indicated in a subset of patients. Delineation of likely spatial extent of epileptogenic zones could in theory be optimised using rigorous localisation methods developed here in conjunction with non-invasive E/MEG recording modalities, especially in the interictal time window. This will provide a more objective and fully automatic method for the localisation of EZ than current methods which often rely on visual detection, manual annotation and operator-dependent analysis of epileptiform patterns85,86.

In conclusion, we showed that powerful classification patterns were embedded within the EEG signal, which could reliably differentiate epileptogenic from non-epileptogenic contacts in every individual. Such patterns could be identified in both interictal and ictal recordings through features such as signal *variance*, *Hjorth mobility* and *complexity* as well as high-frequency power and network features, without taking account of any predetermined figures such as spikes, HFO or known ictal patterns. There were also features that could correctly predict EZ in ictal recordings from interictal recordings. Again, high-frequency *power* and network features were the most contributory features here. Finally, we showed that, while there were differences between epileptiform patterns across patients suggesting subject-specific effects, we could localise the EZ with well above chance precision using interictal and more dominantly ictal activities. The proposed methods and results provide new evidence for generalisability of epileptiform patterns across time and patients and open new avenues for future methods developed for epileptogenic zone localisation. Clues from neural signal changes could also provide new directions for investigating the biological correlates of interictal87 and ictal88 epileptiform activity. Their explainable nature is important for further investigation of pathophysiologic underpinnings of these signal changes. This could help contribute to efforts to develop paradigm-shifting therapeutic possibilities in epilepsy including disease-modifying treatments89, as well as further refining network-based surgical treatments90, seizure forecasting91 and seizure detection92.

## Supporting information

Supplementary material [[supplements/294708_file03.pdf]](pending:yes)

## Data Availability

The dataset used in this study was from previous studies and is available at [https://openneuro.org/datasets/ds004100/versions/1.1.3](https://openneuro.org/datasets/ds004100/versions/1.1.3)

[https://openneuro.org/datasets/ds004100/versions/1.1.3](https://openneuro.org/datasets/ds004100/versions/1.1.3)

## Funding

No funding was received towards this work.

## Competing interests

The authors report no competing interests.

## Acknowledgements

We would like to thank Mater Research Institute for supporting this study and Professor Brian Litt’s group for sharing their dataset online which allowed us to test our methods on their data.

## Footnotes

*   Correction of a few figure captions and addition of a few clarification statements to the manuscript.

*   1 In this manuscript, the word “generalisation” refers to testing machine learning classifiers on data from unseen time windows/patients rather than the conventional epilepsy definition of seizure propagation over the brain.

*   Received August 29, 2023.
*   Revision received December 11, 2023.
*   Accepted December 12, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Kwan P, Brodie MJ. Early Identification of Refractory Epilepsy. New England Journal of Medicine. 2000;342(5):314–319. doi:10.1056/NEJM200002033420503
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM200002033420503&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10660394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085070300003&link_type=ISI) 

2.  2.Chen Z, Brodie MJ, Liew D, Kwan P. Treatment Outcomes in Patients With Newly Diagnosed Epilepsy Treated With Established and New Antiepileptic Drugs. JAMA Neurol. 2018;75(3):279. doi:10.1001/jamaneurol.2017.3949
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamaneurol.2017.3949&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29279892&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

3.  3.Fisher RS, Cross JH, French JA, et al. Operational classification of seizure types by the International League Against Epilepsy: Position Paper of the ILAE Commission for Classification and Terminology. Epilepsia. 2017;58(4):522–530. doi:10.1111/epi.13670
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.13670&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28276060&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

4.  4.Talairach J, Bancaud J, Bonis A, Szikla G, Tournoux P. Functional Stereotaxic Exploration of Epilepsy. Stereotact Funct Neurosurg. 1962;22(3-5):328–331. doi:10.1159/000104378
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000104378&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=13984743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A19629335A00047&link_type=ISI) 

5.  5.Chauvel P. The Epileptogenic Zone: A Critical Reconstruction. In: A Practical Approach to Stereo EEG. Springer Publishing Company; 2020. doi:10.1891/9780826136930.0010
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1891/9780826136930.0010&link_type=DOI) 

6.  6.Najm I, Jehi L, Palmini A, Gonzalez-Martinez J, Paglioli E, Bingaman W. Temporal patterns and mechanisms of epilepsy surgery failure. Epilepsia. 2013;54(5):772–782. doi:10.1111/epi.12152
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.12152&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23586531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000318241000008&link_type=ISI) 

7.  7.de Tisi J, Bell GS, Peacock JL, et al. The long-term outcome of adult epilepsy surgery, patterns of seizure remission, and relapse: a cohort study. The Lancet. 2011;378(9800):1388–1395. doi:10.1016/S0140-6736(11)60890-8
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(11)60890-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22000136&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296269600025&link_type=ISI) 

8.  8.Iasemidis LD, Chris Sackellares J, Zaveri HP, Williams WJ. Phase space topography and the Lyapunov exponent of electrocorticograms in partial seizures. Brain Topogr. 1990;2(3):187–201. doi:10.1007/BF01140588
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF01140588&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2116818&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

9.  9.Bartolomei F, Lagarde S, Wendling F, et al. Defining epileptogenic networks: Contribution of SEEG and signal analysis. Epilepsia. 2017;58(7):1131–1147. doi:10.1111/epi.13791
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.13791&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28543030&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

10. 10.Li A, Huynh C, Fitzgerald Z, et al. Neural fragility as an EEG marker of the seizure onset zone. Nat Neurosci. 2021;24(10):1465–1474. doi:10.1038/s41593-021-00901-w
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-021-00901-w&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34354282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

11. 11.Bernabei JM, Li A, Revell AY, et al. Quantitative approaches to guide epilepsy surgery from intracranial EEG. Brain. 2023;146(6):2248–2258. doi:10.1093/brain/awad007
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awad007&link_type=DOI) 

12. 12.Lagarde S, Buzori S, Trebuchon A, et al. The repertoire of seizure onset patterns in human focal epilepsies: Determinants and prognostic values. Epilepsia. 2019;60(1):85–95. doi:10.1111/epi.14604
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.14604&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30426477&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

13. 13.Bartolomei F, Chauvel P, Wendling F. Epileptogenicity of brain structures in human temporal lobe epilepsy: A quantified study from intracerebral EEG. Brain. 2008;131(7):1818–1830. doi:10.1093/brain/awn111
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awn111&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18556663&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257397100015&link_type=ISI) 

14. 14.Aubert S, Wendling F, Regis J, et al. Local and remote epileptogenicity in focal cortical dysplasias and neurodevelopmental tumours. Brain. 2009;132(11):3072–3086. doi:10.1093/brain/awp242
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awp242&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19770216&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000271389300017&link_type=ISI) 

15. 15.Gnatkovsky V, Francione S, Cardinale F, et al. Identification of reproducible ictal patterns based on quantified frequency analysis of intracranial EEG signals. Epilepsia. 2011;52(3):477–488. doi:10.1111/j.1528-1167.2010.02931.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1528-1167.2010.02931.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21269289&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

16. 16.David O, Blauwblomme T, Job AS, et al. Imaging the seizure onset zone with stereo-electroencephalography. Brain. 2011;134(10):2898–2911. doi:10.1093/brain/awr238
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awr238&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21975587&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295681400011&link_type=ISI) 

17. 17.Grinenko O, Li J, Mosher JC, et al. A fingerprint of the epileptogenic zone in human epilepsies. Brain. 2018;141(1):117–131. doi:10.1093/brain/awx306
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awx306&link_type=DOI) 

18. 18.Li J, Grinenko O, Mosher JC, Gonzalez-Martinez J, Leahy RM, Chauvel P. Learning to define an electrical biomarker of the epileptogenic zone. Hum Brain Mapp. 2020;41(2):429–441. doi:10.1002/hbm.24813
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.24813&link_type=DOI) 

19. 19.Nakatani M, Inouchi M, Daifu-Kobayashi M, et al. Ictal direct current shifts contribute to defining the core ictal focus in epilepsy surgery. Brain Commun. 2022;4(5). doi:10.1093/braincomms/fcac222
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/braincomms/fcac222&link_type=DOI) 

20. 20.De Curtis M, Jefferys JGR, Avoli M. Interictal Epileptiform Discharges in Partial Epilepsy: Complex Neurobiological Mechanisms Based on Experimental and Clinical Evidence.; 2012. [https://www.ncbi.nlm.nih.gov/books/](https://www.ncbi.nlm.nih.gov/books/)
    
    

21. 21.Thomschewski A, Hincapié AS, Frauscher B. Localization of the epileptogenic zone using high frequency oscillations. Front Neurol. 2019;10(FEB). doi:10.3389/fneur.2019.00094
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fneur.2019.00094&link_type=DOI) 

22. 22.Roehri N, Pizzo F, Lagarde S, et al. High-frequency oscillations are not better biomarkers of epileptogenic tissues than spikes. Ann Neurol. 2018;83(1):84–97. doi:10.1002/ana.25124
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ana.25124&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29244226&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

23. 23.Hashemi M, Vattikonda AN, Sip V, et al. The Bayesian Virtual Epileptic Patient: A probabilistic framework designed to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. Neuroimage. 2020;217:116839. doi:10.1016/j.neuroimage.2020.116839
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2020.116839&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32387625&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

24. 24.Avoli M, Biagini G, de Curtis M. Do Interictal Spikes Sustain Seizures and Epileptogenesis? Epilepsy Curr. 2006;6(6):203–207. doi:10.1111/j.1535-7511.2006.00146.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1535-7511.2006.00146.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17260060&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

25. 25.Sklenarova B, Zatloukalova E, Cimbalnik J, et al. Interictal <scp>High-Frequency</scp> Oscillations, Spikes, and Connectivity Profiles: a Fingerprint of Epileptogenic Brain Pathologies. Epilepsia. Published online August 17, 2023. doi:10.1111/epi.17749
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17749&link_type=DOI) 

26. 26.Jasper HH, Arfel-Capdeville G, Rasmussen T. Evaluation of EEG and cortical electrographic studies for prognosis of seizures following surgical excision of epileptogenic lesions. Epilepsia. 1961;2:130–137.
    
    

27. 27.Bancaud J. La stéréo-électroencéphalographie dans l’épilepsie: informations neurophysiopathologiques apportées par l’investigation fonctionnelle stéreotaxique: rapport présenté a la Société d’électroencéphalographie de langue franc̦ aise. Published online October 1965.
    
    

28. 28.Andrzejak RG, Schindler K, Rummel C. Nonrandomness, nonlinear dependence, and nonstationarity of electroencephalographic recordings from epilepsy patients. Phys Rev E Stat Nonlin Soft Matter Phys. 2012;86(4). doi:10.1103/PhysRevE.86.046206
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevE.86.046206&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23214662&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

29. 29.Mooij AH, Frauscher B, Amiri M, Otte WM, Gotman J. Differentiating epileptic from non-epileptic high frequency intracerebral EEG signals with measures of wavelet entropy. Clinical Neurophysiology. 2016;127(12):3529–3536. doi:10.1016/j.clinph.2016.09.011
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2016.09.011&link_type=DOI) 

30. 30.Sato Y, Ochi A, Mizutani T, Otsubo H. Low entropy of interictal gamma oscillations is a biomarker of the seizure onset zone in focal cortical dysplasia type II. Epilepsy and Behavior. 2019;96:155–159. doi:10.1016/j.yebeh.2019.01.030
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.yebeh.2019.01.030&link_type=DOI) 

31. 31.Cimbalnik J, Klimes P, Sladky V, et al. Multi-feature localization of epileptic foci from interictal, intracranial EEG. Clinical Neurophysiology. 2019;130(10):1945–1953. doi:10.1016/j.clinph.2019.07.024
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2019.07.024&link_type=DOI) 

32. 32.Mooij AH, Frauscher B, Gotman J, Huiskamp GJM. A skew-based method for identifying intracranial EEG channels with epileptic activity without detecting spikes, ripples, or fast ripples. Clinical Neurophysiology. 2020;131(1):183–192. doi:10.1016/j.clinph.2019.10.025
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2019.10.025&link_type=DOI) 

33. 33.Taylor PN, Papasavvas CA, Owen TW, et al. Normative brain mapping of interictal intracranial EEG to localize epileptogenic tissue. Brain. 2022;145(3):939–949. doi:10.1093/brain/awab380
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awab380&link_type=DOI) 

34. 34.Travnicek V, Klimes P, Cimbalnik J, et al. Relative entropy is an easy-to-use invasive electroencephalographic biomarker of the epileptogenic zone. Epilepsia. 2023;64(4):962–972. doi:10.1111/epi.17539
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17539&link_type=DOI) 

35. 35.Gallagher R, Sinha N, Pattnaik A, et al. Quantifying interictal intracranial EEG to predict focal epilepsy. ArXiv. 2023;arXiv:2307.15170.
    
    

36. 36.Dauwels J, Eskandar E, Cash S. Localization of seizure onset area from intracranial non-seizure EEG by exploiting locally enhanced synchrony. In: Proceedings of the 31st Annual International Conference of the IEEE Engineering in Medicine and Biology Society: Engineering the Future of Biomedicine, EMBC 2009. IEEE Computer Society; 2009:2180–2183. doi:10.1109/IEMBS.2009.5332447
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/IEMBS.2009.5332447&link_type=DOI) 

37. 37.Gazit T, Doron I, Sagher O, et al. Time-frequency characterization of electrocorticographic recordings of epileptic patients using frequency-entropy similarity: A comparison to other bi-variate measures. J Neurosci Methods. 2011;194(2):358–373. doi:10.1016/j.jneumeth.2010.10.011
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2010.10.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20969891&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

38. 38.Li A, Chennuri B, Subramanian S, et al. Using network analysis to localize the epileptogenic zone from invasive EEG recordings in intractable focal epilepsy. Network Neuroscience. 2018;2(2):218–240. doi:10.1162/netn\_a\_00043
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1162/netn_a_00043&link_type=DOI) 

39. 39.Lagarde S, Roehri N, Lambert I, et al. Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies. Brain-A Journal of Neurology. 2018;(10):2966–2980. doi:10.1093/brain/awy214ï
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awy214&link_type=DOI) 

40. 40.Gunnarsdottir KM, Li A, Smith RJ, et al. Source-sink connectivity: a novel interictal EEG marker for seizure localization. Brain. 2022;145(11):3901–3915. doi:10.1093/brain/awac300
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awac300&link_type=DOI) 

41. 41.Johnson GW, Doss DJ, Morgan VL, et al. The Interictal Suppression Hypothesis in focal epilepsy: network-level supporting evidence. Brain. 2023;146(7):2828–2845. doi:10.1093/brain/awad016
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awad016&link_type=DOI) 

42. 42.Runfola C, Sheheitli H, Bartolomei F, Wang H, Jirsa V. In pursuit of the epileptogenic zone in focal epilepsy:a dynamical network biomarker approach. Commun Nonlinear Sci Numer Simul. 2023;117:106973. doi:10.1016/j.cnsns.2022.106973
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cnsns.2022.106973&link_type=DOI) 

43. 43.Wendling F, Bellanger JJ, Bartolomei F, Chauvel P. Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals. Biol Cybern. 2000;83(4):367–378. doi:10.1007/s004220000160
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s004220000160&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11039701&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089294300005&link_type=ISI) 

44. 44.Spencer SS. Neural Networks in Human Epilepsy: Evidence of and Implications for Treatment. Epilepsia. 2002;43(3):219–227. doi:10.1046/j.1528-1157.2002.26901.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1528-1157.2002.26901.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11906505&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174557800002&link_type=ISI) 

45. 45.Kramer MA, Cash SS. Epilepsy as a disorder of cortical network organization. Neuroscientist. 2012;18(4):360–372. doi:10.1177/1073858411422754
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1073858411422754&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22235060&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306732400010&link_type=ISI) 

46. 46.Kini LG, Bernabei JM, Mikhail F, et al. Virtual resection predicts surgical outcome for drug-resistant epilepsy. Brain. 2019;142(12):3892–3905. doi:10.1093/brain/awz303
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awz303&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31599323&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

47. 47.Balatskaya A, Roehri N, Lagarde S, et al. The “Connectivity Epileptogenicity Index” (cEI), a method for mapping the different seizure onset patterns in StereoElectroEncephalography recorded seizures. Clinical Neurophysiology. 2020;131(8):1947–1955. doi:10.1016/j.clinph.2020.05.029
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2020.05.029&link_type=DOI) 

48. 48.Bernabei JM, Sinha N, Arnold TC, et al. Normative intracranial EEG maps epileptogenic tissues in focal epilepsy. Brain. 2022;145(6):1949–1961. doi:10.1093/brain/awab480
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awab480&link_type=DOI) 

49. 49.Hart P, Stork D, Duda R. Pattern Classification. Hoboken: Wiley; 2000.
    
    

50. 50.Trübutschek D, Yang YF, Gianelli C, et al. EEGManyPipelines: A large-scale, grass-root multi-analyst study of EEG analysis practices in the wild. OSF preprint, doi:10.31222/osf.io/jq342. Published online 2022.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.31222/osf.io/jq342&link_type=DOI) 

51. 51.Kerr MSD, Burns SP, Gale J, Gonzalez-Martinez J, Bulacio J, Sarma S V. Multivariate analysis of SEEG signals during seizure. In: *Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society*, EMBS.; 2011:8279–8282. doi:10.1109/IEMBS.2011.6092041
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/IEMBS.2011.6092041&link_type=DOI) 

52. 52.Wendling F, Chauvel P, Biraben A, Bartolomei F. From intracerebral EEG signals to brain connectivity:Identification of epileptogenic networks in partial epilepsy. Front Syst Neurosci. 2010;4. doi:10.3389/fnsys.2010.00154
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnsys.2010.00154&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21152345&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

53. 53.Schindler K, Leung H, Elger CE, Lehnertz K. Assessing seizure dynamics by analysing the correlation structure of multichannel intracranial EEG. Brain. 2007;130(1):65–77. doi:10.1093/brain/awl304
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awl304&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17082199&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000243061500006&link_type=ISI) 

54. 54.Shah P, Bernabei JM, Kini LG, et al. High interictal connectivity within the resection zone is associated with favorable post-surgical outcomes in focal epilepsy patients. Neuroimage Clin. 2019;23. doi:10.1016/j.nicl.2019.101908
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nicl.2019.101908&link_type=DOI) 

55. 55.Karimi-Rouzbahani H, Shahmohammadi M, Vahab E, Setayeshi S, Carlson T. Temporal Variabilities Provide Additional Category-Related Information in Object Category Decoding: A Systematic Comparison of Informative EEG Features. Neural Comput. Published online August 30, 2021:1–46. doi:10.1162/neco\_a\_01436
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1162/neco_a_01436&link_type=DOI) 

56. 56.Karimi-Rouzbahani H, Woolgar A. When the Whole Is Less Than the Sum of Its Parts: Maximum Object Category Information and Behavioral Prediction in Multiscale Activation Patterns. Front Neurosci. 2022;16. doi:10.3389/fnins.2022.825746
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnins.2022.825746&link_type=DOI) 

57. 57.Karimi-Rouzbahani H. Evidence for multiscale multiplexed representation of visual features in EEG. BioRxiv doi:10.1101/20230709548296. Published online 2023.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/20230709548296&link_type=DOI) 

58. 58.Jiang H, Kokkinos V, Ye S, et al. Interictal SEEG Resting-State Connectivity Localizes the Seizure Onset Zone and Predicts Seizure Outcome. Advanced Science. 2022;9(18). doi:10.1002/advs.202200887
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/advs.202200887&link_type=DOI) 

59. 59.Rouder JN, Morey RD, Speckman PL, Province JM. Default Bayes factors for ANOVA designs. J Math Psychol. 2012;56(5):356–374. doi:10.1016/j.jmp.2012.08.001
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jmp.2012.08.001&link_type=DOI) 

60. 60.Dienes Z. Using Bayes to get the most out of non-significant results. Front Psychol. 2014;5. doi:10.3389/fpsyg.2014.00781
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpsyg.2014.00781&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25120503&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

61. 61.Lee MD, Wagenmakers EJ. Bayesian statistical inference in psychology: Comment on Trafimow (2003). Psychol Rev. 2005;112(3):662–668. doi:10.1037/0033-295X.112.3.662
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/0033-295X.112.3.662&link_type=DOI) 

62. 62.Jeffreys H. The Theory of Probability. OuP Oxford; 1998.
    
    

63. 63.Zellner A, Siow A. Posterior odds ratios for selected regression hypotheses. Trabajos de Estadistica Y de Investigacion Operativa. 1980;31(1):585–603. doi:10.1007/BF02888369
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF02888369&link_type=DOI) 

64. 64.Wang Y, Schroeder GM, Horsley JJ, et al. Temporal stability of intracranial electroencephalographic abnormality maps for localizing epileptogenic tissue. Epilepsia. 2023;64(8):2070–2080. doi:10.1111/epi.17663
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17663&link_type=DOI) 

65. 65.Gliske S V., Irwin ZT, Chestek C, et al. Variability in the location of high frequency oscillations during prolonged intracranial EEG recordings. Nat Commun. 2018;9(1):2155. doi:10.1038/s41467-018-04549-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-04549-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29858570&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

66. 66.Conrad EC, Tomlinson SB, Wong JN, et al. Spatial distribution of interictal spikes fluctuates over time and localizes seizure onset. Brain. 2020;143(2):554–569. doi:10.1093/brain/awz386
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awz386&link_type=DOI) 

67. 67.Castro-Villablanca F, Moeller F, Pujar S, et al. Seizure outcome determinants in children after surgery for single unilateral lesions on magnetic resonance imaging: Role of preoperative ictal and interictal electroencephalography. Epilepsia. 2022;63(12):3168–3179. doi:10.1111/epi.17425
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17425&link_type=DOI) 

68. 68.Gnatkovsky V, De Curtis M, Pastori C, et al. Biomarkers of epileptogenic zone defined by quantified stereo-EEG analysis. Epilepsia. 2014;55(2):296–305. doi:10.1111/epi.12507
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.12507&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24417731&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

69. 69.Wang Y, Yang Y, Li S, et al. Automatic Localization of Seizure Onset Zone Based on Multi-Epileptogenic Biomarkers Analysis of Single-Contact from Interictal SEEG. Bioengineering. 2022;9(12):769. doi:10.3390/bioengineering9120769
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/bioengineering9120769&link_type=DOI) 

70. 70.Shoeibi A, Ghassemi N, Alizadehsani R, et al. A comprehensive comparison of handcrafted features and convolutional autoencoders for epileptic seizures detection in EEG signals. Expert Syst Appl. 2021;163. doi:10.1016/j.eswa.2020.113788
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eswa.2020.113788&link_type=DOI) 

71. 71.Gascoigne SJ, Waldmann L, Schroeder GM, et al. A library of quantitative markers of seizure severity. Epilepsia. 2023;64(4):1074–1086. doi:10.1111/epi.17525
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.17525&link_type=DOI) 

72. 72.Ben-Jacob E, Doron I, Gazit T, Rephaeli E, Sagher O, Towle VL. Mapping and assessment of epileptogenic foci using frequency-entropy templates. Phys Rev E Stat Nonlin Soft Matter Phys. 2007;76(5). doi:10.1103/PhysRevE.76.051903
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevE.76.051903&link_type=DOI) 

73. 73.Smith EH, Liou JY, Merricks EM, et al. Title: Human interictal epileptiform discharges are bidirectional traveling waves echoing ictal discharges. Elife. 2022;11. doi:10.7554/eLife.73541
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.73541&link_type=DOI) 

74. 74.Diamond JM, Withers CP, Chapeton JI, Rahman S, Inati SK, Zaghloul KA. Interictal discharges in the human brain are travelling waves arising from an epileptogenic source. Brain. 2023;146(5):1903–1915. doi:10.1093/brain/awad015
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awad015&link_type=DOI) 

75. 75.Perucca P, Dubeau F, Gotman J. Widespread EEG changes precede focal seizures. PLoS One. 2013;8(11). doi:10.1371/journal.pone.0080972
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0080972&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24260523&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

76. 76.Pearce A, Wulsin D, Blanco JA, Krieger A, Litt B, Stacey WC. Temporal changes of neocortical high-frequency oscillations in epi-lepsy. J Neurophysiol. 2013;110:1167–1179. doi:10.1152/jn.01009.2012.-High-frequency
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/jn.01009.2012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23761699&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

77. 77.Karoly PJ, Freestone DR, Boston R, et al. Interictal spikes and epileptic seizures: Their relationship and underlying rhythmicity. Brain. 2016;139(4):1066–1078. doi:10.1093/brain/aww019
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/aww019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26912639&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

78. 78.Jirsa VK, Stacey WC, Quilichini PP, Ivanov AI, Bernard C. On the nature of seizure dynamics. Brain. 2014;137(8):2210–2230. doi:10.1093/brain/awu133
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awu133&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24919973&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340051100015&link_type=ISI) 

79. 79.Bernabei JM, Arnold TC, Shah P, et al. Electrocorticography and stereo EEG provide distinct measures of brain connectivity: implications for network models. Brain Commun. 2021;3(3). doi:10.1093/braincomms/fcab156
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/braincomms/fcab156&link_type=DOI) 

80. 80.Kundu S. AI in medicine must be explainable. Nat Med. 2021;27(8):1328–1328. doi:10.1038/s41591-021-01461-z
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-021-01461-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34326551&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

81. 81.Saminu S, Xu G, Shuai Z, et al. Application of Deep Learning and WT-SST in Localization of Epileptogenic Zone Using Epileptic EEG Signals. Applied Sciences. 2022;12(10):4879. doi:10.3390/app12104879
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/app12104879&link_type=DOI) 

82. 82.D.K,. T. B.G P, Xiong F. Auto-detection of epileptic seizure events using deep neural network with different feature scaling techniques. Pattern Recognit Lett. 2019;128:544–550. doi:10.1016/j.patrec.2019.10.029
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.patrec.2019.10.029&link_type=DOI) 

83. 83.Boonyakitanont P, Lek-uthai A, Chomtho K, Songsiri J. A review of feature extraction and performance evaluation in epileptic seizure detection using EEG. Biomed Signal Process Control. 2020;57:1–28. doi:10.1016/j.bspc.2019.101702
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bspc.2019.101702&link_type=DOI) 

84. 84.Karimi-Rouzbahani H, Bagheri N, Ebrahimpour R. Average activity, but not variability, is the dominant factor in the representation of object categories in the brain. Neuroscience. 2017;346. doi:10.1016/j.neuroscience.2017.01.002
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroscience.2017.01.002&link_type=DOI) 

85. 85.van Mierlo P, Vorderwülbecke BJ, Staljanssens W, Seeck M, Vulliémoz S. Ictal EEG source localization in focal epilepsy: Review and future perspectives. Clinical Neurophysiology. 2020;131(11):2600–2616. doi:10.1016/j.clinph.2020.08.001
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.clinph.2020.08.001&link_type=DOI) 

86. 86.Plummer C, Vogrin SJ, Woods WP, Murphy MA, Cook MJ, Liley DTJ. Interictal and ictal source localization for epilepsy surgery using high-density EEG with MEG: A prospective long-term study. Brain. 2019;142(4):932–951. doi:10.1093/brain/awz015
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awz015&link_type=DOI) 

87. 87.Wendling F, Bartolomei F, Mina F, Huneau C, Benquet P. Interictal spikes, fast ripples and seizures in partial epilepsies - combining multi-level computational models with experimental data. European Journal of Neuroscience. 2012;36(2):2164–2177. doi:10.1111/j.1460-9568.2012.08039.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1460-9568.2012.08039.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22805062&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

88. 88.Gentiletti D, de Curtis M, Gnatkovsky V, Suffczynski P. Focal seizures are organized by feedback between neural activity and ion concentration changes. Elife. 2022;11. doi:10.7554/eLife.68541
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.68541&link_type=DOI) 

89. 89.Terrone G, Pauletti A, Pascente R, Vezzani A. Preventing epileptogenesis: A realistic goal? Pharmacol Res. 2016;110:96–100. doi:10.1016/j.phrs.2016.05.009
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.phrs.2016.05.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27173399&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F12%2F2023.08.29.23294708.atom) 

90. 90.Piper RJ, Richardson RM, Worrell G, et al. Towards network-guided neuromodulation for epilepsy. Brain. 2022;145(10):3347–3362. doi:10.1093/brain/awac234
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awac234&link_type=DOI) 

91. 91.Karoly PJ, Cook MJ, Maturana M, et al. Forecasting cycles of seizure likelihood. Epilepsia. 2020;61(4):776–786. doi:10.1111/epi.16485
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/epi.16485&link_type=DOI) 

92. 92. González Otárula KA, Mikhaeil-Demo Y, Bachman EM, Balaguera P, Schuele S. Automated seizure detection accuracy for ambulatory EEG recordings. Neurology. 2019;92(14):e1540–e1546. doi:10.1212/WNL.0000000000007237
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0000000000007237&link_type=DOI)

 [1]: /embed/graphic-1.gif