Next Article in Journal
Fisher Information of Free-Electron Landau States
Next Article in Special Issue
Complexity of Body Movements during Sleep in Children with Autism Spectrum Disorder
Previous Article in Journal
Fault Detection Based on Multi-Dimensional KDE and Jensen–Shannon Divergence
Previous Article in Special Issue
A Multi-Class Automatic Sleep Staging Method Based on Photoplethysmography Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Short-Term HRV Analysis Using Nonparametric Sample Entropy for Obstructive Sleep Apnea

1
School of Biomedical Engineering, Sun Yat-Sen University, Guangzhou 510275, China
2
Key Laboratory of Sensing Technology and Biomedical Instruments of Guangdong Province, School of Engineering, Sun Yat-Sen University, Guangzhou 510275, China
3
Guangdong Provincial Engineering and Technology Centre of Advanced and Portable Medical Device, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 8 December 2020 / Revised: 12 January 2021 / Accepted: 13 January 2021 / Published: 24 February 2021
(This article belongs to the Special Issue Entropy and Sleep Disorders II)

Abstract

:
Obstructive sleep apnea (OSA) is associated with reduced heart rate variability (HRV) and autonomic nervous system dysfunction. Sample entropy (SampEn) is commonly used for regularity analysis. However, it has limitations in processing short-term segments of HRV signals due to the extreme dependence of its functional parameters. We used the nonparametric sample entropy (NPSampEn) as a novel index for short-term HRV analysis in the case of OSA. The manuscript included 60 6-h electrocardiogram recordings (20 healthy, 14 mild-moderate OSA, and 26 severe OSA) from the PhysioNet database. The NPSampEn value was compared with the SampEn value and frequency domain indices. The empirical results showed that NPSampEn could better differentiate the three groups (p < 0.01) than the ratio of low frequency power to high frequency power (LF/HF) and SampEn. Moreover, NPSampEn (83.3%) approached a higher OSA screening accuracy than the LF/HF (73.3%) and SampEn (68.3%). Compared with SampEn (|r| = 0.602, p < 0.05), NPSampEn (|r| = 0.756, p < 0.05) had a significantly stronger association with the apnea-hypopnea index (AHI). Hence, NPSampEn can fully overcome the influence of individual differences that are prevalent in biomedical signal processing, and might be useful in processing short-term segments of HRV signal.

1. Introduction

Obstructive sleep apnea (OSA) is the most common type of apnea, which is characterized by repetitive collapse or partial collapse of the pharyngeal airway during sleep [1]. OSA is common in patients with high blood pressure, coronary artery disease, stroke, and atrial fibrillation [2,3]. The primary pathophysiological mechanism is cardiovascular autonomic dysfunction, which mainly manifests as sympathetic overactivation [4]. The prevalence of OSA is 4% in men and 2% in women at the rate of 10 events per hour, but more than 85% of OSA patients never receive a definite diagnosis [5,6,7].
Changes of respiration and blood pressure affect the balance of the autonomic nervous system (ANS), thus affecting heart rate variability (HRV) [8]. The ANS modulates physiological function through the sympathetic nervous system (SNS) and parasympathetic nervous system (PNS) [9]. The parameters of HRV have been proven to be significant indices for analyzing and diagnosing OSA [10]. Traditionally, linear methods such as time and frequency domain analyses have been used to analyze OSA. The ratio of low-frequency to high-frequency power (LF/HF) is the classic frequency domain index, which is considered to reflect the balance of the ANS. A higher LF/HF indicates that the sympathetic nerve is more active, and a lower LF/HF indicates that the parasympathetic nerve is more active [11]. Gula, L.J. et al. found that the value of LF/HF was higher in patients with moderate OSA than that in healthy subjects and severe OSA patients, and Gong, X. et al. found that the value of the apnea-hypopnea index (AHI) was positively correlated with the LF/HF [12,13]. These results suggest that the LF/HF is effective in discriminating individuals with and without OSA.
The autonomic nervous system is nonlinear and nonstationary [14,15]; thus, nonlinear methods are widely used for OSA detection. Entropy is considered as a nonlinear method that is used to analyze the complexity of ANS. Sample entropy (SampEn) analysis has been used to show that the HRV in normal individuals is significantly more complex than that in individuals with OSA (p < 0.005) by Al-Angari, H.M. et al. [16]. Ravelo-Garcia, A. G et al. used permutation entropy to better detect OSA with HRV analysis [17]. Sliding trend fuzzy approximate entropy has been proposed for use in HRV analysis, and demonstrated high accuracy (85%) with a sensitivity of 82.5% and a specificity of 90% for OSA detection [18].
For HRV analysis, most of the previous studies used a 5-min segment of electrocardiogram (ECG) ECG signal to assess the ANS [12,19]. Nevertheless, reference annotations were typically created for every minute with 60–100 R peak to R peak (RR) RR intervals to indicate the presence or absence of apnea during that minute [20,21]. Furthermore, SampEn is highly dependent on parameter r, which is the threshold of distance. SampEn is not suitable for short-term HRV analysis because an improper r value may cause the SampEn value to be undefined [22]. Several studies have also reported that inappropriate r was randomly selected in the SampEn calculation, which led to misleading regular results [23,24]. Therefore, in this study, a method based on non-parametric sample entropy (NPSampEn) is proposed to reflect the influence of self-fluctuation and overcome individual differences, thus improving the performance of OSA screening and quantification.
For NPSampEn, the r sequence was generated based on the unique template distance of the signal. The distribution of the template distance was estimated by r sequence to evaluate the regularity of the signal and the fluctuation of the autonomic nerve. The process was as follows. First, the distance matrix of 1-minute RR segments was calculated as in the SampEn analysis. Then, the extracted sequence containing unique points from the distance matrix was used as the packet interval to calculate the cumulative probability and identify abnormal sympathetic signal variance. Finally, the NPSampEn was used to detect sleep apnea and quantify the severity of OSA disease.

