Next Article in Journal
Key Concepts of Systemological Approach to CPS Adaptive Information Security Monitoring
Previous Article in Journal
Calculating Beamforming Vectors for 5G System Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation on the Prediction of Cardiovascular Events Based on Multi-Scale Time Irreversibility Analysis

1
School of Science, China Pharmaceutical University, Nanjing 211198, China
2
School of Physics & Information Technology, Shaanxi Normal University, Xi’an 710119, China
*
Authors to whom correspondence should be addressed.
These authors contribute equally to this article.
Submission received: 12 October 2021 / Revised: 24 November 2021 / Accepted: 24 November 2021 / Published: 15 December 2021

Abstract

:
Investigation of the risk factors associated with cardiovascular disease (CVD) plays an important part in the prevention and treatment of CVD. This study investigated whether alteration in the multi-scale time irreversibility of sleeping heart rate variability (HRV) was a risk factor for cardiovascular events. The D-value, based on analysis of multi-scale increments in HRV series, was used as the measurement of time irreversibility. Eighty-four subjects from an open-access database (i.e., the Sleep Heart Health Study) were included in this study. None of them had any CVD history at baseline; 42 subjects had cardiovascular events within 1 year after baseline polysomnography and were classed as the CVD group, and the other 42 subjects in the non-CVD group were age matched with those in the CVD group and had no cardiovascular events during the 15-year follow-up period. We compared D-values of sleeping HRV between the CVD and non-CVD groups and found that the D-values of the CVD group were significantly lower than those of the non-CVD group on all 10 scales, even after adjusting for gender and body mass index. Moreover, we investigated the performance of a machine learning model to classify CVD and non-CVD subjects. The model, which was fed with a feature space based on the D-values on 10 scales and trained by a random forest algorithm, achieved an accuracy of 80.8% and a positive prediction rate of 86.7%. These results suggest that the decreased time irreversibility of sleeping HRV is an independent predictor of cardiovascular events that could be used to assist the intelligent prediction of cardiovascular events.

1. Introduction

Cardiovascular disease (CVD) has become the most common chronic disease [1]. The early and accurate prediction of CVD risk is thus of vital importance for the timely prevention and treatment of CVD events. The assessment of CVD risk and prediction of CVD events have been topics of great interest in recent decades. Multiple cardiovascular risk assessment systems, including the Framingham risk score, have been proposed and applied in clinical diagnosis. However, it remains a challenge to improve the accuracy of methods to identify the risk of CVD events automatically. Thus, exploring effective biomarkers of CVD events is of great significance.
As a reflection of the cardiovascular dynamical system, electrocardiograms (ECGs) contain substantial information associated with the physiological and pathological activity of the heart [2]. Heart rate variability (HRV), which is defined as the variation in the continuous interval of heart beats, can be extracted from ECG signals. It is widely accepted that HRV is controlled by a variety of factors, including autonomic modulation, body fluids, and blood pressure. As a non-invasive method, HRV analysis has been used in numerous studies to investigate the alterations in cardiovascular dynamics related to CVD [3,4].
It has been established and is now generally accepted that the cardiovascular system is chaotic and nonlinear [5], as it is coordinated and controlled by many factors; these include the autonomic nervous system (ANS), which is composed of the sympathetic and parasympathetic nervous systems. Time irreversibility is an important nonlinear aspect and fundamental property of the cardiovascular system [6]. If a dynamic system is time irreversible, it can return to its past state when time is reversed. In statistics, a time series is generally regarded as time irreversible only when its statistical characteristics will not change after time is reversed [7].
In recent decades, time irreversibility analysis has been widely applied in many research fields including dynamics, finance, and physiological signals [8,9,10,11,12]. Cammarota et al. proposed that the HRV signal in healthy people was time irreversible [13], whereas Costa et al. reported that the time irreversibility of HRV signals decreased with aging and with the occurrence of CVD [14]. However, despite studies showing that HRV signals in healthy controls were time irreversible and became more symmetrical with the emergence of diseases, especially CVD, it has remained unclear whether the time irreversibility of HRV has predictive value for cardiovascular events.
Almost all previous studies on time irreversibility analysis focused on 24 h HRV signals [13,14] or signals that had been sampled during the daytime when the participants were awake [15,16]. However, about a third of a person’s lifetime is spent in sleep. Numerous studies have shown that sleep has important roles in human attention [17,18], working memory [19], visual learning [20], mood regulation [21], cognitive function [22,23], etc. Furthermore, the sleep–wake cycle also affects the regulation of the cardiovascular system by the ANS [24,25]. As demonstrated by Morris et al. and Neufeld et al., circadian misalignment increases CVD risk factors in humans [26], and sleep deprivation may increase the risk of CVD [27]. In addition, alterations in HRV during sleep were found to be related to the risk of CVD. Nakayama et al. revealed that increased activity during the first hour of sleep in patients with cardiovascular risk factors increased HRV [28]. Eguchi et al. showed that HRV during sleep was closely associated with an increased risk of CVD in people with type 2 diabetes [29]. Ulmer et al. demonstrated that the high-frequency component of HRV in sleep, which reflects the regulatory effect of parasympathetic nerves on the cardiovascular system, is an independent predictor of CVD [30]. Zhang et al. demonstrated that sleep heart rate variability assists the automatic prediction of long-term cardiovascular outcomes [31]. However, most of these studies only measured HRV based on time and frequency domain analyses, which may not be sufficient to determine the characteristics of the complex dynamics of the chaotic and nonlinear cardiovascular system [32].
Therefore, based on a nonlinear analysis method (i.e., time irreversibility analysis), this study evaluated the differences in HRV time irreversibility during sleep in participants who were healthy at baseline but suffered from cardiovascular events in the subsequent year (i.e., the CVD group) and in those who experienced no cardiovascular events during the 15-year follow-up period (i.e., the non-CVD group).We aimed to explore whether the change in HRV time irreversibility during sleep had predictive value for the occurrence of cardiovascular events and could thus contribute to the improvement and development of cardiovascular risk assessment systems.

