Area asymmetry of heart rate variability signal

Prediction of mortality in intensive care unit with short-term heart rate variability: Machine learning-based analysis of the MIMIC-III database

BACKGROUND

Prognosis prediction in the intensive care unit (ICU) traditionally relied on physiological scoring systems based on clinical indicators at admission. Electrocardiogram (ECG) provides easily accessible information, with heart rate variability (HRV) derived from ECG showing prognostic value. However, few studies have conducted a comprehensive analysis of HRV-based prognostic model against established standards, which limits the application of HRV's prognostic value in clinical settings. This study aims to evaluate the utility of HRV in predicting mortality in the ICU. Additionally, we analyzed the applicability and interpretability of the HRV-integrated clinical model and identified the HRV factors that are most significant for patient prognosis.

METHODS

A total of 2838 patients from the MIMIC-III database were retrospectively included in this study. These patients were randomly divided into training and testing sets at a 4:1 ratio. We collected 86 HRV indicators from patients' lead II ECG readings between 0.5h and 2h before the time of death in the ICU of deceased patients or time of discharge from the ICU of alive patients, in addition to 9 clinical parameters upon admission. Subsequently, machine learning models were developed by algorithms including logistic regression (LR), Random Forest (RF), Adaptive Boosting (Adaboost), Gradient Boost (GB), eXtreme Gradient Boosting (XGB), and Light GBM (LGB) algorithms. An ensemble model that integrated these six algorithms, along with a deep neural network model, was also explored. The ten most important variables were identified using the Shapley method. Subsequently, an HRV-modified clinical scoring system was constructed through recursive feature elimination.

RESULTS

The study demonstrated that the integrated model, utilizing both clinical and HRV features, outperformed the model based solely on clinical information in XGB, LGB and LR algorithms (p = 0.005-0.03). The ensemble model exhibited the best performance (AUROC = 0.878), followed closely by XGB algorithm (AUROC = 0.869). Both of these models significantly outperformed the APS III scoring system (AUROC = 0.765). Notably, this improvement is not dependent on a specific disease but rather on the timing of ECG recordings that are closer to clinical endpoints. For parameter analysis, Shapley's method identified MSEn, SD1SD2, DFAα1, and DFAα2 as key HRV features in predicting mortality. These variables also showed significant differences in univariate analysis across patients with different clinical outcomes (p < 0.0001). Additionally, regardless of machine learning, the additive scoring system incorporating HRV showed a significant enhancement in prognostic ability compared to traditional physiological scores APS III (p = 0.02).

CONCLUSIONS

The integration of HRV features into mortality prediction models has been shown to enhance predictive performances in ICU. This enhancement is not limited to specific machine learning models or diseases but is influenced by the timing of HRV measurement relative to clinical endpoints. HRV features, when combined with other clinical parameters, offer high interpretability and significant prognostic value. Furthermore, incorporating HRV into traditional ICU scoring systems can lead to improved predictive performance.

Validity of the Pneumonitor for Analysis of Short-Term Heart Rate Asymmetry Extended with Respiratory Data in Pediatric Cardiac Patients

Background: Wearable technologies have been developed to measure physiological parameters conveniently. To consider the new measurement device valid, the crucial point is to assess its reliability with the gold standard. The study aimed to assess the validity of the Pneumonitor (PM, fs = 250 Hz) for acquisition of 5 min RR intervals (RRi) for analysis of heart rate asymmetry (HRA) in relation to the electrocardiography (ECG, fs = 1000 Hz) in a group of 19 pediatric cardiac patients. Association between HRA and respiratory rate (RespRate) was verified. Methods: The validation comprised Bland-Altman analysis, intraclass correlation coefficient, and Student's t-test. Results: Sufficient agreement between 10 from 16 HRA parameters was observed. Different HRA parameters values calculated based on RRi from both devices were related to different results of correlation analysis between two parameters and RespRate. Conclusions: The PM might be considered valid for recording RRi, which are then processed to calculate selected HRA parameters in a group of pediatric cardiac patients in rest condition. However, RRi recorded using devices with fs < 250 Hz may be not adequate for reliable HRA analysis.

