# Fast monitoring of epileptic seizures using recurrence time statistics of electroencephalography

^{1}Institute of Complexity Science and Big Data Technology, Guangxi University, Nanning, China^{2}PMB Intelligence LLC, Sunnyvale, CA, USA

Epilepsy is a relatively common brain disorder which may be very debilitating. Currently, determination of epileptic seizures often involves tedious, time-consuming visual inspection of electroencephalography (EEG) data by medical experts. To better monitor seizures and make medications more effective, we propose a recurrence time based approach to characterize brain electrical activity. Recurrence times have a number of distinguished properties that make it very effective for forewarning epileptic seizures as well as studying propagation of seizures: (1) recurrence times amount to periods of periodic signals, (2) recurrence times are closely related to information dimension, Lyapunov exponent, and Kolmogorov entropy of chaotic signals, (3) recurrence times embody Shannon and Renyi entropies of random fields, and (4) recurrence times can readily detect bifurcation-like transitions in dynamical systems. In particular, property (4) dictates that unlike many other non-linear methods, recurrence time method does not require the EEG data be chaotic and/or stationary. Moreover, the method only contains a few parameters that are largely signal-independent, and hence, is very easy to use. The method is also very fast—it is fast enough to on-line process multi-channel EEG data with a typical PC. Therefore, it has the potential to be an excellent candidate for real-time monitoring of epileptic seizures in a clinical setting.

## 1. Introduction

Epilepsy is a relatively common brain disorder which may be very debilitating. It affects approximately 1% of the world population (Jallon, 1997) and three million people in the United States alone. It is characterized by intermittent seizures. During a seizure, the normal activity of the central nervous system is disrupted. The concrete symptoms include abnormal running/bouncing fits, clonus of face and forelimbs, or tonic rearing movement as well as simultaneous occurrence of transient EEG signals such as spikes, spike and slow wave complexes or rhythmic slow wave bursts. Clinical effects may include motor, sensory, affective, cognitive, automatic and physical symptomatology. Although epilepsy can be treated effectively in many instances, severe side effects may result from constant medication. Even worse, some patients may become drug-resistant not long after treatment. To make medications more effective, timely detection of seizure is very important.

In the past several decades, considerable efforts have been made to detect/predict seizures through non-linear analysis of EEGs (Kanz and Schreiber, 1997; Gao et al., 2007). Representative non-linear methods proposed for seizure prediction/detection include approaches based on correlation dimension (Lehnertz and Elger, 1995, 1997; Martinerie et al., 1998; Aschenbrenner-Scheibe et al., 2003), Kolmogorov entropy (van Drongelen et al., 2003), permutation entropy (Cao et al., 2004), short time largest Lyapunov exponent (STLmax) (Iasemidis et al., 1990; Lai et al., 2003), dissimilarity measures (Protopopescu et al., 2001; Quyen et al., 2001), long-range-correlation (Hwa and Ferree, 2002; Gao et al., 2006b, 2007, 2011b; Valencia et al., 2008), power-law sensitivity to initial conditions (Gao et al., 2005b), scale-dependent Lyapunov exponent (SDLE) (Gao et al., 2006a, 2012a,b), and synthesis of linear/non-linear methods by using neural networks (Adeli et al., 2007). Readers interested in “what is epilepsy, where, when, and why (how) do seizures occur?” are referred to the April, 2007 issue of *Journal of Clinical Neurophysiology*.

Note that most of the proposed methods assume that EEG signals are chaotic and stationary. As a result, they tend to have performances that are signal- and patient-dependent due to the noisy and non-stationary nature of the EEG within and across patients. In addition, they are computationally expensive. Consequentially, studies of epilepsy still heavily involve visual inspection of multi-channel EEG signals by medical experts. Visual inspection of long (e.g., tens of hours or days) EEG data is, however, tedious, time-consuming, and in-efficient. Therefore, it is important to develop new non-linear seizure monitoring approaches.

In this paper, we explore recurrence time based analysis of EEG (Gao, 1999, 2001; Gao and Cai, 2000; Gao et al., 2003), with the goal of potentially on-line monitoring the occurrence and propagation of seizures. The method does not assume that the underlying dynamics of EEGs be chaotic or stationary. More importantly, it has been tested to be able to readily detect very subtle changes in signals (Gao, 2001; Gao et al., 2003).