2. Materials and Methods

2.1. Participants

The research data for this study were downloaded from a public database (the Sleep Heart Health Study) [33]. All participants underwent baseline polysomnography (PSG) for sleep monitoring between 1995 and 1998. After the baseline PSG monitoring, all the participants were followed up for 15 years, and information about cardiovascular events including coronary heart disease, angina, myocardial infarction, heart failure, and stroke was recorded.
Among the multiple signals recorded in PSG monitoring, there is one ECG channel with a sampling rate of 125 Hz. The R-wave peak of the ECG was identified using the Pan–Tompkins method [34], and the original RR intervals were then calculated. Moreover, to obtain the final HRV series for further analysis, artifacts in the original RR intervals were detected according to the following criteria: (1) less than 300 ms or 0.8 times the median RR intervals, or (2) larger than 1700 ms or 1.2 times the median RR intervals.
In this study, participants who experienced emergence of cardiovascular events during the following year after baseline PSG monitoring were selected as the CVD group. Each subject was required to have no baseline CVD history and a HRV series longer than 10,000 points after sleep onset, with a proportion of artifacts less than 10%. In this way, 42 participants were selected for the CVD group. They comprised 24 males and 18 females, aged 70 ± 9 years (mean ± standard deviation, SD), with a body mass index (BMI) of 29.28 ± 4.69 (mean ± SD) kg/m2. The occurrence time of the first cardiovascular event among these participants was 221 ± 114 days (mean ± SD) after baseline PSG recording. Moreover, 42 aged-matched healthy controls with neither CVD history at baseline nor experience of any cardiovascular events during the 15-year follow-up period were included as the non-CVD group. There were 17 males and 25 females in the non-CVD group, with BMI values of 29.01 ± 4.50 kg/m2 (mean ± SD). In line with the CVD group, for each participant in the non-CVD group, the HRV series during sleep was required to be longer than 10,000 points, with a proportion of artifacts less than 10%.

2.2. Multi-Scale Time Irreversibility Analysis

For each participant, an HRV segment of 10,000 data points (after artifact removal) during sleep was used for time irreversibility analysis. Given an HRV time series { x i | i = 1 , 2 ,   , N } with N data points, a coarse-grained series { y j τ | j = 1 , 2 ,   , N τ } was firstly applied according to Formula (1) [35].
y j τ = 1 τ i = j 1 τ + 1 j τ x i   , 1 j N τ
Here, the integer τ represents the scale factor for coarse graining, and N τ means to round down the ratio of N and τ . When τ equals 1, the coarse-grained time series is actually the original time series x i .
For each coarse-grained series y j τ , differences in every two successive values were then calculated and denoted as { R R k τ |   1 k N τ 1 } [36], as shown in Formula (2).
R R k τ = y k + 1 τ y k τ ,     1 k N τ 1
In { R R k τ } , the components greater than 0 were denoted as { R R + } , while the components less than 0 were denoted as { R R } . According to Hou et al. [36], two measurements regarding the time irreversibility of y j τ could then be calculated according to Formulas (3) and (4).
P τ = N R R N R R + N R R + × 100 %
G τ = i = 1 N R R + R R + i 2 i = 1 N R R + R R + i 2 + i = 1 N R R R R i 2 × 100 %
In Formulas (3) and (4),   N RR +   and N RR represent the numbers of components in RR + and RR , respectively. The time irreversibility of the HRV series could thus be quantified by calculating the Euclidean distance of ( P τ , G τ ) to the symmetric center (50%, 50%), as shown in Formula (5) [36]:
D τ = P τ 50 % 2 + G τ 50 % 2
In this way, we could obtain   D τ when different scale factors were considered, such as τ ranging from 1 to 10 in steps of 1.
In this way, D τ broke through the limitation that the measurement of time irreversibility is limited to a two-dimensional state space and reflects the average condition of the asymmetry of the multi-dimensional vector with respect to the main diagonal on multiple projection planes [37]. A larger value of D τ corresponded to more asymmetric values when reading the HRV series from both forward and backward directions. Thus, the metric D τ was proposed and demonstrated as a measurement of time irreversibility of the HRV series [36].
In this study, for each participant, D τ (τ ranging from 1 to 10 in steps of 1) was calculated on the corresponding HRV segment of data points. The corresponding values of D τ were denoted as D1, D2, …, and D10. Moreover, the average of these D-values was also computed and denoted as Dmean.

