Chaotic behavior in the locomotion ofamoeba proteus

Summary The locomotion ofAmoeba proteus has been investigated by algorithms evaluating correlation dimension and Lyapunov spectrum developed in the field of nonlinear science. It is presumed by these parameters whether the random behavior of the system is stochastic or deterministic. For the analysis of the nonlinear parameters, n-dimensional time-delayed vectors have been reconstructed from a time series of periphery and area ofA. proteus images captured with a charge-coupled-device camera, which characterize its random motion. The correlation dimension analyzed has shown the random motion ofA. proteus is subjected only to 3–4 macrovariables, though the system is a complex system composed of many degrees of freedom. Furthermore, the analysis of the Lyapunov spectrum has shown its largest exponent takes positive values. These results indicate the random behavior ofA. proteus is chaotic and deterministic motion on an attractor with low dimension. It may be important for the elucidation of the cell locomotion to take account of nonlinear interactions among a small number of dynamics such as the sol-gel transformation, the cytoplasmic streaming, and the relating chemical reaction occurring in the cell.     

Nonlinear analysis of continuous ECG during sleep I. Reconstruction

In recent years evidence has accumulated that ECG signals are of a nonlinear nature. It has been recognized that strictly periodic cardiac rhythms are not accompanied by healthy conditions but, on the contrary, by pathological states. Therefore, the application of methods from nonlinear system theory for the analysis of ECG signals has gained increasing interest. Crucial for the application of nonlinear methods is the reconstruction (embedding) of the time series in a phase space with appropriate dimension. In this study continuous ECG signals of 12 healthy subjects recorded during different sleep stages were analysed. Proper embedding dimension was determined by application of two techniques – the false nearest neighbours method and the saturation of the correlation dimension. Results for the ECG signals were compared with findings for simulated data (quasiperiodic dynamics, Lorenz data, white noise) and for phase randomized surrogates. Findings obtained with the two approaches suggest that embedding dimensions from 6 to 8 may be regarded as suitable for the topologically proper reconstruction of ECG signals. Received: 7 June 1999 / Accepted in revised form: 10 December 1999     

Detection of chaos: New approach to atmospheric pollen time-series analysis

Pollen and spores are biological particles that are ubiquitous to the atmosphere and are pathologically significant, causing plant diseases and inhalant allergies. One of the main objectives of aerobiological surveys is forecasting. Prediction models are required in order to apply aerobiological knowledge to medical or agricultural practice; a necessary condition of these models is not to be chaotic. The existence of chaos is detected through the analysis of a time series. The time series comprises hourly counts of atmospheric pollen grains obtained using a Burkard spore trap from 1987 to 1989 at Mar del Plata. Abraham's method to obtain the correlation dimension was applied. A low and fractal dimension shows chaotic dynamics. The predictability of models for atomspheric pollen forecasting is discussed.     

On the correlation dimension of optokinetic nystagmus eye movements: computational parameters, filtering, nonstationarity, and surrogate data

 We discuss the estimation of the correlation dimension of optokinetic nystagmus (OKN), a type of reflexive eye movement. Parameters of the time-delay reconstruction of the attractor are investigated, including the number of data points, the time delay, the window duration, and the duration of the signal being analyzed. Adequate values are recommended. Digital low-pass filtering causes the dimension to increase as the filter cutoff frequency decreases, in accord with a previously published prediction. The stationarity of the correlation dimension is examined; the dimension appears to decrease over the course of 120 s of continuous stimulation. Implications for the reliable estimation of the dimension are considered. Several surrogate data sets are constructed, based on both early (0–30 s) and late (100–130 s) OKN segments. Most of the surrogate data sets randomize some aspect of the original OKN, while maintaining other aspects. Dimensions are found for all surrogates and for the original OKN. Evidence is found that is consistent with some amount of deterministic and nonlinear dynamics in OKN. When this structure is randomized in the surrogate, the dimension changes or the dimension algorithm ceases to converge to a finite value. Implications for further analysis and modeling of OKN are discussed. Received: 30 August 1996/Accepted in revised form: 13 November 1996     