When developing a new method, it is important to compare its performance with that of existing methods. For seizure detection, such a task has been greatly simplified by our recent studies (Gao et al., 2011a, 2012a). By comparing seizure detection using a variety of complexity measures from deterministic chaos theory, random fractal theory, and information theory, we have found that the variations of those complexity measures with time have two patterns—either similar or reciprocal (Gao et al., 2011a). More importantly, we have gained fundamental understanding about the connections among different complexity measures through a new multiscale complexity measure, the SDLE. These results are recapitulated in Figure 1. While we leave the details to our prior works (Gao et al., 2006a, 2007, 2012a,b), these results suggest that it would be sufficient for us to compare the performance of the recurrence time based method for seizure detection with the performance of any of the existing complexity measures. Since some of the EEG data examined here had also been analyzed by the STLmax method and documented results exist, we shall compare our recurrence time method with the STLmax method. We shall show that the recurrence time method is both more accurate and faster than the STLmax method in detecting seizures from EEG.

**Figure 1. The variation with time of (A) λ _{small−ε}, (B) λ_{large−ε}, (C) Lyapunov exponent, (D) correlation entropy, (E) correlation dimension, and (F) the Hurst parameter obtained using DFA.** Adapted from Gao et al. (2011a).

The remainder of the paper is organized as follows. In section 2, we describe the data used here and the recurrence time method and the STLmax method for seizure detection. In section 3, we compare the performance of the recurrence time and STLmax method for seizure detection, as well as study seizure propagation. In section 4, we make a few concluding remarks.

## 2. Materials and Methods

In this section, we first describe EEG data used here, then describe the recurrence time method and the short-time Lyapunov exponent (STLmax) method.

### 2.1. Data

The EEG signals analyzed here are human EEG. They were recorded intracranially with approved clinical equipment by the Shands hospital at the University of Florida, with a sampling frequency of 200 Hz. Figure 2 shows our typical 28 electrode montage used for subdural and depth recordings.

**Figure 2. Schematic diagram of the depth and subdural electrode placement.** This view from the inferior aspect of the brain shows the approximate location of depth electrodes, oriented along the anterior–posterior plane in the hippocampi (RTD, right temporal depth; LTD, left temporal depth), and subdural electrodes located beneath the orbitofrontal and subtemporal cortical surfaces (ROF, right orbitoftrontal; LOF, left orbitofrontal; RST, right subtemporal; LST, left subtemporal).

Intracranial EEG is also called depth EEG, and is considered less contaminated by noise or motion artifacts. However, the clinical equipment used to measure the data has a pre-set, unadjustable maximal amplitude, which is around 5300 μV. This causes clipping of the signals when the signal amplitude is higher than this threshold. This is often the case during seizure episodes, especially for certain electrodes. To a certain extent, clipping complicates seizure detection, since certain seizure signatures may not be captured by the measuring equipment. However, we did not apply any filtering or conditioning methods to preprocess the raw EEG signals when we use our recurrence time method. The good results presented below thus suggest that the method is very reliable.

Altogether we have data of seven patients. The total duration of the measurement for each patient was up to about 3 days, as shown in the 2nd column of Table 1. There were only one or a few seizures for some patients while there were several tens of seizures for some other patients, as shown in the 3rd column of Table 1. Some of the seizures were considered subclinical, i.e., not manifested in the EEG signals. Sometimes the EEG signals may contain signatures distinctly different from background non-seizure signals, due to, for example, the fact that the patient may be eating food, drinking, etc. These non-seizure signatures typically may also be picked up by a seizure monitoring method. In this study, we shall focus on the behavior of the recurrence time and STLmax method in detecting seizures using only three channels EEG data without any preprocessing. As we shall see later, reliable decisions can be made based on single channel EEG data. There appears to be no need to combine multiple channels data.

### 2.2. Recurrence Time Based Method for Seizure Detection

The method involves first partitioning a long EEG signal into (overlapping or non-overlapping) blocks of data sets of short length *k*, and compute the so-called mean recurrence time of the 2nd type, ${\overline{{T}}}_{{2}}{(}{r}{)}$, for each data subset. For non-stationary and transient time series, it has been found (Gao, 1999, 2001; Gao and Cai, 2000; Gao et al., 2003) that ${\overline{{T}}}_{{2}}{(}{r}{)}$ will be different for different blocks of data subsets.