2.3. Conventional HRV Analysis

In this study, we also performed conventional HRV analysis on each HRV segment with 10,000 data points to explore whether the proposed multi-scale time irreversibility analysis outperformed the conventional approach.
Here, we used four well-established HRV indices derived from time domain, frequency domain, and nonlinear complexity analyses, i.e., the standard deviation of all RR intervals (SDNN), the power in the high-frequency (0.15–0.4 Hz) range (HF), the power in the low-frequency (0.04–0.15 Hz) range in normalized units (LFnorm) [38], and the multi-scale sample entropy (MSE) [35]. HF and LFnorm were calculated on each 5 min HRV segment [38] and then averaged for the whole point segment. Similar to the proposed multi-scale time irreversibility analysis, the calculation of MSE was conducted on 10 time scales of each point segment. For each scale, SE was computed with an embedding dimension of 2 and a tolerance of 2 * SD (the SD of the point segment) [35]. Then, the average of these SE values was computed and denoted by SEmean.

2.4. Classification of CVD and Non-CVD Participants by Random Forest (RF) Algorithm

To further investigate the capability of the time irreversibility of sleep HRV for CVD prediction, we used machine learning technology to classify CVD and non-CVD participants based on their D-values over 10 scales. The RF algorithm was chosen as it is widely used in classification tasks [39]. RF is a machine learning method based on the decision tree (DT) and bagging ensemble learning algorithm proposed by Leo Breiman and Adele Cutler in 2001 [40].

2.4.1. The Basic Principle of DT

DT is a process of classifying samples in a training set based on a series of decision rules, including classification DT and regression DT. The core idea behind the construction of a classification DT is to recursively divide the training samples into different subsets from the DT’s root according to the value of a certain feature, until all sample categories in the subset at a node are the same or the features of the samples have been exhausted. For example, Figure 1 shows the decision-making process used to screen participants for the CVD group in this study. Four features were considered, i.e., whether there was a baseline CVD history or not, whether the individual suffered from cardiovascular events during the following year or not, whether the HRV series was longer than 10,000 points or not, and whether the proportion of artifacts was less than 10% or not.

2.4.2. Bagging Ensemble Learning Algorithm

Bagging is a parallel ensemble learning method proposed by Breiman et al. in 1996 [41]. Its basic principle is shown in Figure 2. With the help of the bagging algorithm, a stronger learner can be established by the combination of a series of weaker learners. Each weaker learner is trained independently using a subset of all the available samples based on bootstrap sampling. For classification tasks, a simple voting strategy is usually used to obtain the stronger learner in the combination stage [42].

2.4.3. RF Algorithm

RF is a variant of the bagging ensemble learning algorithm that adopts DTs as the weak learners and introduces random feature selection into the training process of DTs. Compared with a single DT, RF usually shows better generalization ability. The pseudo-code of RF is shown in Algorithm 1.
Algorithm 1: Random forest
Input: original sample set D; number of decision tree T
Output: strong learner G
1  For t = 1, 2, …, T
2  request size-N’ data Dt by bootstrapping with D;
3  obtain tree gt with random feature selection;
4  return G = Uniform({gt });
As a bagging ensemble learning algorithm, the RF algorithm generates T training sets by bootstrapping and then constructs a DT for each training set. Moreover, for each DT, a subset of features should also be randomly selected from all the features. That is, each DT in RF is trained with both independent samples and features.
In this study, the RF model was used to assign each participant to the CVD group or the non-CVD group. For each subject, there were ten features, i.e., the D-values on 10 scales (D1, D2, …, and D10), which were combined to form a feature vector for the participant and then fed into the RF classifiers. To achieve better generalization, five-fold cross-validation was used in the model selection and training processes. Furthermore, two hyperparameters of the RF algorithm, i.e., the number of DTs (denoted as n_estimators) and the maximum depth of the tree (denoted as max_depth), were adjusted by grid search in Python 3.6 with the help of the sklearn library [40].

2.5. Statistical Analysis

To investigate the difference in the time irreversibility of sleep HRV between the CVD and non-CVD groups, we compared the D-values obtained on different scales for both groups. If the D-values followed a normal distribution, a t-test was used; otherwise, the Wilcoxon rank test was used. Moreover, logistic regression analysis was used to explore whether the D-value was an independent predictor of cardiovascular events after adjusting for the effects of gender and BMI, as studies have shown that women have a lower risk of CVD than men [43], and that obesity is a risk factor for CVD [44]. A value of p < 0.05 was considered to indicate a significant difference.
For comparison of the proposed multi-scale time irreversibility analysis with conventional HRV analysis, the variance inflation factor (VIF) was first used to check the multicollinearity of the HRV indices, and the area under the receiver operating characteristic (ROC) curve (AUC) was used as the measurement of the ability to discriminate between the CVD and non-CVD groups.
All the statistical analyses were performed using the SAS9.4 software (SAS Institute Inc., Cary, NC, USA).

