The dimensionality of human's electroencephalogram during sleep

In order to perform an analysis of nonlinear EEG-dynamics we investigated the EEG of ten male probands during sleep. According to Rechtschaffen and Kales (1968) we scored the sleep-EEG and applied an algorithm, proposed by Grassberger and Proccaccia (1983) to compute the correlation dimension of different sleep stages. The correlation dimension characterizes the dynamics of the EEG signal and estimates the degrees of freedom of the signal under study. We could demonstrate, that the EEG of slow wave sleep stages depicts a dimensionality, which is two units smaller than that of light or REM sleep.     

高维混沌降维的投影追踪主分量分析法

一些生理信号,例如脑电是源自于高维混沌系统,因此低维混沌理论和方法不适用于分析这类高维混沌。采用投影追踪主分量分析法（Princiopal Component Analysis based on Projection Pursuit,PP PCA）对高维Lorenz模型系统进行了降维的研究。在用上述方法成功地对线性和非线性噪声－周期模型分别进行了PP PCA分析的基础上,对Lorenz高维混沌系统进行了PPPCA降维的研究。结果表明,正确选用非线性的投影追踪主分量分析法,可以通过简化原系统达到降维的目的,并能保留研究所关心的原系统的主要动态特性。同时也阐明了方法的稳定性和将该方法应用于高维脑电降维的可行性。     

Dual Dimensionality Reduction Reveals Independent Encoding of Motor Features in a Muscle Synergy for Insect Flight Control

What are the features of movement encoded by changing motor commands? Do motor commands encode movement independently or can they be represented in a reduced set of signals (i.e. synergies)? Motor encoding poses a computational and practical challenge because many muscles typically drive movement, and simultaneous electrophysiology recordings of all motor commands are typically not available. Moreover, during a single locomotor period (a stride or wingstroke) the variation in movement may have high dimensionality, even if only a few discrete signals activate the muscles. Here, we apply the method of partial least squares (PLS) to extract the encoded features of movement based on the cross-covariance of motor signals and movement. PLS simultaneously decomposes both datasets and identifies only the variation in movement that relates to the specific muscles of interest. We use this approach to explore how the main downstroke flight muscles of an insect, the hawkmoth Manduca sexta, encode torque during yaw turns. We simultaneously record muscle activity and turning torque in tethered flying moths experiencing wide-field visual stimuli. We ask whether this pair of muscles acts as a muscle synergy (a single linear combination of activity) consistent with their hypothesized function of producing a left-right power differential. Alternatively, each muscle might individually encode variation in movement. We show that PLS feature analysis produces an efficient reduction of dimensionality in torque variation within a wingstroke. At first, the two muscles appear to behave as a synergy when we consider only their wingstroke-averaged torque. However, when we consider the PLS features, the muscles reveal independent encoding of torque. Using these features we can predictably reconstruct the variation in torque corresponding to changes in muscle activation. PLS-based feature analysis provides a general two-sided dimensionality reduction that reveals encoding in high dimensional sensory or motor transformations.     

Deterministic chaos and the first positive Lyapunov exponent: a nonlinear analysis of the human electroencephalogram during sleep

Under selected conditions, nonlinear dynamical systems, which can be described by deterministic models, are able to generate so-called deterministic chaos. In this case the dynamics show a sensitive dependence on initial conditions, which means that different states of a system, being arbitrarily close initially, will become macroscopically separated for sufficiently long times. In this sense, the unpredictability of the EEG might be a basic phenomenon of its chaotic character. Recent investigations of the dimensionality of EEG attractors in phase space have led to the assumption that the EEG can be regarded as a deterministic process which should not be mistaken for simple noise. The calculation of dimensionality estimates the degrees of freedom of a signal. Nevertheless, it is difficult to decide from this kind of analysis whether a process is quasiperiodic or chaotic. Therefore, we performed a new analysis by calculating the first positive Lyapunov exponent L ₁ from sleep EEG data. Lyapunov exponents measure the mean exponential expansion or contraction of a flow in phase space. L ₁ is zero for periodic as well as quasiperiodic processes, but positive in the case of chaotic processes expressing the sensitive dependence on initial conditions. We calculated L ₁ for sleep EEG segments of 15 healthy men corresponding to the sleep stages I, II, III, IV, and REM (according to Rechtschaffen and Kales). Our investigations support the assumption that EEG signals are neither quasiperiodic waves nor a simple noise. Moreover, we found statistically significant differences between the values of L ₁ for different sleep stages. All together, this kind of analysis yields a useful extension of the characterization of EEG signals in terms of nonlinear dynamical system theory.     

