MuTE: a new matlab toolbox for estimating the multivariate transfer entropy in physiological variability series

We present a new time series analysis toolbox, developed in Matlab, for the estimation of the Transfer entropy (TE) between time series taken from a multivariate dataset. The main feature of the toolbox is its fully multivariate implementation, that is made possible by the design of an approach for the non-uniform embedding (NUE) of the observed time series. The toolbox is equipped with parametric (linear) and non-parametric (based on binning or nearest neighbors) entropy estimators. All these estimators, implemented using the NUE approach in comparison with the classical approach based on uniform embedding, are tested on RR interval, systolic pressure and respiration variability series measured from healthy subjects during head-up tilt. The results support the necessity of resorting to NUE for obtaining reliable estimates of the multivariate TE in short-term cardiovascular and cardiorespiratory variability

Efficient transfer entropy analysis of non-stationary neural time series

Information theory allows us to investigate information processing in neural systems in terms of information transfer, storage and modification. Especially the measure of information transfer, transfer entropy, has seen a dramatic surge of interest in neuroscience. Estimating transfer entropy from two processes requires the observation of multiple realizations of these processes to estimate associated probability density functions. To obtain these observations, available estimators assume stationarity of processes to allow pooling of observations over time. This assumption however, is a major obstacle to the application of these estimators in neuroscience as observed processes are often non-stationary. As a solution, Gomez-Herrero and colleagues theoretically showed that the stationarity assumption may be avoided by estimating transfer entropy from an ensemble of realizations. Such an ensemble is often readily available in neuroscience experiments in the form of experimental trials. Thus, in this work we combine the ensemble method with a recently proposed transfer entropy estimator to make transfer entropy estimation applicable to non-stationary time series. We present an efficient implementation of the approach that deals with the increased computational demand of the ensemble method's practical application. In particular, we use a massively parallel implementation for a graphics processing unit to handle the computationally most heavy aspects of the ensemble method. We test the performance and robustness of our implementation on data from simulated stochastic processes and demonstrate the method's applicability to magnetoencephalographic data. While we mainly evaluate the proposed method for neuroscientific data, we expect it to be applicable in a variety of fields that are concerned with the analysis of information transfer in complex biological, social, and artificial systems.Comment: 27 pages, 7 figures, submitted to PLOS ON

Transfer Entropy Analysis of Pulse Arrival Time - Heart Period Interactions during Physiological Stress

Although Heart Period (HP) variability is the most widely used measure to assess cardiovascular oscillations, its evaluation combined with that of Pulse Arrival Time (PAT) variability may provide additional information about cardiac dynamics and cardiovascular interactions. In this study, we computed the transfer entropy from PAT to HP in 76 subjects monitored at rest and during orthostatic and mental stress using both a model-free (k- Nearest Neighbors) and a linear parametric estimator. Our results show how the information flow between these two variables depends on the physiological condition and how the nonlinear measure captures more information than the linear one during orthostatic stress

Information-Theoretic Analysis of Cardio-Respiratory Interactions in Heart Failure Patients: Effects of Arrhythmias and Cardiac Resynchronization Therapy

The properties of cardio-respiratory coupling (CRC) are affected by various pathological conditions related to the cardiovascular and/or respiratory systems. In heart failure, one of the most common cardiac pathological conditions, the degree of CRC changes primarily depend on the type of heart-rhythm alterations. In this work, we investigated CRC in heart-failure patients, applying measures from information theory, i.e., Granger Causality (GC), Transfer Entropy (TE) and Cross Entropy (CE), to quantify the directed coupling and causality between cardiac (RR interval) and respiratory (Resp) time series. Patients were divided into three groups depending on their heart rhythm (sinus rhythm and presence of low/high number of ventricular extrasystoles) and were studied also after cardiac resynchronization therapy (CRT), distinguishing responders and non-responders to the therapy. The information-theoretic analysis of bidirectional cardio-respiratory interactions in HF patients revealed the strong effect of nonlinear components in the RR (high number of ventricular extrasystoles) and in the Resp time series (respiratory sinus arrhythmia) as well as in their causal interactions. We showed that GC as a linear model measure is not sensitive to both nonlinear components and only model free measures as TE and CE may quantify them. CRT responders mainly exhibit unchanged asymmetry in the TE values, with statistically significant dominance of the information flow from Resp to RR over the opposite flow from RR to Resp, before and after CRT. In non-responders this asymmetry was statistically significant only after CRT. Our results indicate that the success of CRT is related to corresponding information transfer between the cardiac and respiratory signal quantified at baseline measurements, which could contribute to a better selection of patients for this type of therapy

