# Nonlinear analysis of EEG signals at different mental states

- Kannathal Natarajan
^{1}Email author, - Rajendra Acharya U
^{1}, - Fadhilah Alias
^{1}, - Thelma Tiboleng
^{1}and - Sadasivan K Puthusserypady
^{2}

**3**:7

**DOI: **10.1186/1475-925X-3-7

© Natarajan et al; licensee BioMed Central Ltd. 2004

**Received: **18 September 2003

**Accepted: **16 March 2004

**Published: **16 March 2004

## Abstract

### Background

The EEG (Electroencephalogram) is a representative signal containing information about the condition of the brain. The shape of the wave may contain useful information about the state of the brain. However, the human observer can not directly monitor these subtle details. Besides, since bio-signals are highly subjective, the symptoms may appear at random in the time scale. Therefore, the EEG signal parameters, extracted and analyzed using computers, are highly useful in diagnostics. This work discusses the effect on the EEG signal due to music and reflexological stimulation.

### Methods

In this work, nonlinear parameters like Correlation Dimension (CD), Largest Lyapunov Exponent (LLE), Hurst Exponent (H) and Approximate Entropy (ApEn) are evaluated from the EEG signals under different mental states.

### Results

The results obtained show that EEG to become less complex relative to the normal state with a confidence level of more than 85% due to stimulation.

### Conclusions

It is found that the measures are significantly lower when the subjects are under sound or reflexologic stimulation as compared to the normal state. The dimension increases with the degree of the cognitive activity. This suggests that when the subjects are under sound or reflexologic stimuli, the number of parallel functional processes active in the brain is less and the brain goes to a more relaxed state

## Background

The electrical activity of a brain measured by Electroencephalogram (EEG) exhibits complex behavior with nonlinear dynamic properties. This behavior takes the form of EEG patterns with different complexities. Considering this, the nonlinear dynamics theory may be a better approach than traditional linear methods in characterizing the intrinsic nature of EEG. The study of nonlinear dynamics and characterization can contribute to the understanding of the EEG dynamics and underlying brain processes and search for its physiological significance. The literature on the study of the application of the nonlinear dynamics theory to analyze physiological signals, shows that nonlinear approaches were used for analysis of heart rate, nerve activity, renal blood flow, arterial pressure, EEG and respiratory signals [1, 2].

The importance of the biological time series analysis, which exhibits typically complex dynamics, has long been recognized in the area of nonlinear analysis. Several features of these approaches have been proposed to detect the hidden important dynamical properties of the physiological phenomenon. The nonlinear dynamical techniques are based on the concept of chaos and it has been applied to many areas including the areas of medicine and biology. The theory of chaos has been used to detect some cardiac arrhythmia such as ventricular fibrillation [3]. Efforts have been made in determining nonlinear parameters for pathological signals and it has been shown that they are useful indicators of pathologies.

Many investigators, for example, Duke et al [4], has proved that complex dynamical evolutions lead to chaotic regimes. In the last thirty years, experimental observations have pointed out that, in fact, chaotic systems are common in nature. A detail of such system is given in by Boccaletti et al. [5]. In theoretical modeling of neural systems, emphasis has been put mainly on either stable or cyclic behaviors. Perhaps chaotic behavior at neural level could be responsible for schizophrenia, insomnia, epilepsy and other disorders [6–8]. In the past wide range of work has been done in understanding the complexities associated with the brain through multiple windows of mathematics, physics, engineering and chemistry, physiology etc [9–11]. During the past, there has been a sustained interest in describing neural processes and brain signals, especially the EEG, within the context of nonlinear dynamics and theory of deterministic chaos [12, 13]. Nonlinear dynamics theory opens new window for understanding behavior of EEG. EEG models were proposed by Freeman et al [14] for neocortical dynamics or by Wright et al[15] for chaotic dynamics to meet requirements in neurobiology. In analysis of EEG data, different chaotic measures such as correlation dimension, Lyapunov exponent and entropy are used in recent literature [16–22]. Approximate Entropy of human respiratory movement during eye-closed walking and different sleep stages was studied by Naoto et al[23].