Nonlinear Heart Rate Variability Analysis for Sleep Stage Classification Using Integration of Ballistocardiogram and Apple Watch

Purpose

Wearable or non-contact, non-intrusive devices present a practical alternative to traditional polysomnography (PSG) for daily assessment of sleep quality. Physiological signals have been known to be nonlinear and nonstationary as the body adapts to states of rest or activity. By integrating more sophisticated nonlinear methodologies, the accuracy of sleep stage identification using such devices can be improved. This advancement enables individuals to monitor and adjust their sleep patterns more effectively without visiting sleep clinics.

Patients and Methods

Six participants slept for three cycles of at least three hours each, wearing PSG as a reference, along with an Apple Watch, an actigraphy device, and a ballistocardiography (BCG) bed sensor. The physiological signals were processed with nonlinear methods and trained with a long short-term memory (LSTM) model to classify sleep stages. Nonlinear methods, such as return maps with advanced techniques to analyze the shape and asymmetry in physiological signals, were used to relate these signals to the autonomic nervous system (ANS). The changing dynamics of cardiac signals in restful or active states, regulated by the ANS, were associated with sleep stages and quality, which were measurable.

Results

Approximately 73% agreement was obtained by comparing the combination of the BCG and Apple Watch signals against a PSG reference system to classify rapid eye movement (REM) and non-REM sleep stages.

Conclusion

Utilizing nonlinear methods to evaluate cardiac dynamics showed an improved sleep quality detection with the non-intrusive devices in this study. A system of non-intrusive devices can provide a comprehensive outlook on health by regularly measuring sleeping patterns and quality over time, offering a relatively accessible method for participants. Additionally, a non-intrusive system can be integrated into a user's or clinic's bedroom environment to measure and evaluate sleep quality without negatively impacting sleep. Devices placed around the bedroom could measure user vitals over longer periods with minimal interaction from the user, representing their natural sleeping trends for more accurate health and sleep disorder diagnosis.

Real-time machine learning model to predict in-hospital cardiac arrest using heart rate variability in ICU

Predicting in-hospital cardiac arrest in patients admitted to an intensive care unit (ICU) allows prompt interventions to improve patient outcomes. We developed and validated a machine learning-based real-time model for in-hospital cardiac arrest predictions using electrocardiogram (ECG)-based heart rate variability (HRV) measures. The HRV measures, including time/frequency domains and nonlinear measures, were calculated from 5 min epochs of ECG signals from ICU patients. A light gradient boosting machine (LGBM) algorithm was used to develop the proposed model for predicting in-hospital cardiac arrest within 0.5-24 h. The LGBM model using 33 HRV measures achieved an area under the receiver operating characteristic curve of 0.881 (95% CI: 0.875-0.887) and an area under the precision-recall curve of 0.104 (95% CI: 0.093-0.116). The most important feature was the baseline width of the triangular interpolation of the RR interval histogram. As our model uses only ECG data, it can be easily applied in clinical practice.

Complexity and Entropy in Physiological Signals (CEPS): Resonance Breathing Rate Assessed Using Measures of Fractal Dimension, Heart Rate Asymmetry and Permutation Entropy

BACKGROUND

As technology becomes more sophisticated, more accessible methods of interpretating Big Data become essential. We have continued to develop Complexity and Entropy in Physiological Signals (CEPS) as an open access MATLAB® GUI (graphical user interface) providing multiple methods for the modification and analysis of physiological data.

METHODS

