Stability, Consistency and Performance of Distribution Entropy in Analysing Short Length Heart Rate Variability (HRV) Signal

Sample, Fuzzy and Distribution Entropies of Heart Rate Variability: What Do They Tell Us on Cardiovascular Complexity?

Distribution Entropy (DistEn) has been introduced as an alternative to Sample Entropy (SampEn) to assess the heart rate variability (HRV) on much shorter series without the arbitrary definition of distance thresholds. However, DistEn, considered a measure of cardiovascular complexity, differs substantially from SampEn or Fuzzy Entropy (FuzzyEn), both measures of HRV randomness. This work aims to compare DistEn, SampEn, and FuzzyEn analyzing postural changes (expected to modify the HRV randomness through a sympatho/vagal shift without affecting the cardiovascular complexity) and low-level spinal cord injuries (SCI, whose impaired integrative regulation may alter the system complexity without affecting the HRV spectrum). We recorded RR intervals in able-bodied (AB) and SCI participants in supine and sitting postures, evaluating DistEn, SampEn, and FuzzyEn over 512 beats. The significance of "case" (AB vs. SCI) and "posture" (supine vs. sitting) was assessed by longitudinal analysis. Multiscale DistEn (mDE), SampEn (mSE), and FuzzyEn (mFE) compared postures and cases at each scale between 2 and 20 beats. Unlike SampEn and FuzzyEn, DistEn is affected by the spinal lesion but not by the postural sympatho/vagal shift. The multiscale approach shows differences between AB and SCI sitting participants at the largest mFE scales and between postures in AB participants at the shortest mSE scales. Thus, our results support the hypothesis that DistEn measures cardiovascular complexity while SampEn/FuzzyEn measure HRV randomness, highlighting that together these methods integrate the information each of them provides.

ComEDA: A new tool for stress assessment based on electrodermal activity

Non-specific sympathetic arousal responses to different stressful elicitations can be easily recognized from the analysis of physiological signals. However, neural patterns of sympathetic arousal during physical and mental fatigue are clearly not unitary. In the context of physiological monitoring through wearable and non-invasive devices, electrodermal activity (EDA) is the most effective and widely used marker of sympathetic activation. This study presents ComEDA, a novel approach for the characterization of complex dynamics of EDA. ComEDA overcomes the methodological limitations related to the application of nonlinear analysis to EDA dynamics, is not parameter-sensitive and is suitable for the analysis of ultra-short time series. We validated the proposed algorithm using synthetic series of white noise and 1/f noise, varying the number of samples from 50 to 5000. By applying our approach, we were able to discriminate a statistically significant increase of complexity in the 1/f noise with respect to white noise, obtaining p-values in the range [4.35 × 10-6, 0.03] after the Mann-Whitney test. Then, we tested ComEDA on both EDA signal and its tonic and phasic components, acquired from healthy subjects during four experimental protocols: two inducing a sympathetic activation through physical efforts and two based on mentally stressful tasks. Results are encouraging and promising, outperforming state of the art metrics such as the Sample Entropy. ComEDA shows good performance not only in discriminating between stressful tasks and resting state (p-value < 0.01 after the Wilcoxon non-parametric statistical test applied to EDA signals of all the four datasets), but also in differentiating different trends of complexity of EDA dynamics when induced by physical and mental stressors. These findings suggest future applications to automatically detect and selectively identify threats due to overwhelming stress impacting both physical and mental health or in the field of telemedicine to monitor autonomic diseases correlated to atypical sympathetic activation. The Matlab code implementing the ComEDA algorithm is available online.

Irregular Autonomic Modulation Predicts Risky Drinking and Altered Ventromedial Prefrontal Cortex Response to Stress in Alcohol Use Disorder

AIMS

Autonomic dysfunction has been associated with risky drinking and alcohol use disorder (AUD). Although autonomic nervous system (ANS) activity has been attributed to the ventromedial prefrontal cortex (VmPFC)-limbic-striatal regions, the specific role of ANS disruption in AUD and its association with these regions remain unclear. Using functional magnetic resonance imaging (fMRI) and concurrent electrocardiogram (ECG), the current study examined neural correlates of ANS activity in AUD and its role in AUD pathology.

METHODS