In this work, EEG is recorded for various mental states such as (1) subjects under normal resting state, (2) subjects listening to classical music, (3) subjects listening to rock music and (4) subjects with foot reflexology stimulation. The effect of the music and reflexological stimulation on the EEG signal is studied by evaluating the nonlinear parameters like Correlation dimension, Approximate Entropy, Largest Lyapunov Exponent and Hurst Exponent.

## Methods

EEG signals are being recorded by using BIOPAC equipment with ACQKNOWLEDGE 3.7.2 as data acquisition software [24]. Silver-Silver chloride electrodes were applied following the 10–20 system. The EEG is recorded at 500 samples/sec with a resolution of 12 bits/sample and for a duration of twenty minutes. Then the data is digitally filtered using 1–50 Hz band pass filter. The EEG is first recorded for the subject under resting condition with eyes closed. In second state, the EEG is recorded with subject listening to classical music and then followed by rock music. According to Voss et al. [25], the classical and rock music exhibit different 1/f distribution. Here the hypothesis is that the since the two music have different 1/f distribution they will have different effects on EEG. Finally the EEGs are recorded with subjects under foot reflexologic stimulation. Reflexology is a science, an art and a skill that uses gentle pressure and manipulation of the feet, hands and ears. It deals with the principle that, there are reflex areas in the feet and hands correspond to various glands, organs and parts of the body. EEGs are recorded from 30 subjects (15 boys and 15 girls) with age ranging from 20 ± 3 years for 20 minutes. The data is recorded in sequential order with subject (1) under normal resting state, (2) listening to classical music (3) listening to rock music and (4) with foot reflexology stimulation. Time gap of 30 minutes is given between any two successive recordings to ensure that a previous state stimulus doesn't affect the subsequent state EEG.

### Analysis

In this work, we have analyzed the EEGs using various characteristic measures like Correlation Dimension, Lyapunov exponent, Hurst Exponent, and Approximate Entropy. A brief description of the measures is given below.

### Correlation Dimension

The dimension of a graph can give much information about the nature of the system from which it is derived. The dimension of a graph determined by experimental data is therefore very useful to the understanding of the system, particularly for deciding if the system is periodic, chaotic, or noisy. Mathematically, the dimension of any graph that is made of a finite set of data points is zero. The correlation dimension however can be used to approximate the dimension of these types of graphs. It is most often computed from the time series diagram. A time series diagram is a simplified phase space diagram constructed from a single data vector.

*EEG signal*X[n] and the Y-axis represents the

*EEG signal after a delay*X[n+delay]. The choice of an appropriate delay is calculated using the minimal mutual information technique [26, 27]. The phase space plot for various mental states is shown in Figure 1.

Correlation Dimension is one of the most widely used measures of Fractal Dimension. Here we adapt the algorithm proposed by Grassberger and Procaccia [28, 29]. The idea is to construct a function *C*(*r*) that is the probability that two arbitrary points on the orbit are closer together than *r*. This is done by calculating the separation between every pair of *N* data points and sorting them into bins of width d*r* proportionate to *r*. A correlation dimension can be calculated using the distances between each pair of points in the set of N number of points, *s*(*i*, *j*) = |*X*
_{
i
}- *X*
_{
j
}|.

A correlation function, C(r), is then calculated using,

Where *X*
_{
i
}, *X*
_{
j
}→ points of the trajectory in the phase space,

*N* → is the number of data points in phase space

*r* → radial distance around each reference point X_{i}

Θ → is the Heaviside function.

Correlation dimension was calculated using the fundamental definition

_{corr}, where D

_{corr}is the correlation Dimension and m is the embedding dimension. But the limitation on the applicability of the theorems is that the measure D and D

_{corr}are not known. For practical applications, it is best to apply the Grassberger and Procaccia algorithm [28] and calculate the D

_{corr}for various embedding dimension. Then the minimum embedding dimension of the attractor for one to one embedding is m+1, where m is the embedding dimension above which the D