Complex nonlinear dynamics of the Hodgkin-Huxley equations induced by time scale changes

 The Hodgkin–Huxley equations with a slight modification are investigated, in which the inactivation process (h) of sodium channels or the activation process of potassium channels (n) is slowed down. We show that the equations produce a variety of action potential waveforms ranging from a plateau potential, such as in heart muscle cells, to chaotic bursting firings. When h is slowed down – differently from the case of n variable being slow – chaotic bursting oscillations are observed for a wide range of parameter values although both variables cause a decrease in the membrane potential. The underlying nonlinear dynamics of various action potentials are analyzed using bifurcation theory and a so-called slow–fast decomposition analysis. It is shown that a simple topological property of the equilibrium curves of slow and fast subsystems is essential to the production of chaotic oscillations, and this is the cause of the large difference in global firing characteristics between the h-slow and n-slow cases. Received: 9 August 2000 / Accepted in revised form: 10 January 2001     

Nonlinear analysis of continuous ECG during sleep II. Dynamical measures

The hypothesis that cardiac rhythms are associated with chaotic dynamics implicating a healthy flexibility has motivated the investigation of continuous ECG with methods of nonlinear system theory. Sleep is known to be associated with modulations of the sympathetic and parasympathetic control of cardiac dynamics. Thus, the differentiation of ECG signals recorded during different sleep stages can serve to determine the usefulness of nonlinear measures in discriminating ECG states in general. For this purpose the following six nonlinear measures were implemented: correlation dimension D2, Lyapunov exponent L1. Kolmogorov entropy K2, as well as three measures derived from the analysis of unstable periodic orbits. Results of this study show that continuous ECG signals can be differentiated from linear stochastic surrogates by each of the nonlinear measures. The most significant finding with respect to the sleep-related differentiation of ECG signals is an increase in dominant chaoticity assessed by L1 and a reduction in the degrees of freedom estimated by D2 during REM sleep compared to slow wave sleep. Our findings suggest that the increase in dominant chaoticity during REM sleep with regard to time-continuous nonlinear analysis is comparable to an increased heart rate variability. The reduction in the correlation dimension may be interpreted as an expression of the withdrawal of respiratory influences during REM sleep. Received: 7 June 1999 / Accepted in revised form: 10 December 1999     

Fractal dimension of electroencephalographic time series and underlying brain processes

 Fractal dimension has been proposed as a useful measure for the characterization of electrophysiological time series. This paper investigates what the pointwise dimension of electroencephalographic (EEG) time series can reveal about underlying neuronal generators. The following theoretical assumptions concerning brain function were made (i) within the cortex, strongly coupled neural assemblies exist which oscillate at certain frequencies when they are active, (ii) several such assemblies can oscillate at a time, and (iii) activity flow between assemblies is minimal. If these assumptions are made, cortical activity can be considered as the weighted sum of a finite number of oscillations (plus noise). It is shown that the correlation dimension of finite time series generated by multiple oscillators increases monotonically with the number of oscillators. Furthermore, it is shown that a reliable estimate of the pointwise dimension of the raw EEG signal can be calculated from a time series as short as a few seconds. These results indicate that (i) The pointwise dimension of the EEG allows conclusions regarding the number of independently oscillating networks in the cortex, and (ii) a reliable estimate of the pointwise dimension of the EEG is possible on the basis of short raw signals. Received: 1 September 1994/Accepted in revised form: 16 May 1995     

Quantification of emotion by nonlinear analysis of the chaotic dynamics of electroencephalograms during perception of 1/f music

The goal of this study is to quantify and determine the way in which the emotional response to music is reflected in the electrical activities of the brain. When the power spectrum of sequences of musical notes is inversely proportional to the frequency on a log-log plot, we call it 1/f music. According to previous research, most listeners agree that 1/f music is much more pleasing than white (1/f ⁰) or brown (1/f ²) music. Based on these studies, we used nonlinear methods to investigate the chaotic dynamics of electroencephalograms (EEGs) elicited by computer-generated 1/f music, white music, and brown music. In this analysis, we used the correlation dimension and the largest Lyapunov exponent as measures of complexity and chaos. We developed a new method that is strikingly faster and more accurate than other algorithms for calculating the nonlinear invariant measures from limited noisy data. At the right temporal lobe, 1/f music elicited lower values of both the correlation dimension and the largest Lyapunov exponent than white or brown music. We observed that brains which feel more pleased show decreased chaotic electrophysiological behavior. By observing that the nonlinear invariant measures for the 1/f distribution of the rhythm with the melody kept constant are lower than those for the 1/f distribution of melody with the rhythm kept constant, we could conclude that the rhythm variations contribute much more to a pleasing response to music than the melody variations do. These results support the assumption that chaos plays an important role in brain function, especially emotion. Received: 30 December 1996 / Accepted in revised form: 18 December 1997     