Let us first define the recurrence time of the 2nd type. Suppose we are given a scalar time series {*x*(*i*), *i* = 1, 2, …}. We first construct vectors of the form: *X*_{i} = [*x*(*i*), *x*(*i* + *L*), …, *x*(*i* + (*m* − 1)*L*)], with *m* being the embedding dimension and *L* the delay time (Packard et al., 1980; Takens, 1981; Sauer et al., 1991). {*X*_{i}, *i* = 1, 2, …, *N*} then represents certain trajectory in a *m*-dimensional space. Next, we arbitrarily choose a reference point *X*_{0} on the reconstructed trajectory, and consider recurrences to its neighborhood of radius *r*: *B*_{r}(*X*_{0}) = {*X*: ∥*X* − *X*_{0}∥ ≤ *r*}. The recurrence points of the 2nd type are defined as the set of points comprised of the first trajectory point getting inside the neighborhood from outside. These are denoted as the dark solid circles in Figure 3. The trajectory may stay inside the neighborhood for a while, thus generating a sequence of points, as designated by open circles in Figure 3. These are called sojourn points (Gao, 1999). It is clear that there will be more such points when the size of the neighborhood gets larger as well as when the trajectory is sampled more densely. The summation of the recurrence points of the second kind and the sojourn points is called the recurrence points of the first kind. These are often called nearest neighbors of the reference point *X*_{0}, and have been used by all other chaos theory-based non-linear methods.

**Figure 3. A schematic showing the recurrence points of the second type (solid circles) and the sojourn points (open circles) in B_{r}(X_{0})**.

Let us be more precise mathematically. We denote the recurrence points of the 1st type by *S*_{1} = {*X*_{t1}, *X*_{t2}, …, *X*_{ti} …}, and the corresponding Poincare recurrence time of the 1st type by {*T*_{1}(*i*) = *t*_{i + 1} − *t*_{i}, *i* = 1, 2, …}. Note the time is computed based on successive returns, not based on the returning points and the reference point. Also note *T*_{1}(*i*) may be 1 (for continuous time systems, this means one unit of the sampling time), for some *i*. This occurs when there are at least one sojourn point. Existence of such points makes further quantitative analysis difficult. Thus, we remove the sojourn points from the set *S*_{1} (which can be easily achieved by monitoring whether the recurrence times of the first type are one or not). Let us denote the remaining set by *S*_{2} = {*X*_{t′1}, *X*_{t′2}, …, *X*_{t′i} …}. *S*_{2} then defines a time sequence {*T*_{2}(*i*) = *t′*_{i + 1} − *t′*_{i}, *i* = 1, 2, …}. These are called the recurrence times of the 2nd type.

*T*_{2}(*i*) has a number of interesting properties: (1) For periodic motions, so long as the size of the neighborhood is not too large, *T*_{2}(*i*) accurately estimates the period of the motion. (2) For discrete sequences, the entire Renyi entropy spectrum can be computed from the moments of *T*_{2} (Gao et al., 2005a). (3) For chaotic motions, *T*_{2}(*i*) is closely related to the Lyapunov exponent, and hence, Kolmogorov entropy (Gao and Cai, 2000). (4) For chaotic motions, *T*_{2}(*i*) is related to the information dimension *d*_{1} by a simple scaling law (Gao, 1999; Gao et al., 2003),

where α takes on value 0 or 1, depending on whether the sojourn points form very few isolated points inside the neighborhood *B*_{r}(*X*_{0}), thus contribute dimension 0, or form a smooth curve inside *B*_{r}(*X*_{0}), thus contribute dimension 1. These properties make the recurrence time based method very versatile and powerful in detecting signal transitions.

We now explain how the mean recurrence time of the 2nd type can be computed. We simply evaluate this quantity for every reference point in a window, then take the mean of those times. Such calculation is carried out for all the data subsets, resulting in a curve which describes how ${\overline{{T}}}_{{2}}{(}{r}{)}$ varies with time. It has been observed (Gao, 1999, 2001; Gao and Cai, 2000; Gao et al., 2003) that the variations of ${\overline{{T}}}_{{2}}{(}{r}{)}$ coincide very well with sudden changes in the signal dynamics, such as bifurcations or transitions from regular motions to chaotic motions in non-stationary data, and vise versa. An example is shown in Figure 4 using the transient logistic map described by

We observe from Figure 4 that the method not only detects all the bifurcations in the signal, but also gives the exact periods of periodic signals. Note that some changes in a signal may be difficult to detect visually (Gao, 2001).

Since there are altogether four parameters involved, namely, the embedding dimension *m* and delay time *L*, the window length *k* for the data subsets, and the neighborhood size *r*, how shall we select them properly? To better illustrate the ideas, we postpone the discussion to section 3.1.1.