2. Materials and Methods

In our study, a new nonlinear analysis method was used for short-term HRV analysis. The framework of the study is shown in Figure 1. First, the ECG signals were downloaded from the database, and the RR intervals were extracted. Then, the corrected 1-min RR segments were obtained and resampled to 2 Hz by interpolation. Then, short-term HRV measures and the NPSampEn were used to determine the indices in the frequency domain and nonlinear analyses. Finally, these short-term indices were validated by statistical analysis and OSA screening.

2.1. Data

To develop and evaluate our HRV-based algorithm, we used the Computers in Cardiology Challenge 2000 of the PhysioNet database [21]. It contains 70 single channel overnight ECG recordings from 32 subjects. The dataset provided a minute-by-minute sequence of reference annotations for the ECG recordings. Based on the apnea-hypopnea index (AHI), which was calculated by determining the sum of apneas and hypopnea events per hour during sleep [25], all recordings were divided into three groups. Recordings with a 1-hour AHI of 10 or more were included in the apnea group. There were 40 recordings in the apnea group with AHI ≥ 5. The subjects were 15 males and one female between 29 and 63 years of age (50 ± 6.6 years). Recordings with five minutes or less of sleep apnea were included in the normal group. There were 20 recordings in the normal group with AHI < 5. The subjects were six males and five females between 27 and 42 years of age (33 ± 5.8 years). In addition, the remaining recordings were included in the borderline group. In this study, we used the normal group and apnea group to develop and evaluate the algorithm.
To analyze the severity of OSA disease, according to the study of Li et al. [26], the apnea group was further divided into a mild-moderate OSA group and a severe OSA group through AHI. The mild-moderate OSA group included 14 recordings with 5 ≤ AHI < 30. The severe OSA group included 26 recordings with AHI ≥ 30.

2.2. Simulation Test

A simulation signal was created to compare the performance of SampEn and NPSampEn in the short-term HRV analysis. The N-point simulation signal MIX (p) was defined as follows [27]:
M I X ( p ) j = ( 1 Z j ) X j + Z j Y j   ( 1 j N ) ,
where X j = 2 sin ( 2 π j 12 ) ;   Y j was a set of real random variables within [ 3 , 3 ] , and Z j was a vector related to p, which consisted of 0 and 1, where Z j = 1 had a probability of p and Z j = 0 had a probability of 1-p. Simply, when p = 0 , Z j fully consisted of 0; when p = 1 , Z j fully consisted of 1. Hence, the value of p reflected the complexity of the M I X ( p ) signals. As the p value increased, the random components and complexity of M I X ( p ) increased. To compare the performance of SampEn and NPSampEn, we calculated the SampEn and NPSampEn of the simulation signals with different lengths in the complexity of MIX (0), MIX (0.3), MIX (0.5).

2.3. Short-Term Heart Rate Variability (HRV) Analysis

2.3.1. Preprocessing

To obtain appropriate RR segments for short-term HRV analysis, the ECG signals were preprocessed. The first six hours of overnight ECG signals were inputted. First, RR intervals were extracted from the ECG signals using the Pan-Tompkins algorithm and corrected using the local median filter proposed by Chen, L et al. [28,29]. Then, the 6-h RR intervals were divided into multiple nonoverlapping 1-min RR segments (RRs) in each recording for short-term HRV analysis. Finally, all 1-min RRs were resampled by 2 Hz with interpolation.

2.3.2. Frequency Domain Analysis

There are three classic indices in short-term HRV analysis [11,12,13]: the power in the low-frequency range (LF, 0.04–0.15 Hz), the power in the high-frequency range (HF, 0.15–0.4 Hz), and the ratio of low-frequency and high-frequency power (LF/HF). The three frequency domain indices of 1-min RRs were calculated by power spectrum density (PSD), and the mean value in 6-h RRs was computed. Due to the short lengths of the RRs, the 6-order autoregressive Burg parametric method was used to calculate the PSD.

2.3.3. Nonlinear Analysis