Demographically matched 20 AUD patients and 20 social drinkers (SD) completed an fMRI task involving repeated exposure to stress, alcohol-cue and neutral-relaxing images in a block design. Based on the known VmPFC-limbic-striatal functions involved in emotions, reward and the ANS, we performed a regions of interest (ROI) analysis to examine the associations between ANS activity and neural responses in the VmPFC, amygdala, and ventral striatum.

RESULTS

Across conditions, AUD patients showed significantly higher levels of overall heart rate (HR) and approximate entropy (ApEn) compared to SD (Ps < 0.05). In all participants, increased HR was associated with greater drinking volume (P < 0.05). In addition, higher ApEn levels were associated with greater drinking volume (P < 0.05) and decreased right VmPFC response to stress (P < 0.05).

DISCUSSION

Our findings demonstrate ANS disruption in AUD indexed by high overall HR and ApEn. The association between ApEn and rVmPFC response suggests that ApEn may play a role in modulating drinking via interactions with neural regions of emotion regulation. These findings provide insight into patterns of ANS disruption and their relevance to AUD pathology.

Heart Rate Variability Analysis for Seizure Detection in Neonatal Intensive Care Units

In Neonatal Intensive Care Units (NICUs), the early detection of neonatal seizures is of utmost importance for a timely clinical intervention. Over the years, several neonatal seizure detection systems were proposed to detect neonatal seizures automatically and speed up seizure diagnosis, most based on the EEG signal analysis. Recently, research has focused on other possible seizure markers, such as electrocardiography (ECG). This work proposes an ECG-based NSD system to investigate the usefulness of heart rate variability (HRV) analysis to detect neonatal seizures in the NICUs. HRV analysis is performed considering time-domain, frequency-domain, entropy and multiscale entropy features. The performance is evaluated on a dataset of ECG signals from 51 full-term babies, 29 seizure-free. The proposed system gives results comparable to those reported in the literature: Area Under the Receiver Operating Characteristic Curve = 62%, Sensitivity = 47%, Specificity = 67%. Moreover, the system’s performance is evaluated in a real clinical environment, inevitably affected by several artefacts. To the best of our knowledge, our study proposes for the first time a multi-feature ECG-based NSD system that also offers a comparative analysis between babies suffering from seizures and seizure-free ones.

Assessment of Physiological Responses During Field Science Task Performance: Feasibility and Future Needs

Objective

By understanding the physiological demands of different types of tasks that will be performed during extravehicular activity (EVA) on Mars, human performance safety risks can be mitigated. In addition, such understanding can assist in planning EVAs with an appropriate balance of human health and safety with scientific mission return.

Background

This paper describes the results of a study of technical feasibility performed within a Mars human research analog, with participants conducting scientifically relevant planetary science sample analysis and return tasks in two distinct field locations.

Methods

The authors collected heart rate, respiration rate, and heart rate variability (HRV) data, using commercial off-the-shelf hardware and software from study participants as they performed field science tasks within a concept of operations for a Mars science return human expedition mission. These data were remotely monitored, shared in real time, and later analyzed to identify different responses to different tasks in order to determine if there were any predictable or consistent patterns among participants.

Results

It was ultimately determined that, while differences exist between responses to tasks, they are highly subject to multiple sources of individual variability, dynamics of evolving field science tasks, and demands of a demanding physical environment. Further, distributional analyses of participants do not support parametric statistical analysis techniques.

Conclusion

The authors conclude that the physiology of individual astronauts should be extensively studied and modeled to support individualized automated monitoring tools for each crew member that is sent to Mars. Application: Physiological monitoring for specialized populations will require significant individual-level analysis, baselining, and bootstrap statistical methods to enable appropriate human performance determinations.

Ischemic Stroke Risk Assessment by Multiscale Entropy Analysis of Heart Rate Variability in Patients with Persistent Atrial Fibrillation