Instantaneous Transfer Entropy for the Study of Cardiovascular and Cardio-Respiratory Nonstationary Dynamics

Objective: Measures of Transfer Entropy (TE) quantify the direction and strength of coupling between two complex systems. Standard approaches assume stationarity of the observations, and therefore are unable to track time-varying changes in nonlinear information transfer with high temporal resolution. In this study, we aim to define and validate novel instantaneous measures of transfer entropy to provide an im- proved assessment of complex non-stationary cardio-respiratory interactions

Identifying changes in EEG information transfer during drowsy driving by transfer entropy

© 2015 Huang, Pal, Chuang and Lin. Drowsy driving is a major cause of automobile accidents. Previous studies used neuroimaging based approaches such as analysis of electroencephalogram (EEG) activities to understand the brain dynamics of different cortical regions during drowsy driving. However, the coupling between brain regions responding to this vigilance change is still unclear. To have a comprehensive understanding of neural mechanisms underlying drowsy driving, in this study we use transfer entropy, a model-free measure of effective connectivity based on information theory. We investigate the pattern of information transfer between brain regions when the vigilance level, which is derived from the driving performance, changes from alertness to drowsiness. Results show that the couplings between pairs of frontal, central, and parietal areas increased at the intermediate level of vigilance, which suggests that an enhancement of the cortico-cortical interaction is necessary to maintain the task performance and prevent behavioral lapses. Additionally, the occipital-related connectivity magnitudes monotonically decreases as the vigilance level declines, which further supports the cortical gating of sensory stimuli during drowsiness. Neurophysiological evidence of mutual relationships between brain regions measured by transfer entropy might enhance the understanding of cortico-cortical communication during drowsy driving

A mathematical model for breath gas analysis of volatile organic compounds with special emphasis on acetone

Recommended standardized procedures for determining exhaled lower respiratory nitric oxide and nasal nitric oxide have been developed by task forces of the European Respiratory Society and the American Thoracic Society. These recommendations have paved the way for the measurement of nitric oxide to become a diagnostic tool for specific clinical applications. It would be desirable to develop similar guidelines for the sampling of other trace gases in exhaled breath, especially volatile organic compounds (VOCs) which reflect ongoing metabolism. The concentrations of water-soluble, blood-borne substances in exhaled breath are influenced by: (i) breathing patterns affecting gas exchange in the conducting airways; (ii) the concentrations in the tracheo-bronchial lining fluid; (iii) the alveolar and systemic concentrations of the compound. The classical Farhi equation takes only the alveolar concentrations into account. Real-time measurements of acetone in end-tidal breath under an ergometer challenge show characteristics which cannot be explained within the Farhi setting. Here we develop a compartment model that reliably captures these profiles and is capable of relating breath to the systemic concentrations of acetone. By comparison with experimental data it is inferred that the major part of variability in breath acetone concentrations (e.g., in response to moderate exercise or altered breathing patterns) can be attributed to airway gas exchange, with minimal changes of the underlying blood and tissue concentrations. Moreover, it is deduced that measured end-tidal breath concentrations of acetone determined during resting conditions and free breathing will be rather poor indicators for endogenous levels. Particularly, the current formulation includes the classical Farhi and the Scheid series inhomogeneity model as special limiting cases.Comment: 38 page