SampEn is a nonlinear method used to analyze the regularity of sequences by computing the similarity between template vectors. The signal vector was divided into a number of template vectors at the embedding dimensions m and m + 1, and the distance between the template vector and every other vector was computed. Finally, the tolerance parameter r (usually 0.1 to 0.25 times the standard deviation of the data [30]) was used as the threshold, and the probability of the template distance being less than the parameter ‘r’ was calculated to determine the similarity between template vectors. However, SampEn depends heavily on parameter ‘r’, and an incorrect r choice also leads to the opposite result. Hence, a SampEn measure independent of parameter ‘r’ was proposed by Udhayakumar, R.K. et al. [31].
The core process of SampEn is to compare the template distance with the threshold ‘r’ to compute the probability, and NPSampEn is a measure that maintains this core process with improvements. For NPSampEn, the artificially selected parameter ‘r’ was not used as the threshold, instead, a tolerant vector from the signal itself was generated, and then, the probability was calculated based on the tolerant vector. The schemes of SampEn and NPSampEn are shown in Figure 2, and the algorithms of NPSampEn and SampEn are as follows:
For a time series with N data points, denoted as { x ( n ) | n = 1 , 2 , , N } .
Step 1: RR reconstruction. The time series x(n) was divided into a number of template vectors at the embedding dimensions m and m + 1. When 1 i N m :
At the embedding dimension m,
X i m = { x ( i ) , , x ( i + m 1 ) }
Similarly, at the embedding dimension m + 1,
X i m + 1 = { x ( i ) , , x ( i + m ) }
Step 2: Distance calculation. At dimensions m and m + 1, the distance between the template vector and every other vector was calculated.
d i j m = m a x | x ( i + k ) x ( j + k ) |
0 k m 1 , 1 j N m   a n d   j i
Similarly,
d i j m + 1 = m a x | x ( i + k ) x ( j + k ) |
0 k m , 1 j N m   a n d   j i
Self-matches were excluded in both SampEn and NPSampEn, and two distance matrices, d m and d m + 1 , were obtained with N m dimensions without d i j m and d i j m + 1 , when i = j .
Step 3: In the probability calculation, the NPSampEn value is different from the SampEn values.
For NPSampEn, let scope be the vector containing all unique elements of d i j m and d i j m + 1 , sorted in ascending order with the length of nbin. The pr is the probability, scope(q) is the element of the scope vector ( 1 q n b i n ), and the cumulative distribution function p d i m is calculated as follows:
p d i m ( q ) = p r ( d i m s c o p e ( q ) ) , 1 q n b i n
Similarly, at the embedding dimension m + 1,
p d i m + 1 ( q ) = p r ( d i m + 1 s c o p e ( q ) ) , 1 q n b i n
The cumulative probability matrix is as follows:
p d m = [ p d 1 m ( 1 ) p d 1 m ( n b i n ) p d N m m ( 1 ) p d N m m ( n b i n ) ]
p d m + 1 = [ p d 1 m + 1 ( 1 ) p d 1 m + 1 ( n b i n ) p d N m m + 1 ( 1 ) p d N m m + 1 ( n b i n ) ]
For SampEn, the tolerance parameter r (usually 0.1 to 0.25 times the standard deviation of the data) was used as the threshold.
p d i m = p r ( d i m r ) , 1 i N m
p d i m + 1 = p r ( d i m + 1 r ) , 1 i N m
The probability vectors were obtained as follows:
p d m = [ p d 1 m , , p d N m m ]
p d m + 1 = [ p d 1 m + 1 , , p d N m m + 1 ]
Step 4: The values of NPSampEn and SampEn were calculated. We obtained the probability matrices for NPSampEn and the probability vectors for SampEn from step 3.
For NPSampEn, the vector Φ m containing the column-average of probability was obtained.
Φ m ( q ) = 1 N m i = 1 N m p d i m ( q )
Φ m + 1 ( q ) = 1 N m i = 1 N m p d i m + 1 ( q )
The NPSampEn value was finally calculated
N P S a m p E n = 1 n b i n q = 1 n b i n l n Φ m ( q ) Φ m + 1 ( q )
For SampEn, the mean of probability vector Φ m was calculated.
Φ m = 1 N m i = 1 N m p d i m
Φ m + 1 = 1 N m i = 1 N m p d i m + 1
The SampEn value was finally calculated
S a m p E n = l n Φ m Φ m + 1

2.4. Validation

To evaluate the performance of a single index, statistical one-way ANOVA, followed by post-hoc analysis with the least significant difference (LSD) test was used to analyze the statistical significance between the normal, mild-moderate OSA, and severe OSA groups. p < 0.05 was considered statistically significant, and decreasing p values represented a tendency toward statistical significance. The data are expressed as the mean ± SD. Furthermore, support vector machine (SVM) in scikit-learn support of NuSVC class in Python 3 was used for classification, and the accuracy, sensitivity, and specificity were obtained using 2-fold cross validation [32]. Moreover, the receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) of each index were computed using MATLAB. Finally, the Pearson coefficient of the association of SampEn and NPSampEn with AHI was computed using SPSS. These tests were performed using SPSS version 22.0.0.0 (SPSS Inc., Chicago, IL, USA) and MATLAB R2019a, Python 3.

3. Results

3.1. Comparisons between the Simulation of Two Entropies

To compare the performance between NPSampEn and SampEn, three simulated signals MIX (p) were constructed. Monotonicity, consistency, and continuity are three important factors that reflect the performance of entropies in complexity analysis [33]. Figure 3 indicates that entropies change with the length of signals. As shown in Figure 3d, the NPSampEn values increased with an increase in the complexity of the simulated signals, and decreased with the increase in the lengths of the simulated signals. There was no cross between the NPSampEn values for the three simulated signals. As illustrated in Figure 3e, there were crosses among the SampEn values between MIX (0.3) and MIX (0.5) when N was small (nearly 50), and there were fluctuations in the SampEn of MIX (0.3) and MIX (0.5) when N was small (N < 100). This result indicates that the value of SampEn was significantly affected by the length of the signal when the complexity of the signal was high. In Figure 3f, the increasing r value resulted in a tendency toward a lower SampEn value and represented high system regularity. The SampEn value of MIX (p) theoretically increased with the p value of the MIX (p) signal. However, the SampEn value of MIX (0.3) was sometimes higher than MIX (0.5). These results indicate that an improper r value can lead to the opposite conclusion.

3.2. HRV Analysis among the Three Groups

