Next Article in Journal
A Novel Adaptive Feature Fusion Strategy for Image Retrieval
Previous Article in Journal
Update of Prior Probabilities by Minimal Divergence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Similarity Changes Analysis for Heart Rate Fluctuation Regularity as a New Screening Method for Congestive Heart Failure

1
School of Biomedical Engineering, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
2
School of Science, Hua Zhong Agricultural University, Wuhan 430070, China
*
Authors to whom correspondence should be addressed.
Entropy 2021, 23(12), 1669; https://doi.org/10.3390/e23121669
Submission received: 10 November 2021 / Revised: 26 November 2021 / Accepted: 30 November 2021 / Published: 11 December 2021
(This article belongs to the Section Signal and Data Analysis)

Abstract

:
Congestive heart failure (CHF) is a chronic cardiovascular condition associated with dysfunction of the autonomic nervous system (ANS). Heart rate variability (HRV) has been widely used to assess ANS. This paper proposes a new HRV analysis method, which uses information-based similarity (IBS) transformation and fuzzy approximate entropy (fApEn) algorithm to obtain the fApEn_IBS index, which is used to observe the complexity of autonomic fluctuations in CHF within 24 h. We used 98 ECG records (54 health records and 44 CHF records) from the PhysioNet database. The fApEn_IBS index was statistically significant between the control and CHF groups (p < 0.001). Compared with the classical indices low-to-high frequency power ratio (LF/HF) and IBS, the fApEn_IBS index further utilizes the changes in the rhythm of heart rate (HR) fluctuations between RR intervals to fully extract relevant information between adjacent time intervals and significantly improves the performance of CHF screening. The CHF classification accuracy of fApEn_IBS was 84.69%, higher than LF/HF (77.55%) and IBS (83.67%). Moreover, the combination of IBS, fApEn_IBS, and LF/HF reached the highest CHF screening accuracy (98.98%) with the random forest (RF) classifier, indicating that the IBS and LF/HF had good complementarity. Therefore, fApEn_IBS effusively reflects the complexity of autonomic nerves in CHF and is a valuable CHF assessment tool.

1. Introduction

Congestive heart failure (CHF) is a clinical condition with inadequate ventricle filling or manifestation of inadequate myocardial contraction (myocardial failure), caused by changes in cardiac structure and function [1]. It has been identified as a major public health concern and extensively studied within the past two decades. There are about 3 million patients in the United States, nearly 1.5% of the adult population suffering from CHF [2]. The amount of CHF patients is rapidly growing because of the population aging and the increase of survival after myocardial infarction [3]. In addition, CHF may also cause some underlying heart disease [4], or bring structural and functional derangements to the liver with congestion [5]. Thus, the diagnosis of CHF patients is critical.
A previous study proved that cardiac arrhythmias in CHF is associated with the instability of autonomic nervous system (ANS) [6]. To this end, scholars usually extract heart rate variability (HRV) indicators from ECG signals, which can reflect the function of ANS, such as the vagal and sympathetic function, to estimate the activity of ANS in CHF [7,8]. Traditionally, HRV analysis is usually based on frequency and time domain. Takase et al. [9] reported that the HRV index on time-domain could successfully distinguish patients and normal individuals. Rovere et al. [10] found that power in low-frequency LF (0.04–0.15 Hz) can effectively predict the probability of sudden death in patients with CHF. However, these methods are effective, HR fluctuations are considered complex behaviors derived from nonlinear and non-stationary processes [11,12,13]. So, the traditional linear methods may not suit CHF studies because they would be disturbed by fluctuations and cause experimental errors.
Therefore, many nonlinear HRV analysis methods have been extensively studied in recent years, including the fractal dimension [14], Poincare diagram [15], complexity analysis [16] and other methods, all of which have achieved certain results. However, Li et al. [17] pointed out that few scholars have studied the dynamic changes of short-term HR fluctuations, and most of them only evaluated the overall level of HRV. Nevertheless, information-based similarity (IBS) was proved to be a valuable way of nonlinear methods to study HRV [18,19], many scientists demonstrated that some underlying dynamics existed in the complex fluctuations in the output signals of the physiological system [12,20]. IBS method considers nonlinear dynamic characteristics of a physiological signal. It studies the acceleration/deceleration mode of heart rate fluctuations based on the sequencing frequency analysis. By transforming HR fluctuations into symbols, it can effectively correct the noise made by the antagonism effect of sympathetic nerve system (SNS) and parasympathetic nervous system (PNS) effectively [21]. In this way, the IBS method can be applied to get the instantaneous HR fluctuation differences between CHF patients and normal. Cui et al. [21] found that IBS index could correctly classify Atrial fibrillation (AF) and normal heart rate. Yang et al. [22] used IBS to study β-AR polymorphism, which reflected the underlying pattern of ANS and the difference between ANS changes and dynamic heart rate changes.
Entropy is also a feasible method for CHF study for the complexity of heart rate fluctuations [23,24]. Many studies proved the effectiveness of entropy in CHF screening [25], which supported that the HRV analysis was helpful in CHF diagnosis. Moreover, introducing the concept of “fuzzy set”, fuzzy approximate entropy (fApEn) is an ideal tool to measure the complexity of time series [26]. Li et al. [27] indirectly processed ECG signals to expose ANS imbalance in OSA patients by introducing fuzzy approximate entropy to calculate the complexity of slope series, rather than directly calculating the entropy value of the ECG data. In this way, the fApEn could indirectly estimate the complexity of HR fluctuations during a short time period. Sokunbi et al. [28] found that fApEn not only depends heavily on data length, but also has a better performance in monotonicity, relative consistency, and robustness to noise than other entropy methods when processing different complexities signals.
Therefore, the combination of IBS and fApEn cannot only reveal the law of HR fluctuations but also better study sympathetic disorders by measuring the complexity of different wave patterns, and also reveal the mechanism of CHF patients. First, the IBS method was used to convert the HRV signal into an IBS value. Secondly, fApEn_IBS was extracted to study the complexity of short-term heart rate fluctuations. In addition, LF/HF was also calculated as the contrast, the significance analysis of the three indicators was conducted. Finally, the Fisher discriminant function was used for single index discrimination. Since machine learning methods have been proven to be able to successfully detect and classify heart failures [29,30,31], we decided to use machine learning classifiers for multi-feature discrimination to improve screening accuracy. The samples were divided into normal and CHF groups, and the accuracy, specificity and sensitivity were calculated.

2. System and Method

In this study, HRV characteristics of normal and CHF subjects were analyzed. The framework of the system is shown in Figure 1. Firstly, the 24-h RR-interval signals were extracted after collecting ECG signals. Then, the corrected 1-min RR-interval signals (RRI) were obtained through preprocessing. Furthermore, the traditional HRV research methods and the IBS method were used to calculate the indices. Finally, LF/HF ratio and IBS index were used for statistical analysis and severe CHF detection.

2.1. Data Collection