_{corr}saturates. In this work, we have calculated D

_{corr}with an embedding dimension varying from 1 to 10 for all the subjects. The graph of D

_{corr}Vs Embedding Dimension is shown in Figure 2. From the graph it can be seen that D

_{corr}saturates after the embedding dimension of 9. So we have chosen the embedding dimension of m+1 = 10 for analysis. The software used for analysis is CDA Pro Data analyzer [33].

The data points used in calculating the autocorrelation come from a time series, which has too small or too large time resolution may introduce spurious effects. Here, the time resolution is too small, the data may contain multiple copies of essentially the same measurements, which leads to multiple-counting. The correlation dimension is then artificially low because all the points are temporally close to each other. This effect occurs when the time resolution Δ*t* of the analyzed data is much smaller than the autocorrelation time τ_{
ac
}(or any characteristic time scale).

The correction proposed by Theiler [34] is, for each reference point *X*
_{
i
}to include measurements *X*
_{
j
}which are at least *w* = τ_{
ac
}steps away from *X*
_{
i
}

### Approximate Entropy

Entropy is a thermodynamic quantity describing the amount of disorder in the system. From an information theory perspective, the above concept of entropy is generalized as the amount of information stored in a more general probability distribution. Recently a number of different entropy estimators [35–40] have been applied to EEG data to quantify the complexity of the EEG signal. These techniques do not measure the shape of the distribution of the EEG voltages per se, but instead describe how the EEG signals changes with time – either in frequency-space or phase-space. Information about how the EEG signal fluctuates with time is obtained by comparing the time series with itself, but lagged by a specified time interval. This practice is referred as embedding the one-dimensional signal in a 'phase-space'. Intuitively it seems sensible that, if an EEG signal is irregular, the position of a particular point will not be easily predicted using knowledge of its previous points; whereas in a regular signal the position of the point will be more reliably predicted. The number of previous (lagged) points used in making the prediction is the embedding dimension (*m*). For a process whose underlying dimension is *n* (i.e. which can be described uniquely in terms of *n* parameters), the required embedding dimension is; *m* ≥ *2n+1*, and the required minimum date size to extract these *n* parameters is 10^{
m
}. It is impractical to properly embed the EEG signal. Thus these techniques are NOT able to fulfill their theoretical promise and extract high-dimensional information from the univariate EEG data-stream. Using these embeddings, the theoretical measure of the rate of "information" generation by a system is the Kolmogorov-Sinai entropy [29]. However this measure diverges to a value of infinity when the signal is contaminated by the slightest noise. Pincus, Gladstone et al. proposed Approximate Entropy (ApEn) as a solution to these problems and successfully applied it to relatively short and noisy data [41, 42]. The ApEn was used by Bruhn [43] to analyze EEG signals in patients under general anesthesia. According to Steyn-Ross et al [44] the approximate entropy when applied to EEG signals reflects the intra-cortical information flow in the brain. In the calculation ApEn, two parameters *m* and *r* must be chosen prior to the computation of ApEn where, *m* specifies the pattern length, and *r* is the noise threshold. The ApEn measure is obtained by

For this study, *m* is set to 2 and *r* is set to 15% of the standard deviation of each time series. These values are selected on the basis of previous studies indicating good statistical validity for ApEn within these variable ranges [42].

### Hurst Exponent (H)

The Hurst Exponent is a measure that has been widely used to evaluate the self-similarity and correlation properties of fractional Brownian noise, the time series produced by a fractional (fractal) Gaussian process. Hurst exponent is used to evaluate the presence or absence of long-range dependence and its degree in a time series. However, local trends (nonstationarities) is often present in physiological data and may compromise the ability of some methods to measure self-similarity. Hurst Exponent is the measure of the smoothness of a fractal time series based on the asymptotic behavior of the rescaled range of the process. In time series analysis of EEG, Hurst Exponent H is used by [45] for characterize the non-stationary behavior of the sleep EEG episodes. The hurst exponent H is defined as:

H = log(R/S)/log (T)

where T is the duration of the sample of data and R/S the corresponding value of rescaled range. The above expression is obtained from the Hurst's generalized equation of time series that is also valid for Brownian motion. If H = 0.5, the behavior of the time-series is similar to a random walk. If H < 0.5, the time-series cover less "distance" than a random walk. But if H > 0.5, the time-series covers more "distance" than a random walk. H is related to the fractal dimension D given by

H = E+1-D

where E is the Euclidean dimension.

### Largest Lyapunov Exponent (LLE)

Lyapunov Exponent (λ) is a quantitative measure of the sensitive dependence on the initial conditions. It defines the average rate of divergence of two neighboring trajectories. An exponential divergence of initially nearby trajectories in phase space coupled with folding of trajectories, to ensure that the solutions will remain finite, is the general mechanism for generating deterministic randomness and unpredictability. Therefore, the existence of a positive λ for almost all initial conditions in a bounded dynamical system is widely used definition of deterministic chaos. To discriminate between chaotic dynamics and periodic signals Lyapunov exponent (λ) are often used. It is a measure of the rate at which the trajectories separate one from other. The trajectories of chaotic signals in phase space follow typical patterns. Closely spaced trajectories converge and diverge exponentially, relative to each other. For dynamical systems, sensitivity to initial conditions is quantified by the Lyapunov exponent (λ). They characterize the average rate of divergence of these neighboring trajectories. A negative exponent implies that the orbits approach a common fixed point. A zero exponent means the orbits maintain their relative positions; they are on a stable attractor. Finally, a positive exponent implies the orbits are on a chaotic attractor.

The algorithm proposed by Wolf et al. [46] is used to Largest LE (LLE) from EEG data. For Given the time series x(t) for m dimensional phase space with delay coordinate t, that is a point on the attractor is given by