To demonstrate the functionality of the software, data were collected from 44 healthy adults for a study investigating the effects on vagal tone of breathing paced at five different rates, as well as self-paced and un-paced. Five-minute 15-s recordings were used. Results were also compared with those from shorter segments of the data. Electrocardiogram (ECG), electrodermal activity (EDA) and Respiration (RSP) data were recorded. Particular attention was paid to COVID risk mitigation, and to parameter tuning for the CEPS measures. For comparison, data were processed using Kubios HRV, RR-APET and DynamicalSystems.jl software. We also compared findings for ECG RR interval (RRi) data resampled at 4 Hz (4R) or 10 Hz (10R), and non-resampled (noR). In total, we used around 190-220 measures from CEPS at various scales, depending on the analysis undertaken, with our investigation focused on three families of measures: 22 fractal dimension (FD) measures, 40 heart rate asymmetries or measures derived from Poincaré plots (HRA), and 8 measures based on permutation entropy (PE).

RESULTS

FDs for the RRi data differentiated strongly between breathing rates, whether data were resampled or not, increasing between 5 and 7 breaths per minute (BrPM). Largest effect sizes for RRi (4R and noR) differentiation between breathing rates were found for the PE-based measures. Measures that both differentiated well between breathing rates and were consistent across different RRi data lengths (1-5 min) included five PE-based (noR) and three FDs (4R). Of the top 12 measures with short-data values consistently within ± 5% of their values for the 5-min data, five were FDs, one was PE-based, and none were HRAs. Effect sizes were usually greater for CEPS measures than for those implemented in DynamicalSystems.jl.

CONCLUSION

The updated CEPS software enables visualisation and analysis of multichannel physiological data using a variety of established and recently introduced complexity entropy measures. Although equal resampling is theoretically important for FD estimation, it appears that FD measures may also be usefully applied to non-resampled data.

Comprehensive breathing variability indices enhance the prediction of extubation failure in patients on mechanical ventilation

BACKGROUND AND OBJECTIVE

Despite the numerous studies on extubation readiness assessment for patients who are invasively ventilated in the intensive care unit, a 10-15% extubation failure rate persists. Although breathing variability has been proposed as a potential predictor of extubation failure, it is mainly assessed using simple statistical metrics applied to basic respiratory parameters. Therefore, the complex pattern of breathing variability conveyed by continuous ventilation waveforms may be underexplored.

METHODS

Here, we aimed to develop novel breathing variability indices to predict extubation failure among invasively ventilated patients. First, breath-to-breath basic and comprehensive respiratory parameters were computed from continuous ventilation waveforms 1 h before extubation. Subsequently, the basic and advanced variability methods were applied to the respiratory parameter sequences to derive comprehensive breathing variability indices, and their role in predicting extubation failure was assessed. Finally, after reducing the feature dimensionality using the forward search method, the combined effect of the indices was evaluated by inputting them into the machine learning models, including logistic regression, random forest, support vector machine, and eXtreme Gradient Boosting (XGBoost).

RESULTS

The coefficient of variation of the dynamic mechanical power per breath (CV-MPd[J/breath]) exhibited the highest area under the receiver operating characteristic curve (AUC) of 0.777 among the individual indices. Furthermore, the XGBoost model obtained the best AUC (0.902) by combining multiple selected variability indices.

CONCLUSIONS

These results suggest that the proposed novel breathing variability indices can improve extubation failure prediction in invasively ventilated patients.

An assessment of heart rate and blood pressure asymmetry in the diagnosis of vasovagal syncope in females