In this study, the 24-h RR interval signals of 54 healthy subjects (31 males and 23 females, aged 61.38 ± 11.63 years) were collected from the Normal Sinus Rhythms RR Interval database, and 44 CHF subjects (19 males and 6 females, 19 subjects’ gender were unknown, aged 55.51 ± 11.44) were acquired from the Beth Israel Deaconess Medical Center (BIDMC) Congestive Heart Failure database (15 subjects) and the Congestive Heart Failure RR Interval database (29 subjects). All these databases and ECG signals are open source from the PhysioNet open database [32,33,34].

2.2. Preprocess

To reduce the affect produced by noise, the 24-h HRV signals were preprocessed. First, we remove the first and last RR intervals of each 24-h record, as well as the RR intervals, which are longer than 3 s, to remove singular value interference [35]. Then, we did not choose to delete lower outliers because the RR intervals of healthy groups are usually higher than CHF groups [9]. Many researchers considered 5-min series as a standard interval of spectral HRV studies [36]. However, based on the situation the 1-min segments would have a good real-time performance for efficiently capturing changes in heart rate. In this way, RRI may also be suitable for CHF detection [37] and we separated the 24-h HRV signals into non-overlapping 1-min RRIs.

2.3. Indices Calculation

2.3.1. Frequency Domain Analysis

In frequency domain, the low-frequency power (0.04–0.15 Hz, LF) is closely related to the SNS activity, and the high-frequency power (0.15–0.4 Hz, HF) is usually associated with the activity of the PNS. Then, the LF to HF ratio (LF/HF) can reflect the tension balance between SNS and PNS [38]. In this study, a fast Fourier transform would be used to compute the power spectral density for each RRI. The LF/HF ratio was calculated by the LF and HF components index, which is widely accepted to reflect the ANS balance [39]. The formula is shown as:
LH   Ratio = LF   power LH   power

2.3.2. Similarity Change Analysis

In the previous studies, information-based similarity (IBS) of RR-interval signals can successfully analyze the similarity of the regularity in short-term HR fluctuation [21,40]. Then, fuzzy approximate entropy (fApEn) can reflect the complexity of short-term HR fluctuations’ similarity. Based on similarity analysis and entropy measurement, the fApEn_IBS was proposed as a novel HRV analysis method to determine the dynamic rhythm changes of the ANS by observing the complexity changes of fluctuation patterns. The scheme of the fApEn_IBS method is shown in Figure 2, and the details of the fApEn_IBS are as follows.
Step 1 (RRI reconstruction): the raw RRI is defined as follows:
RRI = ( r 1 , r 2 , r 3 , , r k )   k = 1 , 2 , , n
We assume that the coarse-grained is s, that is, the data of every s RR interval will be classified into a set for non-overlapping mean coarse-grained processing. This means that you will get a new variable yj, whose value is the average of the data in the corresponding set:
y j = k = ( j 1 ) s + 1 j s r k s .   j = 1 , 2 , ,   n .
Through non-overlapping mean coarse-grained processing, the original RRI sequence is reconstructed into a new sequence Y, which is formed by yj:
Y = ( y 1 , y 2 , , y   n )   n = [ n s ] .
Figure 3 illustrates the details of RRI reconstruction. As shown in Figure 3, the first s RR interval data r1 to rs are classified into the first set, and y1 is obtained after averaging these s data. By analogy, the new sequence Y composed of y1 to yn replaces the original RRI sequence.
Step 2 (Construction of m-bit words sequence): The increase and decrease between the yi and yi+1 were defined as 1 and 0, respectively. Based on this standard, the Y sequence with n′ length would be converted into a 0–1 binary sequence of length n′-1. Then, a binary sequence of length m is represented as a word, called a m-bit word. Each different m-bit word could reflect a specific pattern of HR fluctuations [40]. In this way, by applying the sliding-window method to Y, it became a series of m-bit words.
Step 3 (Distance calculation): By quantifying the number of each word in the sequence, the frequency of occurrence of different words can be obtained. Moreover, finally, every word’s position in the sequence would depend on its frequency after sorting. To calculate the IBS value between sequences of adjacent words, the formulas are as follows [21]:
w ( x i ) =   [ p 1 ( x i ) log   p 1 ( x i ) p 2 ( x i ) logp 2 ( x i ) ] i = 1 L [ p 1 ( x i ) logp 1 ( x i ) p 2 ( x i ) logp 2 ( x i ) ]
where w is the weight of a decimal m-bit word xi, L represents the total number of different m-bit words, and p denotes the possibility of each word.
D ( R 1 , R 2 ) = i = 1 L | K 1 ( x i ) K 2 ( x i ) | w ( x i ) L
D(R1, R2) denotes the distance between adjacent words sequences, as well as the IBS value between adjacent raw RRIs. K denotes the rank of word.
Step 4 (Average calculation): After calculating the average of all IBS values between adjacent RRIs in one recording, the IBS index would be obtained.
Step 5 (Space reconstruction & Space vectors distance calculation): From the time series containing n IBS values in one recording u(t) = u(1), u(2), ..., u(N), the following m-dimensional vector is defined as:
XX i m = { u ( i ) , , u ( i + m 1 ) } u 0 ( i )
i = 1 , 2 , , n m + 1
where u 0 ( i ) is calculated as:
u 0 ( i ) = k = 0 i + m 1 u ( k ) m
Distance d ij m between vectors XX i m and XX j m is based on maximum absolute difference calculation. For a given similarity tolerance of function r, the similarity degree D ij m is determined by a fuzzy membership function FZ( d ij m , r), which employs the Gaussian function.
D ij m = FZ ( d ij m , r )
Step 6 (Approximation calculation): define the possibility of vectors XX i m and XX j m will match as function φ i m ,
φ i m = j = 1 N m D ij m N m + 1
Moreover, define the function Φ m ( r ) as follows:
Φ m ( r ) = i = 1 N m + 1 ln ( φ i m ) N m + 1
Because the input IBS parameter is a finite time series [41],
fApEn ( m , r , N ) = Φ m ( r ) Φ m + 1 ( r ) = ln ( N m + 1 ) 1 i 1 N m + 1 φ i m + 1 ( N m ) 1 i 1 N m φ i m
Based on these, fApEn_IBS value could be calculated.

2.4. Indice Validation

To evaluate the HRV indices, first we calculate LF/HF ratio and IBS index as well as fApEn_IBS index based on MATLAB (R2017b, MathWorks, Natick, MA, USA). Then, a t-test was used to analyze the significant difference between the control and CHF groups based on SPSS version 22.0.0.0 (SPSS Inc., Chicago, IL, USA), the results were expressed as the mean ± SD. Moreover, p < 0.05 was considered statistically significant. Finally, to verify the performance of these HRV indicators, Fisher linear discriminant of SPSS was used for CHF screening. In addition, based on the Python 3.6.5 environment, we use support vector machine (SVM) [42], k-nearest neighbor (KNN) [43], and Random forests (RF) [44] to perform five cross-validation for multiple features. The results were expressed as accuracy, sensitivity, and specificity. The receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) of the algorithms were calculated using MATLAB (R2017b, MathWorks, Natick, MA, USA).

3. Results

3.1. HRV Analysis of the Difference between the Control and CHF Group