复杂度脑电地形图研究

脑电地形图是近年脑电分析的热点之一。通过对各种复杂度算法的分析得出,近似熵由于所需要的时间序列长度较短,大大减少了脑电非平稳性所带来的困难,且无需粗粒化,在对生物医学信号的复杂度分析中有其一定的优点,采用近似熵对多道脑电信号的复杂度运算结果,通过空间插值,构建复杂性动态脑地形图,以便于观察大脑各部ＥＥＧ信号复杂度在同一时刻的相对强弱关系和这种关系随时间的变化。并通过对一些脑疾病患者脑电数据的分析,     

对脑电相关维数计算中有关参数的探讨

应用Ｇｒａｓｓｂｅｒｇｅｒ－Ｐｒｏｃａｃｉａ ａｌｇｏｒｉｔｈｍ算法,通过点列长度、时间延迟、嵌入维数等参数的变化,对人的８导脑电时间序列的相关维数值的影响进行了研究。结果表明不仅参数的选用是否合理对计算值影响很大,而且脑电信号不像是一个稳定的混沌吸引子,因此拟用相关维数作用一项客观指标来准确地表示大脑的不同功能状态,几乎是不可能的。     

SPHARA - A Generalized Spatial Fourier Analysis for Multi-Sensor Systems with Non-Uniformly Arranged Sensors: Application to EEG

Important requirements for the analysis of multichannel EEG data are efficient techniques for signal enhancement, signal decomposition, feature extraction, and dimensionality reduction. We propose a new approach for spatial harmonic analysis (SPHARA) that extends the classical spatial Fourier analysis to EEG sensors positioned non-uniformly on the surface of the head. The proposed method is based on the eigenanalysis of the discrete Laplace-Beltrami operator defined on a triangular mesh. We present several ways to discretize the continuous Laplace-Beltrami operator and compare the properties of the resulting basis functions computed using these discretization methods. We apply SPHARA to somatosensory evoked potential data from eleven volunteers and demonstrate the ability of the method for spatial data decomposition, dimensionality reduction and noise suppression. When employing SPHARA for dimensionality reduction, a significantly more compact representation can be achieved using the FEM approach, compared to the other discretization methods. Using FEM, to recover 95% and 99% of the total energy of the EEG data, on average only 35% and 58% of the coefficients are necessary. The capability of SPHARA for noise suppression is shown using artificial data. We conclude that SPHARA can be used for spatial harmonic analysis of multi-sensor data at arbitrary positions and can be utilized in a variety of other applications.     

Nonlinear analysis of EEG signals: Surrogate data analysis

ObjectivesThe electroencephalogram (EEG) signal contains information about the state and condition of the brain. The aim of the study is to conduct a nonlinear analysis of the EEG signals and to compare the differences in the nonlinear characteristics of the EEG during normal state and the epileptic state.DataThe EEG data used for this study – which consisted of epileptic EEG and normal EEG – were obtained from the EEG database available with the Bonn University, Germany.ResultsThe attractors seen in normal and epileptic human brain dynamics were studied and compared. Surrogate data analyses were conducted on two nonlinear measures, namely the largest Lyapunov exponent and the correlation dimension, to test the hypothesis whether EEG signals were in accordance with linear stochastic models.DiscussionsThe existence of deterministic chaos in brain activity is confirmed by the existence of a chaotic attractor; also, saturation of the correlation dimension towards a definite value is the manifestation of a deterministic dynamics. Also a reduction is observed between the dimensionalities of the brain attractors from normal state to the epileptic state. The evaluation of the largest Lyapunov exponent also confirms the lowering of complexity during an episode of seizure.ConclusionIn case of Lyapunov exponent of EEG data, the change due to surrogating is small suggesting that it is not representing the system complexity properly but there is a marked change in the case of correlation dimension value due to surrogating.     

睡眠脑电的分形维数分析

为了研究是否用维数来描述脑功能的复杂程序,我们计算了５例正常人整夜睡眠时不同阶段的ＥＥＧ资料维数。结果显示,随着睡眠深度的增加而维数降低。这一事实和生理概念及其他作者所报道的结果相一致,故我们考虑维数分析在表示脑功能状态上有一定的理论意义。如果我们用简化法来计算维数,此法就有可能实时地应用于临床诊断。     

Fractal dimensions of pacing and grip force in drawing and handwriting production