Introduction: Heart Rate Asymmetry (HRA) describes a phenomenon of differences between accelerations and decelerations in human heart rate. Methods used for HRA assessment can be further implemented in the evaluation of asymmetry in blood pressure variations (Blood Pressure Asymmetry-BPA). Methods: We have analyzed retrospectively the series of heartbeat intervals extracted from ECG and beat-to-beat blood pressure signals from 16 vasovagal patients (age: 32.1 ± 13.3; BMI: 21.6 ± 3.8; all female) and 19 healthy subjects (age: 34.6 ± 7.6; BMI: 22.1 ± 3.4; all female) who have undergone tilt test (70°). Asymmetry was evaluated with Poincaré plot-based methods for 5 min recordings from supine and tilt stages of the test. The analyzed biosignals were heart rate (RR), diastolic (dBP) and systolic Blood Pressure (sBP) and Pulse Pressure (PP). In the paper we explored the differences between healthy and vasovagal women. Results: The changes of HRA indicators between supine and tilt were observed only in the control group (Porta Index p = 0.026 and Guzik Index p = 0.005). No significant differences in beat-to-beat variability (i.e. spread of points across the line of identity in Poincaré plot-SD1) of dBP was noted between supine and tilt in the vasovagal group (p = 0.433 in comparison to p = 0.014 in healthy females). Moreover, in vasovagal patients the PP was significantly different (supine: 41.47; tilt: 39.27 mmHg) comparing to healthy subjects (supine: 35.87; tilt: 33.50 mmHg) in supine (p = 0.019) and in tilt (p = 0.014). Discussion: Analysis of HRA and BPA represents a promising method for the evaluation of cardiovascular response to orthostatic stressors, however currently it is difficult to determine a subject's underlying health condition based only on these parameters.

Comprehensive HRV estimation pipeline in Python using Neurokit2: Application to sleep physiology

NeuroKit2 is a Python Toolbox for Neurophysiological Signal Processing. The presented method is an adaptation of NeuroKit2 to simplify and automate computation of the various mathematical estimates of heart rate variability (HRV) or similar time series. By default, the present approach accepts as input electrocardiogram's R-R intervals (RRIs) or peak times, i.e., timestamp of each consecutive R peak, but the RRIs or peak times can also stem from other biosensors such as photoplethysmography (PPGs) or represent more general kinds of biological or non-biological time series oscillations. The data may be derived from a single or several sources such as conventional univariate heart rate time series or intermittently weakly coupled fetal and maternal heart rate data. The method describes preprocessing and computation of an output of 124 HRV measures including measures with a dynamic, time-series-specific optimal time delay-based complexity estimation with a user-definable time window length. I also provide an additional layer of HRV estimation looking at the temporal fluctuations of the HRV estimates themselves, an approach not yet widely used in the field, yet showing promise (doi: 10.3389/fphys.2017.01112). To demonstrate the application of the methodology, I present an approach to studying the dynamic relationships between sleep state architecture and multi-dimensional HRV metrics in 31 subjects. NeuroKit2's documentation is extensive. Here, I attempted to simplify things summarizing all you need to produce the most extensive HRV estimation output available to date as open source and all in one place. The presented Jupyter notebooks allow the user to run HRV analyses quickly and at scale on univariate or multivariate time-series data. I gratefully acknowledge the excellent support from the NeuroKit team.•Univariate or multivariate time series input; ingestion, preprocessing, and computation of 124 HRV metrics.•Estimation of intra- and inter-individual higher order temporal fluctuations of HRV metrics.•Application to a sleep dataset recorded using Apple Watch and expert sleep labeling.

Analysis of Short-Term Heart Rate Asymmetry in High-Performance Athletes and Non-Athletes

Heart rate asymmetry (HRA) refers to how asymmetrically the acceleration and deceleration patterns in heartbeat fluctuations are distributed. There is limited evidence regarding HRA changes in athletes and their association with autonomic regulation. This study aimed to compare the short-term HRA of high-performance athletes and non-athletes during an autonomic function test by calculating relevant HRA measures. This exploratory study obtained beat-to-beat RR interval time series from 15 high-performance athletes and 12 non-athletes during a standardized autonomic function test. This test includes rest, postural change, controlled respiration, prolonged orthostatism, exercise, and recovery phases. The following HRA parameters were computed from the RR time series for both groups: asymmetric spread index (ASI), slope index (SI), Porta’s index (PI), Guzik’s index (GI), and Ehlers’ index (EI). We found significant differences (p < 0.01) in the mean value of several HRA parameters between athletes and non-athletes and across the autonomic function test phases, mainly in postural change and recovery phases. Our results indicate that high-performance athletes manifest a higher number and magnitude of cardiac decelerations than non-athletes after an orthostatic challenge, as indicated by GI and EI. In addition, lower HRA was found in athletes in the recovery phase than in non-athletes, as indicated by ASI.