3. Results

For each subject, time irreversibility analysis and conventional analysis were applied to its HRV time series with data points, which was obtained after the participant fell asleep in the night and preprocessed with artifact removal. D1, D2, …, and D10 were calculated and compared between the CVD and non-CVD groups. As shown in Table 1 and Figure 3, compared with the non-CVD group, all the D-values significantly decreased in the CVD group (p < 0.05), indicating a weakening of time irreversibility with a higher risk of CVD.
Moreover, logistic regression models were used to test whether the differences in D-values between the two groups were associated with gender or BMI. As shown in Table 2, a significant reduction in the D-value remained an independent predictor of cardiovascular events after correction for gender and BMI (p < 0.05).
Furthermore, we used the RF algorithm to classify CVD and non-CVD subjects automatically using a feature space constructed based on the values of D1, D2, …, and D10. The results of the grid search and five-fold cross-validation show that the highest performance (i.e., accuracy of 80.8%, specificity of 81.3%, sensitivity of 80.0%, negative prediction rate of 72.7%, and positive prediction rate of 86.7%) was achieved when max_depth was set to 2 and n_estimators was set to 20.
Finally, we compared the ability of the various HRV indices, i.e., Dmean, SDNN, HF, LFnorm, and SEmean, to discriminate between CVD and non-CVD participants (Figure 4). The multicollinearity test showed no collinearity among these indices (VIF less than 10), and the highest AUC (0.7948) was achieved by Dmean. These results suggest that the proposed method could supplement conventional HRV analysis and might be more sensitive to the pathological dynamics of CVD than conventional HRV indices. Moreover, we calculated the AUC for D1 to separate CVD and non-CVD participants. The result (0.744) was also lower than that achieved with Dmean.

4. Conclusions and Discussion

In this study, we investigated the role of the time irreversibility of sleep HRV in the prediction of CVD events. Forty-two participants who were healthy at baseline but suffered from cardiovascular events in the subsequent year were included as potential CVD patients. Compared with those participants who were healthy at baseline and had no cardiovascular events in the following 15 years, the results show the following: (1) the time irreversibility of the cardiovascular system during sleep was significantly reduced in people at risk of CVD, regardless of the scale used; (2) the significant decrease in time irreversibility remained after adjusting for gender and BMI; (3) the machine learning model based on the RF algorithm showed an accuracy of 80.8% in predicting whether the participant would have CVD events within 1 year by using features derived from the proposed time irreversibility analysis, i.e., D1, D2, …, and D10; (4) Dmean, i.e., the averaged time irreversibility on all 10 scales, outperformed several conventional HRV indices as well as D1 in distinguishing potential CVD patients. These findings suggest that the proposed method might be a more sensitive means of capturing alterations in cardiovascular dynamics in CVD patients than conventional HRV analysis. Moreover, the measurements of time irreversibility used here (i.e., D1, D2, …, and D10) were potential predictors of cardiovascular events.
The human body is a complex physiological system. In order to maintain the stability of the body’s internal environment, the ANS regulates the heart through various cardiovascular reflexes [45]. These reflexes often act on the heart’s dynamics system with different time delays, leading to complex fluctuations in the heart rate on multiple time scales. This is the main reason why multi-scale analysis has been widely used in the study of physiological signals [35]. In this work, we calculated the D-values of 1–10 scales to explore the predictive value of time irreversibility analysis for cardiovascular events. Our results show that the AUC of Dmean was higher than that of D1, indicating that multi-scale analysis was more informative than only focusing on the original time scale in cardiac dynamics research.
Since the 1990s, nonlinear concepts and research methods have been widely used in the analysis of HRV time series. Time irreversibility analysis, a nonlinear method, has also been used successfully [13,14,46,47]. Previous studies have shown that time irreversibility was reduced among CVD patients, especially in those suffering with heart failure [14]. Our work extends the results of previous studies to demonstrate that the decrease in the time irreversibility of sleep HRV may be a potential biomarker for CVD at a very early stage. Furthermore, our results indicate that time irreversibility analysis might be associated with some aspect of the underlying mechanism of the cardiovascular system that can be captured by nonlinear complexity analysis as well as conventional time and frequency domain analyses.
Owing to the limited sample size, we only used the RF algorithm to construct a machine learning model to distinguish CVD and non-CVD participants. The model’s classification was solely based on the D-values. Although the model was simple, it achieved an accuracy of 80.8% and a positive prediction rate of 86.7%, indicating its promise for applications in the early prediction of cardiovascular events to enable timely intervention. For example, it would be informative for clinical staff to include the proposed biomarker in cardiovascular risk assessment systems. Moreover, people engaged in the construction of CVD prediction models could use this biomarker as one characteristic input of their model. In addition, with the development of wearable devices, e.g., smart watches, HRV data obtained during sleep are becoming increasingly accessible. These wearable devices could use this biomarker to monitor the user’s cardiovascular health.
In conclusion, in this study, we demonstrated the predictive value of the multi-scale time irreversibility of sleep HRV for CVD events. However, there were still some limitations. Considering the influence of sample quantity on the accuracy of models, it would be necessary to confirm the results of this study in a wider range of subjects. Moreover, adverse events such as snoring and sleep apnea could interfere with ECG, but information about such events was not available in the database. Therefore, future efforts could also focus on exploring the effects of these adverse events on the time irreversibility of sleep HRV.