The recordings were divided into normal (N), mild-moderate OSA (MOSA), and severe OSA (SOSA) group. The mean ± SD values of the entropies and frequency domain indices are shown in Table 1. As shown in Table 1, in the frequency domain indices, there was no statistically significant difference in the LF or HF between any two of the three groups. The LF/HF was significantly different between the N and MOSA groups (p < 0.05) and N and SOSA groups (p < 0.001). Thus, the LF/HF can distinguish OSA patients and healthy persons, but cannot distinguish mild-moderate OSA patients and severe OSA patients. In the nonlinear indices, there were statistically significant differences in the N and SOSA groups and MOSA and SOSA groups. No significant differences were observed for SampEn between the N and MOSA groups, while NPSampEn was significantly different between the N and MOSA groups (p < 0.01). Figure 4 indicates that with the aggravation of pathogenetic conditions (from the N to MOSA to SOSA groups), the SampEn and NPSampEn values decreased, and the LF/HF value in the frequency domain increased. From a physiological point of view, the more severe the pathogenetic condition, the more complex HRV is, and the more active the sympathetic nerve is.
As shown in Figure 5, the NPSampEn values of all RRs for three recordings from the three groups were rearranged into three matrices in order. Each number in the NPSampEn matrix represents the NPSampEn value computed from the RRs time series from a given recording. For example, the first row of the NPSampEn matrix includes the first to the 18th NPSampEn values, the second row includes the 19th to the 36th NPSampEn values, and so on. These matrices were made into graphs using MATLAB (R2019a, Mathworks, Natick, MA, USA) MATLAB. The color of each rectangle reflects the NPSampEn value of one minute. The normal group is represented by warm colors, while the severe OSA group is represented by cool colors. Therefore, based on three typical cases, the color became cooler with the severity of OSA. In short, NPSampEn is believed to be a more appropriate index for OSA detection and quantization.

3.3. Relevance and Obstructive Sleep Apnea (OSA) Screening

Figure 6 indicates the correlation between AHI and NPSampEn and the correlation between AHI and SampEn. When it is significantly different (p < 0.05), the greater the absolute Pearson correlation (R) and the stronger the correlation. Negative R represents a negative correlation. The correlation between AHI and NPSampEn was strong (|r| = 0.756, p < 0.001), while the correlation between AHI and SampEn was weak (|r| = 0.600, p < 0.001).
Samples used during the OSA recording screening were obtained by calculating the mean of the segments for each recording. Fisher’s discriminant function (SPSS) was used to classify each sample as normal or apneic. The accuracy (ACC), sensitivity (SEN), specificity (SPE), and area under the receiver operating characteristic curve (AUC) values of SampEn, NPSampEn, and LF/HF were computed and are displayed in Table 2. The ACC was defined as the percentage of correctly classified samples. The SEN was defined as the percentage of correctly classified OSA samples. The SPE was defined as the percentage of correctly classified healthy samples. The AUC reflected the recall ratio of the indices. As shown in Table 2, for 60 samples, NPSampEn achieved the highest ACC (83.3%), SEN (77.5%), and SPE (95%) among these indices, which were higher than those of SampEn and LF/HF. Moreover, Figure 7 shows the receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) of SampEn, NPSampEn, and LF/HF. The closer the curve is to the upper left corner, the larger the area under the curve, and the higher the recall ratio of the index. As shown in Figure 7, NPSampEn achieved the highest AUC value (0.795), and SampEn and LF/HF yielded lower AUC values (0.581 and 0.628, respectively). Therefore, NPSampEn had the highest recall ratio.
Taking these results into consideration, it was indisputable that NPSampEn improved the OSA screening accuracy, and the classification results showed that NPSampEn is a more reliable index for ruling out OSA compared with SampEn and LF/HF.

4. Discussion

4.1. Comparison and Summary

The purpose of the study was to identify a useful single index that has interpretability and can improve the recognition performance of a classifier. Among the classic indices, the LF/HF had great OSA detection performance (Table 2), which is consistent with the results in other studies [12,13,26]. Moreover, NPSampEn yielded high screening performance and enabled the quantification of OSA severity (Figure 4, Table 2). The simulation results also proved that NPSampEn is more suitable than the other indices for short-term complexity analysis (Figure 3). Therefore, NPSampEn is more suitable for the short-term assessment of autonomic nerves and can eliminate the influence of individual differences and short-term fluctuations.
There are some previous works using the same dataset to detect OSA, although the number of features and lengths of RR segment were not completely consistent [16,26,34,35,36]. As can be seen in Table 3, some of the studies needing longer RR segments showed the performance of their introduced methods [26,34,35]. Reference notes are usually created every minute to indicate whether apnea occurred in that minute. The 1-min RR segments had better real-time performance and are more suitable for OSA detection compared with 1000-second and 5-min RR segments [37,38]. In a previous study, Al-Angari et al. analyzed the 1-min RR segment and calculated the SampEn index and frequency domain index of each minute segment. The results showed that the HRV analysis of 1-min RR segments was effective to detect OSA [16]. Some of the other studies also used 1-min RR segments to detect OSA, but these were based on multi-feature and achieved lower accuracy than that of our proposed method [16,36]. Our method not only guarantees the real-time performance, but also achieved high accuracy.

4.2. Method Motivation Analysis