Longitudinal Changes and Recovery in Heart Rate Variability of Young Healthy Subjects When Exposure to a Hypobaric Hypoxic Environment

Background: The autonomic nervous system (ANS) is crucial for acclimatization. Investigating the responses of acute exposure to a hypoxic environment may provide some knowledge of the cardiopulmonary system’s adjustment mechanism.

Objective: The present study investigates the longitudinal changes and recovery in heart rate variability (HRV) in a young healthy population when exposed to a simulated plateau environment.

Methods: The study followed a strict experimental paradigm in which physiological signals were collected from 33 healthy college students (26 ± 2 years, 171 cm ± 7 cm, 64 ± 11 kg) using a medical-grade wearable device. The subjects were asked to sit in normoxic (approximately 101 kPa) and hypoxic (4,000 m above sea level, about 62 kPa) environments. The whole experimental process was divided into four stable resting measurement segments in chronological order to analyze the longitudinal changes of physical stress and recovery phases. Seventy-six time-domain, frequency-domain, and non-linear indicators characterizing rhythm variability were analyzed in the four groups.

Results: Compared to normobaric normoxia, participants in hypobaric hypoxia had significantly lower HRV time-domain metrics, such as RMSSD, MeanNN, and MedianNN (p &lt; 0.01), substantially higher frequency domain metrics such as LF/HF ratio (p &lt; 0.05), significantly lower Poincaré plot parameters such as SD1/SD2 ratio and other Poincaré plot parameters are reduced considerably (p &lt; 0.01), and Refined Composite Multi-Scale Entropy (RCMSE) curves are reduced significantly (p &lt; 0.01).

Conclusion: The present study shows that elevated heart rates, sympathetic activation, and reduced overall complexity were observed in healthy subjects exposed to a hypobaric and hypoxic environment. Moreover, the results indicated that Multiscale Entropy (MSE) analysis of RR interval series could characterize the degree of minor physiological changes. This novel index of HRV can better explain changes in the human ANS.

The effect of passive lower limb training on heart rate asymmetry

OBJECTIVE

Heart rate asymmetry (HRA) is an approach for quantitatively assessing the uneven distribution of heart rate accelerations and decelerations for sinus rhythm. We aimed to investigate whether automatic regulation led to HRA alternation during passive lower limb training.

METHODS

Thirty healthy participants were recruited in this study. The protocol included a baseline (Pre-E) and three passive lower limb training trials (E1, E2 and E3) with a randomized order. Several variance-based HRA variables were established. Heart rate variability (HRV) parameters, i.e., mean RR, SDNN, RMSSD, LF (n.u.), HF (n.u.) and VLF (ms2), and HRA variables, i.e., SD1a, SD1d, SD2a, SD2d, SDNNa and SDNNd, were calculated by using 5-min RR time series, as well as the normalized HRA variables, i.e., C1a, C1d, C2a, C2d, Ca and Cd.

RESULTS

Our results showed that the performance of HRA was distinguished. The normalized HRA was observed with significant changes in E1, E2 and E3 compared to Pre -E. Moreover, parts of non-normalized HRA variables correlated with HRV parameters, which indicated that HRA might benefit in assessing cardiovascular modulation in passive lower limb training.

CONCLUSIONS

In summary, this study suggested that passive training led to significant HRA alternation and the application of HRA gave us the possibility for autonomic assessment.

EntropyHub: An open-source toolkit for entropic time series analysis