{x(t), x(t + t),..., x(t + (m-1)t}

We locate nearest neighbor to initial point

{x(t0), x(t0+ t),..., x(t0 + (m-1)t}

And denote the distance between these two points as L(t0). At a later time t1, initial length will evolve to length L'(t1). The mean exponential rate of divergence of two initially close orbits is characterized by

In implementation of this program, the following set of numerical parameters has to be chosen:

P = {m, t, T, Smax, Smin, thmax}

where m is the embedding dimension, t is delay, T being evaluation time (= t k+1 - t k-1) and S_{max}, S_{min} are the maximum and minimum separations of replacement point respectively and th_{max} is the maximum orientation error. According to Das et al [47] an embedding dimension between 5 to 20 and a delay of 1 should be chosen when calculating LE for EEG data. In our analysis we have chosen an embedding dimension of 10 and delay of 1.

### Surrogate data

To test if the attractor geometry and correlation dimension are truly due to chaotic dynamics one must examine these characteristics for surrogate data sets.

The purpose of surrogate data is to test for any nonlinearity in the original data. Nonlinear indexes such as ApEn are computed for several surrogate data series. Their values are compared with that assumed by the nonlinear index computed for the original index [48]. The demonstration of statistically significant difference in ApEn between the original and surrogate data are in keeping with the presence of nonlinear dynamics in the original data.

Surrogate data have Fourier decomposition with the same amplitudes as the empirical data decomposition but with random phase components. This is obtained from the Chaos Data Analyzer.

To test for a statistical significance of difference in original ApEn and the surrogate data, 10 surrogate data series were generated to match each original signal. Then it is subjected to Student t test distribution. We found that, the surrogate data ApEn and Original data ApEn, are close to each other with confidence level of 95%. The similar procedure is repeated for Correlation Dimension. The surrogate data Correlation Dimension and original data Correlation Dimension are different from each other with confidence level of 93%. Hence the original data contain nonlinear features.

## Results

_{corr}is estimated for all the experimental data sets and result is shown in Figure 3. It quantifies the variability in a time series. The signal becomes less complex (less random) when the person is subjected music (of his choice) or under reflexological stimulation. Similarly Largest Lyapunov Exponent, Approximate Entropy and Hurst Exponent are calculated for all the data sets are presented in Figure 4, 5, 6. The Figure 4, depicts the predictability of the signal increases and the signal becomes less chaotic under the influence of music and reflexologic stimulation. Largest Lyapunov Exponent behaves as indicator of long-term behavior. The LLE value closer to 1 indicates the chaotic behavior of the series. This value falls due to the influence of the music and reflexological stimulation. Hence the randomness falls due this external influence. From the Figure 5, it can be seen that the same trend is obtained for ApEn. ApEn is the measure of dynamic changes of the EEG signal in time domain. A decrease in entropy indicates higher predictability and a reduced stochastic behavior. The self-similarity parameter, Hurst Exponent H is obtained for all data sets is shown in Figure 6. If H = 0.5 indicates the random behavior of the signal and this value falls under the influence of external stimulus. This means that the randomness decreases due to music/reflexology. The interesting finding is that the randomness or the disorderliness of the EEG is reduced after the subject hears music or reflexologic stimulation. Having computed the values of nonlinear parameters of all the subjects under normal resting condition and under music or reflexologic stimulation, we have performed the statistical analysis on these numeral characteristics. All the results are presented as mean ± SD with "p" values in Table 1. The results indicate after the music/ reflexological stimulus, the brain goes relaxed state. In neurophysiology the most cited indicator of relaxation is the rise in alpha wave frequencies will be predominant.

Range of nonlinear parameter variation for different mental states

Parameter | Normal resting state | With classical music | With rock music | With foot reflexologic stimulation | p-value |
---|---|---|---|---|---|

Correlation Dimension | 5.0426 ± 0.286 | 4.921 ± 0.258 | 4.846 ± 0.284 | 4.828 ± 0.269 | 0.0956 |

Lyapunov Exponent | 0.6185 ± 0.0107 | 0.6133 ± 0.0112 | 0.6072 ± 0.0122 | 0.5297 ± 0.0231 | 0.0701 |

Approximate Entropy | 2.8214 ± 0.146 | 2.762 ± 0.157 | 2.6833 ± 0.181 | 2.615 ± 0.191 | 0.1346 |

Hurst Exponent | 0.2887 ± 0.0568 | 0.3133 ± 0.0432 | 0.3372 ± 0.0516 | 0.3497 ± 0.0523 | 0.064 |

## Discussion

The statistical analysis is performed by comparing the results obtained from various mental states with respect to the normal resting state. Table 1 depicts the statistical outcome of the nonlinear analysis between the normal and other states. It is found that the measures are significantly lower when the subjects are under sound or reflexologic stimulation as compared to the normal state. It is well known that the dimension of EEG time series is closely related to the cognitive activity of the brain [49]. The dimension increases with the degree of the cognitive activity. Our results show that dimension reduces during reflexological stimulation indicating that the brain is not involved in cognitive tasks or thinking rigorously, meaning that the brain is in a stoic, stolid, passive state of relaxation. From the result it can be seen that there is a distinct difference in the correlation dimension in the different mental states with a confidence level of more than 90%. The correlation dimension for the normal group (mean = 5.0426 and SD = 0.286) is much higher compared to the other groups. A decreased value indicates that the randomness of the brain activity is reduced. This suggests that when the subjects are under sound or reflexologic stimuli, the number of parallel functional processes active in the brain is less and the brain goes to a more relaxed state [50–55]. This gives rise to the increase in alpha frequencies in the brain waves. The same trend is seen for parameters like Lyapunov exponent and ApEn (LLE (p < 0.07), ApEn (p < 0.13) and H(p < 0.064)). From Table 1 it is evident that the nonlinear complexity of the EEG is reduced when subjected to music and reflexolgical stimulus. Hurst exponent (p < 0.06) also shows the same trend because of the very definition of the term Hurst exponent [H = E+1-D, where 'H' is the Hurst exponent, 'E' is the Euclidean dimension and 'D' is the correlation dimension]. Hence as 'D' decreases 'H' has to increase [45]. These results show that the effect of sound and reflexologic stimulation causes EEG to become less complex relative to the normal state.

## Conclusion

EEG signal can be used as a reliable indicator of the state of the mind. Ever since the birth of 'nonlinear science' chaoticians of physiology, biomedical engineering and theoretical biology are searching for meaningful chaotic parameters in physiological processes. Here we have evaluated linear parameters and nonlinear parameters: CD, LLE, ApEn and Hurst Exponent. The results of these were subjected to ANOVA test with more than 90% confidence interval giving excellent 'p' values in all cases. The importance of ApEn lies in the fact that it is measure of the disorder in the EEG signal. This ApEn falls indicating the decrease in the disorder of the EEG signal due to rock/classical music and reflexology. This indicates that, the ∝ wave tends to become more predominant due to the influence of the music and reflexology. Similarly, Correlation Dimension, LLE and Hurst Exponent decreases due to the music and reflexology. In this work, we have proposed a set of ranges for these nonlinear parameters for various mental states.

## Declarations

## Authors’ Affiliations

## References

- Stein KM, Lippman N, Kligfield P:
**Fractal rhythms of the heart.***Journal of Electrocardiology*1992,**24:**72–76.View Article - Hoyer D, Schmidt K, Bauer R, Zwiener U, Kohler M, Luthke B, Eiselt M:
**Nonlinear analysis of heart rate and respiratory dynamics.***IEEE Engineering in Medicine and Biology Magazine*1997,**16**(1):31–39. 10.1109/51.566150View Article - Kaplan DK, Cohen JR:
**Searching for Chaos in fibrillation.***Annals NY Academic Science*1991, 367–374. - Duke W, Pritchard WS, Krieble KK:
**Dimensional analysis of resting human EEG II: Surrogate data testing indicates nonlinearity but not low-dimensional chaos.***Psychophysiology*1995,**32:**486–491.View Article - Boccaletti S, Grebogi C, Lai YC, Mancini H, Mazaet D:
**The control of chaos: Theory and applications.***Physics Reports*2000,**329:**108–109. 10.1016/S0370-1573(99)00096-4View Article - Glass L, Michel RG, Mackey M, Shrier A:
**Chaos in neurobiology.***IEEE Trans Sys Man & Cybernatics*1983,**13**(5):790–798. - Jaeseung J, Jeong-Ho C, Kim SY, Seol-Heui H:
**Nonlinear dynamical analysis of the EEG in patients with Alzheimer's disease and vacular dementia.***Clin Neurophysiol*2001,**18**(1):58–67.View Article - Philippe F, Henri K:
**Is there chaos in the brain? Concepts of nonlinear dynamics and methods of investigation.***Life Sciences*2001,**324:**773–793. - Babloyantz A:
**Evidence of chaotic dynamics during sleep cycle. In Dimensions and entropies in Chaotic system.***(Edited by: Mayer-Kress G).*Berlin: Springer-Verlag 1986. - Bessar E:
**Biophysical and physiological systems analysis.**London: Addison-Wesley 1960. - Mayer G, Layne SC:
**Dimensionality of the human encephalogram. In perspectives in Biological Dynamics and Theoretical Medicine.***Annals NY Academic Science**(Edited by: Koslow SH, Madell AJ, Shlessinger MF).*1987, 504–507. - Rapp PE, Bashore T, Martinerie J, Albano A, Zimmerman I, Mess A:
**Dynamics of brain electrical activity.***Brain Topography*1989,**2:**99–118.View Article - Rapp PE:
**Chaos in the neurosciences: cautionary tales from the frontier.***Biologist*1993,**40:**89–94. - Freeman WJ:
**A proposed name for aperiodic brain activity: stochastic chaos.***Neural Networks*2000,**13:**11–13. 10.1016/S0893-6080(99)00093-3MathSciNetView Article - Wright JJ, Liley DTJ:
**Dynamics of the brain at global and microscopic scales. Neural networks and the EEG.***Behavioral and Brain Sciences*1996,**19:**285–320.View Article - Pijn JPM:
**Quantitative evaluation of EEG signals in epilepsy: nonlinear association time delays and nonlinear dynamics.***Ph.D. thesis*University of Amsterdam 1990. - Pijn JP, Velis DN, van der Heyden MJ, DeGoede J, van Veelen CW, Lopes da Silva FH:
**Nonlinear dynamics of epileptic seizures on basis of intracranial EEG recordings.***Brain Topography*1997,**9:**249–270.View Article - Rombouts SARB, Keunen RWM, Stam CJ:
**Investigation of nonlinear structure in multichannel EEG.***Phys Lett A*1995,**202:**352–358. 10.1016/0375-9601(95)00335-ZView Article - Lamberts J, Van den Broek PLC, Bener J, Van Egmond J, Dirksen R, Cohen AML:
**Correlation dimension of the human electroencephalogram corresponding to cognitive load.***Neuropsychobiology*2000,**41**(3):149–153. 10.1159/000026647View Article - Bradley E:
**Time-series analysis. Intelligent data analysis: An introduction.***(Edited by: Berthold M, Hand D).*Berlin: Springer-Verlag 1999. - Theiler J:
**Spurious dimension from correlation algorithms applied to limited time-series data.***Phys Rev A*1986,**34:**2427–2432. 10.1103/PhysRevA.34.2427View Article - Theiler J, Eubank S, Longtin A, Galdrikian B, Farmer JD:
**Testing for nonlinearity in time series: the method of surrogate data.***Physica D*1992,**58:**77–94. 10.1016/0167-2789(92)90102-SView Article - Burioka N, Cornelissen G, Halberg F, Kaplan DT, Suyama H, Sako T, Shimizu E:
**Approximate entropy of human respiratory movement during eye-closed waking and different sleep stages.***Chest*2003,**123:**80–86. 10.1378/chest.123.1.80View Article -
**Bipoac Owners Guide**2001. - Voss RF, Clarke J:
**1/f noise in music and speech.***Nature*1975,**258:**317–318.View Article - Fraser AM, Swinney HL:
**Independent coordinates for strange attractors from mutual information.***Phys Rev A*1986,**33:**1134–1140. 10.1103/PhysRevA.33.1134MathSciNetView Article - Fraser AM:
**Information and entropy in strange attractors.***IEEE transactions on Information Theory*1989,**35:**245–262. 10.1109/18.32121MathSciNetView Article - Grassberger P, Procaccia I:
**Characterization of strange attractors.***Phys Rev Lett*1983,**50:**346–349. 10.1103/PhysRevLett.50.346MathSciNetView Article - Grassberger P, Procaccia I:
**Measuring the strangeness of strange attractors.***Physica D*1983,**9:**189–208. 10.1016/0167-2789(83)90298-1MathSciNetView Article - Packard NH, Crutchfield JP, Farmer JD, Shaw RS:
**Geometry from a time series.***Phys Rev Lett*1980,**45:**712–716. 10.1103/PhysRevLett.45.712View Article - Takens F:
**Detecting Strange Attractors in turbulence.***In Dynamical Systems and Turbulence*Berlin: Springer 1981. - Sauer T, Yorke J, Casdagli M:
**Embedology.***Journal of Statistical Physics*1994,**65:**579–616.MathSciNetView Article -
**CDA User Manual**2001. - Theiler J:
**Spurious dimensions from correlation algorithms applied to limited time-series data.***Phys Rev A*1986,**34:**2427–2432. 10.1103/PhysRevA.34.2427View Article - Inouye T, Shinosaki K:
**Quantification of EEG irregularity by use of the entropy of the power spectrum.***Electroencephalogr Clin Neurophysiol*1991,**79**(3):204–210. 10.1016/0013-4694(91)90138-TView Article - Fell J, Roschke J:
**A comparison between spectral and nonlinear EEG measures.***Electroencephalogr Clin Neurophysiol*1996,**98**(5):401–410. 10.1016/0013-4694(96)95636-9View Article - Amari S:
**Differential-Geometrical methods in statistics.**New York: Springer-Verlag 1985.View Article - Grassberger P, Schrieber T:
**Nonlinear time sequence analysis.***Int J Bifurcation and Chaos*1991,**1**(3):512–547.View Article - Gonzalez Andino SL, de Peralta Menendez Grave:
**Measuring the complexity of time series: an application to neurophysiological signals.***Hum Brain Mapp*2000,**11**(1):46–57. Publisher Full Text 10.1002/1097-0193(200009)11:1<46::AID-HBM40>3.3.CO;2-XView Article - Quiroga RQ, Arnhold J, Lehnertz K, Grassberger P:
**Kulback-Leibler and renormalized entropies: applications to electroencephalograms of epilepsy patients.***Stat Phys Plasmas Fluids Relat Interdiscip Topics*2000,**62:**8380–8386. 10.1103/PhysRevE.62.8380 - Pincus SM:
**Approximate entropy as a measure of system complexity.***Proc Natl Acad Sci USA*1991,**88:**2297–2301.MathSciNetView Article - Pincus SM, Goldberger AL:
**Physiological time-series analysis: what does regularity quantify?***Am J Physiol*1994,**266:**H1643-H1656. - Bruhn J, Ropcke H:
**Approximate entropy as an electroencephalographic measure of anesthetic drug effect during desflurane anesthesia.***Anesthesiology*2000,**92**(3):715–26.View Article - Steyn-Ross ML, Steyn-Ross DA, Sleigh JW, Liley DT:
**Theoretical electroencephalogram stationary spectrum for a white-noise-driven cortex: evidence for a general anesthetic-induced phase transition.***Phys Rev E*1999,**60**(6):7299–7311. 10.1103/PhysRevE.60.7299View Article - Dangel S, Meier PF, Moser HR, Plibersek S, Shen Y:
**Time series analysis of sleep EEG.***Computer assisted Physics*1999, 93–95. - Wolf A, Swift JB, Swinney LH, Vastano JA:
**Determining Lyapunov exponent from a time series.***Physica D*1985,**16:**285–317. 10.1016/0167-2789(85)90011-9MathSciNetView Article - Das A, Das P, Roy AB:
**Applicability of Lyapunov Exponent in EEG Data Analysis.***Complexity International*2002. - Theiler J, Eubank S, Longtin A, Galdrikian B, Farmer JD:
**Testing for nonlinearity in time series: the method of surrogate data.***Physica D*1992,**58:**77–94. 10.1016/0167-2789(92)90102-SView Article - Bruce JW:
**Fractal Physiology and Chaos in Medicine.**World Scientific Press 1990. - Battacharya J:
**Complexity analysis of spontaneous EEG.***Acta Neurobiol Exp*2000,**60:**495–501. - Klonowski W, Jernajczyk W, Niedzielska K, Rydz A, Stepien R:
**Quantitative measure of complexity of EEG signal dynamics.***Acta Neurobiol Exp (Wars)*1999,**59:**315–321. - Palus M, Dvorak I, David I:
**Remarks on spatial and temporal dynamics of EEG.***In: Mathematical approaches to brain functioning diagnostics**(Edited by: Dvorak I, Holden AV).*Manchester University Press 1991, 369–385. - Soong AC, Stuart CI:
**Evidence of chaotic dynamics underlying the human alpha-rhythm electroencephalogram.***Biol Cybern*1989,**62:**55–62.View Article - Stam CJ, Pijn JPM, Suffczynski P, Lopes da Silva FH:
**Dynamics of the human alpha rhythm: evidence for non-linearity?***Clin Neurophysiol*1999,**110:**1801–1813. 10.1016/S1388-2457(99)00099-1View Article - Szelenberger W, Wackermann J, Skalski M, Niemcewicz S, Drojewski J:
**Analysis of complexity of EEG during sleep.***Acta Neurobiol Exp (Wars)*1996,**56:**165–169.

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.