In previous studies, the LF/HF has been considered as a useful linear tool for OSA screening; however, it is not suitable for analyzing nonlinear and nonstationary systems [14,15]. Entropy methods have been widely used to identify subtle changes and complexities in HRV analysis. Classic nonlinear analysis methods including SampEn and fuzzy entropy have some advantages, but there are some limitations. On one hand, reference annotations are usually made for every minute with 60–100 RR intervals to determine the presence or absence of apnea during the corresponding minute [21]. Nevertheless, SampEn and fuzzy entropy have limitations in short-term HRV analysis for two reasons. The values of SampEn and fuzzy entropy are undefined when inappropriate parameters are selected such as an inappropriate r parameter. In addition, the values of SampEn and fuzzy entropy may not accurately reflect the complexity of the signal (Figure 3f). On the other hand, the SampEn and fuzzy entropy values depend on the selection of the parameter r, which is closely related to the standard deviation of the sequence. Nevertheless, the standard deviation can only reflect the fluctuation of the sequence integrally; it cannot fully reflect the subtle fluctuations when quantizing the template distance. In addition, r values that are too similar will lead to unnecessary redundancy in the values, and increase the computational expense, while this problem occurs with fuzzy entropy due to distance quantization. Therefore, NPSampEn was proposed for use in short-term data analysis.
The results showed that NPSampEn (83.3%) significantly improved the OSA screening performance compared to SampEn (68.3%) and LF/HF (73.3%), and successfully quantified the severity of OSA (p < 0.05). Instead of the constant r value, the r solution of NPSampEn is defined based on the distances between the template vectors. The cumulative probability of the template distance determined by the set of r values can more accurately and precisely quantify the distribution of template distances in each region, and describe the signal fluctuation more specifically. The advantages of NPSampEn include its sensitivity to subtle fluctuations in signals, as the set of r values are defined by the unique template distance. As they constitute a set of adaptive values, the set of r values depends more on the dynamics of the signal and is less affected by the signal length. Therefore, NPSampEn is an adaptive method that is suitable for short-term sequences and fully reflects the fluctuations of sequences; furthermore, it can be used to overcome the influence of individual differences and identify the subtle fluctuations in physiological signals.

4.3. Physiological Significance

OSA affects the balance of the autonomic nervous system (ANS), thus affecting HRV [8]. In this study, increased LF/HF values were observed in OSA patients (Table 1, Figure 4a). This result indicates that the sympathetic nerve is more active in OSA patients. A number of mechanisms may be responsible for this finding. Relative hypoxemia may act on central chemoreceptors to increase the sympathetic tone [12].
The regulation of the autonomic nerve to the heart rate is a nonlinear and nonstationary process, and short-term sequence analysis has a considerable advantage in overcoming nonstationary interference [39]. Therefore, NPSampEn, a short-term data analysis method, is suitable for analyzing the complexity of nonstationary systems, which reflect disruptions of the ANS.
The ANS disorder caused by OSA events is related to individual differences [40]; therefore, it is not appropriate to use the classical time-frequency domain indices for OSA screening and quantification. In the calculation of NPSampEn, the adaptive parameter r is determined according to the physiological signals from different individuals, which can effectively avoid the influence of individual differences to better estimate the regularity of the physiological signals. Therefore, NPSampEn provides a new method for the analysis of HRV in the ANS of OSA patients.

4.4. Limitations

There are some limitations in the present study. First, the approach could face the problem in a more robust and reliable way and should be done on a larger set of recordings. More databases should be used to validate the effect of the method in future study. Second, it should define a training, validation, and testing framework in order to have comparable performance measures with the literature in future study.

5. Conclusions

In this study, nonparametric sample entropy (NPSampEn) was proposed as a nonlinear index to evaluate the heart rate variability and complexity of the autonomic nervous system (ANS) in OSA patients. The results showed that the ANS changes were significantly different between the normal and OSA groups. The results also showed that the complexity of ANS decreased significantly as the severity of OSA increased. Therefore, the study in this paper basically proves that NPSampEn is extremely suitable for OSA short-term detection and performs well in autonomic nerve complexity analysis.

Author Contributions

