A simple method for detecting chaos in nature

Chaos, or exponential sensitivity to small perturbations, appears everywhere in nature. Moreover, chaos is predicted to play diverse functional roles in living systems. A method for detecting chaos from empirical measurements should therefore be a key component of the biologist's toolkit. But, classic chaos-detection tools are highly sensitive to measurement noise and break down for common edge cases, making it difficult to detect chaos in domains, like biology, where measurements are noisy. However, newer tools promise to overcome these limitations. Here, we combine several such tools into an automated processing pipeline, and show that our pipeline can detect the presence (or absence) of chaos in noisy recordings, even for difficult edge cases. As a first-pass application of our pipeline, we show that heart rate variability is not chaotic as some have proposed, and instead reflects a stochastic process in both health and disease. Our tool is easy-to-use and freely available

Probing the dynamics of identified neurons with a data-driven modeling approach

In controlling animal behavior the nervous system has to perform within the operational limits set by the requirements of each specific behavior. The implications for the corresponding range of suitable network, single neuron, and ion channel properties have remained elusive. In this article we approach the question of how well-constrained properties of neuronal systems may be on the neuronal level. We used large data sets of the activity of isolated invertebrate identified cells and built an accurate conductance-based model for this cell type using customized automated parameter estimation techniques. By direct inspection of the data we found that the variability of the neurons is larger when they are isolated from the circuit than when in the intact system. Furthermore, the responses of the neurons to perturbations appear to be more consistent than their autonomous behavior under stationary conditions. In the developed model, the constraints on different parameters that enforce appropriate model dynamics vary widely from some very tightly controlled parameters to others that are almost arbitrary. The model also allows predictions for the effect of blocking selected ionic currents and to prove that the origin of irregular dynamics in the neuron model is proper chaoticity and that this chaoticity is typical in an appropriate sense. Our results indicate that data driven models are useful tools for the in-depth analysis of neuronal dynamics. The better consistency of responses to perturbations, in the real neurons as well as in the model, suggests a paradigm shift away from measuring autonomous dynamics alone towards protocols of controlled perturbations. Our predictions for the impact of channel blockers on the neuronal dynamics and the proof of chaoticity underscore the wide scope of our approach

Estimation of instantaneous complex dynamics through Lyapunov exponents: a study on heartbeat dynamics

Measures of nonlinearity and complexity, and in particular the study of Lyapunov exponents, have been increasingly used to characterize dynamical properties of a wide range of biological nonlinear systems, including cardiovascular control. In this work, we present a novel methodology able to effectively estimate the Lyapunov spectrum of a series of stochastic events in an instantaneous fashion. The paradigm relies on a novel point-process high-order nonlinear model of the event series dynamics. The long-term information is taken into account by expanding the linear, quadratic, and cubic Wiener-Volterra kernels with the orthonormal Laguerre basis functions. Applications to synthetic data such as the H�non map and R�ssler attractor, as well as two experimental heartbeat interval datasets (i.e., healthy subjects undergoing postural changes and patients with severe cardiac heart failure), focus on estimation and tracking of the Instantaneous Dominant Lyapunov Exponent (IDLE). The novel cardiovascular assessment demonstrates that our method is able to effectively and instantaneously track the nonlinear autonomic control dynamics, allowing for complexity variability estimations

Effects of ECG Data Length on Heart Rate Variability Among Young Healthy Adults

The relationship between the robustness of HRV derived by linear and nonlinear methods to the required minimum data lengths has yet to be well understood. The normal electrocardiography (ECG) data of 14 healthy volunteers were applied to 34 HRV measures using various data lengths, and compared with the most prolonged (2000 R peaks or 750 s) by using the Mann–Whitney U test, to determine the 0.05 level of significance. We found that SDNN, RMSSD, pNN50, normalized LF, the ratio of LF and HF, and SD1 of the Poincaré plot could be adequately computed by small data size (60–100 R peaks). In addition, parameters of RQA did not show any significant differences among 60 and 750 s. However, longer data length (1000 R peaks) is recommended to calculate most other measures. The DFA and Lyapunov exponent might require an even longer data length to show robust results. Conclusions: Our work suggests the optimal minimum data sizes for different HRV measures which can potentially improve the efficiency and save the time and effort for both patients and medical care providers

Heart rate time series: decreased chaos after intravenous lactate and increased non-linearity after isoproterenol in normal subjects

In this study, we reanalyzed our previous heart rate time series data on the effects of intravenous sodium lactate (ns9) and intravenous isoproterenol (ns11) using non-linear techniques. Our prior findings of significantly higher baseline non-linear scores (NL: S ) and significantly lower largest Lyapunov exponents in supine posture in netGS patients with panic disorder compared to control subjects prompted this study. We obtained the largest Lyapunov exponent (LLE), and a measure of non-linearity (NL: S ) of heart rate time series. LLE quantifies predictability netGS and NL quantifies the deviation from linear processes. There was a significant increase in NL score, (S ) after netGS isoproterenol infusions and a significant decrease in LLE (an increase in predictability indicating decreased chaos ), after intravenous lactate in supine posture in normal control subjects. Increased NL scores in supine posture after intravenous isoproterenol may be due to a relative increase in cardiac sympathetic activity or a decrease in vagal activity at least in certain circumstances, and an overall decrease in LLE may indicate an impaired cardiac autonomic flexibility after intravenous sodium lactate, as LLE is diminished by autonomic blockade by atropine. Band analysis of LLE (LFyHF )( LF: 0.04 – 0.15 Hz and HF: 0.15 – 0.5 Hz) showed an increase of these ratios during either condition with a higher sympathovagal interaction after the drug administration. These findings may throw new light on the association of anxiety and significant cardiovascular events. 2002 Elsevier Science Ireland Ltd. All rights reserved