It has been recognized that heart rate variability (HRV), defined as the fluctuation of ventricular response intervals in atrial fibrillation (AFib) patients, is not completely random, and its nonlinear characteristics, such as multiscale entropy (MSE), contain clinically significant information. We investigated the relationship between ischemic stroke risk and HRV with a large number of stroke-naïve AFib patients (628 patients), focusing on those who had never developed an ischemic/hemorrhagic stroke before the heart rate measurement. The CHA2DS2−VASc score was calculated from the baseline clinical characteristics, while the HRV analysis was made from the recording of morning, afternoon, and evening. Subsequently, we performed Kaplan–Meier method and cumulative incidence function with mortality as a competing risk to estimate the survival time function. We found that patients with sample entropy (SE(s)) ≥ 0.68 at 210 s had a significantly higher risk of an ischemic stroke occurrence in the morning recording. Meanwhile, the afternoon recording showed that those with SE(s) ≥ 0.76 at 240 s and SE(s) ≥ 0.78 at 270 s had a significantly lower risk of ischemic stroke occurrence. Therefore, SE(s) at 210 s (morning) and 240 s ≤ s ≤ 270 s (afternoon) demonstrated a statistically significant predictive value for ischemic stroke in stroke-naïve AFib patients.

Resting Heartbeat Complexity Predicts All-Cause and Cardiorespiratory Mortality in Middle- to Older-Aged Adults From the UK Biobank

Background

Spontaneous heart rate fluctuations contain rich information related to health and illness in terms of physiological complexity, an accepted indicator of plasticity and adaptability. However, it is challenging to make inferences on complexity from shorter, more practical epochs of data. Distribution entropy (DistEn) is a recently introduced complexity measure that is designed specifically for shorter duration heartbeat recordings. We hypothesized that reduced DistEn predicted increased mortality in a large population cohort.

Method and Results

The prognostic value of DistEn was examined in 7631 middle‐older–aged UK Biobank participants who had 2‐minute resting ECGs conducted (mean age, 59.5 years; 60.4% women). During a median follow‐up period of 7.8 years, 451 (5.9%) participants died. In Cox proportional hazards models with adjustment for demographics, lifestyle factors, physical activity, cardiovascular risks, and comorbidities, for each 1‐SD decrease in DistEn, the risk increased by 36%, 56%, and 73% for all‐cause, cardiovascular, and respiratory disease–related mortality, respectively. These effect sizes were equivalent to the risk of death from being >5 years older, having been a former smoker, or having diabetes mellitus. Lower DistEn was most predictive of death in those <55 years with a prior myocardial infarction, representing an additional 56% risk for mortality compared with older participants without prior myocardial infarction. These observations remained after controlling for traditional mortality predictors, resting heart rate, and heart rate variability.

Conclusions

Resting heartbeat complexity from short, resting ECGs was independently associated with mortality in middle‐ to older‐aged adults. These risks appear most pronounced in middle‐aged participants with prior MI, and may uniquely contribute to mortality risk screening.

Modified Distribution Entropy as a Complexity Measure of Heart Rate Variability (HRV) Signal

The complexity of a heart rate variability (HRV) signal is considered an important nonlinear feature to detect cardiac abnormalities. This work aims at explaining the physiological meaning of a recently developed complexity measurement method, namely, distribution entropy (DistEn), in the context of HRV signal analysis. We thereby propose modified distribution entropy (mDistEn) to remove the physiological discrepancy involved in the computation of DistEn. The proposed method generates a distance matrix that is devoid of over-exerted multi-lag signal changes. Restricted element selection in the distance matrix makes "mDistEn" a computationally inexpensive and physiologically more relevant complexity measure in comparison to DistEn.

Permutation Entropy: Enhancing Discriminating Power by Using Relative Frequencies Vector of Ordinal Patterns Instead of Their Shannon Entropy

Many measures to quantify the nonlinear dynamics of a time series are based on estimating the probability of certain features from their relative frequencies. Once a normalised histogram of events is computed, a single result is usually derived. This process can be broadly viewed as a nonlinear IRn mapping into IR, where n is the number of bins in the histogram. However, this mapping might entail a loss of information that could be critical for time series classification purposes. In this respect, the present study assessed such impact using permutation entropy (PE) and a diverse set of time series. We first devised a method of generating synthetic sequences of ordinal patterns using hidden Markov models. This way, it was possible to control the histogram distribution and quantify its influence on classification results. Next, real body temperature records are also used to illustrate the same phenomenon. The experiments results confirmed the improved classification accuracy achieved using raw histogram data instead of the PE final values. Thus, this study can provide a very valuable guidance for the improvement of the discriminating capability not only of PE, but of many similar histogram-based measures.

Comparison of QT interval variability of coronary patients without myocardial infarction with that of patients with old myocardial infarction