基于混沌降噪的神经元放电峰峰间期序列分析

围绕如何来消除神经元峰峰间期序列中随机噪声影响从而提取出决定不规则性的确定性动力学关系这个问题,本文首先简要介绍峰峰间期序列样本的制备,然后着重讨论一个简单可行的混沌时间序列降噪方法的原理和算法实现,最终将该方法运用到神经元放电活动数值模拟和实验记录到的峰峰间期时间序列样本分析中。本文分析结果再次证明神经放电活动中确实存在着不规则混沌运动,而且降噪结果进一步揭示了神经电生理实验中决定混沌放电的不连续但分段光滑的单峰函数关系     

基于复杂性度量的心率变异信号非线性分析

假设心率变异信号是累积－发放模型（Ｉｎｔｅｇｒａｔｅ－ｆｉｒｅ）与非线性动力学系统耦合产生的峰电位链（ＳｐｉｋｅＴｒａｉｎ）。以符号动力学为基础,提出利用峰电位间隔（ｉｎｔｅｒｓｐｉｋｅｉｎｔｅｒｖａｌ,ＩＳＩ）及其随机替代数据的Ｃ１、Ｃ２复杂度来定量刻划非线性动力学系统特性。结果表明：确定性驱动产生的峰电位间隔序列可以与随机性驱动产生的峰电位间隔序列区分开。因此,在噪声干扰较强的生理信号中,尤其是在不清楚非线性动力系统变量和峰电位间隔序列之间是否存在微分同胚的情况下,以复杂性度量来代替以Ｔａｋｅｎｓ嵌入定理为基础的关联维数、Ｌｙａｐｎｏｖ指数等描述动力系统特征的方法是合适的。最后通过２类共３７个个体,每个个体的心电数据为１０００个Ｒ－Ｒ间期的微分序列检验心率变异信号的非线性结构。     

Analysis of cedar pollen time series: no evidence of low-dimensional chaotic behavior

Much of the current interest in pollen time series analysis is motivated by the possibility that pollen series arise from low-dimensional chaotic systems. If this is the case, short-range prediction using nonlinear modeling is justified and would produce high-quality forecasts that could be useful in providing pollen alerts to allergy sufferers. To date, contradictory reports about the characterization of the dynamics of pollen series can be found in the literature. Pollen series have been alternatively described as featuring and not featuring deterministic chaotic behavior. We showed that the choice of test for detection of deterministic chaos in pollen series is difficult because pollen series exhibit power spectra. This is a characteristic that is also produced by colored noise series, which mimic deterministic chaos in most tests. We proposed to apply the Ikeguchi–Aihara test to properly detect the presence of deterministic chaos in pollen series. We examined the dynamics of cedar (Cryptomeria japonica) hourly pollen series by means of the Ikeguchi–Aihara test and concluded that these pollen series cannot be described as low-dimensional deterministic chaos. Therefore, the application of low-dimensional chaotic deterministic models to the prediction of short-range pollen concentration will not result in high-accuracy pollen forecasts even though these models may provide useful forecasts for certain applications. We believe that our conclusion can be generalized to pollen series from other wind-pollinated plant species, as wind speed, the forcing parameter of the pollen emission and transport, is best described as a nondeterministic series that originates in the high dimensionality of the atmosphere.     

Automatic estimation of the correlation dimension for the analysis of electrocardiograms