We performed a repeated measures experiment to show that the pacing of a cyclic, ballistic drawing task has a fractal dimension. We also estimated the dimensionality of the force used to grip the drawing implement. Finally, we present an analysis of pediatric data to show that grip force has a fractal dimension in an actual handwriting task. In our experiment, subjects drew circles of varying sizes and at varying rates on a digitizing tablet, using a pen instrumented to measure radial force applied to its barrel. Subjects also drew circles in synchrony with a metronome. We found strong evidence for fractal scaling of both drawing period and grip force in the circle-drawing study. The dimensionality ranged from fractal Gaussian noise (fGn) to fractal Brownian motion, with Hurst coefficients clustering around the value for 1/f noise. When the subjects were required to synchronize their drawing with a metronome, the Hurst coefficient for the drawing period decreased, while the coefficient for grip force did not. This result indicates that independent processes control the variations in pacing and grip force. Grip force in the handwriting study also displayed fractal properties, with Hurst coefficients in the range of correlated fGn. We draw parallels between our handwriting measurements and studies of human gait.     

Coupling of cardiac and locomotor rhythms

The pressure within exercising skeletal muscle rises and falls rhythmically during normal human locomotion, the peak pressure reaching levels that intermittently impede blood flow to the exercising muscle. Speculating that a reciprocal relationship between the timing of peak intramuscular and pulsatile arterial pressures should optimize blood flow through muscle and minimize cardiac load, we tested the hypothesis that heart rate becomes entrained with walking and running cadence at some locomotion speeds, by means of electrocardiography and an accelerometer to provide signals reflecting heart rate and cadence, respectively. In 18 of 25 subjects, 1:1 coupling of heart and step rates was present at one or more speeds on a motorized treadmill, generally at moderate to high exercise intensities. To determine how exercise specific this phenomenon is, and to refute the competing hypothesis that coupling is due to vertical accelerations of the heart during locomotion, we had 12 other subjects cycle on an electronically braked bicycle ergometer. Coupling was found between heart rate and pedaling frequency in 10 of them. Cardiac-locomotor coupling appears to be a normal physiological phenomenon, and its identification provides a fresh perspective from which to study endurance.     

Rates of diffusion controlled reactions in one, two and three dimensions

The dimensionality of diffusion may markedly affect the rate and economy of diffusion controlled reactions. Moreover, the degree of dependence of the steady state rate of these reactions on the concentration of each of the two reacting species is also dictated by the dimensionality and it ranges from linear dependence in the three dimensional case to a nearly square dependence in the one dimensional case. These theoretical observations emerge from a direct analysis of the steady state diffusion controlled rates which are derived here using a simple straightforward approach. This approach is based on the conjecture that in the steady state the rate of diffusional encounters between the two reaction partners equals to the sum of the encounter rates of two independent processes which are obtained by alternately immobilizing one of the reaction partners while the other partner diffuses freely. Unlike Smoluchowski's classical approach, the presented point of view permits to obtain in a unified fashion reaction rates for all dimensionalities.     

Investigation of the dynamics underlying periodic complexes in the EEG

Periodic complexes (PC), occurring lateralised or diffuse, are relatively rare EEG phenomena which reflect acute severe brain disease. The pathophysiology is still incompletely understood. One hypothesis suggested by the alpha rhythm model of Lopes da Silva is that periodic complexes reflect limit cycle dynamics of cortical networks caused by excessive excitatory feedback. We examined this hypothesis by applying a recently developed technique to EEGs displaying periodic complexes and to periodic complexes generated by the model. The technique, non-linear cross prediction, characterises how well a time series can be predicted, and how much amplitude and time asymmetry is present. Amplitude and time asymmetry are indications of non-linearity. In accordance with the model, most EEG channels with PC showed clear evidence of amplitude and time asymmetry, pointing to non-linear dynamics. However, the non-linear predictability of true PC was substantially lower than that of PC generated by the model. Furthermore, no finite value for the correlation dimension could be obtained for the real EEG data, whereas the model time series had a dimension slighter higher than one, consistent with a limit cycle attractor. Thus we can conclude that PC reflect non-linear dynamics, but a limit cycle attractor is too simple an explanation. The possibility of more complex (high dimensional and spatio-temporal) non-linear dynamics should be investigated. Received: 26 February 1998 / Accepted in revised form: 24 August 1998     

Nonlinearity in normal human EEG: cycles,temporal asymmetry,nonstationarity and randomness,not chaos