The analysis results of the frequency domain, IBS and fApEn_IBS indices are shown in Table 1 as mean ± SD values. There were significant differences between the healthy and CHF groups (p < 0.001) in the three indices, so LF/HF and IBS, fApEn_IBS indexes are all considered as useful indices to distinguish the control and CHF groups.
Figure 4 shows the representative adjacent raw RRIs between the normal and CHF group. Moreover, it also reveals that the regularity of HR fluctuations between adjacent RRIs in the two groups was different, the regularity of HR fluctuations in CHF group is more similar than normal people. Moreover, as shown in Figure 5, two records were selected from the normal and CHF groups as typical examples to reflect the distribution of IBS index in 24 h. The complexity of color distribution indicated that the degree of disturbance of HR fluctuation was irregular between the two groups.

3.2. CHF Sreening

Using the Fisher linear discriminant function, each sample was classified as normal and CHF patient to access the performance of each index. The accuracy, sensitivity and specificity are defined as the percentage of correctly classified normal or CHF patients, the percentage of correctly classified CHF patients, and the percentage of correctly classified normal samples, respectively. Table 2 shows that fApEn_IBS is the best screening indicator with the highest accuracy (84.69%) and specificity (98.15%) for 98 samples. Therefore, fApEn_IBS improves detection accuracy and is a more useful index than LF/HF and IBS. As shown in Figure 6, two groups could be separated from three-dimensional observations, proving the rationality of screening by three indicators.
When the samples have multiple features and a large amount of data, machine learning classifiers will have significant advantages over other methods [45]. We selected the following three machine learning methods for multi-index screening:
(1)
K Nearest Neighbor (KNN): KNN is a type of instance-based learning. New cases are classified based on a similarity measure in the vector space model. Here the neighbor number is defined as 5 [46].
(2)
Random forest classifier (RF): RF is an ensemble learning method that consists of a multitude of decision trees. The result is determined by the output model of a single tree [47].
(3)
Support Vector Machine (SVM): SVM is a supervised clustering method that maps data points to high-dimensional space through kernel functions for classification. In this paper, we select polynomial kernel function to classify data points [48].
Combining the above five-fold cross-validation analysis, the screening effect is better. With multiple features, including IBS, LF/HF, and fApEn_IBS index, the random forest classifier achieved the best classifier result, as shown in Table 3. Compared with other classifiers, RF achieved a better classification performance under the same input conditions, with the highest accuracy and more balanced sensitivity and specificity (98.98% accuracy, 98.15% sensitivity and 100% specificity). The confusion matrix also points out that only one CHF patient was wrongly identified, which proved that the RF classifier is an effective method for CHF detection. Furthermore, Figure 7 also supports this conclusion through the ROC curve that is the closest to the top with the largest area (AUC = 0.9996).

3.3. Parameter Selection

When calculating the fApEn_IBS index, the influences of different bits of constructed words and the coarse-grained scale s on the CHF screening performance are shown in Figure 8. When s and m are 2 at the same time, the screening performance is optimal. Therefore, s = 2 and m = 2 are selected as the average coarse-granulation scale and the word length of the construction. When m = 2, it can better capture the relationship between adjacent letters, and different numbers are more likely to appear. When s = 2, the information of the original RR interval sequence can be preserved more finely and noise would be eliminated. In conclusion, when s = 2 and m = 2, CHF detection performance is the best and screening accuracy is improved. So, m = 2 and s = 2 are selected as the number of bits of the words and coarse-grained scale to calculate the fApEn_IBS index for HRV analysis.

4. Discussion

4.1. Comparison and Summary

Previous studies have shown that LF/HF is a useful frequency domain indicator for HRV analysis in CHF patients [49,50], and our results also confirmed that there is a statistical difference in LF/HF between healthy people and CHF patients (p < 0.001). However, as the study progressed, the screening accuracy of the LF/HF indicator was not high enough (less than 80%), the single index screening accuracy of IBS and fApEn_IBS was higher.
Furthermore, fApEn_IBS has the most significant performance in distinguishing normal people and CHF patients. Compared with the IBS index, fApEn_IBS can not only analyze the similarity of HR fluctuation, but also analyze the change of fluctuation regularity. Table 2 shows that the accuracy of fApEn_IBS in CHF screening is 84.69%, while the screening accuracy of IBS, LF/HF is 83.67% and only 77.55%, respectively. In addition, the specificity and sensitivity of fApEn_IBS were 98.15% and 68.18%, respectively. The high specificity indicates that fApEn_IBS has strong screening ability for patients with congestive heart failure, and its low sensitivity may be due to the smaller number of positive patients than normal people, so small fluctuations in the number have a greater impact on the calculation of the index. In summary, the fApEn_IBS index is a valid index for CHF screening.

4.2. Comparison with Previous Studies

To expand the sample size, we collected data from three databases. Focusing on study methods and screening results, we compared them with related studies that analyzed ECG data [25,51,52,53,54], as shown in Table 4.
Some of the previous studies only focused on the screening by nonlinear methods [25,51] and did not combine machine learning with improving the screening accuracy. Other studies had smaller sample sizes or inadequate screening accuracy [52,53]. Compared with the results of other studies, the classification accuracy of our method is higher than most with the accuracy (99%), sensitivity (97.8%) and specificity (100%) all above 95%. These results indicate that fApEn_IBS has a certain clinical reference value in evaluating ANS complexity for CHF patients.

4.3. Method Propsed and Parameter Selection

A previous study has reported that 5-min series is the standard time interval for HRV studies [36], but considering that RRI could be keenly aware of heart rate changes, the 1-min segment has good real-time performance. Although the traditional classical frequency-domain method [10] is feasible, the regularity and complexity of HR fluctuations are considered more comprehensively. Therefore, the fApEn_IBS index was constructed by a 1-min RR segment for CHF diagnosis in this paper.
It is found that the interval time series of heartbeats are correlated [55], and it can be changed with some diseases [56]. As shown in Figure 6, CHF patient has more similar HR fluctuations and a higher correlation between adjacent RRIs than normal subject. Therefore, correlation analysis can be used to detect CHF. However, most studies did not quantify the degree of correlation. As a method to measure the similarity between symbolic sequences, IBS had been used effectively to quantify the correlation between interbeat interval time series [18,40].
As a nonlinear system, the function of the cardiopulmonary coupling system is affected by CHF, and the regulation of the system to HR through ANS is also nonlinear [57]. Considering that HRV also contains nonlinear components in signal generation and regulation, compared with the linear method, more information can be obtained by processing HRV signal with the nonlinear method [58]. IBS as a novel nonlinear method, before and after using the numbers reflect the change of HR fluctuations, thus weakening the influence of volatility. Moreover, after converting the 0–1 binary sequences into m-bit words sequences, IBS is only associated with the relative proportion of letters that appeared instead of the absolute number and time node of occurrences of letters. Therefore, the IBS method can analyze the similarity of nonlinear system signals by comparing the variation rule, and this is a reasonable strategy to detect CHF by evaluating the similarity of HR fluctuation regularity.
Complexity is a feature of normal cardiovascular regulation, and the nonlinear regulatory pattern of the cardiovascular system leads to the complexity and irregularity of RRI [59]. As a complexity analysis method, fuzzy approximate entropy is not strongly dependent on data length when processing different complex signals, and has better performance in terms of robustness to noise [28]. Thus, after studying the HR fluctuation rule by IBS sequences, the complexity of IBS sequence would be studied by combining fuzzy approximate entropy algorithm, and the fApEn_IBS index was constructed to observe the complexity difference HR fluctuation between normal people and CHF patients.
To improve index performance and detection accuracy, the parameters of different bits of constructed words and the coarse-grained scale and fApEn were also analyzed. Wu et al. [40] revealed that if the constructed words in sequences for RRI are too long, there would be too few m-bit words. Therefore, we calculated the IBS index when M was 2–6. If the degree of coarse granitization is too large, NAN value will be generated and noise interference cannot be effectively eliminated. In addition, considering the result of retaining original sample data when s value is 1, the value range of S is 1–7. Three parameters of fApEn must be determined: parameter N, which represents the number of sequence points; M is the embedding dimension of sequence reconstruction; r is the similarity tolerance of the exponential function. Sun, Zheng et al. [60,61] proved the feasibility of the method when the parameters N, M and r were set to 2, 2 and 0.25, respectively, and the results were consistent with other results.