Author Contributions

Conceptualization, X.W. and J.L.; data curation, Q.Y.; formal analysis, Q.Y.; investigation, Q.Y.; methodology, X.W. and Q.Y.; project administration, X.W.; resources, F.H.; software, Q.Y.; supervision, F.H.; validation, X.W. and Q.Y.; visualization, X.W. and Q.Y.; writing—original draft, X.W.; writing—review and editing, F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 11974231 and No. 11774212) and the Natural Science Foundation of Shaanxi Province, China (Grant No. 2020SF-134).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

In this study, we employed data from an open-access database, i.e., the Sleep-EDF database available at Sleep Heart Health Study. The study protocol of the used datasets was approved by the institutional review board of each participating center, and each participant signed informed consent. All methods were carried out in accordance with relevant guidelines and regulations. The current study only analyzed de-identified data from those databases and did not involve a research protocol requiring approval by the relevant institutional review board or ethics committee.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Sleep-EDF Database Expanded: https://sleepdata.org/datasets/shhs (accessed on 22 September 2021).

Acknowledgments

The authors would like to acknowledge a graduated student, Lulu Zhang, for her useful suggestions on the experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Cardiovascular Diseases (CVDs). Available online:https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds) (accessed on 20 June 2021).
  2. Harris, P.R.E. The Normal Electrocardiogram: Resting 12-Lead and Electrocardiogram Monitoring in the Hospital. Crit. Care Nurs. Clin. N. Am. 2016, 28, 281–296. [Google Scholar] [CrossRef] [PubMed]
  3. María, E.R.; Lluís, C.; Ramos, C.J.; Cleofé, Z.M.; Sara, M.; José, A.; Jesús, C.-M. Reduced heart rate variability predicts fatigue severity in individuals with chronic fatigue syndrome/myalgic encephalomyelitis. J. Transl. Med. 2020, 18, 4. [Google Scholar] [CrossRef]
  4. Yan, X.; Zhang, L.; Li, J.; Du, D.; Hou, F. Entropy-Based Measures of Hypnopompic Heart Rate Variability Contribute to the Automatic Prediction of Cardiovascular Events. Entropy 2020, 22, 241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Karavaev, A.S.; Ishbulatov, Y.M.; Ponomarenko, V.I.; Bezruchko, B.P.; Kiselev, A.R.; Prokhorov, M.D. Autonomic control is a source of dynamical chaos in the cardiovascular system. Chaos 2019, 29, 121101. [Google Scholar] [CrossRef]
  6. Humeau-Heurtier, A.; Mahé, G.; Chapeau-Blondeau, F.; Rousseau, D.; Abraham, P. Study of time reversibility/irreversibility of cardiovascular data: Theoretical results and application to laser Doppler flowmetry and heart rate variability signals. Phys. Med. Biol. 2012, 57, 4335–4351. [Google Scholar] [CrossRef] [PubMed]
  7. Kennel, M.B. Testing time symmetry in time series using data compression dictionaries. Phys. Rev. E Stat. Nonlin. Soft. Matter. Phys. 2004, 69, 056208. [Google Scholar] [CrossRef] [Green Version]
  8. Huang, Y.; Yang, D.; Wang, L.; Wang, K. Classifying of welding time series based on multi-scale time irreversibility analysis and extreme learning machine. Chaos Solitons Fractals 2020, 139, 110040. [Google Scholar] [CrossRef]
  9. Lucia, U.; Grisolia, G.; Kuzemsky, A.L. Time, Irreversibility and Entropy Production in Nonequilibrium Systems. Entropy 2020, 22, 887. [Google Scholar] [CrossRef]
  10. Xia, J.; Shang, P.; Wang, J.; Shi, W. Classifying of financial time series based on multiscale entropy and multiscale time irreversibility. Phys. A Stat. Mech. Its Appl. 2014, 400, 151–158. [Google Scholar] [CrossRef]
  11. Yao, W.; Dai, J.; Perc, M.; Wang, J.; Yao, D.; Guo, D. Permutation-based time irreversibility in epileptic electroencephalograms. Nonlinear Dyn. 2020, 100, 907–919. [Google Scholar] [CrossRef]
  12. Zanin, M.; Güntekin, B.; Aktürk, T.; Hanoğlu, L.; Papo, D. Time Irreversibility of Resting-State Activity in the Healthy Brain and Pathology. Front. Physiol. 2019, 10, 1619. [Google Scholar] [CrossRef]
  13. Cammarota, C.; Rogora, E. Time reversal, symbolic series and irreversibility of human heartbeat. Chaos Solitons Fractals 2006, 32, 1649–1654. [Google Scholar] [CrossRef]
  14. Costa, M.; Goldberger, A.L.; Peng, C.K. Broken asymmetry of the human heartbeat: Loss of time irreversibility in aging and disease. Phys. Rev. Lett. 2005, 95, 198102. [Google Scholar] [CrossRef] [PubMed]
  15. Chladekova, L.; Czippelova, B.; Turianikova, Z.; Tonhajzerova, I.; Calkovska, A.; Baumert, M.; Javorka, M. Multiscale time irreversibility of heart rate and blood pressure variability during orthostasis. Physiol. Meas. 2012, 33, 1747–1756. [Google Scholar] [CrossRef]
  16. Tonhajzerova, I.; Ondrejka, I.; Chladekova, L.; Farsky, I.; Visnovcova, Z.; Calkovska, A.; Jurko, A.; Javorka, M. Heart rate time irreversibility is impaired in adolescent major depression. Prog. Neuropsychopharmacol. Biol. Psychiatry 2012, 39, 212–217. [Google Scholar] [CrossRef] [PubMed]
  17. Spiegel, K.; Tasali, E.; Leproult, R.; Van Cauter, E. Effects of poor and short sleep on glucose metabolism and obesity risk. Nat. Rev. Endocrinol. 2009, 5, 253–261. [Google Scholar] [CrossRef] [PubMed]
  18. Faraut, B.; Boudjeltia, K.Z.; Vanhamme, L.; Kerkhofs, M. Immune, inflammatory and cardiovascular consequences of sleep restriction and recovery. Sleep Med. Rev. 2012, 16, 137–149. [Google Scholar] [CrossRef]
  19. Born, J.; Wilhelm, I. System consolidation of memory during sleep. Psychol. Res. 2012, 76, 192–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Ma, Y.; Hou, F.; Yang, A.C.; Ahn, A.C.; Fan, L.; Peng, C.-K. Symbolic dynamics of electroencephalography is associated with the sleep depth and overall sleep quality in healthy adults. Phys. A Stat. Mech. Its Appl. 2018, 513, 22–31. [Google Scholar] [CrossRef]
  21. Ford, D.E.; Cooper-Patrick, L. Sleep disturbances and mood disorders: An epidemiologic perspective. Depress. Anxiety 2001, 14, 3–6. [Google Scholar] [CrossRef]
  22. Song, Y.; Blackwell, T.; Yaffe, K.; Ancoli-Israel, S.; Redline, S.; Stone, K.L. Relationships between sleep stages and changes in cognitive function in older men: The MrOS Sleep Study. Sleep 2015, 38, 411–421. [Google Scholar] [CrossRef] [Green Version]
  23. Blackwell, T.; Yaffe, K.; Ancoli-Israel, S.; Schneider, J.L.; Cauley, J.A.; Hillier, T.A.; Fink, H.A.; Stone, K.L. Poor sleep is associated with impaired cognitive function in older women: The study of osteoporotic fractures. J. Gerontol. A Biol. Sci. Med. Sci. 2006, 61, 405–410. [Google Scholar] [CrossRef] [Green Version]
  24. Aslan, S.; Erbil, N.; Tezer, F.I. Heart Rate Variability During Nocturnal Sleep and Daytime Naps in Patients With Narcolepsy Type 1 and Type 2. J. Clin. Neurophysiol. 2019, 36, 104–111. [Google Scholar] [CrossRef]
  25. Bazilio, D.S.; Bonagamba, L.G.H.; Moraes, D.J.A.; Machado, B.H. Cardiovascular and respiratory profiles during the sleep-wake cycle of rats previously submitted to chronic intermittent hypoxia. Exp. Physiol. 2019, 104, 1408–1419. [Google Scholar] [CrossRef] [PubMed]
  26. Morris, C.J.; Purvis, T.E.; Hu, K.; Scheer, F.A. Circadian misalignment increases cardiovascular disease risk factors in humans. Proc. Natl. Acad. Sci. USA 2016, 113, E1402–E1411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Neufeld, E.V.; Carney, J.J.; Dolezal, B.A.; Boland, D.M.; Cooper, C.B. Exploratory Study of Heart Rate Variability and Sleep among Emergency Medical Services Shift Workers. Prehosp. Emerg. Care 2017, 21, 18–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Nakayama, N.; Miyachi, M.; Tamakoshi, K.; Negi, K.; Watanabe, K.; Hirai, M. Increased Activity in Patients with Cardiovascular Risk Factors Increases Heart Rate Variability. West J. Nurs. Res. 2020, 42, 431–436. [Google Scholar] [CrossRef]
  29. Eguchi, K.; Schwartz, J.E.; Pickering, T.G.; Hoshide, S.; Ishikawa, J.; Shimada, K.; Kario, K. Increased heart rate variability during sleep is a predictor for future cardiovascular events in patients with type 2 diabetes. Hypertens. Res. 2010, 33, 737–742. [Google Scholar] [CrossRef]
  30. Ulmer, C.S.; Hall, M.H.; Dennis, P.A.; Beckham, J.C.; Germain, A. Posttraumatic stress disorder diagnosis is associated with reduced parasympathetic activity during sleep in US veterans and military service members of the Iraq and Afghanistan wars. Sleep 2018, 41, zsy174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Zhang, L.; Wu, H.; Zhang, X.; Wei, X.; Hou, F.; Ma, Y. Sleep heart rate variability assists the automatic prediction of long-term cardiovascular outcomes. Sleep Med. 2020, 67, 217–224. [Google Scholar] [CrossRef]
  32. Fang, S.C.; Wu, Y.L.; Tsai, P.S. Heart Rate Variability and Risk of All-Cause Death and Cardiovascular Events in Patients With Cardiovascular Disease: A Meta-Analysis of Cohort Studies. Biol. Res. Nurs. 2020, 22, 45–56. [Google Scholar] [CrossRef]
  33. Dean, D.A.; Goldberger, A.L.; Mueller, R.; Kim, M.; Rueschman, M.; Mobley, D.; Sahoo, S.S.; Jayapandian, C.P.; Cui, L.; Morrical, M.G.; et al. Scaling Up Scientific Discovery in Sleep Medicine: The National Sleep Research Resource. Sleep 2016, 39, 1151–1164. [Google Scholar] [CrossRef] [Green Version]
  34. Pan, J.; Tompkins, W.J. A real-time QRS detection algorithm. IEEE Trans. Biomed. Eng. 1985, 32, 230–236. [Google Scholar] [CrossRef]
  35. Costa, M.; Goldberger, A.L.; Peng, C.K. Multiscale entropy analysis of complex physiologic time series. Phys. Rev. Lett. 2002, 89, 068102. [Google Scholar] [CrossRef] [Green Version]
  36. Hou, F.; Zhuang, J.; Bian, C.; Tong, T.; Chen, Y.; Yin, J.; Qiu, X.; Ning, X. Analysis of heartbeat asymmetry based on multi-scale time irreversibility test. Phys. A Stat. Mech. Its Appl. 2009, 389, 754–760. [Google Scholar] [CrossRef]
  37. Hou, F.; Ning, X.; Zhuang, J.; Huang, X.; Fu, M.; Bian, C. High-dimensional time irreversibility analysis of human interbeat intervals. Med. Eng. Phys. 2011, 33, 633–637. [Google Scholar] [CrossRef]
  38. Camm, A.J.; Malik, M.; Bigger, J.T.; Breithardt, G.; Cerutti, S.; Cohen, R.J.; Coumel, P.; Fallen, E.L.; Kennedy, H.L.; Kleiger, R.E.; et al. Heart rate variability: Standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Circulation 1996, 93, 1043–1065. [Google Scholar]
  39. Fernández-Delgado, M.; Cernadas, E.; Barro, S.; Amorim, D. Do we need hundreds of classifiers to solve real world classification problems? J. Mach. Learn. Res. 2014, 15, 3133–3181. [Google Scholar] [CrossRef]
  40. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  41. Breiman, L. Bagging predictors. Mach. Learn. 1996, 24, 123–140. [Google Scholar] [CrossRef] [Green Version]
  42. Bayar, B.; Stamm, M.C. Constrained Convolutional Neural Networks: A New Approach Towards General Purpose Image Manipulation Detection. IEEE Trans. Inf. Forensics Secur. 2018, 13, 2691–2706. [Google Scholar] [CrossRef]
  43. Stramba-Badiale, M.; Fox, K.M.; Priori, S.G.; Collins, P.; Daly, C.; Graham, I.; Jonsson, B.; Schenck-Gustafsson, K.; Tendera, M. Cardiovascular diseases in women: A statement from the policy conference of the European Society of Cardiology. Eur. Heart J. 2006, 27, 994–1005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hubert, H.B.; Feinleib, M.; McNamara, P.M.; Castelli, W.P. Obesity as an independent risk factor for cardiovascular disease: A 26-year follow-up of participants in the Framingham Heart Study. Circulation 1983, 67, 968–977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Goldberger, A.L.; Peng, C.K.; Lipsitz, L.A. What is physiologic complexity and how does it change with aging and disease? Neurobiol. Aging 2002, 23, 23–26. [Google Scholar] [CrossRef]
  46. Porta, A.; Guzzetti, S.; Montano, N.; Gnecchi-Ruscone, T.; Furlan, R.; Malliani, A. Time Reversibility in Short-Term Heart Period Variability. In Proceedings of the 2006 Computers in Cardiology, Valencia, Spain, 17–20 September 2006; IEEE: Piscataway, NJ, USA, 2008; pp. 77–80. [Google Scholar]
  47. Porta, A.; D’Addio, G.; Bassani, T.; Maestri, R.; Pinna, G.D. Assessment of cardiovascular regulation through irreversibility analysis of heart period variability: A 24 hours Holter study in healthy and chronic heart failure populations. Philos. Trans. A Math Phys. Eng. Sci. 2009, 367, 1359–1375. [Google Scholar] [CrossRef] [PubMed]
