 Original article
 Open Access
 Published:
Symbolic time series analysis of electroencephalographic (EEG) epileptic seizure and brain dynamics with eyeopen and eyeclosed subjects during resting states
Journal of Physiological Anthropology volume 36, Article number: 21 (2017)
Abstract
Objective
Epilepsy is a neuronal disorder for which the electrical discharge in the brain is synchronized, abnormal and excessive. To detect the epileptic seizures and to analyse brain activities during different mental states, various methods in nonlinear dynamics have been proposed. This study is an attempt to quantify the complexity of control and epileptic subject with and without seizure as well as to distinguish eyeopen (EO) and eyeclosed (EC) conditions using thresholdbased symbolic entropy.
Methods
The thresholddependent symbolic entropy was applied to distinguish the healthy and epileptic subjects with seizure and seizurefree intervals (i.e. interictal and ictal) as well as to distinguish EO and EC conditions. The original time series data was converted into symbol sequences using quantization level, and word series of symbol sequences was generated using a word length of three or more. Then, normalized corrected Shannon entropy (NCSE) was computed to quantify the complexity. The NCSE values were not following the normal distribution, and the nonparametric Mann–Whitney–Wilcoxon (MWW) test was used to find significant differences among various groups at 0.05 significance level. The values of NCSE were presented in a form of topographic maps to show significant brain regions during EC and EO conditions. The results of the study were compared to those of the multiscale entropy (MSE).
Results
The results indicated that the dynamics of healthy subjects are more complex compared to epileptic subjects (during seizure and seizurefree intervals) in both EO and EC conditions. The comparison of the dynamics of epileptic subjects revealed that seizurefree intervals are more complex than seizure intervals. The dynamics of healthy subjects during EO conditions are more complex compared to those during EC conditions. Further, the results clearly demonstrated that thresholddependent symbolic entropy outperform MSE in distinguishing different physiological and pathological conditions.
Conclusion
The threshold symbolic entropy has provided improved accuracy in quantifying the dynamics of healthy and epileptic subjects during EC an EO conditions for each electrode compared to the MSE.
Introduction
Epilepsy is monitored using electroencephalography signals (EEGs) and epileptic seizure detection algorithms [1]. About 50 million people across the globe are suffering from epilepsy [2], including patients of all age groups from newborns [3] to senior adults [4]. The behaviour, cognitive functions and mood of epileptic patients affect the epileptic activities within the brain. Moreover, patient’s psychological and social adaptation can be modified according to their epileptic experience. Due to the interactions between these aspects, people who suffered from epilepsy may face many psychological and cultural problems [1, 5,6,7,8,9].
Various techniques have been developed for understanding the mechanism of epileptic disorders and epileptic seizure detection [10,11,12] based on time–frequency decomposition [13] and waveletbased spare functional linear model [14]. Alkan and Kiymik [15] used AR and Welch methods for detection of epileptic seizure and by examining the power spectra and power spectral densities. Buteneers et al. [16] used reservoir computing (RC) to detect epileptic seizure on intercranial rate data. Bogaarts et al. [17] employed a support vector machine (SVM) to classify and optimize neonatal EEG seizure detection by first filtering the EEG features and data using Kalman filter (KF) in order to increase the temporal precision. Fergus et al. [18] used an advanced machine learning approach on generalized epileptic seizure detection of CHBMIT database. Recently, researchers have employed DWTbased ApEn and artificial neural network [19], probability distribution based on equal frequency discretization [20], and best basis wavelet functions in temporal lobe mimetic [21] for detection and analysis of EEG epileptic seizures.
The nonlinear dynamics in normal restingstate EEG are primarily concerned with studying the dynamics in normal EEG particularly in alpha rhythm. Generally, alpha activity in EEG is dominant in normal individuals during an eyeclosed resting condition and suppresses as visual stimulation [22,23,24,25]. Alpha activity decreased in occipital regions and also in posterior regions when the individuals opened their eyes [26,27,28,29,30]. These studies suggest that alpha desynchronization is reflecting the increased visual system functioning due to visual stimulation being mediated by the reticular activating system [23, 25]. Alpha rhythm biofeedback has gotten some successes in humans for seizure suppression and for depression treatment [26]. Aich [31] examined the relationship between epilepsy, seizure activity and alpha activity in EEG. The findings revealed that the absence of alpha rhythm activity is significantly correlated with the presence of seizure activity. Sherman et al. [32] tracked thalamocortical association during pentylenetetrazol (PTZ) seizures in rats with and without prior treatment with antiepileptic drug (AED) that raises the threshold for seizure. The findings revealed the formation of nonlinearities at specific frequencies in the recorded EEG due to the increase in low alpha wave harmonic AED therapy.
Previous research evidences reported that physiological systems operate across multiple temporal and spatial scales [33,34,35,36]. Traditional entropybased methods such as approximate entropy [37] and sample entropy [38] which are singlescalebased models quantify the complexity of physiological systems by computing the repetitive patterns of a time series to quantify the degree of its regularity. To address this issue, multiscale entropy (MSE) was proposed [33] which incorporates multiple time scale accurate entropy estimate. Kang et al. [39] used MSE to track and differentiate dynamical changes in complexity of each subband under hypothermia and normothermia conditions. Park et al. [40] demonstrated that dynamical alternations owing to Alzheimer disease can be effectively described by MSE curves. In a study conducted by Ouyang et al. [41], dynamical characteristics during seizurefree, preseizure and seizure states in epileptic rates were investigated. The results revealed that the dynamics of seizurefree state are more complex than preseizure and during seizure state dynamics.
The methods derived from symbolic dynamics provide another framework to deal with the underlying multiscale character of physiological systems and could improve the accuracy in quantifying the dynamics of health and epileptic subjects. Aziz and Arif [35] proposed thresholddependent symbolic entropy for distinguishing control (healthy) and neurodegenerative disease patients (ALS, Huntington and Parkinson). The symbolic time series was also used to analyse the temporal gait dynamics of human locomotor system during constrained and metronomically walking protocols and observed that the output of locomotor system during unconstrained normal walking are more complex than slow, fast or metronomically paced walking [42].
The presented study was conducted to determine the dynamics of healthy and epileptic seizure subjects as well as to discriminate eyeclosed (EC) and eyeopen (EO) states during resting conditions, using thresholddependent symbolic entropy. The thresholdbased symbolic entropy was employed to distinguish EEG healthy (both with the eyes open and closed) and EEG epileptic seizure subjects (epileptic subjects, i.e. interictal during seizurefree interval; epileptic subjects, i.e. ictal during seizure interval). The normalized corrected Shannon entropy (NCSE) [31] was used to quantify dynamics, the complexity of healthy (set O and set Z) subjects and epileptic (set F, set N and set S) subjects at different threshold values as well as to discriminate EC from EO during resting states at various thresholds. The results of the study were compared to those of the MSE [33].
Methods
EEG datasets
The performance of thresholddependent symbolic entropy and MSE was evaluated using two different datasets. The first dataset was taken from a publicly available database (http://epileptologiebonn.de/cms/front_content.php?idcat=193&lang=3&changelang=3), made available by the Department of Epileptology, Bonn University, and its detailed description is provided by Andrzejak et al. [43]. The EEG data was recorded with 128channel amplifier system, using an average common reference. The data were digitized at a sampling frequency of 173.61 Hz using a 12b resolution. The spectral bandwidth of the data acquisition system varied from 0.5 to 85 Hz. The whole EEG data comprised of five sets (denoted as Z, O, N, F and S), each containing 23.6sduration 100 singlechannel EEG segments. The sets O and Z were recorded from five healthy volunteers during awake state with the eyes closed (set O) and the eyes open (set Z) using a standardized electrode placement scheme. Sets N, F and S originated from an EEG archive of presurgical diagnosis. Segments in set F were acquired from the epileptogenic zone and those in set N were acquired from the hippocampal formation of the opposite hemisphere of the brain. The segments in sets N and F contained EEG recordings acquired during seizurefree intervals, and segments in set S contained seizure activity.
The second dataset comprised of EEG recordings acquired from 16 healthy (6 males and 10 females) subjects during EC and EO conditions in resting state, and this data has been used in various studies [44,45,46,47]. The EEG data is publically available at (https://www.nbtwiki.net/) and has been used as demo dataset for Neurophysiological Biomarker Toolbox. Two versions, raw and clean EEGs, are available, and in the present study, the clean version of data is used. All volunteers refrained from consuming caffeine, alcohol, and nicotine for at least 6 h before the experiment. The EEG was recorded outside the MRI scanner using two BrainAmp amplifiers (Brain Products, Munich, Germany) connected to the EEG monitor via optical fibres [47]. The EEG was recorded in DC mode at a resolution of 0.1 V, with a 250Hz lowpass filter [47]. The standard 10–20 positions of the electrodes were used as reported by Halder et al. [44], Maurer et al. [45], and Brem et al. [46]. A 2.5min EC restingstate EEG recordings at a resolution of 3.3 mV with a sampling frequency of 5000 Hz were acquired. These recording were downsampled to 500 Hz offline [47]. The subjects were instructed to relax and stay on the scanner while remaining calm and awake. For both resting conditions, i.e. EC and EO, the subjects were additionally instructed accordingly. The EEG datasets with EC and EO were preprocessed offline using Brain Vision Analyzer software (version 1.05, brain Products, Munich, Germany) and the Neurophysiological Biomarker Toolbox (http://www.nbtwiki.net/). Major large muscles or motion artefacts (typical duration 2.5 s, range 0.4–31.6 s) were removed from the EEG before ICA. Moreover, data were digitally highpass filtered for ICA cleaning (0.5Hz finite impulse response filter, 4s Hanning window).
Thresholddependent symbolic entropy
The symbolic time analysis involves the transformation of the original time series into a series of discretized symbols that are processed to extract useful information about the state of the system generating the process. There have been many recent applications of symbolic analysis for biological systems. The NCSE is calculated for discretized symbols. This method requires the following steps: while using the symbolic time series analysis, most of the microscopic detail of dynamics may be lost; however, the temporal dynamics are embedded in the structure according to length—L distribution [20,49,, 48–50]. The symbolic time analysis is efficient in terms of computer time and storage and noise robustness.
Symbolization process
The EEG time series of both epileptic seizure and eyeopen and eyeclosed conditions in resting states are transformed into symbol sequence and word sequence by following the procedure detailed below [35, 51].
Step I:
Consider an EEG time series T = {T _{ i }, i = 1, 2, … … …. N}. The time series is transformed into symbol sequence as t ^{Ψ} = {t _{ i } ^{Ψ} , i = 1, 2, …. N} having fixed values of quantization level (psi Ψ) labelled with number zero to Ψ − 1. Each sample of time series is normalized, and we obtained a time series t = {t _{ i }, i = 1, 2, … … …. N}.
We used the quantization level 2 (symbol 0 and 1) according to the following symbolization criteria:
where θ is the threshold and is computed as the mean of the EEG time series.
Consider, for example, EEG time series of an EC resting state taking only 10 data points to illustrate the symbolization process such as the following:
T = {0.5378, 2.0403, 5.1684, 7.8292, 7.3310, 3.4433, 0.4669, −0.3348, −0.7980, −2.2799}, where threshold (θ) is \( \overline{t}\ \left(\mathrm{mean}\ \mathrm{of}\ \mathrm{time}\ \mathrm{series}\right)=2.3404 \)
Symbol series using Eq. (1) was formed as {0 0 1 1 1 0 0 0 0 0 0 } as shown in Fig. 1a, b.
Step II:
The word sequence is formed by dividing the symbol sequence of a word length L of three or more words. To construct the symbol sequence (words) from symbol series in step I, the group of symbols are collected together in temporal order by moving the symbol series one step at a time, where each step reveals a new sequence using:
The following word series is generated from the symbol series generated using Eq. (2) taking word length L = 3.
Step III:
Each possible sequence generated by step II is represented in terms of a unique identifier; this new time series is referred to as a code series [13] using the following procedure:
Finally, the code series are as follows:
To unfold the dynamics of time series, word length plays a vital role as described by [32]. They also pointed out that the optimal word length correlates with the number of data points in the time series and by increasing the word length, NCSE decreases. We also observed that NCSE decreases by increasing the word length [35], and the optimal results are obtained using world length L = 3 in this study.
Symbolic sequence analysis
Data symbolization can be done with several methods (e.g. heart rate variability using two bins), the value ‘1’ or ‘0’ is assigned to each bin according to its occurrence [50] as described in the below example [52, 53]. This is usually done by replacing the original data with relative symbols according to the threshold criterion as described in [35]. Moreover, Cysarz et al. [51] quantified the dynamics in a heart rate using different symbolic dynamics binning approaches. Symbolic dynamics have been widely used to quantify the nonlinear dynamics of physiological time series. There are several approaches to symbolize the time series, e.g. one approach is to symbolize the time series using deviation from an average heart rate as used by Kurths et al. [54] and Wessel et al. [55]. After obtaining the symbol patterns, further complexity can be analysed using Shannon or Renyi entropy [35, 51]. The nonlinear dynamics have also been investigated using time–frequency methods such as correlation, coherence, and interaction and synchronization measures. The basic symbolization approach is also extended to the timedelayed coupling methods in HRV and EEG as described in [56].
Symbol sequence analysis depends on quantitative measures of symbol sequence frequencies. Classical theory (containing examples from Euclidean norm and chisquare statistics) and information theory (containing examples from Shannon and generalized Renyi entropy) are the two general divisions of symbolic sequence analysis. In this study, we have followed the information theoretic approach using Shannon entropy.

1.
Shannon entropy
Shannon entropy of order L and quantization level Ψ can be defined as
$$ \mathrm{S}\mathrm{E}\left( L,\Psi \right)={\displaystyle \sum } p\left({t}_L^{\Psi}\right) \cdot { \log}_2 p\left({t}_L^{\Psi}\right) $$(4)where p(t _{ L } ^{Ψ} ) is the probability of the time series t _{ L } ^{Ψ} being the pattern and Ψ is the quantization level labelled with number from zero to Ψ − 1. Eguia et al. [57] reported that systematic error or bias and random errors affect the estimates computed by using Shannon entropy and thus reported a leading correction of entropy.

2.
Corrected Shannon entropy
The corrected Shannon entropy is defined as
$$ \mathrm{C}\mathrm{S}\mathrm{E}\left( L,\Psi \right)=\mathrm{S}\mathrm{E}\left( L,\Psi \right)+\frac{C_R1}{2 M \ln 2} $$(5)where M is the total number of words and C _{ R } is the number of words which occurred among the possible outcomes of words. For certain quantization level Ψ and word length L, the value of CSE is the maximum, and when all words M occur uniformly, thus, max CSE will be:
$$ {\mathrm{CSE}}^{\max}\left( L,\Psi \right)={ \log}_2\left(\frac{1}{M}\right) + \frac{M1}{2 M \ln 2} $$(6)From Eq. (6), it can be observed clearly that the maximum value of CSE cannot be the same for two different word lengths and when L increases, M also increases thereby increasing the maximum value of CSE. Likewise, it is not possible to compare two values of CSE having different word lengths at the same threshold and quantization level. To overcome this problem, NCSE is proposed and used by Aziz and Arif [35].

3.
Normalized corrected Shannon entropy
Thus, NCSE can be defined as
$$ \mathrm{NCSE}\left( L,\Psi \right)=\frac{\mathrm{CSE}\left( L,\Psi \right)}{{\mathrm{CSE}}^{\max}\left( L,\Psi \right)}=\frac{\mathrm{CSE}\left( L,\Psi \right)}{{ \log}_2\left(\frac{1}{M}\right) + \frac{M1}{2 M \ln 2}} $$(7)Using Eq. (7) for any quantization level (Ψ) and word length (L), the value of NCSE will always vary from 0 to 1.
Statistical analysis
The NCSE values were not following the normal distribution, and Mann–Whitney–Wilcoxon (MWW) test was used to find significant differences among various groups at 0.05 significance level. The MWW is a nonparametric analogue of t test, which does not assume any properties regarding the distribution of the independent variable. All the codes for the topographic maps and statistical and symbolic time series analysis were developed in Matlab.
Results
NCSE was employed to distinguish the EEG healthy (set O and set Z) with epileptic subjects (interictal and ictal intervals) and eyeclosed resting (ECR) and eyeopen resting (EOR) states. The mean ± SD against each group are shown in Tables 1 and 2. The results show that the NCSE value of the healthy EEG subjects is higher than that of the epileptic subjects (both with seizure and seizurefree intervals) as shown in Table 1 as well as the NCSE value of the EO condition than that of the EC condition during resting states in Table 2 at a smaller threshold range. The results also reveal that an epileptic seizurefree interval, i.e. interictal (F and N), has higher complexity than the epileptic with seizure (set S) at the same smaller threshold values. The maximum NCSE value for the healthy subjects (with eyeopen set O) was found at a threshold of 30 ms, whereas the maximum threshold for the healthy subject (with eyeclosed set Z), epileptic seizurefree intervals (interictal intervals) was found at a threshold of 25 ms; however, the NCSE value for epileptic seizure (ictal interval) increases until a threshold of 80 ms.
The NCSE values are computed from a threshold of 1 to 30 ms in which the significance results are obtained overall from both eyeclosed and eyeopen conditions during resting states for 19 channels C3, C4, Cz, F3, F4, F7, F8, Fp1, Fp2, Fz, O1, O2, P3, P4, Pz, T7, T8, P7 and P8 according to 10–20 system. MSE was also estimated with m = 1 and r = 0.25 times the SD of the original time series. The results are summarized in Table 2. The values of EO for NCSE and MSE at all electrodes are higher than those of EC. These electrodes are chosen according to 10–20 international standard system and are representing the overall 129 other electrodes as usually used in the research to investigate the nonlinear dynamics such as correlation, coherence, complexity or any spectral measure.
Table 2 reflects the median NCSE values of EC and EO during resting states with maximum separation, i.e. significance values in the threshold range 1 to 30 using nonparametric Wilcoxon rank sum test, and results are compared to MSE with a time scale from 1 to 20. The maximum medium NCSE values with the highest separation (minimum P value) significance value for EC and EO during resting states are reflected in Table 2 with a corresponding threshold also reflected where these highest significance levels are obtained. The column Significance (θ, 1–30 ms) in Table 2 reflects the threshold ranges in which the statistically significance was obtained from threshold from 1 to 30 ms. The results show that 18 electrodes out of 19 selected standard electrodes exhibit the significant results in the threshold range 1 to 30, and only Fp2 electrodes did not show any significant result on any threshold. From the results, it is also evident that the highest significance level at all electrodes was obtained at thresholds 1.2 and 1.4 in most of the cases and at thresholds 1 and 6 in few cases. The results also reflect that NCSE gives significant results at frontal electrodes in threshold range 1–2 and 4–10 in most of the cases whereas the occipital and parietal region electrodes exhibit significance in lower range 1–2 and upper range 4–30 thresholds. The central electrodes C3, C4 and Cz give significant results at smaller threshold values in the range 1–2 only whereas central, occipital and parietal probes exhibit these significances even at higher thresholds from 4 to 15/30. Overall, NCSE gives the highest significant results in all the regions such as frontal, central, parietal, occipital and temporal regions using 19 selected electrodes as per standard in the threshold range 1–30. The maximum median values of the NCSE values and corresponding significance values for both EO and EC are reflected in the form of topographic maps in Fig. 2a–d at thresholds 1.2 and 1.4 and a world length of 3. The MSE values for EC and EO are also quantified for scale 1 to 20, and the corresponding P values are computed using Wilcoxon rank sum test. Table 2 also reflect the maximum median values for each 19 electrodes using MSE, and the minimum P values are reflected from scale 1 to 20.
The performance of thresholdbased symbolic entropy was also evaluated for distinguishing EC from EO during restingstate conditions on 14 selected electrode, and the results were compared to MSE. The results revealed that symbolic time series gives higher significant results than MSE at all electrodes except Fp2 which does not show any significant result as depicted in Table 2 and Figs. 2 and 3 in the form of topographic maps, whereas MSE shows the significance results in few electrodes only. The significant results using symbolic time series analysis were obtained at frontal electrodes F4 (P value 0.00043) and F7 (P value 0.00049); occipital region O1 (P value 0.00013) and O2 (P value 0.00043); and parietal region P7 (P value 0.00032) and P8 (P value 0.00085). While high significant results were obtained at electrodes C3 (P value 0.00704), F3 (P value 0.00143), F8 (P value 0.007), Fz (P value 0.00273), P3 and P4 (P value 0.0043) and T8 (P value 0.00983). Moreover, the electrodes C4, Cz, Pz and T8 gives only significant results. In comparison, the MSE gives only significant result to distinguish the EC from EO conditions during resting states at few electrodes such as C3, Cz, F3, F4, F7, O1, O2, P4, Pz, T8 and P8 only, and the other electrodes do not show any significant result.
In Table 2, we have compared the symbolic time series and MSE for which the minimum P values are reflected against each 19 electrodes. Symbolic entropy showed statistically significant results for 18 electrodes except Fp2 whereas MSE has not showed significant results at electrodes C4, F8, Fp1, Fp2, Fz, P3, T7 and P7. Moreover, the results showed that P values using symbolic time series analysis are lower than those using the MSE analysis method which is consistent with previous studies [31, 38, 47].
The NCSE strength is reflected on topographic maps which show the regions where the maximum NCSE values are obtained in the form of colourmaps. The colour scales reflect the strength of the NCSE values on each figure. Table 2 reflects only the maximum median values whereas the topographic maps reflect the median NCSE values at thresholds 1.2 and 1.4 where the maximum separation was obtained. The topographic maps also help to understand visually where the maximum strength of NCSE is obtained in the brain region. Thus, entropic measures can be best reflected in the form of topographic maps to investigate the overall underlying dynamics in a system.
Thus from the topographic maps in Fig. 3 (e), it can be seen that at threshold 1.2 and word length 3, the electrodes F7, F3 and P8 exhibits strictly significant results; electrodes Fz, F4, F8, T7, T8, C3, C4, Cz, P3, P4, P7, O1 and O2 shows very significant results where electrodes Fp1 and Pz only shows just significant results. Figure 3 (f) shows the statistically significant results in the form of topographic maps at threshold 1.4 and word length 3. In this case the electrodes F7, F3, F4, F8, T7, P7, O1, P4, and P8 exhibit strictly statistically significant results. The electrodes Fz, C3, Cz, T8, P3 and O2 also shows very significant results and Fp1, C4 and Pz shows just significant results.
In the topographic maps in Fig. 4e, the overall MSE values of EC are found greater than those of EO and only Cz electrodes gives very significant result, whereas the electrodes F7, F3, F4, Fp2, C3, P3 and Pz were just found to be statistically significant. Moreover, the electrodes Fp1, Fz, P4, P8 and T8 are found almost significant, and T7, C4, P7, O1 and O2 exhibit no significance at all. In the topographic maps in Fig. 4f, again, only Cz gives very significant statistical results. Electrodes F7, F3 and F4 exhibit just significant results. and C3, Fp2, F8 and T8 are almost found significant. A larger number of electrodes as Fp1, Fz, T7, C4, P7, P8, P3, P4, Pz, O1 and O2 did not show any significance at all.
Figure 4a, b shows the comparisons of the NCSE values of the healthy subjects (with eyeopen set O and eyeclose set Z) with epileptic seizure subject (ictal interval) with an error bar. From the results, it is evident that the healthy subjects in both cases exhibit higher complexity than the epileptic subjects at certain threshold values. Figure 5a shows the comparison of the healthy subjects (eyeclosed condition) to the epileptic seizurefree (interictal interval) subject and also depicts the higher NCSE values of the healthy subject than those of the interictal interval at a certain threshold. Figures 4 and 5 and Tables 1 also reflect that the epileptic seizurefree (interictal) intervals (both focal and nonfocal) have also higher complexity values than the epileptic seizure (ictal) subject. In summary, it is clearly observed that the healthy subjects in both conditions, i.e. eyeopen and eyeclosed condition, exhibits higher complexity than the epileptic seizure and seizurefree intervals, whereas the seizurefree groups show higher complexity than the epileptic seizure (ict intervals) subjects at a certain threshold value. Figure 5b also shows the comparisons of the NCSE values of eyeclosed and eyeopen subjects during resting states at occipital electrode (O1). From Fig. 5b, it can be clearly depicted that the EO subjects shows greater NCSE values than the EC during resting states at certain smaller thresholds. Thus, EO is more complex than EC at these smaller threshold ranges.
Figure 4a, b shows the mean NCSE distribution of the healthy subjects with the eyes open (set O) vs the epileptic with seizures (ictal interval), the healthy with the eyes closed (set Z) vs the epileptic with seizures (ictal interval), the healthy with the eyes closed (set Z) and epileptic seizurefree interval (interictal set focal) and EC and EO during resting states. In all the cases, the symbols are used to denote the means, and bars to denote the standard error. The standard error is calculated using the formula Standard error = SD/√ n where n is the number of subjects in each case and SD denote the standard deviation. The most significant differences are seen at threshold 30 ms for set O with mean ± SD = 0.84 ± 0.00, sets Z, F and N at a threshold of 25 ms with mean ± SD of 0.83 ± 0.00, 0.68 ± 0.01 and 0.69 ± 0.01, respectively, whereas the mean ± SD for the epileptic seizure subject increased till threshold of 80 ms.
Discussions and conclusions
The physiological system is composed of various subsystems, which are controlling the regularity mechanism of that system [33,34,35,36]. If all the subsystems are functional and working properly, the controlling mechanism of the overall system will be appropriate and the system is healthy. The healthy systems evolve with time and their adaptive capability is higher, resulting in higher complexity. The alternations in the structural components and/or decreased functional capability of the subsystem cause dysfunction in the regularity mechanism of the overall system, which results in the loss of complexity [33,34,35,36]. Various entropybased methods have been proposed by researchers [33,34,35,36,37,38], each having its advantages and disadvantages in quantifying the dynamics of physiological system in healthy and disease. In this study, we have used thresholdbased symbolic entropy to distinguish epileptic and healthy subjects as well as eyeopen and eyeclosed conditions. The results of the study were compared to those of the MSE [33], which has been used in several neurological studies to differentiate hypothermia and normothermia conditions [39], describe dynamical alternation owing to Alzheimer disease and dynamics of seizurefree, preseizure and during seizure states [41].
The NCSE values of the healthy subjects were higher than the epileptic subjects (including both ictal and interictal intervals) at wide range threshold values. Higher NCSE values for the healthy subjects at a wide range of thresholds manifest that probability of distribution of patterns (symbols sequence) is more uniform (almost equal chance of occurrence of all possible patterns). The higher probability of occurrence of all patterns (higher NCSE) reflects that the complexity of the healthy subjects higher. The loss of structural components and/or decreased functional capability of the subsystem due to ageing or disease perturbs the normal functioning of the overall system [36]. This perturbation may affect the probability of occurrence of all possible patterns which results in nonuniform distribution of pattern [35]. The decrease in uniformity of patterns due to dominancy of some patterns and/or nonoccurrence of some pattern will result in low NCSE values. Thus, reduction in the NCSE values for epileptic subjects reflects the loss complexity, which is in line with the hypothesis ‘loss of complexity is a generic feature of pathological systems’ [33,34,35,36]. The MSE results also showed that the dynamics of healthy subjects are more complex compared to epileptic subjects; however, the results of NCSE are statistically more significant in distinguishing healthy and epileptic subjects. The one more advantage of using symbolic entropy measure is that it is robust short time series data, whereas MSE provides dynamically incorrect information at large temporal scales and may induce undefined entropy estimates.
The human brain displays five different waves to cope different situations, and each have a specific frequency band. The δ (0.5−4 Hz) wave activity is predominant in infants and occurs during deep sleep. θ (4−7.5 Hz) is found in normal awake adults but more prevalent during drowsiness. α (7.5−14 Hz) activity is attenuated by attention and found in the posterior head region in individuals who are relaxed, are awake and relaxed, and have their eyes closed. β (14−22 Hz) activity is associated with active attention and thinking by focusing outside the world or solving special types of problems usually the waking rhythm whereas γ (22−80 Hz) activities have very low amplitude and occur often. Alpha activity is dominant in normal individuals during resting states and is suppressed by visual stimulation [22,23,24,25,26]. Any brain wave either overproduced and/or underproduced can cause a problem. In this study, we also investigated changes in the dynamics of EO and EC conditions during resting states. The results indicated EO dynamics are more complex compared to EC dynamics. The resting state involves the different brain oscillations and particularly the alpha rhythm. The EO condition reflects increased visual system activity due to visual stimulation, probably mediated by the reticular activating system [27, 29, 58]. The increased visual activity may involve more structural components and hence escalate the coupling between additional functional resulting in more complex dynamics.
Limitations of study and future recommendations
The research reported in this manuscript is focused on quantifying the dynamics of EEG signals during ictal and interictal intervals and to distinguish eyeclosed condition from eyeopen condition during resting states. The study revealed very interesting results; however, there are several limitations of this study such as the number of subjects is small and the lack of ageing and gender base analysis. Furthermore, dataset were taken from publicly available databases, and clinical profile of data was not available. Future studies can be designed to quantify the dynamics for the larger groups based on gender and age as well as to assess the suitability threshold symbolic for specific type of patients with epilepsy.
References
 1.
Oğulata SN, Şahin C, Erol R. Neural networkbased computeraided diagnosis in classification of primary generalized epilepsy by EEG signals. J Med Syst. 2009;33(2):107–12.
 2.
WHO (2009) http://www.who.int/mediacentre/factsheets/fs999/en/. Accessed 26 May 2015.
 3.
Rankine L, Mesbah M, Boashash B. A matching pursuitbased signal complexity measure for the analysis of newborn EEG. Med Biol Eng Comput. 2007;45(3):251–60.
 4.
Jenssen S, Schere D. Treatment and management of epilepsy in the elderly demented patient. Am J Alzheimers Dis Other Demen. 2010;25(1):18–26.
 5.
Aslan K, Bozdemir H, Şahin C, Oğulata SN, Erol R. A radial basis function neural network model for classification of epilepsy using EEG signals. J Med Syst. 2008;32(5):403–8.
 6.
Betts TA, Merskey H, Pond DA. Psychiatry. In: Laidlaw J, Richens A, editors. Textbook of epilepsy. Edinburgh: Churchill Livingstone; 1976. pp. 145–84. ISBN 0443014264.
 7.
Buchanan N. Understanding epilepsy: what it is and how it can affect your life. Sydney: Simon & Schuster; 2002.
 8.
Bishop M, Allen CA. The impact of epilepsy on quality of life: a qualitative analysis. Epilepsy Behav. 2003;4(3):226–33.
 9.
Gumnit RJ. Living well with epilepsy. New York: Demos Medical Publishing; 1997. pp. 173–81.
 10.
Srinivasan V, Eswaran C, Sriraam AN. Artificial neural network based epileptic detection using timedomain and frequencydomain features. J Med Syst. 2005;29(6):647–60.
 11.
Tzallas AT, Tsipouras MG, Fotiadis DI. Automatic seizure detection based on timefrequency analysis and artificial neural networks. Comput Intell Neurosci. 2007;2007.
 12.
GuerreroMosquera C, Trigueros AM, Franco JI, NaviaVázquez Á. New feature extraction approach for epileptic EEG signal detection using timefrequency distributions. Med Biol Eng Comput. 2010;48(4):321–30.
 13.
Nason GP, Von Sachs R, Kroisandt G. Wavelet processes and adaptive estimation of the evolutionary wavelet spectrum. J R Stat Soc Series B Stat Methodol. 2000;62(2):271–92.
 14.
Xie S, Krishnan S. Waveletbased sparse functional linear model with applications to EEGs seizure detection and epilepsy diagnosis. Med Biol Eng Comput. 2013;51(12):49–60.
 15.
Alkan A, Kiymik MK. Comparison of AR and Welch methods in epileptic seizure detection. J Med Syst. 2006;30(6):413–9.
 16.
Buteneers P, Schrauwen B, Verstraeten D, Stroobandt D. Realtime epileptic seizure detection on intracranial rat data using reservoir computing. In: International Conference on Neural Information Processing Nov. 25–28, 2008, Auckland, New Zealand. Germany: Springer Berlin Heidelberg; 2008. pp. 56–63.
 17.
Bogaarts JG, Gommer ED, Hilkman DM, van KranenMastenbroek VH, Reulen JP. EEG feature preprocessing for neonatal epileptic seizure detection. Ann Biomed Eng. 2014;42(11):2360–8.
 18.
Fergus P, Hignett D, Hussain AJ, AlJumeily D. An advanced machine learning approach to generalised epileptic seizure detection. In: International Conference on Intelligent Computing Aug. 3–6, 2014, Taiyuan, China. Switzerland: Springer International Publishing; 2014. pp. 112–18.
 19.
Kumar Y, Dewal ML, Anand RS. Epileptic seizures detection in EEG using DWTbased ApEn and artificial neural network. Signal Image Video P. 2014;8(7):1323–34.
 20.
Palazzolo JA, Estafanous FG, Murray PA. Entropy measures of heart rate variation in conscious dogs. Am J Physiol Heart Circ Physiol. 1998;274(4):H1099–105.
 21.
Abibullaev B, Kim MS, Seo HD. Seizure detection in temporal lobe epileptic EEGs using the best basis wavelet functions. J Med Syst. 2010;34(4):755–65.
 22.
Berger H. Über das elektrenkephalogramm des menschen. Eur Arch Psychiatry Clin Neurosci. 1933;98(1):231–54.
 23.
Adrian ED. The bebgee ehythm: potential changes feom the occipital lobes in man. Physiol Lab Camb. 1934;57(4):357–85.
 24.
Jasper HH, Solomon P, Bradley C. Electroencephalographic analyses of behavior problem children. Am J Psychiatry. 1938;95(3):641–58.
 25.
Smith HF. An empirical law describing heterogeneity in the yields of agricultural crops. J Agr Sc. 1938;28(01):1–23.
 26.
Chapman RM, Armington JC, Bragdon HR. A quantitative survey of kappa and alpha EEG activity. Electroencephalogr Clin Neurophysiol. 1962;14(6):858–68.
 27.
Volavka J, Matoušek M, Roubíček J. Mental arithmetic and eye opening. An EEG frequency analysis and GSR study. Electroencephalogr Clin Neurophysiol. 1967;22(2):174–6.
 28.
Glass A, Kwiatkowski AW. Power spectral density changes in the EEG during mental arithmetic and eyeopening. Psychol Forsch. 1970;33(2):85–99.
 29.
Gale A, Coles M, Boyd E. Variation in visual input and the occipital EEG: II. Psychon Sci. 1971;23(1):99–100.
 30.
Kraft U. Train your brain. Sci Am Mind. 2006;1(17):58–63.
 31.
Aich TK. Absent posterior alpha rhythm: an indirect indicator of seizure disorder? Indian J Psychiatry. 2014;56(1):61.
 32.
Sherman D, Zhang N, Garg S, Thakor NV, Mirski MA, White MA, Hinich MJ. Detection of nonlinear interactions of EEG alpha waves in the brain by a new coherence measure and its application to epilepsy and antiepileptic drug therapy. Int J Neural Syst. 2011;21(02):115–26.
 33.
Costa M, Goldberger AL, Peng CK. Multiscale entropy analysis of complex physiologic time series. Phys Rev Lett. 2002;89(6):068102.
 34.
Costa M, Goldberger AL, Peng CK. Multiscale entropy analysis of biological signals. Phys Rev E. 2005;71(2):021906.
 35.
Aziz W, Arif M. Complexity analysis of stride interval time series by threshold dependent symbolic entropy. Eur J Appl Physiol. 2006;98(1):30–40.
 36.
Peng CK, Goldberger AL. Adaptive data analysis of complex fluctuations in physiological time series. Adv Adap Data Analy. 2009;1(1):61–70.
 37.
Pincus SM. Approximate entropy as a measure of system complexity. Proc Natl Acad Sci. 1991;88(6):2297–301.
 38.
Richman JS, Moorman JR. Physiological timeseries analysis using approximate entropy and sample entropy. Am J Physiol Heart Circ Physiol. 2000;278(6):H2039–49.
 39.
Kang X, Jia X, Geocadin RG, Thakor NV, Maybhate A. Multiscale entropy analysis of EEG for assessment of postcardiac arrest neurological recovery under hypothermia in rats. IEEE Trans Biomed Eng. 2009;56(4):1023–31.
 40.
Park JH, Kim S, Kim CH, Cichocki A, Kim K. Multiscale entropy analysis of EEG from patients under different pathological conditions. Fractals. 2007;15(04):399–404.
 41.
Ouyang G, Li X, Dang C, Richards DA. Using recurrence plot for determinism analysis of EEG recordings in genetic absence epilepsy rats. Clin Neurophysiol. 2008;119(8):1747–55.
 42.
Abbasi AQ, Loun WA. Symbolic time series analysis of temporal gait dynamics. J Sig Proc Systs. 2014;74(3):417–22.
 43.
Andrzejak RG, Lehnertz K, Rieke C, Mormann F, David P, Elger CE. Indications of nonlinear deterministic and finite dimensional structures in time series of brain electrical activity: dependence on recording region and brain state. Phys Rev E. 2001;64:061907.
 44.
Halder P, Sterr A, Brem S, Bucher K, Kollias S, Brandeis D. Electrophysiological evidence for cortical plasticity with movement repetition. Eur J Neurosc. 2005;21(8):2271–7.
 45.
Maurer U, Brem S, Bucher K, Kranz F, Benz R, Steinhausen HC, Brandeis D. Impaired tuning of a fast occipitotemporal response for print in dyslexic children learning to read. Brain. 2007;130(12):3200–10.
 46.
Brem S, Bach S, Kucian K, Kujala JV, Guttorm TK, Martin E, Lyytinen H, Brandeis D, Richardson U. Brain sensitivity to print emerges when children learn letter–speech sound correspondences. Proc Natl Acad Sci. 2010;107(17):7939–44.
 47.
Gorman RLO, Poli SS, Brandeis D, Klaver P, Bollman S, Ghisleni C, Luchinger R, Martin E, Shankerranarayanan A, Aslop DC, Michels L. Coupling between resting cerebral perfusion and EEG. Brain Topogr. 2013;26:442–57.
 48.
Lehrman M, Rechester AB, White RB. Symbolic analysis of chaotic signals and turbulent fluctuations. Phys Rev Lett. 1997;78(1):54.
 49.
Porta A, GnecchiRuscone T, Tobaldini E, Guzzetti S, Furlan R, Montano N. Progressive decrease of heart period variability entropybased complexity during graded headup tilt. J Appl Physiol. 2007;103(4):1143–9.
 50.
Tong S, Bezerianos A, Malhotra A, Zhu Y, Thakor N. Parameterized entropy analysis of EEG following hypoxic–ischemic brain injury. Phys Lett A. 2003;314(5):354–61.
 51.
Cysarz D, Porta A, Montano N, Leeuwen PV, Kurths J, Wessel N. Quantifying heart rate dynamics using different approaches of symbolic dynamics. Eur Phys J Spec Top. 2013;222(2):487–500.
 52.
Voss A, Hnatkova K, Wessel N, Kurths J, Sander A, Schirdewan A, Camm AJ, Malik M. Multiparametric analysis of heart rate variability used for risk stratification among survivors of acute myocardial infarction. Pacing Clin Electrophysiol. 1998;21(1):186–96.
 53.
Voss A, Kurths J, Kleiner HJ, Witt A, Wessel N. Improved analysis of heart rate variability by methods of nonlinear dynamics. J Electrocardiol. 1995;28:81–8.
 54.
Kurths J, Voss A, Saparin P, Witt A, Kleiner HJ, Wessel N. Quantitative analysis of heart rate variability. Chaos. 1995;5(1):88–94.
 55.
Wessel N, Malberg H, Bauernschmitt R, Kurths J. Nonlinear methods of cardiovascular physics and their clinical applicability. Int J Bifurcar Chaos. 2007;17(10):3325–71.
 56.
Suhrbier A, Riedl M, Malberg H, Penzel T, Bretthauer G, Kurths J, Wessel N. Cardiovascular regulation during sleep quantified by symbolic coupling traces. Chaos. 2010;20(4):045124.
 57.
Eguia MC, Rabinovich MI, Abarbanel HD. Information transmission and recovery in neural communications channels. Phys Rev E. 2000;62(5):7111.
 58.
Härdle W, Gasser T, Bächer P. EEGresponsiveness to eye opening and closing in mildly retarded children compared to a control group. Biol Psychol. 1984;18(3):185–99.
Acknowledgements
The work is supported as part of Ph.D. for 6 months at the UK by the Higher Education Commission, Government of Pakistan, under the International Research Support Initiative Program (IRSIP), Govt. of Pakistan vide grant no. 18/HEC/HRD/20/43467, PIN: IRSIP 27 Engg 05.
Authors’ contributions
LH had written the Matlab codes to obtain the results, drafted the manuscript, analysed and interpreted the data and reviewed the manuscript. WA, JSA, NH, MR, SS and SZHK had critically reviewed the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Symbolic dynamics
 Physiological complexity
 Normalized corrected Shannon entropy
 Epileptic seizure
 Interictal and ictal seizure states
 Resting states
 Brain activities