BACKGROUND

The significant association of myocardial ischemia with elevated QT interval variability (QTV) has been reported in myocardial infarction (MI) patients. However, the influence of the time course of MI on QTV has not been investigated systematically.

METHOD

Short-term QT and RR interval time series were constructed from the 5 min electrocardiograms of 49 coronary patients without MI and 26 patients with old MI (OMI). The QTV, heart rate variability (HRV), and QT-RR coupling of the two groups were analyzed using various time series analysis tools in the time- and frequency-domains, as well as nonlinear dynamics.

RESULTS

Nearly all of the tested QTV indices for coronary patients with OMI were higher than those for patients without MI. However, no significant differences were found between the two groups in any of the variables employed to assess the HRV and QT-RR coupling. All of the markers that showed statistical significances in univariate analyses still possessed the capabilities of distinguishing between the two groups even after adjusting for studied baseline characteristics, including the coronary atherosclerotic burden.

CONCLUSIONS

The results suggested that the QTV increased in coronary patients with OMI compared to those without MI, which might reflect the influence of post-MI remodeling on the beat-to-beat temporal variability of ventricular repolarization. The non-significant differences in the HRV and QT-RR couplings could indicate that there were no differences in the modulation of the autonomic nervous system and interaction of QT with the RR intervals between the two groups.

Renyi Distribution Entropy Analysis of Short-Term Heart Rate Variability Signals and Its Application in Coronary Artery Disease Detection

Coronary artery disease (CAD) is a life-threatening condition that, unless treated at an early stage, can lead to congestive heart failure, ischemic heart disease, and myocardial infarction. Early detection of diagnostic features underlying electrocardiography signals is crucial for the identification and treatment of CAD. In the present work, we proposed novel entropy called Renyi Distribution Entropy (RdisEn) for the analysis of short-term heart rate variability (HRV) signals and the detection of CAD. Our simulation experiment with synthetic, physiological, and pathological signals demonstrated that RdisEn could distinguish effectively among different subject groups. Compared to the values of sample entropy or approximation entropy, the RdisEn value was less affected by the parameter choice, and it remained stable even in short-term HRV. We have developed a combined CAD detection scheme with RdisEn and wavelet packet decomposition (WPD): (1) Normal and CAD HRV beats obtained were divided into two equal parts. (2) Feature acquisition: RdisEn and WPD-based statistical features were calculated from one part of HRV beats, and student’s t-test was performed to select clinically significant features. (3) Classification: selected features were computed from the remaining part of HRV beats and fed into K-nearest neighbor and support vector machine, to separate CAD from normal subjects. The proposed scheme automatically detected CAD with 97.5% accuracy, 100% sensitivity and 95% specificity and performed better than most of the existing schemes.

Embedded Dimension and Time Series Length. Practical Influence on Permutation Entropy and Its Applications

Permutation Entropy (PE) is a time series complexity measure commonly used in a variety of contexts, with medicine being the prime example. In its general form, it requires three input parameters for its calculation: time series length N, embedded dimension m, and embedded delay τ. Inappropriate choices of these parameters may potentially lead to incorrect interpretations. However, there are no specific guidelines for an optimal selection of N, m, or τ, only general recommendations such as N>>m!, τ=1, or m=3,…,7. This paper deals specifically with the study of the practical implications of N>>m!, since long time series are often not available, or non-stationary, and other preliminary results suggest that low N values do not necessarily invalidate PE usefulness. Our study analyses the PE variation as a function of the series length N and embedded dimension m in the context of a diverse experimental set, both synthetic (random, spikes, or logistic model time series) and real–world (climatology, seismic, financial, or biomedical time series), and the classification performance achieved with varying N and m. The results seem to indicate that shorter lengths than those suggested by N>>m! are sufficient for a stable PE calculation, and even very short time series can be robustly classified based on PE measurements before the stability point is reached. This may be due to the fact that there are forbidden patterns in chaotic time series, not all the patterns are equally informative, and differences among classes are already apparent at very short lengths.

EZ Entropy: a software application for the entropy analysis of physiological time-series

Background