4.4. Physiological Significance

CHF is a chronic cardiovascular syndrome associated with ANS dysfunction [6], HRV is a useful tool for evaluating ANS [7,8], and previous studies on the changes of HRV signal energy have attempted to screen CHF patients. In this paper, three indexes of LF/HF, IBS and fApEN_IBS were constructed to analyze the changes of HRV signal by evaluating the balance of ANS function, the similarity and complexity of heart rate fluctuation rhythm. Finally, according to p-value in the Table 1, it was concluded that there were significant differences in three aspects between patients and healthy people.

4.4.1. Why Does ANS Balance Differ Significantly between the Two Groups?

The traditional HRV frequency domain method, LF/HF, reflects the balance of ANS; that is, the level of sympathetic and parasympathetic nervous system activity, demonstrating a significant difference between the normal group and CHF group in our study. Valenza et al. [62] indicated that vagal activity was significantly reduced in patients with CHF. Kishi et al. [63] found an imbalance between the sympathetic nerve and the vagus nerve during a CHF attack. Hasking et al. [64] proposed that to compensate for the pumping function of the heart, patients with CHF increase the level of norepinephrine through a vicious cycle, thereby increasing the sympathetic outflow of the heart and leading to more sympathetic nerve activity. These are the reasons for the statistically significant difference in LF/HF between the two groups, which show a significant difference in ANS balance.

4.4.2. Why Does the Similarity of HR Fluctuations Change Significantly between the Two Groups?

Rahko et al. [65] found that systolic cardiac function was obstructed in patients with CHF, characterized by a significant abnormality in the total systolic interval. Moreover, increased sympathetic activity is accompanied by increased heart rate and energy expenditure, resulting in decreased myocardial oxygen delivery, gradual delay and deceleration of diastole [66,67], which is closely associated with malignant arrhythmias [68]. All these contribute to decreased ability of the autonomic nerve to regulate the change of heart rate and the slow fluctuation of heart rate, leading to the significant difference in the similarity of HR fluctuation rule between the two groups. Figure 6 also confirms that patients’ adjacent RR interval fluctuation is more similar.
Understanding from each point in time, IBS reflects the transient fluctuation of heart rate; understanding from a period of time, and IBS reflects the similarity of the HR fluctuation rule of the sample. CHF patients had lower IBS values, indicating that the adjacent RR intervals fluctuated less and the fluctuation pattern was more similar. With the decrease of oxygen intake, patients cannot change cardiac output by heart rate regulation. The heart cannot fully discharge the venous return of blood from the heart, leading to blood stasis in the venous system and blood perfusion in the arterial system [69]. This increased atrial pressure, which in turn leads to increased capillary pressure, causes cardiac circulation disorders such as congestion in the lungs, and thus patients with CHF are less able to adapt to their environment [70,71].

4.4.3. Why Does the Complexity of HR Fluctuations Differ between the Two Groups?

Saul et al. [72] found that during heart failure, vagal activity decreased. However, sympathetic modulation relatively increased, and these findings reflect that the regulation of HR in CHF patients was significantly disorganized. Naturally, there is a difference in the complexity of HR fluctuations between patients and normal, but dysregulated HR regulation does not mean great complexity. Zhao et al. [73] pointed out that CHF subjects would lose complexity due to pathology. Yin et al. [74] found that the nonlinear complexity of the healthy group was higher than that of the CHF group because the complexity of the cardiovascular system of patients would be reduced, and the entropy value of the complexity of the ECG description would be correspondingly reduced. All of these conclusions are consistent with our experimental results.
The fApEn_IBS of CHF patients was smaller, indicating that the complexity of short-term HR fluctuations in the 24 h was reduced, reflecting that patients are accompanied by impairment of cardiac function and decline of ventricular systolic and diastolic functions [75], and that ANS injury is severe and imbalanced, leading to a relatively simple pattern of short-term heart rate fluctuations. Autonomic nerves do not effectively regulate HR to increase cardiac output and restore systolic/diastolic function [69].
Limitation:
There are some limitations to this study. Firstly, more databases should be used to expand the sample size to verify the effect of the method. Secondly, we only studied the validity of 1-min clips, and further studies on the length of time should be conducted in the future. Thirdly, gender and age inconsistencies between the CHF and control groups may also lead to differences in experimental results. Finally, the effect of drug use on CHF patients should be valued and quantified, excluding the interference of irrelevant variables, so this limitation should be considered in future studies.

5. Conclusions

In this study, the fApEn_IBS method was used to extract the information of different RR intervals. This method not only evaluated the similarity of the HR fluctuation regularity between adjacent RR intervals, but also reflected the complexity of the fluctuation rule, thus significantly improving the CHF screening performance. Our results showed that the similarity and complexity of HR fluctuation law between adjacent RRI in the CHF group were significantly reduced, and the accuracy of fApEn_IBS in CHF detection was 84.69%, and the accuracy of fApEn_IBS, IBS and LF/HF combined with random forest classifier was 98.98%. Therefore, the fApEn_IBS method can be used to analyze the changes of short-term HR fluctuation regularity and detect the patients with CHF.

Author Contributions

Conceptualization, G.L.; methodology, Z.L. and G.L.; software, Z.L. and T.C.; validation, Z.L. and K.W.; formal analysis, Z.L.; investigation, B.L. and Z.L.; resources, G.L.; data curation, T.C. and K.W.; writing—original draft preparation, Z.L.; writing—review and editing, B.L. and G.L.; visualization, Z.L. and T.C.; supervision, B.L.; project administration, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Shenzhen Science and Technology Plan for Fundamental Research, grant number JCYJ20180307153213863 and JCY20190807162003696, JCYJ20170818163426597.

Institutional Review Board Statement