Methodology, Software, Validation, Writing—original draft, D.L. and S.W.; Writing—review & editing, L.T.; Data curation, K.F.; Funding acquisition, Project administration, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Guangdong Basic and Applied Basic Research Foundation (No. 2020A1515010701), and the Shenzhen Science and Technology Plan for fundamental research (No. JCYJ20180307153213863 and JCY20190807162003696).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ballard, R.D.; Tan, W.C.; Kelly, P.L.; Pak, J.; Pandey, R.; Martin, R.J. Effect of Sleep and Sleep Deprivation on Ventilatory Response to Bronchoconstriction. J. Appl. Physiol. 1990, 69, 490–497. [Google Scholar] [CrossRef]
  2. Malhotra, A.; White, D.P. Obstructive Sleep Apnoea. Lancet 2002, 360, 237–245. [Google Scholar] [CrossRef]
  3. Caples, S.M.; Gami, A.S.; Somers, V.K. Obstructive Sleep Apnea. Ann. Intern. Med. 2005, 142, 187–197. [Google Scholar] [CrossRef] [PubMed]
  4. Khoo, M.C.K.; Kim, T.S.; Berry, R.B. Spectral Indices of Cardiac Autonomic Function in Obstructive Sleep Apnea. Sleep 1999, 22, 443–451. [Google Scholar] [CrossRef] [PubMed]
  5. Cartwright, R. Obstructive Sleep Apnea: A Sleep Disorder with Major Effects on Health. Dm Dis. Mon. 2001, 47, 109–147. [Google Scholar] [CrossRef] [PubMed]
  6. Kapur, V.; Strohl, K.P.; Redline, S.; Iber, C.; O′Connor, G.; Nieto, J. Underdiagnosis of Sleep Apnea Syndrome in U.S. Communities. Sleep Breath. Schlaf Atm. 2002, 6, 49–54. [Google Scholar] [CrossRef]
  7. Young, T.; Evans, L.; Finn, L.; Palta, M. Estimation of the Clinically Diagnosed Proportion of Sleep Apnea Syndrome in Middle-Aged Men and Women. Sleep 1997, 20, 705–706. [Google Scholar] [CrossRef]
  8. Babaeizadeh, S.; White, D.P.; Pittman, S.D.; Zhou, S.H. Automatic Detection and Quantification of Sleep Apnea Using Heart Rate Variability. J. Electrocardiol. 2010, 43, 535–541. [Google Scholar] [CrossRef]
  9. Griffiths, K.R.; Quintana, D.S.; Hermens, D.F.; Spooner, C.; Tsang, T.W.; Clarke, S.; Kohn, M.R. Sustained Attention and Heart Rate Variability in Children and Adolescents with Adhd. Biol. Psychol. 2017, 124, 11–20. [Google Scholar] [CrossRef]
  10. Song, M.-K.; Ha, J.H.; Ryu, S.-H.; Yu, J.; Park, D.-H. The Effect of Aging and Severity of Sleep Apnea on Heart Rate Variability Indices in Obstructive Sleep Apnea Syndrome. Psychiatry Investig. 2012, 9, 65–72. [Google Scholar] [CrossRef] [Green Version]
  11. Tarvainen, M.P.; Niskanen, J.-P.; Lipponen, J.A.; Ranta-aho, P.O.; Karjalainen, P.A. Kubios Hrv-Heart Rate Variability Analysis Software. Comput. Methods Programs Biomed. 2014, 113, 210–220. [Google Scholar] [CrossRef] [PubMed]
  12. Gula, L.J.; Krahn, A.D.; Skanes, A.; Ferguson, K.A.; George, C.; Yee, R.; Klein, G.J. Heart Rate Variability in Obstructive Sleep Apnea: A Prospective Study and Frequency Domain Analysis. Ann. Noninvasive Electrocardiol. 2003, 8, 144–149. [Google Scholar] [CrossRef]
  13. Gong, X.; Huang, L.; Liu, X.; Li, C.; Mao, X.; Liu, W.; Huang, X.; Chu, H.; Wang, Y.; Wu, W.; et al. Correlation Analysis between Polysomnography Diagnostic Indices and Heart Rate Variability Parameters among Patients with Obstructive Sleep Apnea Hypopnea Syndrome. PLoS ONE 2016, 11, e0156628. [Google Scholar] [CrossRef] [Green Version]
  14. Makikallio, T.H.; Tapanainen, J.M.; Tulppo, M.P.; Huikuri, H.V. Clinical Applicability of Heart Rate Variability Analysis by Methods Based on Nonlinear Dynamics. Card. Electrophysiol. Rev. 2002, 6, 250–255. [Google Scholar] [CrossRef] [PubMed]
  15. Jons, C.; Raatikainen, P.; Gang, U.J.; Huikuri, H.V.; Joergensen, R.M.; Johannesen, A.; Dixen, U.; Messier, M.; McNitt, S.; Thomsen, P.E.B.; et al. Autonomic Dysfunction and New-Onset Atrial Fibrillation in Patients with Left Ventricular Systolic Dysfunction after Acute Myocardial Infarction: A Carisma Substudy. J. Cardiovasc. Electrophysiol. 2010, 21, 983–990. [Google Scholar] [CrossRef]
  16. Al-Angari, H.M.; Sahakian, A.V. Use of Sample Entropy Approach to Study Heart Rate Variability in Obstructive Sleep Apnea Syndrome. IEEE Trans. Biomed. Eng. 2007, 54, 1900–2004. [Google Scholar] [CrossRef]
  17. Ravelo-Garcia, A.G.; Navarro-Mesa, J.L.; Casanova-Blancas, U.; Martin-Gonzalez, S.; Quintana-Morales, P.; Guerra-Moreno, I.; Canino-Rodriguez, J.M.; Hernandez-Perez, E. Application of the Permutation Entropy over the Heart Rate Variability for the Improvement of Electrocardiogram-Based Sleep Breathing Pause Detection. Entropy 2015, 17, 914–927. [Google Scholar] [CrossRef]
  18. Li, Y.; Pan, W.; Li, K.; Jiang, Q.; Liu, G. Sliding Trend Fuzzy Approximate Entropy as a Novel Descriptor of Heart Rate Variability in Obstructive Sleep Apnea. IEEE J. Biomed. Health Inform. 2019, 23, 175–183. [Google Scholar] [CrossRef] [PubMed]
  19. Voss, A.; Heitmann, A.; Schroeder, R.; Peters, A.; Perz, S. Short-Term Heart Rate Variability-Age Dependence in Healthy Subjects. Physiol. Meas. 2012, 33, 1289–1311. [Google Scholar] [CrossRef] [PubMed]
  20. Khandoker, A.H.; Karmakar, C.K.; Palaniswami, M. Comparison of Pulse Rate Variability with Heart Rate Variability during Obstructive Sleep Apnea. Med. Eng. Phys. 2011, 33, 204–209. [Google Scholar] [CrossRef] [PubMed]
  21. Penzel, T.; Moody, G.B.; Mark, R.G.; Goldberger, A.L.; Peter, J.H.; Ieee, I. The Apnea-Ecg Database. Comput. Cardiol. 2000, 27, 255–258. [Google Scholar]
  22. Li, P.; Liu, C.; Li, K.; Zheng, D.; Liu, C.; Hou, Y. Assessing the Complexity of Short-Term Heartbeat Interval Series by Distribution Entropy. Med Biol. Eng. Comput. 2015, 53, 77–87. [Google Scholar] [CrossRef]
  23. Mayer, C.C.; Bachler, M.; Hörtenhuber, M.; Stocker, C.; Holzinger, A.; Wassertheurer, S. Selection of entropy-measure parameters for knowledge discovery in heart rate variability data. BMC Bioinform. 2014, 15, S2. [Google Scholar] [CrossRef] [Green Version]
  24. Castiglioni, P.; Di Rienzo, M. How the threshold r influences approximate entropy analysis of heart-rate variability. Comput. Cardiol. 2008, 561–564. [Google Scholar] [CrossRef]
  25. Shahar, E.; Whitney, C.W.; Redline, S.; Lee, E.T.; Newman, A.B.; Nieto, F.J.; O′Connor, G.T.; Boland, L.L.; Schwartz, J.E.; Samet, J.M. Sleep-Disordered Breathing and Cardiovascular Disease: Cross-Sectional Results of the Sleep Heart Health Study. Am. J. Respir. Crit. Care Med. 2001, 163, 19–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Li, Y.F.; Wu, S.; Yang, Q.A.N.; Liu, G.Z.; Ge, L.J. 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]
  27. Pincus, S.M. Approximate Entropy as a Measure of System Complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Penzel, T.; Kantelhardt, J.W.; Grote, L.; Peter, J.H.; Bunde, A. Comparison of Detrended Fluctuation Analysis and Spectral Analysis for Heart Rate Variability in Sleep and Sleep Apnea. IEEE Trans. Biomed. Eng. 2003, 50, 1143–1151. [Google Scholar] [CrossRef] [Green Version]
  29. Chen, L.; Zhang, X.; Song, C. An Automatic Screening Approach for Obstructive Sleep Apnea Diagnosis Based on Single-Lead Electrocardiogram. IEEE Trans. Autom. Sci. Eng. 2015, 12, 106–115. [Google Scholar] [CrossRef]
  30. Lake, D.E.; Richman, J.S.; Griffin, M.P.; Moorman, J.R. Sample Entropy Analysis of Neonatal Heart Rate Variability. Am. J. Physiol. -Regul. Integr. Comp. Physiol. 2002, 283, R789–R797. [Google Scholar] [CrossRef] [Green Version]
  31. Udhayakumar, R.K.; Karmakar, C.; Palaniswami, M. Understanding Irregularity Characteristics of Short-Term Hrv Signals Using Sample Entropy Profile. IEEE Trans. Biomed. Eng. 2018, 65, 2569–2579. [Google Scholar] [CrossRef] [PubMed]
  32. Hamid, N.A.; Sjarif, N.N.A. Handwritten recognition using SVM, KNN and neural network. arXiv 2017, arXiv:1702.00723. [Google Scholar]
  33. Ao, D.; Sun, R.; Tong, K.-Y.; Song, R. Characterization of Stroke- and Aging-Related Changes in the Complexity of Emg Signals During Tracking Tasks. Ann. Biomed. Eng. 2015, 43, 990–1002. [Google Scholar] [CrossRef] [PubMed]
  34. Pietrzak, A.; Cybulski, G. New Indices for Sleep Apnea Detection from Long-Time Ecg Recordings. In Proceedings of the 2015 Computing in Cardiology Conference, Nice, France, 6–9 September 2015. [Google Scholar]
  35. Ravelo-García, A.G.; Casanova-Blancas, U.; Martín-González, S.; Hernández-Pérez, E.; Morales, P.Q.; Navarro-Mesa, J.L. An Approach to the Improvement of Electrocardiogram-Based Sleep Breathing Pauses Detection by Means of Permutation Entropy of the Heart Rate Variability. In Proceedings of the 3rd IEEE International Work-Conference on Bioinspired Intelligence, Liberia, Costa Rica, 16–18 July 2014. [Google Scholar]
  36. Varon, C.; Caicedo, A.; Testelmans, D.; Buyse, B.; Van Huffel, S. A Novel Algorithm for the Automatic Detection of Sleep Apnea from Single-Lead Ecg. IEEE Trans. Biomed. Eng. 2015, 62, 2269–2278. [Google Scholar] [CrossRef] [PubMed]
  37. Li, K.; Pan, W.; Li, Y.; Jiang, Q.; Liu, G. A method to detect sleep apnea based on deep neural network and hidden Markov model using single-lead ECG signal. Neurocomputing 2018, 294, 94–101. [Google Scholar] [CrossRef]
  38. Feng, K.; Qin, H.; Wu, S.; Pan, W.; Liu, G. A sleep apnea detection method based on unsupervised feature learning and single-lead electrocardiogram. IEEE Trans. Instrum. Meas. 2020, 70. [Google Scholar] [CrossRef]
  39. 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 2021, 65, 102370. [Google Scholar] [CrossRef]
  40. Liu, G.; Wang, Q.; Chen, S.; Zhou, G.; Chen, W.; Wu, Y. Robustness Evaluation of Heart Rate Variability Measures for Age Gender Related Autonomic Changes in Healthy Volunteers. Australas. Phys. Eng. Sci. Med. 2014, 37, 567–574. [Google Scholar] [CrossRef]