The main purpose of the present work is the definition of a fully automatic procedure for correlation dimension (D₂) estimation. In the first part, the procedure for the estimation of the correlation dimension (D₂) is proposed and tested on various types of mathematical models: chaotic (Lorenz and Henon models), periodical (sinusoidal waves) and stochastic (Gaussian and uniform noise). In all cases, accurate D₂ estimates were obtained. The procedure can detect the presence of multiple scaling regions in the correlation integral function. The connection between the presence of multiple scaling regions and multiple dynamic activities cooperating in a system is investigated through the study of composite time series. In the second part of the paper, the proposed algorithm is applied to the study of cardiac electrical activity through the analysis of electrocardiographic signals (ECG) obtained from the commercially available MIT-BIH ECG arrhythmia database. Three groups of ECG signals have been considered: the ECGs of normal subjects and ECGs of subjects with atrial fibrillation and with premature ventricular contraction. D₂ estimates are computed on single ECG intervals (static analysis) of appropriate duration, striking a balance between stationarity requisites and accurate computation requirements. In addition, D₂ temporal variability is studied by analyzing consecutive intervals of ECG tracings (dynamic analysis). The procedure reveals the presence of multiple scaling regions in many ECG signals, and the D₂ temporal variability differs in the three ECG groups considered; it is greater in the case of atrial fibrillation than in normal sinus rhythms. This study points out the importance of considering both the static and dynamic D₂ analysis for a more complete study of the system under analysis. While the static analysis visualizes the underlying heart activity, dynamic D₂ analysis insights the time evolution of the underlying system. Received: 11 April 1997 / Accepted in revised form: 19 March 1999     

Non-linear and linear forecasting of the EEG time series

The method of non-linear forecasting of time series was applied to different simulated signals and EEG in order to check its ability of distinguishing chaotic from noisy time series. The goodness of prediction was estimated, in terms of the correlation coefficient between forecasted and real time series, for non-linear and autoregressive (AR) methods. For the EEG signal both methods gave similar results. It seems that the EEG signal, in spite of its chaotic character, is well described by the AR model.     

A nonlinear dynamics model for simulating long range correlations of cognitive bistability

Simulation results of bistable perception due to ambiguous visual stimuli are presented which are obtained with a behavioral nonlinear dynamics model using perception–attention–memory coupling. This model provides an explanation of recent experimental results of Gao et al. (Cogn Process 7:105–112, 2006a) and it supports their speculation that the fractal character of perceptual dominance time series may be understood in terms of nonlinear and reentrant dynamics of brain processing. Percept reversals are induced by attention fatigue and noise, with an attention bias which balances the relative percept duration. Dynamical coupling of the attention bias to the perception state introduces memory effects leading to significant long range correlations of perceptual duration times as quantified by the Hurst parameter H > 0.5 (Mandelbrot, The fractal geometry of nature, 1991), in agreement with Gao et al. (Cogn Process 7:105–112, 2006a).     

Correlation analysis of Normalized Different Vegetation Index (NDVI) difference series and climate variables in the Xilingole steppe, China from 1983 to 1999

There is a crucial need in the study of global change to understand how terrestrial ecosystems respond to the climate system. It has been demonstrated by many researches that Normalized Different Vegetation Index (NDVI) time series from remotely sensed data, which provide effective information of vegetation conditions on a large scale with highly temporal resolution, have a good relation with meteorological factors. However, few of these studies have taken the cumulative property of NDVI time series into account. In this study, NDVI difference series were proposed to replace the original NDVI time series with NDVI difference series to reappraise the relationship between NDVI and meteorological factors. As a proxy of the vegetation growing process, NDVI difference represents net primary productivity of vegetation at a certain time interval under an environment controlled by certain climatic conditions and other factors. This data replacement is helpful to eliminate the cumulative effect that exist in original NDVI time series, and thus is more appropriate to understand how climate system affects vegetation growth in a short time scale. By using the correlation analysis method, we studied the relationship between NOAA/AVHRR ten-day NDVI difference series and corresponding meteorological data from 1983 to 1999 from 11 meteorological stations located in the Xilingole steppe in Inner Mongolia. The results show that: (1) meteorological factors are found to be more significantly correlation with NDVI difference at the biomass-rising phase than that at the falling phase; (2) the relationship between NDVI difference and climate variables varies with vegetation types and vegetation communities. In a typical steppe dominated by Leymus chinensis, temperature has higher correlation with NDVI difference than precipitation does, and in a typical steppe dominated by Stipa krylovii, the correlation between temperature and NDVI difference is lower than that between precipitation and NDVI difference. In a typical steppe dominated by Stipa grandis, there is no significant difference between the two correlations. Precipitation is the key factor influencing vegetation growth in a desert steppe, and temperature has poor correlation with NDVI difference; (3) the response of NDVI difference to precipitation is fast and almost simultaneous both in a typical steppe and desert steppe, however, mean temperature exhibits a time-lag effect especially in the desert steppe and some typical steppe dominated by Stipa krylovii; (4) the relationship between NDVI difference and temperature is becoming stronger with global warming. __________ Translated from Acta Phytoecologica Sinica, 2005, 29(5): 753–765 [译自: 植物生态学报]     