### 2.3. STLmax Method for Seizure Detection

The basic idea is to compute the largest positive Lyapunov exponent for each window's EEG signal using the Wolf et al.'s algorithm (Wolf et al., 1985) or its simple variants. Therefore, it is sufficient to describe the Wolf et al.'s algorithm (Wolf et al., 1985) and point out how it can be modified.

To apply the Wolf et al.'s algorithm (Wolf et al., 1985), one selects a reference trajectory and follows the divergence of its neighboring trajectory from it. Denote the reference and the neighboring trajectories by *X*_{i} = [*x*(*i*), *x*(*i* + *L*), …, *x*(*i* + (*m* − 1)*L*)], *X*_{j} = [*x*(*j*), *x*(*j* + *L*), …, *x*(*j* + (*m* − 1)*L*)], *i*, = 1, 2, …, *j* = *K*, *K* + 1, …, respectively. At the start of the time (which corresponds to *i* = 1), *X*_{K} is usually taken as the nearest neighbor of *X*_{1}. That is, *j* = *K* minimizes the distance between *X*_{j} and *X*_{1}. When time evolves, the distance between *X*_{i} and *X*_{j} also changes. Let the spacing between the two trajectories at time *t*_{i} and *t*_{i + 1} be *d*′_{i} and *d*_{i + 1}, respectively. Assuming *d*_{i + 1} ~ *d*′_{i}*e*^{λ1 (ti + 1 − ti)}, the rate of divergence of the trajectory, λ_{1}, over a time interval of *t*_{i + 1} − *t*_{i} is then

To ensure that the separation between the two trajectories is always small, when *d*_{i + 1} exceeds certain threshold value, it has to be renormalized: a new point in the direction of the vector of *d*_{i + 1} is picked up so that *d*′_{i + 1} is very small compared to the size of the attractor. After *n* repetitions of stretching and renormalizing the spacing, one obtains the following formula:

Note that this algorithm assumes but does not verify exponential divergence. In fact, the algorithm can yield a positive value of λ_{1} for any type of noisy process so long as all the distances involved are small. The reason for this is that when *d*′_{i} is small, evolution would move *d*′_{i} to the most probable spacing, which is typically much larger than *d*′_{i}. Then, *d*_{i + 1}, being in the middle step of this evolution, will also be larger than *d*′_{i}; therefore, a quantity calculated based on Equation (3) will be positive. This argument makes it clear that the algorithm cannot distinguish chaos from noise. In other words, even if the algorithm returns a positive λ_{1} from EEG data, one cannot conclude that the data are chaotic.

It is worth noting that in practice, to simplify implementation of the algorithm, one may replace the renormalization procedure described above by requiring that *d*′_{i + 1} is constructed whenever *t*_{i + 1} = *t*_{i} + *T*, where *T* is a small time interval. Such a procedure may be called periodic renormalization. In contrast, the original version of the algorithm is an aperiodic renormalization.

## 3. Results

### 3.1. Seizure Detection Using Recurrence Time Method

As we pointed out earlier, the method contains four parameters: the embedding dimension *m* and delay time *L*, the window length *k* for the data subsets, and the neighborhood size *r*. In this subsection, we first discuss how to choose these four parameters properly. Then we evaluate the effectiveness of the method for detecting epileptic seizures. For ease of presentation, we assume that the data have been normalized to the unit interval [0, 1] before further analysis.

#### 3.1.1. Parameter selection

First, we consider the window length *k* for data subsets. Since our purpose is to find transitions in the signal dynamics, the data subset has to be small. In order to estimate the interesting statistics reliably, a rule of thumb is that so long as a data subset contains several periods of “oscillations”, it would be fine (assuming the motion defines certain periodicity-like time scales). For our EEG sampled with a frequency of 200 Hz, we have found that *k* in the range of 500–2000 are all fine. Figures 5A,B show two examples, for *k* = 1000 and 2000, respectively. Clearly, in both cases, the two seizures have been detected correctly.

**Figure 5. Dependence of ${\overline{{T}}}_{{2}}$ on the parameters of the algorithm. (A–F)** Correspond to (*k*, *m*, *L*, *r*) = (1000, 4, 4, 2^{−4}), (2000, 4, 4, 2^{−4}), (2000, 4, 4, 2^{−3}), (2000, 3, 4, 2^{−4}), (2000, 4, 2, 2^{−4}), and (2000, 4, 6, 2^{−4}), respectively.