Figure 1. An example of a decision-making process based on a decision tree. The decision “included” or “excluded” was made according to whether the conditions (in the rectangles) were met.
Figure 1. An example of a decision-making process based on a decision tree. The decision “included” or “excluded” was made according to whether the conditions (in the rectangles) were met.
Symmetry 13 02424 g001
Figure 2. Schematic diagram of the bagging ensemble learning algorithm. Bagging extracts m samples from the original sample set but ensures different content each time by bootstrapping and carries out T rounds of extraction to obtain T independent training sets.
Figure 2. Schematic diagram of the bagging ensemble learning algorithm. Bagging extracts m samples from the original sample set but ensures different content each time by bootstrapping and carries out T rounds of extraction to obtain T independent training sets.
Symmetry 13 02424 g002
Figure 3. Comparison of D-values (mean ± standard error) between non-CVD group and CVD group on 1–10 scale. * Significant difference in D-value between non-CVD group and CVD group (paired t-test, p < 0.05) on the corresponding scale.
Figure 3. Comparison of D-values (mean ± standard error) between non-CVD group and CVD group on 1–10 scale. * Significant difference in D-value between non-CVD group and CVD group (paired t-test, p < 0.05) on the corresponding scale.
Symmetry 13 02424 g003
Figure 4. ROC curves of the HRV indices used to discriminate between CVD and non-CVD participants. The AUC of the ROC curves of these HRV indices reflects their sensitivity in capturing the pathological dynamics of CVD.
Figure 4. ROC curves of the HRV indices used to discriminate between CVD and non-CVD participants. The AUC of the ROC curves of these HRV indices reflects their sensitivity in capturing the pathological dynamics of CVD.
Symmetry 13 02424 g004
Table 1. Comparison of D-value between non-CVD group and CVD group on 1–10 scale.
Table 1. Comparison of D-value between non-CVD group and CVD group on 1–10 scale.
VariableNon-CVD GroupCVD Groupp
D 1, M [Q1, Q3]9.52 [7.92, 11.65]6.53 [2.6, 8.83]<0.001
D 2, M [Q1, Q3]7.13 [4.99, 8.59]3.46 [1.92, 6.12]<0.001
D 3, Mean ± SD7.43 ± 3.084.15 ± 2.54<0.001
D 4, Mean ± SD7.6 ± 3.324.12 ± 2.83<0.001
D 5, M [Q1, Q3]7.1 [5.75, 9.07]3.17 [2.13, 5.8]<0.001
D 6, Mean ± SD7.02 ± 3.474.07 ± 2.73<0.001
D 7, Mean ± SD6.61 ± 3.254.21 ± 2.47<0.001
D 8, Mean ± SD6.29 ± 3.093.99 ± 2.53<0.001
D 9, M [Q1, Q3]5.93 [3.65, 8.01]3.01 [1.84, 4.87]<0.001
D 10, M [Q1, Q3]5.84 [3.72, 7.22]3.2 [2.32, 5.78]0.002
Table 2. Predictive value of time irreversibility to cardiovascular events by logistic regression analysis.
Table 2. Predictive value of time irreversibility to cardiovascular events by logistic regression analysis.
VariableβSEOR (95% CI)p
D 1−0.2970.0790.743 (0.636–0.868)<0.001
D 2−0.3160.0900.729 (0.611–0.870)<0.001
D 3−0.4510.1110.637 (0.513–0.791)<0.001
D 4−0.3790.0960.685 (0.567–0.826)<0.001
D 5−0.3650.0930.694 (0.578–0.833)<0.001
D 6−0.3200.0910.726 (0.607–0.868)<0.001
D 7−0.2910.0930.747 (0.623–0.896)0.002
D 8−0.2850.0910.752 (0.629–0.899)0.002
D 9−0.3160.1000.729 (0.600–0.886)0.002
D 10−0.2820.1010.754 (0.619–0.919)0.005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, X.; Yang, Q.; Li, J.; Hou, F. Investigation on the Prediction of Cardiovascular Events Based on Multi-Scale Time Irreversibility Analysis. Symmetry 2021, 13, 2424. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13122424

AMA Style

Wu X, Yang Q, Li J, Hou F. Investigation on the Prediction of Cardiovascular Events Based on Multi-Scale Time Irreversibility Analysis. Symmetry. 2021; 13(12):2424. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13122424

Chicago/Turabian Style

Wu, Xiaochuan, Qianru Yang, Jin Li, and Fengzhen Hou. 2021. "Investigation on the Prediction of Cardiovascular Events Based on Multi-Scale Time Irreversibility Analysis" Symmetry 13, no. 12: 2424. https://0-doi-org.brum.beds.ac.uk/10.3390/sym13122424

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop