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 biosignals 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 [68]. 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 [911]. 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 [1622]. Approximate Entropy of human respiratory movement during eyeclosed 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]. SilverSilver 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.
In the phasespace plot shown in Figure 1, Xaxis represents the EEG signal X[n] and the Yaxis 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.
Figure 1. Phase space plots for various mental states
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 dr 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
The accuracy of the nonlinear time series analysis lies in the selection of optimum embedding dimension [30]. Embedding theorems defined by Takens and Sauer [31,32] states that for a strange attractor of Fractal Dimension D, the embedding using time delay coordinates is one to one if m ≥ 2D+1 or m ≥ D_{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].
Figure 2. Correlation Dimension Vs Embedding Dimension
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 multiplecounting. The correlation dimension is then artificially low because all the points are temporally close to each other. This effect occurs when the time resolution Δtof 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 [3540] 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 frequencyspace or phasespace. 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 onedimensional signal in a 'phasespace'. 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 highdimensional information from the univariate EEG datastream. Using these embeddings, the theoretical measure of the rate of "information" generation by a system is the KolmogorovSinai 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 SteynRoss et al [44] the approximate entropy when applied to EEG signals reflects the intracortical 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
where is the correlation integral with embedding dimension m and time lag of 1.
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 selfsimilarity 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 longrange 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 selfsimilarity. 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 nonstationary 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 timeseries is similar to a random walk. If H < 0.5, the timeseries cover less "distance" than a random walk. But if H > 0.5, the timeseries covers more "distance" than a random walk. H is related to the fractal dimension D given by
H = E+1D
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 + (m1)t}
We locate nearest neighbor to initial point
{x(t0), x(t0+ t),..., x(t0 + (m1)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 k1) 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
The phase space plot show the unique pattern for each mental state (Figure 1(a),1(b),1(c) and 1(d)). The correlation dimension parameter D_{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 longterm 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 selfsimilarity 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.
Figure 3. The variation of CD for various mental states
Figure 4. The variation of LLE for various mental states
Figure 5. The variation of ApEn for various mental states
Figure 6. The variation of HE for various mental states
Table 1. Range of nonlinear parameter variation for different mental states
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 [5055]. 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+1D, 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.
Authors' Contributions
Thelma T and Fadilah A carried out the software implementation. Kannathal N and Rajendra U Acharya carried out the design, analysis and coordination of the work. Sadasivan K Puthusserypady conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.
References

Stein KM, Lippman N, Kligfield P: Fractal rhythms of the heart.
Journal of Electrocardiology 1992, 24:7276. PubMed Abstract

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):3139. PubMed Abstract  Publisher Full Text

Duke W, Pritchard WS, Krieble KK: Dimensional analysis of resting human EEG II: Surrogate data testing indicates nonlinearity but not lowdimensional chaos.
Psychophysiology 1995, 32:486491. PubMed Abstract

Boccaletti S, Grebogi C, Lai YC, Mancini H, Mazaet D: The control of chaos: Theory and applications.
Physics Reports 2000, 329:108109. Publisher Full Text

Glass L, Michel RG, Mackey M, Shrier A: Chaos in neurobiology.

Jaeseung J, JeongHo C, Kim SY, SeolHeui H: Nonlinear dynamical analysis of the EEG in patients with Alzheimer's disease and vacular dementia.

Philippe F, Henri K: Is there chaos in the brain? Concepts of nonlinear dynamics and methods of investigation.
Life Sciences 2001, 324:773793. PubMed Abstract  Publisher Full Text

Babloyantz A: Evidence of chaotic dynamics during sleep cycle. In Dimensions and entropies in Chaotic system. Edited by MayerKress G. Berlin: SpringerVerlag; 1986.

Bessar E: Biophysical and physiological systems analysis. London: AddisonWesley; 1960.

Mayer G, Layne SC: Dimensionality of the human encephalogram. In perspectives in Biological Dynamics and Theoretical Medicine.
In Annals NY Academic Science Edited by Koslow SH, Madell AJ, Shlessinger MF. 1987, 504507.

Rapp PE, Bashore T, Martinerie J, Albano A, Zimmerman I, Mess A: Dynamics of brain electrical activity.
Brain Topography 1989, 2:99118. PubMed Abstract

Rapp PE: Chaos in the neurosciences: cautionary tales from the frontier.

Freeman WJ: A proposed name for aperiodic brain activity: stochastic chaos.
Neural Networks 2000, 13:1113. PubMed Abstract  Publisher Full Text

Wright JJ, Liley DTJ: Dynamics of the brain at global and microscopic scales. Neural networks and the EEG.

Pijn JPM: Quantitative evaluation of EEG signals in epilepsy: nonlinear association time delays and nonlinear dynamics. In 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:249270. PubMed Abstract

Rombouts SARB, Keunen RWM, Stam CJ: Investigation of nonlinear structure in multichannel EEG.
Phys Lett A 1995, 202:352358. Publisher Full Text

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):149153. PubMed Abstract  Publisher Full Text

Bradley E: Timeseries analysis. Intelligent data analysis: An introduction. Edited by Berthold M, Hand D. Berlin: SpringerVerlag; 1999.

Theiler J: Spurious dimension from correlation algorithms applied to limited timeseries data.
Phys Rev A 1986, 34:24272432. PubMed Abstract  Publisher Full Text

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:7794. Publisher Full Text

Burioka N, Cornelissen G, Halberg F, Kaplan DT, Suyama H, Sako T, Shimizu E: Approximate entropy of human respiratory movement during eyeclosed waking and different sleep stages.
Chest 2003, 123:8086. PubMed Abstract  Publisher Full Text

2001.

Fraser AM, Swinney HL: Independent coordinates for strange attractors from mutual information.
Phys Rev A 1986, 33:11341140. PubMed Abstract  Publisher Full Text

Fraser AM: Information and entropy in strange attractors.
IEEE transactions on Information Theory 1989, 35:245262. Publisher Full Text

Grassberger P, Procaccia I: Characterization of strange attractors.
Phys Rev Lett 1983, 50:346349. Publisher Full Text

Grassberger P, Procaccia I: Measuring the strangeness of strange attractors.
Physica D 1983, 9:189208. Publisher Full Text

Packard NH, Crutchfield JP, Farmer JD, Shaw RS: Geometry from a time series.
Phys Rev Lett 1980, 45:712716. Publisher Full Text

Takens F: Detecting Strange Attractors in turbulence. In In Dynamical Systems and Turbulence. Berlin: Springer; 1981.

2001.

Theiler J: Spurious dimensions from correlation algorithms applied to limited timeseries data.
Phys Rev A 1986, 34:24272432. PubMed Abstract  Publisher Full Text

Inouye T, Shinosaki K: Quantification of EEG irregularity by use of the entropy of the power spectrum.
Electroencephalogr Clin Neurophysiol 1991, 79(3):204210. PubMed Abstract  Publisher Full Text

Fell J, Roschke J: A comparison between spectral and nonlinear EEG measures.
Electroencephalogr Clin Neurophysiol 1996, 98(5):401410. PubMed Abstract  Publisher Full Text

Amari S: DifferentialGeometrical methods in statistics. New York: SpringerVerlag; 1985.

Grassberger P, Schrieber T: Nonlinear time sequence analysis.

Gonzalez Andino SL, de Peralta Menendez Grave: Measuring the complexity of time series: an application to neurophysiological signals.
Hum Brain Mapp 2000, 11(1):4657. PubMed Abstract  Publisher Full Text

Quiroga RQ, Arnhold J, Lehnertz K, Grassberger P: KulbackLeibler and renormalized entropies: applications to electroencephalograms of epilepsy patients. In Stat Phys Plasmas Fluids Relat Interdiscip Topics. Volume 62. ; 2000::83808386. Publisher Full Text

Pincus SM: Approximate entropy as a measure of system complexity.
Proc Natl Acad Sci USA 1991, 88:22972301. PubMed Abstract  Publisher Full Text

Pincus SM, Goldberger AL: Physiological timeseries analysis: what does regularity quantify?
Am J Physiol 1994, 266:H1643H1656. PubMed Abstract

Bruhn J, Ropcke H: Approximate entropy as an electroencephalographic measure of anesthetic drug effect during desflurane anesthesia.
Anesthesiology 2000, 92(3):71526. PubMed Abstract  Publisher Full Text

SteynRoss ML, SteynRoss DA, Sleigh JW, Liley DT: Theoretical electroencephalogram stationary spectrum for a whitenoisedriven cortex: evidence for a general anestheticinduced phase transition.
Phys Rev E 1999, 60(6):72997311. Publisher Full Text

Dangel S, Meier PF, Moser HR, Plibersek S, Shen Y: Time series analysis of sleep EEG.

Wolf A, Swift JB, Swinney LH, Vastano JA: Determining Lyapunov exponent from a time series.
Physica D 1985, 16:285317. Publisher Full Text

Das A, Das P, Roy AB: Applicability of Lyapunov Exponent in EEG Data Analysis.

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:7794. Publisher Full Text

Bruce JW: Fractal Physiology and Chaos in Medicine. World Scientific Press; 1990.

Klonowski W, Jernajczyk W, Niedzielska K, Rydz A, Stepien R: Quantitative measure of complexity of EEG signal dynamics.
Acta Neurobiol Exp (Wars) 1999, 59:315321. PubMed Abstract

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:369385.

Soong AC, Stuart CI: Evidence of chaotic dynamics underlying the human alpharhythm electroencephalogram.
Biol Cybern 1989, 62:5562. PubMed Abstract

Stam CJ, Pijn JPM, Suffczynski P, Lopes da Silva FH: Dynamics of the human alpha rhythm: evidence for nonlinearity?
Clin Neurophysiol 1999, 110:18011813. PubMed Abstract  Publisher Full Text

Szelenberger W, Wackermann J, Skalski M, Niemcewicz S, Drojewski J: Analysis of complexity of EEG during sleep.
Acta Neurobiol Exp (Wars) 1996, 56:165169. PubMed Abstract