Low-dimensional chaotic attractors in the rat brain

The existence of chaotic attractors for discrete time series, derived from the occurrences of spikes during electrophysiological recordings, was investigated. The time series included between 800 and 5200 points per analyzed record. The spike trains were recorded in the substantia nigra pars reticulata (n=13) and in the auditory thalamus (n=14). The experiments were performed on anesthetized rats during spontaneous activity and during auditory stimulation. According to standard methods of dynamical systems theory, an embedding space was constructed using delay coordinates. The embedding and correlation dimensions were computed by means of the correlation integrals. For 7 of 27 samples, a deterministic structure with a low embedding dimension (ranging between 2 and 6) and a correlation dimension between 0.14 and 3.3 could be determined. Evidence was found that the sensory stimulation may affect the chaotic behavior. Single units recorded simultaneously from the same electrode tip may display different chaotic dynamics, even with a similar time-locked response to the stimulus onset.     

Is breathing in infants chaotic? Dimension estimates for respiratory patterns during quiet sleep

Wedescribe an analysis of dynamic behavior apparent in times-seriesrecordings of infant breathing during sleep. Three principal techniqueswere used: estimation of correlation dimension, surrogate dataanalysis, and reduced linear (autoregressive) modeling (RARM). Correlation dimension can be used to quantify the complexity of timeseries and has been applied to a variety of physiological andbiological measurements. However, the methods most commonly used toestimate correlation dimension suffer from some technical problems thatcan produce misleading results if not correctly applied. We used a newtechnique of estimating correlation dimension that has fewer problems.We tested the significance of dimension estimates by comparingestimates with artificial data sets (surrogate data). On the basis ofthe analysis, we conclude that the dynamics of infant breathing duringquiet sleep can best be described as a nonlinear dynamic system withlarge-scale, low-dimensional and small-scale, high-dimensionalbehavior; more specifically, a noise-driven nonlinear system with atwo-dimensional periodic orbit. Using our RARM technique, we identifiedthe second period as cyclic amplitude modulation of the same period asperiodic breathing. We conclude that our data are consistent withrespiration being chaotic.

     

Surrogate data analysis of sleep electroencephalograms reveals evidence for nonlinearity

 We tested the hypothesis of whether sleep electroencephalographic (EEG) signals of different time windows (164 s, 82 s, 41 s and 20.5 s) are in accordance with linear stochastic models. For this purpose we analyzed the all-night sleep electroencephalogram of a healthy subject and corresponding Gaussian-rescaled phase randomized surrogates with a battery of five nonlinear measures. The following nonlinear measures were implemented: largest Lyapunov exponent L1, correlation dimension D2, and the Green-Savit measures δ2, δ4 and δ6. The hypothesis of linear stochastic data was rejected with high statistical significance. L1 and D2 yielded the most pronounced effects, while the Green-Savit measures were only partially successful in differentiating EEG epochs from the phase randomized surrogates. For L1 and D2 the efficiency of distinguishing EEG signals from linear stochastic data decreased with shortening of the time window. Altogether, our results indicate that EEG signals exhibit nonlinear elements and cannot completely be described by linear stochastic models. Received: 21 December 1995/Accepted in revised form: 19 March 1996     

The presence of both stochastic and deterministic dynamics in electroretinograms measured under identical experimental conditions

We have studied whether living things investigated under the same measuring conditions can generate signals with different types of dynamics. We also wanted to detect the possible effects of low-intensity microwaves using parameters of deterministic chaos. For this purpose, two sets of electroretinograms were analysed by methods aimed at recognizing different types of dynamics. Both sets included the time series recorded from objects exposed to low-intensity microwaves and those that were not exposed. The analytical methods are based on nonlinear forecasting and a “surrogate data” technique. Although the experimental conditions were identical for the two sets, we have shown that both have time series with deterministic and stochastic dynamics. We also found that the use of parameters of deterministic dynamics is insufficient to distinguish between the sets. Received: 18 January 1994 / Accepted: 13 June 1997