Two-hour vigilance and sleep electroencephalogram (EEG) recordings from five healthy volunteers were analyzed using a method for identifying nonlinearity and chaos which combines the redundancy–linear redundancy approach with the surrogate data technique. A nonlinear component in the EEG was detected, however, inconsistent with the hypothesis of low-dimensional chaos. A possibility that a temporally asymmetric process may underlie or influence the EEG dynamics was indicated. A process that merges nonstationary nonlinear deterministic oscillations with randomness is proposed for an explanation of observed properties of the analyzed EEG signals. Taking these results into consideration, the use of dimensional and related chaos-based algorithms in quantitative EEG analysis is critically discussed. Received: 25 September 1994 / Accepted in revised form: 10 July 1996     

Correlation dimension based lossless compression of EEG signals

Transmission of long duration EEG signals without loss of information is essential for telemedicine based applications. In this work, a lossless compression scheme for EEG signals based on neural network predictors using the concept of correlation dimension (CD) is proposed. EEG signals which are considered as irregular time series of chaotic processes can be characterized by the non-linear dynamic parameter CD which is a measure of the correlation among the EEG samples. The EEG samples are first divided into segments of 1 s duration and for each segment, the value of CD is calculated. Blocks of EEG samples are then constructed such that each block contains segments with closer CD values. By arranging the EEG samples in this fashion, the accuracy of the predictor is improved as it makes use of highly correlated samples. As a result, the magnitude of the prediction error decreases leading to less number of bits for transmission. Experiments are conducted using EEG signals recorded under different physiological conditions. Different neural network predictors as well as classical predictors are considered. Experimental results show that the proposed CD based preprocessing scheme improves the compression performance of the predictors significantly.     

Multivariate and multiorgan analysis of cardiorespiratory variability signals: the CAP sleep case.

Signals from different systems are analyzed during sleep on a beat-to-beat basis to provide a quantitative measure of synchronization with the heart rate variability (HRV) signal, oscillations of which reflect the action of the autonomic nervous system. Beat-to-beat variability signals synchronized to QRS occurrence on ECG signals were extracted from respiration, electroencephalogram (EEG) and electromyogram (EMG) traces. The analysis was restricted to sleep stage 2. Cyclic alternating pattern (CAP) periods were detected from EEG signals and the following conditions were identified: stage 2 non-CAP (2 NCAP), stage 2 CAP (2 CAP) and stage 2 CAP with myoclonus (2 CAP MC). The coupling relationships between pairs of variability signals were studied in both the time and frequency domains. Passing from 2 NCAP to 2 CAP, sympathetic activation is indicated by tachycardia and reduced respiratory arrhythmia in the heart rate signal. At the same time, we observed a marked link between EEG and HRV at the CAP frequency. During 2 CAP MC, the increased synchronization involved myoclonus and respiration. The underlying mechanism seems to be related to a global control system at the central level that involves the different systems.     

A New Method for Nonlinear Dynamic Analysis of the Multi-kinetics Neural Mass Model

ObjectiveThe aim was to use a new method to analyze the nonlinear dynamic characteristics of the multi-kinetics neural mass model. We hope that this new method can be as an auxiliary judgment tool for the diagnosis of brain diseases and the identification of brain activity states.MethodsWe apply the Lorenz plot to analyze the nonlinear dynamic characteristics of electroencephalogram (EEG) signals from the multi-kinetics neural mass models. The standard deviations in two orthogonal directions of the Lorenz plot are further used to quantify the nonlinear dynamic characteristics of EEG signals.ResultsThe results show that the normalized signal frequency power spectrum may not be able to distinguish normal EEG signals and epileptiform spikes, but the Lorenz plot can distinguish the normal EEG signals and epileptiform spikes effectively. For EEG signals with multi-rhythms, the Lorenz plot of all the simulated signals are oval, but the value of SD1/SD2 increases monotonically when the multi-rhythm EEG signals change from low frequency to high frequency.ConclusionThe Lorenz plot of EEG signals with different rhythms presents different distribution. It is an effective nonlinear analysis method for EEG signals.     

Fundamental patterns of bilateral muscle activity in human locomotion

Human gait is characterized by smooth, regular and repeating movements but the control system is complex: there are many more actuators (i.e. muscles) than degrees of freedom in the system. Statistical pattern-recognition techniques have been applied to examine muscle activity signals, but these have all concentrated exclusively on unilateral gait. We report here the application of factor analysis to the electromyographic patterns of 16 muscles (eight bilateral pairs) in ten normal subjects. Consistent with our prior work, we have established two factors, named loading response and propulsion, which correspond with important phases in the gait cycle. In addition, we have also discovered a third factor, which we have named the coordinating factor, that maintains the phase shift between the left and right sides. These findings suggest that the central nervous system solves the problem of high dimensionality by generating a few fundamental signals which control the major muscle groups in both legs.