Entropy analysis has been attracting increasing attentions in the recent two or three decades. It assesses complexity, or irregularity, of time-series which is extraordinarily relevant to physiology and diseases as demonstrated by tremendous studies. However, the complexity can hardly be appreciated by traditional methods including time-, frequency-domain analysis, and time-frequency analysis that are the common built-in options in commercialized measurement and statistical software. To facilitate the entropy analysis of physiological time-series, a new software application, namely EZ Entropy, was developed and introduced in this article.

Results

EZ Entropy was developed in MATLAB^® environment. It was programmed in an object-oriented style and was constructed with a graphical user interface. EZ Entropy is easy to operate through its compact graphical interface, thus allowing researchers without knowledge of programming like clinicians and physiologists to perform such kind of analysis. Besides, it offers various settings to meet different analysis needs including (1) processing single data recording, (2) batch processing multiple data files, (3) sliding window calculations, (4) recall, (5) displaying intermediate data and final results, (6) adjusting input parameters, and (7) exporting calculation results after the run or in real-time during the analysis. The analysis results could be exported, either manually or automatically, to comma-separated ASCII files, thus being compatible to and easily imported into the common statistical analysis software. Code-wise, EZ Entropy is object-oriented, thus being quite easy to maintain and extend.

Conclusions

EZ Entropy is a user-friendly software application to perform the entropy analysis of time-series, as well as to simplify and to speed up this useful analysis.

Electronic supplementary material

The online version of this article (10.1186/s12938-019-0650-5) contains supplementary material, which is available to authorized users.

Cardiovascular assessment of supportive doctor-patient communication using multi-scale and multi-lag analysis of heartbeat dynamics

Emphatic doctor-patient communication has been associated with improved psycho-physiological well-being involving cardiovascular and neuroendocrine responses. Nevertheless, a comprehensive assessment of heartbeat linear and nonlinear dynamics throughout the communication of a life-threatening disease has not been performed yet. To this extent, we studied linear heartbeat dynamics through the extraction of time-frequency domain measurements, as well as heartbeat nonlinear and complex dynamics through novel approaches to compute multi-scale and multi-lag series analyses: namely, the multi-scale distribution entropy and lagged Poincaré plot symbolic analysis. Heart rate variability series were recorded from 54 healthy female subjects who were blind to the aim of the experiment. Participants were randomly assigned into two groups: 27 subjects watched a video where an oncologist discloses the diagnosis of a cancer metastasis to a patient, whereas the remaining 27 watched the same video including four additional supportive comments by the clinician. Considering differences between the beginning and the end of each communication video, results from non-parametric Wilcoxon tests demonstrated that, at a group level, significant differences occurred in heartbeat linear and nonlinear dynamics, with lower complexity during nonsupportive communication. Furthermore, a support vector machine algorithm, validated using a leave-one-subject-out procedure, was able to discern the supportive experience at a single-subject level with an accuracy of 83.33% when nonlinear features were considered, dropping to 51.85% when using standard HRV features only. In conclusion, heartbeat nonlinear and complex dynamics can be a viable tool for the psycho-physiological evaluation of supportive doctor-patient communication. Graphical Abstract Scheme of the three main stages of the study: signal acquisition during doctor-patient communication, ECG signal processing and pattern recognition results.

Estimation of Complexity of Sampled Biomedical Continuous Time Signals Using Approximate Entropy

Non-linear analysis found many applications in biomedicine. Approximate entropy (ApEn) is a popular index of complexity often applied to biomedical data, as it provides quite stable indications when processing short and noisy epochs. However, ApEn strongly depends on parameters, which were chosen in the literature in wide ranges. This paper points out that ApEn depends on sampling rate of continuous time signals, embedding dimension, tolerance (under which a match is identified), epoch duration and low frequency trends. Moreover, contradicting results can be obtained changing parameters. This was found both in simulations and in experimental EEG. These limitations of ApEn suggest the introduction of an alternative index, here called modified ApEn, which is based on the following principles: oversampling is compensated, self-recurrences are ignored, a fixed percentage of recurrences is selected and low frequency trends are removed. The modified index allows to get more stable measurements of the complexity of the tested simulated data and EEG. The final conclusions are that, in order to get a reliable estimation of complexity using ApEn, parameters should be chosen within specific ranges, data must be sampled close to the Nyquist limit and low frequency trends should be removed. Following these indications, different studies could be more easily compared, interpreted and replicated. Moreover, the modified ApEn can be an interesting alternative, which extends the range of parameters for which stable indications can be achieved.