This research is approved by the Key Laboratory of Sensing Technology and Biomedical Instruments of Guangdong Province, School of Biomedical Engineering, Sun Yat-sen University, Guangzhou 510275, China, and with the Guangdong Provincial Engineering and Technology Centre of Advanced and Portable Medical Device, Guangzhou 510275, China.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dougherty, A.H.; Naccarelli, G.V.; Gray, E.L.; Hicks, C.H.; Goldstein, R.A. Congestive heart failure with normal systolic function. Am. J. Cardiol. 1984, 54, 778–782. [Google Scholar] [CrossRef]
  2. Garg, R.; Packer, M.; Pitt, B.; Yusuf, S. Heart failure in the 1990s: Evolutions of a major public health problem in cardio-vascular medicine. J. Am. Coll. Cardiol. 1993, 22, A3–A5. [Google Scholar] [CrossRef] [Green Version]
  3. Mårtensson, J.; Karlsson, J.E.; Fridlund, B. Male patients with congestive heart failure and their conception of the life situation. J. Adv. Nurs. 2010, 25, 579–586. [Google Scholar] [CrossRef]
  4. Poppas, A.; Rounds, S. Congestive Heart Failure. Am. J. Respir. Crit. Care Med. 2002, 165, 4–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Dunn, G.D.; Hayes, P.; Breen, K.J.; Schenker, S. The liver in congestive heart failure: A review. Am. J. Med. Sci. 1973, 265, 174–189. [Google Scholar] [CrossRef]
  6. Ogawa, M.; Zhou, S.; Tan, A.Y.; Song, J.; Gholmieh, G.; Fishbein, M.C.; Luo, H.; Siegel, R.J.; Karagueuzian, H.S.; Chen, L.S.; et al. Left Stellate Ganglion and Vagal Nerve Activity and Cardiac Arrhythmias in Ambulatory Dogs With Pacing-Induced Congestive Heart Failure. J. Am. Coll. Cardiol. 2007, 50, 335–343. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, W.; Zheng, L.; Li, K.; Wang, Q.; Liu, G.; Jiang, Q. A Novel and Effective Method for Congestive Heart Failure Detection and Quantification Using Dynamic Heart Rate Variability Measurement. PLoS ONE 2016, 11, e0165304. [Google Scholar] [CrossRef] [PubMed]
  8. Ravenswaaij-Arts, V.; Conny, M.A. Heart rate variability. Ann. Intern. Med. 1993, 118, 436. [Google Scholar] [CrossRef]
  9. Takase, B.; Kurita, A.; Noritake, M.; Uehata, A.; Maruyama, T.; Nagayoshi, H.; Nishioka, T.; Mizuno, K.; Nakamura, H. Heart rate variability in patients with diabetes mellitus, ischemic heart disease, and congestive heart failure. J. Electrocardiol. 1992, 25, 79–88. [Google Scholar] [CrossRef]
  10. La Rovere, M.T.; Pinna, G.D.; Maestri, R.; Mortara, A.; Capomolla, S.; Febo, O.; Ferrari, R.; Franchini, M.; Gnemmi, M.; Opasich, C.; et al. Short-Term Heart Rate Variability Strongly Predicts Sudden Cardiac Death in Chronic Heart Failure Patients. Circulation 2003, 107, 565–570. [Google Scholar] [CrossRef] [PubMed]
  11. Lo, M.T.; Tsai, P.H.; Lin, P.F.; Lin, C.; Hsin, Y.L. The nonlinear and nonstationary properties in EEG signals: Probing the complex fluctuations by Hilbert–Huang transform. Adv. Adapt. Data Anal. 2011, 1, 461–482. [Google Scholar] [CrossRef]
  12. Peng, C.-K.; Costa, M.; Goldberger, A.L. Adaptive data analysis of complex fluctuations in physiologic time series. Adv. Adapt. Data Anal. 2009, 1, 61–70. [Google Scholar] [CrossRef] [PubMed]
  13. Goldberger, A.L.; Amaral, L.; Hausdorff, J.M.; Ivanov, P.C.; Peng, C.K.; Stanley, H.E. Fractal dynamics in physiology: Alteraions with disease and aging. Proc. Natl. Acad. Sci. USA 2002, 99 (Suppl. 1), 2466–2472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Beckers, F.; Verheyden, B.; Couckuyt, K.; Aubert, A.E. Fractal dimension in health and heart failure. Biomed. Tech. Eng. 2006, 51, 194–197. [Google Scholar] [CrossRef] [PubMed]
  15. Brennan, M.; Palaniswami, M.; Kamen, P. Do Existing Measures of Poincare Plot Geometry Reflect Nonlinear Features of Heart Rate Variability? IEEE Trans. Biomed. Eng. 2001, 48, 1342–1347. [Google Scholar] [CrossRef]
  16. Khandoker, A.H.; Jelinek, H.F.; Palaniswami, M. Identifying diabetic patients with cardiac autonomic neuropathy by heart rate complexity analysis. Biomed. Eng. Online 2009, 8, 3–12. [Google Scholar] [CrossRef] [Green Version]
  17. Li, Y.; Wu, S.; Yang, Q.; Liu, G.; Ge, L. Application of the Variance Delay Fuzzy Approximate Entropy for Autonomic Nervous System Fluctuation Analysis in Obstructive Sleep Apnea Patients. Entropy 2020, 22, 915. [Google Scholar] [CrossRef]
  18. Peng, C.K.; Yang, C.C.; Goldberger, A.L. Statistical physics approach to categorize biologic signals: From heart rate dynamics to DNA sequences. Chaos 2007, 17, 88–433. [Google Scholar] [CrossRef]
  19. Yang, C.C.; Hseu, S.S.; Yien, H.W.; Goldberger, A.L.; Peng, C.K. Linguistic Analysis of the Human Heartbeat Using Fre-quency and Rank Order Statistics. Phys. Rev. Lett. 2003, 90, 108103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. 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]
  21. Cui, X.; Chang, E.; Yang, W.-H.; Jiang, B.C.; Yang, A.C.; Peng, C.-K. Automated Detection of Paroxysmal Atrial Fibrillation Using an Information-Based Similarity Approach. Entropy 2017, 19, 677. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, A.C.; Tsai, S.-J.; Hong, C.-J.; Wang, C.; Chen, T.-J.; Liou, Y.-J.; Peng, C.-K. Clustering Heart Rate Dynamics Is Associated with β-Adrenergic Receptor Polymorphisms: Analysis by Information-Based Similarity Index. PLoS ONE 2011, 6, e19232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pincus, S.M.; Goldberger, A.L. Physiological time-series analysis: What does regularity quantify? Am. J. Physiol. Circ. Physiol. 1994, 266, H1643–H1656. [Google Scholar] [CrossRef] [PubMed]
  24. Richman, S.; Joshua, M.; Randall, J. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Pan, W.; He, A.; Feng, K.; Li, Y.; Wu, D.; Liu, G. Multi-Frequency Components Entropy as Novel Heart Rate Variability Indices in Congestive Heart Failure Assessment. IEEE Access 2019, 7, 37708–37717. [Google Scholar] [CrossRef]
  26. Chen, W.; Wang, Z.; Xie, H.; Yu, W. Characterization of Surface EMG Signal Based on Fuzzy Entropy. IEEE Trans. Neural Syst. Rehabil. Eng. 2007, 15, 266–272. [Google Scholar] [CrossRef]
  27. Li, Y.; Pan, W.; Li, K.; Jiang, Q.; Liu, G. Sliding Trend Fuzzy Approximate Entropy as a Novel De-scriptor of Heart Rate Variability in Obstructive Sleep Apnea. IEEE J. Biomed. Health Inform. 2019, 23, 175–183. [Google Scholar] [CrossRef] [PubMed]
  28. Sokunbi, M.O.; Cameron, G.G.; Ahearn, T.S.; Murray, A.D.; Staff, R.T. Fuzzy approximate entropy analysis of resting state fMRI signal complexity across the adult life span. Med. Eng. Phys. 2015, 37, 1082–1090. [Google Scholar] [CrossRef]
  29. Melillo, P.; Fusco, R.; Sansone, M.; Bracale, M.; Pecchia, L. Discrimination power of long-term heart rate variability measures for chronic heart failure detection. Med. Biol. Eng. Comput. 2011, 49, 67–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Leandro, P.; Paolo, M.; Marcello, B. Remote health monitoring of heart failure with data mining via CART method on HRV features. IEEE Trans. Biomed. Eng. 2010, 58, 800–804. [Google Scholar]
  31. İşler, Y.; Kuntalp, M. Combining classical HRV indices with wavelet entropy measures improves to performance in diagnosing congestive heart failure. Comput. Biol. Med. 2007, 37, 1502–1510. [Google Scholar] [CrossRef]
  32. Krum, H.; Bigger, T.; Goldsmith, R.L.; Packer, M. Effect of long-term digoxin therapy on autonomic function in patients with chronic heart failure. J. Am. Coll. Cardiol. 1995, 25, 289–294. [Google Scholar] [CrossRef] [Green Version]
  33. Goldsmith, R.L.; Bigger, J.T.; Bloomfield, D.M.; Krum, H.; Steinman, R.C.; Sackner-Bernstein, J.; Packer, M. Long-term carve-dilol therapy increases parasympathetic nervous system activity in chronic congestive heart failure. Am. J. Cardiol. 1997, 80, 1101–1104. [Google Scholar] [CrossRef]
  34. E Mietus, J.; Peng, C.-K.; Henry, I.; Goldsmith, R.L.; Goldberger, A.L. The pNNx files: Re-examining a widely used heart rate variability measure. Heart 2002, 88, 378–380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. 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]
  36. Singh, J.; Sharma, R.K.; Gupta, A.K. A method of REM-NREM sleep distinction using ECG signal for unobtrusive personal monitoring. Comput. Biol. Med. 2016, 78, 138–143. [Google Scholar] [CrossRef] [PubMed]
  37. Michael, S.; Jay, O.; Halaki, M.; Graham, K.; Davis, G.M. Submaximal exercise intensity modulates acute post-exercise heart rate variability. Graefe’s Arch. Clin. Exp. Ophthalmol. 2016, 116, 697–706. [Google Scholar] [CrossRef] [PubMed]
  38. Berntson, G.G.; Thomas Bigger, J., Jr.; Eckberg, D.L.; Grossman, P.; Kaufmann, P.G.; Malik, M.; Nagaraja, H.N.; Porges, S.W.; Saul, J.P.; Stone, P.H.; et al. Heart rate variability: Origins, methods, and interpretive caveats. Psychophysiology 1997, 34, 623–648. [Google Scholar] [CrossRef] [PubMed]
  39. Saboul, D.; Pialoux, V.; Hautier, C. The breathing effect of the LF/HF ratio in the heart rate variability measurements of athletes. Eur. J. Sport Sci. 2011, 14, S282–S288. [Google Scholar] [CrossRef] [PubMed]
  40. Wu, S.; Liang, D.; Yang, Q.; Liu, G. Regularity of heart rate fluctuations analysis in obstructive sleep apnea patients using information-based similarity. Biomed. Signal Process. Control. 2020, 65, 102370. [Google Scholar] [CrossRef]
  41. Kumar, Y.; Dewal, M.; Anand, R. Epileptic seizure detection using DWT based fuzzy approximate entropy and support vector machine. Neurocomputing 2014, 133, 271–279. [Google Scholar] [CrossRef]
  42. Ramírez, J.; Monasterio, V.; Mincholé, A.; Llamedo, M.; Lenis, G.; Cygankiewicz, I.; de Luna, A.B.; Malik, M.; Martínez, J.P.; Laguna, P.; et al. Automatic SVM classification of sudden cardiac death and pump failure death from autonomic and repolarization ECG markers. J. Electrocardiol. 2015, 48, 551–557. [Google Scholar] [CrossRef] [PubMed]
  43. Liu, G.; Wang, L.; Wang, Q.; Zhou, G.; Wang, Y.; Jiang, Q. A New Approach to Detect Congestive Heart Failure Using Short-Term Heart Rate Variability Measures. PLoS ONE 2014, 9, e93399. [Google Scholar] [CrossRef]
  44. Masetic, Z.; Subasi, A. Congestive heart failure detection using random forest classifier. Comput. Methods Programs Biomed. 2016, 130, 54–64. [Google Scholar] [CrossRef] [PubMed]
  45. Ngiam, K.Y.; Khor, W. Big data and machine learning algorithms for health-care delivery. Lancet Oncol. 2019, 20, e262–e273. [Google Scholar] [CrossRef]
  46. Poddar, M.G.; Birajdar, A.C.; Virmani, J. Automated Classification of Hypertension and Coronary Artery Disease Patients by PNN, KNN, and SVM Classifiers Using HRV Analysis; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  47. Alarsan, F.I.; Younes, M. Mamoon, Analysis and classification of heart diseases using heartbeat features and machine learning algorithms. J. Big Data 2019, 6, 1–15. [Google Scholar] [CrossRef] [Green Version]
  48. Lameski, P.; Zdravevski, E.; Koceski, S.; Kulakov, A.; Trajkovik, V. Vladimir, Suppression of intensive care unit false alarms based on the arterial blood pressure signal. IEEE Access 2017, 5, 5829–5836. [Google Scholar]
  49. Wang, Y.; Wei, S.; Zhang, S.; Zhang, Y.; Zhao, L.; Liu, C.; Murray, A. Comparison of time-domain, frequency-domain and non-linear analysis for distinguishing congestive heart failure patients from normal sinus rhythm subjects. Biomed. Signal Process. Control. 2018, 42, 30–36. [Google Scholar] [CrossRef] [Green Version]
  50. Pecchia, L.; Melillo, P.; Sansone, M.; Bracale, M. Discrimination power of short-term heart rate variability measures for CHF assessment. IEEE Trans. Inf. Technol. Biomed. 2011, 15, 40–46. [Google Scholar] [CrossRef]
  51. Luo, D.; Pan, W.; Li, Y.; Feng, K.; Liu, G. The Interaction Analysis between the Sympathetic and Parasympathetic Systems in CHF by Using Transfer Entropy Method. Entropy 2018, 20, 795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Cornforth, D.; Jelinek, H.F. Detection of Congestive Heart Failure using Renyi Entropy. In 2016 Computing in Cardiology Conference (CinC); IEEE: Manhattan, NY, USA, 2016. [Google Scholar] [CrossRef]
  53. Chen, W.; Liu, G.; Su, S.; Jiang, Q.; Nguyen, H. A CHF detection method based on deep learning with RR intervals. In Proceedings of the 2017 39th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Jeju, Korea, 11–15 July 2017; pp. 3369–3372. [Google Scholar]
  54. Wang, L.; Zhou, W.; Liu, N.; Xing, Y.; Zhou, X. CHF Detection with LSTM Neural Network. In Proceedings of the 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Honolulu, HI, USA, 18–21 July 2018. [Google Scholar]
  55. Iyengar, N.; Peng, C.-K.; Morin, R.; Goldberger, A.L.; Lipsitz, L.A. Age-related alterations in the fractal scaling of cardiac interbeat interval dynamics. Am. J. Physiol. Integr. Comp. Physiol. 1996, 271, R1078–R1084. [Google Scholar] [CrossRef] [Green Version]
  56. Peng, C.K.; Hausdorff, J.M.; Havlin, S.; Mietus, J.E.; Stanley, H.E.; Goldberger, A.L. Multiple-time scales analysis of physiological time series under neural control. Phys. A Stat. Mech. Appl. 1998, 249, 491–500. [Google Scholar] [CrossRef]
  57. Hoyer, D.; Schmidt, K.; Bauer, R.; Zwiener, U.; Kohler, M.; Luthke, B.; Eiselt, M. Nonlinear analysis of heart rate and respiratory dynamics. IEEE Eng. Med. Biol. Mag. 1997, 16, 31–39. [Google Scholar] [CrossRef] [PubMed]
  58. Joo, S.; Choi, K.-J.; Huh, S.-J. Prediction of spontaneous ventricular tachyarrhythmia by an artificial neural network using parameters gleaned from short-term heart rate variability. Expert Syst. Appl. 2011, 39, 3862–3866. [Google Scholar] [CrossRef]
  59. Batchinsky, A.I.; Cancio, L.C.; Salinas, J.; Kuusela, T.; Cooke, W.H.; Wang, J.J.; Boehme, M.; Convertino, V.A.; Holcomb, J.B. Prehospital Loss of R-to-R Interval Complexity is Associated with Mortality in Trauma Patients. J. Trauma: Inj. Infect. Crit. Care 2007, 63, 512–518. [Google Scholar] [CrossRef] [Green Version]
  60. Ao, D.; Sun, R.; Song, R. Comparison of complexity of EMG signals between a normal subject and a patient after stroke—A case study. In Proceedings of the 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 3–7 July 2013; pp. 4965–4968. [Google Scholar] [CrossRef]
  61. Zheng, L.; Pan, W.; Li, Y.; Luo, D.; Wang, Q.; Liu, G. Use of mutual information and transfer entropy to assess interaction between parasympathetic and sympathetic activities of nervous system from HRV. Entropy 2017, 19, 489. [Google Scholar] [CrossRef] [Green Version]
  62. Valenza, G.; Citi, L.; Saul, P.; Barbieri, R. ECG-Derived Sympathetic and Parasympathetic Nervous System Dynamics: A Congestive Heart Failure Study. In Proceedings of the 2018 Computing in Cardiology Conference (CinC), Maastricht, The Netherlands, 23–26 September 2018. [Google Scholar] [CrossRef]
  63. Kishi, T. Heart failure as an autonomic nervous system dysfunction. J. Cardiol. 2012, 59, 117–122. [Google Scholar] [CrossRef] [Green Version]
  64. Hasking, G.J.; Esler, M.D.; Jennings, G.L.; Burton, D.; Johns, J.A.; Korner, P.I. Norepinephrine spillover to plasma in patients with congestive heart failure: Evidence of increased overall and cardiorenal sympathetic nervous activity. Circulation 1986, 73, 615–621. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Rahko, P.S.; Shaver, J.A.; Salerni, R.; Uretsky, B.F. Noninvasive evaluation of systolic and diastolic function in severe congestive heart failure secondary to coronary artery disease or idiopathic dilated cardiomyopathy. Am. J. Cardiol. 1986, 57, 1315–1322. [Google Scholar] [CrossRef]
  66. Böhm, M.; Perez, A.-C.; Jhund, P.S.; Reil, J.C.; Komajda, M.; Zile, M.R.; McKelvie, R.S.; Anand, I.S.; Massie, B.M.; Carson, P.E.; et al. Relationship between heart rate and mortality and morbidity in the irbesartan patients with heart failure and preserved systolic function trial (I-Preserve). Eur. J. Heart Fail. 2014, 16, 778–787. [Google Scholar] [CrossRef] [PubMed]
  67. A Mulieri, L.; Hasenfuss, G.; Leavitt, B.; Allen, P.D.; Alpert, N.R. Altered myocardial force-frequency relation in human heart failure. Circulation 1992, 85, 1743–1750. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Houle, M.S.; Billman, G.E. Low-frequency component of the heart rate variability spectrum: A poor marker of sympathetic activity. Am. J. Physiol. Heart Circ. Physiol. 1999, 276, H215–H223. [Google Scholar] [CrossRef] [PubMed]
  69. Tanabe, Y.; Takahashi, M.; Hosaka, Y.; Ito, M.; Ito, E.; Suzuki, K. Prolonged recovery of cardiac output after maximal exercise in patients with chronic heart failure. J. Am. Coll. Cardiol. 2000, 35, 1228–1236. [Google Scholar] [CrossRef] [Green Version]
  70. Kemp, C.D.; Conte, J.V. The pathophysiology of heart failure. Cardiovasc. Pathol. 2012, 21, 365–371. [Google Scholar] [CrossRef] [PubMed]
  71. Picano, E.; Gargani, L.; Gheorghiade, M. Why, when, and how to assess pulmonary congestion in heart failure: Patho-physiological, clinical, and methodological implications. Heart Fail. Rev. 2010, 15, 63–72. [Google Scholar] [CrossRef] [PubMed]
  72. Saul, J.P.; Arai, Y.; Berger, R.D.; Lilly, L.S.; Colucci, W.S.; Cohen, R.J. Assessment of autonomic regulation in chronic congestive heart failure by heart rate spectral analysis. Am. J. Cardiol. 1988, 61, 1292–1299. [Google Scholar] [CrossRef]
  73. Zhao, L.; Wei, S.; Zhang, C.; Zhang, Y.; Jiang, X.; Liu, F.; Liu, C. Determination of Sample Entropy and Fuzzy Measure Entropy Parameters for Distinguishing Congestive Heart Failure from Normal Sinus Rhythm Subjects. Entropy 2015, 17, 6270–6288. [Google Scholar] [CrossRef] [Green Version]
  74. Yin, J.; Xiao, P.; Li, J.; Liu, Y.; Yan, C.; Zhang, Y. Parameters Analysis of Sample Entropy, Permutation Entropy and Permutation Ratio Entropy for RR Interval Time Series. Inf. Process. Manag. 2020, 57, 102283. [Google Scholar] [CrossRef]
  75. Qin, F.; Shite, J.; Liang, C.-S. Antioxidants attenuate myocyte apoptosis and improve cardiac function in CHF: Association with changes in MAPK pathways. Am. J. Physiol. Circ. Physiol. 2003, 285, H822–H832. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Framework of the proposed HRV analysis system.