Next, we consider the size *r* of the neighborhood. It can be readily appreciated that when *r* is large, there will be a lot of recurrences, while when *r* is small, recurrences will be rather rare. This means ${\overline{{T}}}_{{2}}{(}{r}{)}$ will be large for small *r* but small for large *r*. Such expectations have been extensively observed in practice. For EEG signals, we have found that although the values of ${\overline{{T}}}_{{2}}{(}{r}{)}$ may vary with *r*, the pattern of the variation basically remains the same for a wide range of *r*. Two examples are shown in Figures 5B,C, where *r* differs by a factor of 2. Our experience is that choice of this parameter is not very critical, in so far as seizure monitoring is concerned.

Finally, we consider the embedding parameters. As is well known, the embedding parameters critically control the geometrical structure formed by the constructed vectors. Because of this feature, optimal embedding is a critical issue, especially when geometrical or dynamical quantities of the dynamics are concerned, such as the fractal dimension, Lyapunov exponents, and Kolmogorov entropy. For an in-depth discussion of this issue, we refer to Gao et al. (2007). Here, we wish to point out that the time scales associated with the motion are typically much less sensitive to the embedding parameters than the quantities such as the fractal dimension, Lyapunov exponents, and Kolmogorov entropy. To appreciate this feature, we have schematically shown in Figure 6 two different sets of embeddings. It is clear that the reconstructed trajectory shown in Figure 6A is fairly uniform, while that in Figure 6B is less so. One can readily conceive that when Figure 6B is further squeezed, the embedding quality is even worse. Judged by most optimal embedding criteria, the embedding shown in Figure 6A is considered a much better one than that shown in Figure 6B. However, it can be readily seen that ${\overline{{T}}}_{{2}}{(}{r}{)}$ for both Figures 6A,B are more or less the same. This means that the selection of *m* and *L* for computing ${\overline{{T}}}_{{2}}{(}{r}{)}$ is much less critical than that for computing other dynamical quantities. One good rule of thumb is that as long as the geometrical structure formed by the vectors are reasonably space-filling, the embedding is considered fine. Our experience with computing ${\overline{{T}}}_{{2}}{(}{r}{)}$ from EEG is that 3 ≤ *m* ≤ 6 are all fine, and with a sampling frequency of 200 Hz, *L* may be chosen 2–6. This discussion may be better appreciated by comparing Figures 5B,D–F, where four sets of (*m*, *L*) are illustrated. Clearly, all the parameter combinations have detected the two seizures accurately.

**Figure 6. A schematic showing the effect of embedding on the recurrence times of the second type. (A)** and **(B)** show examples of the reconstructed trajectory that is fairly uniform or less uniform, respectively.

To summarize, the recurrence method is much less sensitive to the parameters when compared with other non-linear methods, where embedding and other parameters have to be chosen carefully, and have to be specifically adapted to each patient for good results. For our recurrence time method, however, we have used the same parameter combination (*k*, *m*, *L*, *r*) = (2000, 4, 4, 2^{−4}) for all seven patients' data.

#### 3.1.2. Performance evaluation of the method

To illustrate the idea, we shall arbitrarily pick up three channels of EEG data,^{1} from one patient, and compare the patterns of variation of ${\overline{{T}}}_{{2}}{(}{r}{)}$ with that of STLmax. One typical result is shown in Figure 7. Vertical dotted lines indicate the seizure occurrence time determined by medical experts by viewing videotapes as well as the EEG signals. There are three seizures in Figure 7 during the period of time plotted. We observe that ${\overline{{T}}}_{{2}}{(}{r}{)}$ curves very cleanly and accurately detect all the seizures occurred. In fact, if one ignores the propagation-related slight timing difference (on the order of a few seconds up to 1 min; this will be further discussed later) among different electrodes, then most of the channels can be considered equivalent. In other words, decision can be based on single channel EEG data. This feature makes automatic detection of seizure by thresholding almost trivial. In contrast, the STLmax curves are much noisier than the ${\overline{{T}}}_{{2}}{(}{r}{)}$ curves. Although STLmax curves can be further post-processed to better reveal seizure information (Iasemidis et al., 2003), those features are still much weaker than those revealed by the recurrence time method.

**Figure 7. ${\overline{{T}}}_{{2}}{(}{r}{)}$ and STLmax vs. time curves for the EEG signals measured from three electrodes for a human patient. (A)** and **(B)** are for electrode LTD1, **(C)** and **(D)** for electrode LTD2, and **(E)** and **(F)** for electrode LTD4.

To more systematically compare the performance of the two methods in detecting seizures, we have computed positive detection (or equivalently, sensitivity) and false alarm per hour for the two methods. Positive detection is defined as the ratio between the number of seizures correctly detected and the total number of seizures. The false alarm per hour is simply the number of falsely detected seizures divided by the total time period. Table 1 summarizes the results. Clearly, the recurrence time method is more accurate than the STLmax method. This accuracy becomes even more attractive if one notices that the recurrence time method only involves simple thresholding, while the STLmax method involves a lot of further analysis (Iasemidis et al., 2003).

#### 3.1.3. Computational cost

The recurrence time method is very fast. With an ordinary PC (CPU speed less than 2 GHz), computation of ${\overline{{T}}}_{{2}}{(}{r}{)}$ from one channel EEG data of duration 1 h with sampling frequency of 200 Hz takes about 1 min CPU time. Computation of STLmax, on the other hand, takes more than 10 min. Hence, the recurrence time based method is much faster than the STLmax method. In fact, even with an ordinary PC, one is able to process all 28 channels of 1-h EEG data in about half an hour, therefore, faster than the data being continuously collected. With a more powerful PC, of course, the speed becomes even faster. Such a speed implies that the method can be used to real-time on-line process continuously collected all channels of EEG data. From an engineering perspective, the fast computation of recurrence time statistics can be considered overwhelming.

### 3.2. Propagation of Epileptic Seizures in the Brain

Formation and propagation of epileptic seizures in the brain is an outstanding example of complex spatial-temporal pattern formations. One of the most desirable ways of studying these problems is to understand how and when information flows from one region of the system to other regions. To resolve this issue, it is critical to accurately providing timing information for interesting events occurring in the system. With the exact timing information, one can then use concepts such as cross correlation and cross spectrum, mutual information, or measures from chaos theory, such as related to cross recurrence plots, to more fully characterize the spatial–temporal patterns. Recurrence time method can effectively provide such a timing information. To illustrate this point, we have shown in Figure 8 an example of analysis of multi-channel EEG signals using the recurrence time method. For the specific seizure studied, it was known that the seizure occurred around 200 s, and lasted about 2 min. While the recurrence time method has accurately detected the seizure, we note that the seizure activity recorded by electrode LTD3 and LTD5 was about 10 and 40 s later than that indicated by electrode LTD1, respectively. Hence, the recurrence time method not only accurately detects the seizure, but also provides invaluable timing information for the development of the seizure.

**Figure 8. ${\overline{{T}}}_{{2}}{(}{r}{)}$ curves for EEG signals measured by three electrodes.** The dashed vertical line indicate the seizure starting position around 200 s. The seizure lasted for about 2 min. Note that from **(A)** LTD1 to **(B)** LTD3, the seizure activity is delayed about 10 s, while from **(A)** LTD1 to **(C)** LTD5, the seizure activity is delayed about 30–40 s.

## 4. Conclusions

Motivated by developing a non-linear method without the limitations of assuming that EEG signals are chaotic and stationary, we have proposed a recurrence time based approach to characterize brain electrical activity. The method is very easy to use, as it only contains a few parameters that are largely signal-independent. It very accurately detects epileptic seizures from EEG signals. Most critically, the method is very fast—it is fast enough to real-time on-line process multi-channel EEG data with a typical PC. Therefore, it has the potential to be an excellent candidate for real-time monitoring of epileptic seizures in a clinical setting.

The recurrence time method is also able to accurately give the timing information critical for understanding seizure propagation. Therefore, it may help characterize epilepsy type, lateralization and seizure classification (Holmes, 2008; Napolitano and Orriols, 2008; Plummer et al., 2008). To more thoroughly understand the capabilities of recurrence time method in characterizing seizure propagation, it would be desirable to combine recurrence time analysis of EEG with studies based on MEG and MRI exams.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Footnotes

- ^In fact, the three chosen channels EEG data may not correspond to where a seizure was localized. This further indicates the robustness of our method.

## References

Adeli, H., Ghosh-Dastidar, S., and Dadmehr, N. (2007). A wavelet-chaos methodology for analysis of eegs and eeg subbands to detect seizure and epilepsy. *IEEE Trans. Biomed. Eng*. 54, 205–211. doi: 10.1109/TBME.2006.886855

Aschenbrenner-Scheibe, R., Maiwald, T., Winterhalder, M., Voss, H., Timmer, J., and Schulze-Bonhage, A. (2003). How well can epileptic seizures be predicted? an evaluation of a nonlinear method. *Brain* 126, 2616. doi: 10.1093/brain/awg265

Cao, Y., Tung, W., Gao, J., Protopopescu, V., and Hively, L. (2004). Detecting dynamical changes in time series using the permutation entropy. *Phys. Rev. E* 70:046217. doi: 10.1103/PhysRevE.70.046217

Gao, J. (1999). Recurrence time statistics for chaotic systems and their applicaitons. *Phys. Rev. Lett*. 83, 3178–3181. doi: 10.1103/PhysRevLett.83.3178

Gao, J. (2001). Detecting nonstationarity and state transitions in a time series. *Phys. Rev. E* 63:066202. doi: 10.1103/PhysRevE.63.066202

Gao, J., and Cai, H. (2000). On the structures and quantification of recurrence plots. *Phys. Lett. A* 270, 75–87. doi: 10.1016/S0375-9601(00)00304-2

Gao, J., Cao, Y., Gu, L., Harris, J., and Principe, J. (2003). Detection of weak transitions in signal dynamics using recurrence time statistics. *Phys. Lett. A* 317, 64–72. doi: 10.1016/j.physleta.2003.08.018

Gao, J., Cao, Y., Qi, Y., and Hu, J. (2005a). Building innovative representations of dna sequences to facilitate gene finding. *IEEE Intell. Syst*. 20, 34–39. doi: 10.1109/MIS.2005.100

Gao, J., Tung, W., Cao, Y., Hu, J., and Qi, Y. (2005b). Power-law sensitivity to initial conditions in a time series with applications to epileptic seizure detection. *Physica A* 353, 613–624. doi: 10.1016/j.physa.2005.01.027

Gao, J., Cao, Y., Tung, W., and Hu, J. (2007). Multiscale Analysis of Complex Time Series – Integration of Chaos and Random Fractal Theory, and Beyond. Hoboken, NJ: Wiley. doi: 10.1002/9780470191651

Gao, J., Hu, J., and Tung, W. (2011a). Complexity measures of brain wave dynamics. *Cogn. Neurodyn*. 5, 171–182. doi: 10.1007/s11571-011-9151-3

Gao, J., Hu, J., and Tung, W. (2011b). Facilitating joint chaos and fractal analysis of biosignals through nonlinear adaptive filtering. *PLoS ONE* 6:e24331. doi: 10.1371/journal.pone.0024331

Gao, J., Hu, J., and Tung, W. (2012a). Entropy measures for biological signal analysis. *Nonlin. Dyn*. 68, 431–444. doi: 10.1007/s11071-011-0281-2

Gao, J., Hu, J., Tung, W., and Blasch, E. (2012b). Multiscale analysis of physiological data by scale-dependent lyapunov exponent. *Front. Fractal Physiol*. 2:110. doi: 10.3389/fphys.2011.00110

Gao, J., Hu, J., Tung, W., and Cao, Y. (2006a). Distinguishing chaos from noise by scale-dependent lyapunov exponent. *Phys. Rev. E* 74:066204. doi: 10.1103/PhysRevE.74.066204

Gao, J., Hu, J., Tung, W., Cao, Y., Sarshar, N., and Roychowdhury, V. (2006b). Assessment of long range correlation in time series: how to avoid pitfalls. *Phys. Rev. E* 73:016117. doi: 10.1103/PhysRevE.73.016117

Holmes, M. (2008). Dense array eeg: methodology and new hypothesis on epilepsy syndromes. *Epilepsia* 49(Suppl. 3), 3–14. doi: 10.1111/j.1528-1167.2008.01505.x

Hwa, R., and Ferree, T. (2002). Scaling properties of fluctuations in the human electroencephalogram. *Phys. Rev. E* 66:021901. doi: 10.1103/PhysRevE.66.021901

Iasemidis, L., Pardalos, P., Shiau, D., Chaovalitwongse, W., Narayanan, K., Kumar, S., et al. (2003). Prediction of human epileptic seizures based on optimization and phase changes of brain electrical activity. *Optim. Method Softw*. 18, 81–104. doi: 10.1080/1055678021000054998

Iasemidis, L., Sackellares, J., Zaveri, H., and Williams, W. (1990). hase space topography and the lyapunov exponent of electrocorticograms in partial seizures. *Brain Topogr*. 2, 187–201. doi: 10.1007/BF01140588