Figure 1. Framework of the proposed method of analyzing heart rate variability (HRV).
Figure 1. Framework of the proposed method of analyzing heart rate variability (HRV).
Entropy 23 00267 g001
Figure 2. Nonlinear analysis method of NPSampEn and SampEn.
Figure 2. Nonlinear analysis method of NPSampEn and SampEn.
Entropy 23 00267 g002
Figure 3. Fifty points simulated for the MIX (0), MIX (0.3), and MIX (0.5) signals (ac). N corresponds to the data length, N increases from 50 to 200 by intervals of 2. The changes in NPSampEn (d), and SampEn (e) (r = 0.25). The change in SampEn (f) with an increase in r from 0.1 to 1 by 0.02 (N = 100) and the corresponding changes in the complexity of MIX (0), MIX (0.3), MIX (0.5).
Figure 3. Fifty points simulated for the MIX (0), MIX (0.3), and MIX (0.5) signals (ac). N corresponds to the data length, N increases from 50 to 200 by intervals of 2. The changes in NPSampEn (d), and SampEn (e) (r = 0.25). The change in SampEn (f) with an increase in r from 0.1 to 1 by 0.02 (N = 100) and the corresponding changes in the complexity of MIX (0), MIX (0.3), MIX (0.5).
Entropy 23 00267 g003
Figure 4. Comparison of the ratio of low frequency power to high frequency power (LF/HF) (a), SampEn (b), and NPSampEn (c) for the normal, mild-moderate OSA, and severe OSA groups. *, **, *** represent p < 0.05, p < 0.01, and p < 0.001, respectively
Figure 4. Comparison of the ratio of low frequency power to high frequency power (LF/HF) (a), SampEn (b), and NPSampEn (c) for the normal, mild-moderate OSA, and severe OSA groups. *, **, *** represent p < 0.05, p < 0.01, and p < 0.001, respectively
Entropy 23 00267 g004
Figure 5. The NPSampEn values for three recordings. (a) NPSampEn values of a recording from the N group, (b) NPSampEn values of a recording from the mild-moderate OSA (MOSA) group, and (c) NPSampEn values of a recording from the severe OSA (SOSA) group.
Figure 5. The NPSampEn values for three recordings. (a) NPSampEn values of a recording from the N group, (b) NPSampEn values of a recording from the mild-moderate OSA (MOSA) group, and (c) NPSampEn values of a recording from the severe OSA (SOSA) group.
Entropy 23 00267 g005
Figure 6. The relevance of NPSampEn (a) and SampEn (b) with apnea-hypopnea index (AHI). r: Pearson’s correlation coefficient; p: p-value.0.5 < |r| < 0.8 represent significant correlation.
Figure 6. The relevance of NPSampEn (a) and SampEn (b) with apnea-hypopnea index (AHI). r: Pearson’s correlation coefficient; p: p-value.0.5 < |r| < 0.8 represent significant correlation.
Entropy 23 00267 g006
Figure 7. The receiver operating characteristic (ROC) curves of LF/HF, NPSampEn, and SampEn.
Figure 7. The receiver operating characteristic (ROC) curves of LF/HF, NPSampEn, and SampEn.
Entropy 23 00267 g007
Table 1. Entropy/frequency domain indices for the normal, mild-moderate, and severe obstructive sleep apnea (OSA) groups.
Table 1. Entropy/frequency domain indices for the normal, mild-moderate, and severe obstructive sleep apnea (OSA) groups.
IndicesN
(Mean ± SD)
MOSA
(Mean ± SD)
SOSA
(Mean ± SD)
p Value
N & MOSAN & SOSAMOSA & SOSA
EntropySampEn0.587 ± 0.0580.564 ± 0.0940.443 ± 0.0970.4260 ***0 ***
NPSampEn0.270 ± 0.0140.243 ± 0.0390.199 ± 0.0290.006 **0 ***0 ***
Frequency domainLF0.0060 ± 0.00750.0010 ± 0.00070.022 ± 0.0540.6870.1500.088
HF0.0050 ± 0.00780.00047 ± 0.000420.0074 ± 0.0220.4040.6080.184
LF/HF2.049 ± 0.7173.785 ± 2.0004.975 ± 2.4370.012 *0 ***0.067
*, **, *** represent p < 0.05, p < 0.01, and p < 0.001, respectively
Table 2. Performance comparisons of different indices for obstructive sleep apnea (OSA) screening.
Table 2. Performance comparisons of different indices for obstructive sleep apnea (OSA) screening.
IndicesTPTNFPFNACCSENSPEAUC
Nonlinear
indices
SampEn291281168.3%72.5%60.0%0.581
NPSampEn31191983.3%77.5%95.0%0.795
Frequency domainLF/HF291551173.3%72.5%75%0.628
NPSampEn: non-parametric sample entropy; LF/HF: ratio of low frequency power to high frequency power. TP: true positive. TN: true negative; FP: false positive; FN: false negative; ACC: accuracy; SEN: sensitivity; SPE: specificity; AUC: area under the receiver operating character curve.
Table 3. Comparison of the classification results between our proposed method and previous works on the Physionet database.
Table 3. Comparison of the classification results between our proposed method and previous works on the Physionet database.
ReferenceMethodFeaturesLength of RR SegmentResults
Pietrzak et al. [34]Standard deviation of successive differencesingle feature1000 sACC = 88.5%
SEN = 96.0%
SPE = 70.0%
Ravelo-García et al. [35]Permutation entropymulti-feature5 minACC = 78.0%
SEN = 64.3%
SPE = 86.5%
Li et al. [26]Variance delay fuzzy approximate entropysingle feature5 minACC = 90.0%
SEN = 87.5%
SPE = 95.0%
Al-Angari et al. [16]Sample entropy m = 1, 2, 3multi-feature1 minACC = 70.3%
SEN = 69.5%
SPE = 70.8%
Varon et al. [36]support vector machinemulti-feature1 minACC = 84.7%
SEN = 84.7%
SPE = 84.7%
Proposed methodNonparametric sample entropysingle feature1 minACC = 83.3%
SEN = 77.5%
SPE = 95.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

Liang, D.; Wu, S.; Tang, L.; Feng, K.; Liu, G. Short-Term HRV Analysis Using Nonparametric Sample Entropy for Obstructive Sleep Apnea. Entropy 2021, 23, 267. https://0-doi-org.brum.beds.ac.uk/10.3390/e23030267

AMA Style

Liang D, Wu S, Tang L, Feng K, Liu G. Short-Term HRV Analysis Using Nonparametric Sample Entropy for Obstructive Sleep Apnea. Entropy. 2021; 23(3):267. https://0-doi-org.brum.beds.ac.uk/10.3390/e23030267

Chicago/Turabian Style

Liang, Duan, Shan Wu, Lan Tang, Kaicheng Feng, and Guanzheng Liu. 2021. "Short-Term HRV Analysis Using Nonparametric Sample Entropy for Obstructive Sleep Apnea" Entropy 23, no. 3: 267. https://0-doi-org.brum.beds.ac.uk/10.3390/e23030267

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