Figure 1. Framework of the proposed HRV analysis system.
Entropy 23 01669 g001
Figure 2. Scheme of similarity change analysis.
Figure 2. Scheme of similarity change analysis.
Entropy 23 01669 g002
Figure 3. An illustration of RRI reconstruction. (a) A 1-min RRI = ( r 1 , r 2 , r 3 , , r n ) . (b) Multiple non-overlapping rks were divided from RRI. (c) The reconstructed sequence Y is composed of a time duration of 1/s minute.
Figure 3. An illustration of RRI reconstruction. (a) A 1-min RRI = ( r 1 , r 2 , r 3 , , r n ) . (b) Multiple non-overlapping rks were divided from RRI. (c) The reconstructed sequence Y is composed of a time duration of 1/s minute.
Entropy 23 01669 g003
Figure 4. Adjacent raw RRIs of a recording from the Normal group (a), the CHF group (b). The blue line represents the RR(i−1), the gray line represents the RR(i), (2 ≤ i ≤ n, and n is the total amount of RRI in the selected recordings).
Figure 4. Adjacent raw RRIs of a recording from the Normal group (a), the CHF group (b). The blue line represents the RR(i−1), the gray line represents the RR(i), (2 ≤ i ≤ n, and n is the total amount of RRI in the selected recordings).
Entropy 23 01669 g004
Figure 5. Distributions of 24-h period IBS index of a typical subject from the Normal group (a), CHF group (b).
Figure 5. Distributions of 24-h period IBS index of a typical subject from the Normal group (a), CHF group (b).
Entropy 23 01669 g005
Figure 6. A three-dimensional image for the normal and CHF groups with LF/HF, fApEn_IBS, IBS index as the x, y, z-axis.
Figure 6. A three-dimensional image for the normal and CHF groups with LF/HF, fApEn_IBS, IBS index as the x, y, z-axis.
Entropy 23 01669 g006
Figure 7. ROC curves of KNN, SVM and RF classifier for the Normal and CHF groups.
Figure 7. ROC curves of KNN, SVM and RF classifier for the Normal and CHF groups.
Entropy 23 01669 g007
Figure 8. Comparison of CHF screening accuracy by the fApEn_IBS index with different m, s. fApEn_IBS: fuzzy approximate entropy of IBS index. m: the number of bits of the constructed words; s: coarse-grained scale.
Figure 8. Comparison of CHF screening accuracy by the fApEn_IBS index with different m, s. fApEn_IBS: fuzzy approximate entropy of IBS index. m: the number of bits of the constructed words; s: coarse-grained scale.
Entropy 23 01669 g008
Table 1. LF/HF, IBS and changes in the similarity of HR fluctuations indices for the CHF and normal groups.
Table 1. LF/HF, IBS and changes in the similarity of HR fluctuations indices for the CHF and normal groups.
When m = 2; s = 2
IndicesCHFNormalp-Value
LF/HF1.2408 ± 0.579422.6092 ± 1.069020.000 ***
IBS0.1770 ± 0.037440.2273 ± 0.018740.000 ***
fApEn_IBS1.5710 ± 0.142791.7229 ± 0.035110.000 ***
LF/HF: low frequency (LF)/high frequency (HF) ratio; IBS: information-based similarity index; fApEn_IBS: fuzzy approximate entropy of IBS index. CHF: CHF group; Normal: normal group. *** represents p < 0.001.
Table 2. Performance comparisons of CHF screening indicators by Fisher linear discriminant function.
Table 2. Performance comparisons of CHF screening indicators by Fisher linear discriminant function.
IndicesTNTPFNFPAcc (%)Sen (%)Spe (%)
LF/HF383861677.5586.3670.37
IBS493311583.677590.74
fApEn_IBS533014184.6968.1898.15
LF/HF: low frequency (LF)/high frequency (HF) ratio; IBS: information-based similarity index; fApEn_IBS: fuzzy approximate entropy of IBS index. TP: true positive; TN: true negative; FP: false positive; FN: false negative; Acc: accuracy; Sen: sensitivity; Spe: specificity.
Table 3. Performance comparisons of CHF screening indicators by Fisher linear discriminant function.
Table 3. Performance comparisons of CHF screening indicators by Fisher linear discriminant function.
IndicesClassifierTNTPFNFPAcc (%)Sen (%)Spe (%)
LF/HFRF54431098.9897.73100
IBSKNN48431692.8697.7388.89
fApEn_IBSSVM48368685.7181.8288.89
LF/HF: low frequency (LF)/high frequency (HF) ratio; IBS: information-based similarity index; fApEn_IBS: fuzzy approximate entropy of IBS.
Table 4. Comparison of classification results between previous studies and our method.
Table 4. Comparison of classification results between previous studies and our method.
ReferenceFeatureNumber of RecordingsLength of RR SegmentClassifierClassification Result
Pan et al. [25]Single-feature—MFCN (Multi-Frequency Components Entropy)985 minFisher DiscriminantAcc = 86.7%
Sen = 79.5%
Spe = 92.6%
Luo et al. [51]Multi-feature—LF/HF+TE(LF→HF)+TE(HF→LF) Transfer Entropy (TE)985 minFisher DiscriminantAcc = 83.7%
Sen = 86.4%
Spe = 81.5%
David et al. [52]Multi-feature—Renyi Entropy+SDNN+RMSSD3313 minNearest NeighborAcc = 87.9%
Sen = 80.0%
Spe = 94.4%
Chen et al. [53]Multi-feature—50 features1165 minA two-layer deep neural network model based on an SAE-based DL algorithmAcc = 72.4%
Sen = -
Spe = -
Wang et al. [54]No need10173.9 sLong Short-Term MemoryAcc = 85.1%
Sen = 73.6%
Spe = 91.8%
Other methodMulti-feature—fApEn_IBS+IBS+LF/HF981 minRandom ForcesAcc = 99.0%
Sen = 97.8%
Spe = 100.0%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Z.; Chen, T.; Wei, K.; Liu, G.; Liu, B. Similarity Changes Analysis for Heart Rate Fluctuation Regularity as a New Screening Method for Congestive Heart Failure. Entropy 2021, 23, 1669. https://doi.org/10.3390/e23121669

AMA Style

Liu Z, Chen T, Wei K, Liu G, Liu B. Similarity Changes Analysis for Heart Rate Fluctuation Regularity as a New Screening Method for Congestive Heart Failure. Entropy. 2021; 23(12):1669. https://doi.org/10.3390/e23121669

Chicago/Turabian Style

Liu, Zeming, Tian Chen, Keming Wei, Guanzheng Liu, and Bin Liu. 2021. "Similarity Changes Analysis for Heart Rate Fluctuation Regularity as a New Screening Method for Congestive Heart Failure" Entropy 23, no. 12: 1669. https://doi.org/10.3390/e23121669

APA Style

Liu, Z., Chen, T., Wei, K., Liu, G., & Liu, B. (2021). Similarity Changes Analysis for Heart Rate Fluctuation Regularity as a New Screening Method for Congestive Heart Failure. Entropy, 23(12), 1669. https://doi.org/10.3390/e23121669

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