Jallon, P. (1997). Epilepsy in developing countries. *Epilepsia* 38, 1143–1151. doi: 10.1111/j.1528-1157.1997.tb01205.x

Kanz, H., and Schreiber, T. (1997). Nonlinear Time Series Analysis. Cambridge, NY: Cambridge University Press.

Lai, Y., Harrison, M., Frei, M., and Osorio, I. (2003). Inability of lyapunov exponents to predict epileptic seizures. *Phys. Rev. Lett*. 91:068102. doi: 10.1103/PhysRevLett.91.068102

Lehnertz, K., and Elger, C. (1995). Spatiotemporal dynamics of the primary epileptogenic area in temporal-lobe epilepsy characterized by neuronal complexity loss. *Electroencepphalogr. Clin. Neurophysiol*. 95, 108–117. doi: 10.1016/0013-4694(95)00071-6

Lehnertz, K., and Elger, C. (1997). Neuronal complexity loss in temporal lobe epilepsy: effects of carbamazepine on the dynamics of the epileptogenic focus. *Electroencepphalogr. Clin. Neurophysiol*. 103, 376–380. doi: 10.1016/S0013-4694(97)00027-1

Martinerie, J., Adam, C., Quyen, M. L. V., Baulac, M., Clemenceau, S., Renault, B., et al. (1998). Epileptic seizures can be anticipated by non-linear analysis. *Nat. Med*. 4, 1173–1176. doi: 10.1038/2667

Napolitano, C., and Orriols, M. (2008). Two types of remote propagation in mesial temporal epilepsy: analysis with scalp ictal EEG. *J. Clin. Neurophysiol*. 25, 69–76. doi: 10.1097/WNP.0b013e31816a8f09

Packard, N., Crutchfield, J., Farmer, J., and Shaw, R. (1980). Gemomtry from time-series. *Phys. Rev. Lett*. 45, 712–716. doi: 10.1103/PhysRevLett.45.712

Plummer, C., Harvey, S., and Cook, M. (2008). Eeg source localization in focal epilepsy: where are we now? *Epilepsia* 49, 201–218. doi: 10.1111/j.1528-1167.2007.01381.x

Protopopescu, V., Hively, L., and Gailey, P. (2001). Epileptic event forewarning from scalp EEG. *J. Clin. Neurophysiol*. 18, 223–245. doi: 10.1097/00004691-200105000-00003

Quyen, M. L. V., Martinerie, J., Navarro, V., Boon, P., D'Have, M., Adam, C., et al. (2001). Anticipation of epileptic seizures from standard EEG recordings. *Lancet* 357, 183. doi: 10.1016/S0140-6736(00)03591-1

Sauer, T., Yorke, J., and Casdagli, M. (1991). Embedology. *J. Stat. Phys*. 65, 579–616. doi: 10.1007/BF01053745

Takens, F. (1981). “Detecting strange attractors in turbulence,” in *Dynamical Systems and Turbulence*. Lecture Notes in Mathematics, Vol. 898, eds D. A. Rand and L. S. Young (Berlin: Springer-Verlag), 366.

Valencia, M., Artieda, J., Alegre, M., and Maza, D. (2008). Influence of filters in the detrended fluctuation analysis of digital electroencephalographic data. *J. Neurosci. Methods* 170, 310–316. doi: 10.1016/j.jneumeth.2008.01.010

van Drongelen, W., Nayak, S., Frim, D., Kohrman, M., Towle, V., Lee, H., et al. (2003). Seizure anticipation in pediatric epilepsy: use of kolmogorov entropy. *Pediatr. Neurol*. 29, 207–213. doi: 10.1016/S0887-8994(03)00145-0

Keywords: EEG, recurrence time, seizure detection, seizure propagation, brain complexity

Citation: Gao J and Hu J (2013) Fast monitoring of epileptic seizures using recurrence time statistics of electroencephalography. *Front. Comput. Neurosci*. **7**:122. doi: 10.3389/fncom.2013.00122

Received: 28 July 2013; Accepted: 04 September 2013;

Published online: 01 October 2013.

Edited by:

Tobias A. Mattei, Ohio State University, USAReviewed by:

Tobias A. Mattei, Ohio State University, USAPaul Rapp, Uniformed Services University of the Health Sciences, USA

Copyright © 2013 Gao and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jianbo Gao, Institute of Complexity Science and Big Data Technology, Guangxi University, 100 Daxue Road, Nanning 530005, China e-mail: jbgao.pmb@gmail.com