An increasing number of studies across many research fields from biomedical engineering to finance are employing measures of entropy to quantify the regularity, variability or randomness of time series and image data. Entropy, as it relates to information theory and dynamical systems theory, can be estimated in many ways, with newly developed methods being continuously introduced in the scientific literature. Despite the growing interest in entropic time series and image analysis, there is a shortage of validated, open-source software tools that enable researchers to apply these methods. To date, packages for performing entropy analysis are often run using graphical user interfaces, lack the necessary supporting documentation, or do not include functions for more advanced entropy methods, such as cross-entropy, multiscale cross-entropy or bidimensional entropy. In light of this, this paper introduces EntropyHub, an open-source toolkit for performing entropic time series analysis in MATLAB, Python and Julia. EntropyHub (version 0.1) provides an extensive range of more than forty functions for estimating cross-, multiscale, multiscale cross-, and bidimensional entropy, each including a number of keyword arguments that allows the user to specify multiple parameters in the entropy calculation. Instructions for installation, descriptions of function syntax, and examples of use are fully detailed in the supporting documentation, available on the EntropyHub website- www.EntropyHub.xyz. Compatible with Windows, Mac and Linux operating systems, EntropyHub is hosted on GitHub, as well as the native package repository for MATLAB, Python and Julia, respectively. The goal of EntropyHub is to integrate the many established entropy methods into one complete resource, providing tools that make advanced entropic time series analysis straightforward and reproducible.

Preoperative Heart Rate Variability During Sleep Predicts Vagus Nerve Stimulation Outcome Better in Patients With Drug-Resistant Epilepsy

Objective: Vagus nerve stimulation (VNS) is an adjunctive and well-established treatment for patients with drug-resistant epilepsy (DRE). However, it is still difficult to identify patients who may benefit from VNS surgery. Our study aims to propose a VNS outcome prediction model based on machine learning with multidimensional preoperative heart rate variability (HRV) indices. Methods: The preoperative electrocardiography (ECG) of 59 patients with DRE and of 50 healthy controls were analyzed. Responders were defined as having at least 50% average monthly seizure frequency reduction at 1-year follow-up. Time domain, frequency domain, and non-linear indices of HRV were compared between 30 responders and 29 non-responders in awake and sleep states, respectively. For feature selection, univariate filter and recursive feature elimination (RFE) algorithms were performed to assess the importance of different HRV indices to VNS outcome prediction and improve the classification performance. Random forest (RF) was used to train the classifier, and leave-one-out (LOO) cross-validation was performed to evaluate the prediction model. Results: Among 52 HRV indices, 49 showed significant differences between DRE patients and healthy controls. In sleep state, 35 HRV indices of responders were significantly higher than those of non-responders, while 16 of them showed the same differences in awake state. Low-frequency power (LF) ranked first in the importance ranking results by univariate filter and RFE methods, respectively. With HRV indices in sleep state, our model achieved 74.6% accuracy, 80% precision, 70.6% recall, and 75% F1 for VNS outcome prediction, which was better than the optimal performance in awake state (65.3% accuracy, 66.4% precision, 70.5% recall, and 68.4% F1). Significance: With the ECG during sleep state and machine learning techniques, the statistical model based on preoperative HRV could achieve a better performance of VNS outcome prediction and, therefore, help patients who are not suitable for VNS to avoid the high cost of surgery and possible risks of long-term stimulation.

Novel Measures of Similarity and Asymmetry in Upper Limb Activities for Identifying Hemiparetic Severity in Stroke Survivors

Stroke survivors are often characterized by hemiparesis, i.e., paralysis in one half of the body, severely affecting upper limb movements. Monitoring the progression of hemiparesis requires manual observation of limb movements at regular intervals, and hence is a labour intensive process. In this work, we use wrist-worn accelerometers for automated assessment of hemiparesis in acute stroke. We propose novel measures of similarity and asymmetry in hand activities through bivariate Poincaré analysis between two-hand accelerometer data for quantifying hemiparetic severity. The proposed descriptors characterize the distribution of activity surrogates derived from acceleration of the two hands, on a 2D bivariate Poincaré Plot. Experiments show that while the descriptors CSD1 and CSD2 can identify hemiparetic patients from control subjects, their normalized difference CSDR and the descriptors Complex Cross-Correlation Measure ( C3M) and Activity Asymmetry Index ( AAI) can distinguish between mild, moderate and severe hemiparesis. These measures are compared with traditional measures of cross-correlation and evaluated against the National Institutes of Health Stroke Scale (NIHSS), the clinical gold standard for hemiparetic severity estimation. This study, undertaken on 40 acute stroke patients with varying levels of hemiparesis and 15 healthy controls, validates the use of short length ( 5 minutes) wearable accelerometry data for identifying hemiparesis with greater clinical sensitivity. Results show that the proposed descriptors with a hierarchical classification model outperform state-of-the-art methods with overall accuracy of 0.78 and 0.85 for 4-class and 3-class hemiparesis identification respectively.

Short-Term Effect of Percutaneous Coronary Intervention on Heart Rate Variability in Patients with Coronary Artery Disease

Myocardial ischemia in patients with coronary artery disease (CAD) leads to imbalanced autonomic control that increases the risk of morbidity and mortality. To systematically examine how autonomic function responds to percutaneous coronary intervention (PCI) treatment, we analyzed data of 27 CAD patients who had admitted for PCI in this pilot study. For each patient, five-minute resting electrocardiogram (ECG) signals were collected before and after the PCI procedure. The time intervals between ECG collection and PCI were both within 24 h. To assess autonomic function, normal sinus RR intervals were extracted and were analyzed quantitatively using traditional linear time- and frequency-domain measures [i.e., standard deviation of the normal-normal intervals (SDNN), the root mean square of successive differences (RMSSD), powers of low frequency (LF) and high frequency (HF) components, LF/HF] and nonlinear entropy measures [i.e., sample entropy (SampEn), distribution entropy (DistEn), and conditional entropy (CE)], as well as graphical metrics derived from Poincaré plot [i.e., Porta’s index (PI), Guzik’s index (GI), slope index (SI) and area index (AI)]. Results showed that after PCI, AI and PI decreased significantly (p < 0.002 and 0.015, respectively) with effect sizes of 0.88 and 0.70 as measured by Cohen’s d static. These changes were independent of sex. The results suggest that graphical AI and PI metrics derived from Poincaré plot of short-term ECG may be potential for sensing the beneficial effect of PCI on cardiovascular autonomic control. Further studies with bigger sample sizes are warranted to verify these observations.

Heart Rate Asymmetry Analysis During Head-Up Tilt Test in Healthy Men

The purpose of this study is to assess the cardiovascular system response to orthostatic stress in a group of 133 healthy men using heart rate asymmetry (HRA) methods. HRA is a feature of variability in human heart rate which is dependent upon external and internal body conditions. The initial phases of head-up tilt test (HUTT), namely, supine and tilt, were chosen as the external body affecting factors. Various calculation methods of HRA, such as Porta's index (PI), Guzik's index (GI), and its variance based components, were used to assess the heart rate variability (HRV) and its asymmetry. We compared 5-min ECG recordings from both supine and tilt phases of HUT test. Short-term HRA was observed in 54.1% of men in supine phase and 65.4% of men in tilt phase. The study revealed significant increase of GI (from 0.50 to 0.52, p < 0.001) in the tilt phase as well as significant changes in HRV descriptors between HUTT phases. Our results showed that the variability of human heart rate and its asymmetry are sensitive to orthostatic stress. The study of short-term HRA is a potential additional tool to increase sensitivity in conditions where HUTT is a diagnostic tool, such as vasovagal syncope.

Asymmetry of lagged Poincare plot in heart rate signals during meditation

Background and aim

Heart rate variability (HRV) quantifies the variability in the heart’s beat-to-beat intervals. This signal is a potential marker of cardiac function in normal, pathological, and psychological states. Signal asymmetry refers to an unequal distribution in the signal, which can be found by a two-dimensional Poincare plot. Earlier, heart rate asymmetry (HRA) was assessed using a conventional Poincare plot (lag of 1). In this study, we have investigated the effect of delay on the phase space asymmetry using lagged Poincare’s plot.

Experimental procedure

This study compared the presence/lack of asymmetries in the HRV data of 12 meditators (four Kundalini yoga (Yoga) at an advanced level of meditation, eight Chinese Chi meditators (Chi) ∼1–3 months) to 25 non-meditators (11 spontaneous nocturnal breathing (Normal) and 14 metronomic breathing (Metron)). Poincare’s plots were constructed with six different lags, and HRA was calculated. The analysis was conducted using HRV data provided in the Physionet database.

Results

The results showed that using conventional Poincare’s plot (lag of 1), the lowest HRA was observed in the Metron group. In addition, the HRA index was different between meditators and non-meditator groups. Moreover, as the most significant difference between groups was observed in a delay of 6, the role of the delay selection on the signal asymmetry was revealed.

Conclusion

The difference between lagged HRA responses on Yoga in comparison with other groups can be an emphasis on the importance of choosing the type of meditation technique and its effects on the cardiovascular system.

•

Asymmetries in HRV was assessed in different meditator and non-meditator groups.

•

The role of delay selection was explored on the phase space asymmetry using lagged Poincare plot.

•

A weaker asymmetry was observed in the metronomic breathing group.

•

The most significant difference between groups was perceived in a delay of six.

Phase entropy: a new complexity measure for heart rate variability

OBJECTIVE

Information entropy is generally employed for analysing the complexity of physiological signals. However, most definitions of entropy estimate the degree of compressibility and thus quantify the randomness. Physiological signals are very complex because of nonlinear relationships and interactions between various systems and subsystems of the body. Therefore, analysis of randomness may not be sufficient to describe this complexity. To analyse the complexity of heart rate variability (HRV), a new entropy method, phase entropy (PhEn), has been proposed as a quantification of two-dimensional phase space.

APPROACH

The second-order difference plot (SODP), a two-dimensional phase space, provides a visual summary of the rate of variability. The distribution of scatter points in a SODP provides information about the dynamics of the underlying system. PhEn estimates the Shannon entropy of the weighted distribution in a coarse-grained SODP.

MAIN RESULTS

The performance of PhEn has been evaluated using simulated signals, synthetic HRV signals and real HRV signals. PhEn shows a better discriminating power and stability than other entropy measures. It is computationally efficient. Moreover, it has the ability to assess temporal asymmetry of physiological signals.

SIGNIFICANCE

PhEn quantifies the multiplicity and rate of variability associated with physiological signals. It is sensitive to time irreversibility. Therefore, it appears to be a promising tool for analysing physiological signals such as HRV.

Does the Temporal Asymmetry of Short-Term Heart Rate Variability Change during Regular Walking? A Pilot Study of Healthy Young Subjects

The acceleration and deceleration patterns in heartbeat fluctuations distribute asymmetrically, which is known as heart rate asymmetry (HRA). It is hypothesized that HRA reflects the balancing regulation of the sympathetic and parasympathetic nervous systems. This study was designed to examine whether altered autonomic balance during exercise can lead to HRA changes. Sixteen healthy college students were enrolled, and each student undertook two 5-min ECG measurements: one in a resting seated position and another while walking on a treadmill at a regular speed of 5 km/h. The two measurements were conducted in a randomized order, and a 30-min rest was required between them. RR interval time series were extracted from the 5-min ECG data, and HRA (short-term) was estimated using four established metrics, that is, Porta's index (PI), Guzik's index (GI), slope index (SI), and area index (AI), from both raw RR interval time series and the time series after wavelet detrending that removes the low-frequency component of <~0.03 Hz. Our pilot data showed a reduced PI but unchanged GI, SI, and AI during walking compared to resting seated position based on the raw data. Based on the wavelet-detrended data, reduced PI, SI, and AI were observed while GI still showed no significant changes. The reduced PI during walking based on both raw and detrended data which suggests less short-term HRA may underline the belief that vagal tone is withdrawn during low-intensity exercise. GI may not be sensitive to short-term HRA. The reduced SI and AI based on detrended data suggest that they may capture both short- and long-term HRA features and that the expected change in short-term HRA is amplified after removing the trend that is supposed to link to long-term component. Further studies with more subjects and longer measurements are warranted to validate our observations and to examine these additional hypotheses.