In this paper the permutation entropy (PE) obtained from heart rate variability (HRV) is analyzed in a statistical model. In this model we also integrate other feature extraction techniques, the cepstrum coefficients derived from the same HRV and a set of band powers obtained from the electrocardiogram derived respiratory (EDR) signal. The aim of the model is detecting obstructive sleep apnea (OSA) events. For this purpose, we apply two statistical classification methods: Logistic Regression (LR) and Quadratic Discriminant Analysis (QDA). For testing the models we use seventy ECG recordings from the Physionet database which are divided into equal-size learning and testing sets. Both sets consist of 35 recordings, each containing a single ECG signal. In our experiments we have found that the features extracted from the EDR signal present a sensitivity of 65.6% and specificity of 87.7% (auc = 85) in the LR classifier, and sensitivity of 59.4% and specificity of 90.3% (auc = 83.9) in the QDA classifier. The HRV-based cepstrum coefficients present a sensitivity of 63.8% and specificity of 89.2% (auc = 86) in the LR classifier, and sensitivity of 67.2% and specificity of 86.8% (auc = 86.9) in the QDA. Subsequent tests show that the contribution of the permutation entropy increases the performance of the classifiers, implying that the complexity of RR interval time series play an important role in the breathing pauses detection. Particularly, when all features are jointly used, the quantification task reaches a sensitivity of 71.9% and specificity of 92.1% (auc = 90.3) for LR. Similarly, for QDA the sensitivity is 75.1% and the specificity is 90.5% (auc = 91.7).

Our contribution focuses on the characterization of sleep apnea from a cardiac rate point of view, using Recurrence Quantification Analysis (RQA), based on a Heart Rate Variability (HRV) feature selection process. Three parameters are crucial in RQA: those related to the embedding process (dimension and delay) and the threshold distance. There are no overall accepted parameters for the study of HRV using RQA in sleep apnea. We focus on finding an overall acceptable combination, sweeping a range of values for each of them simultaneously. Together with the commonly used RQA measures, we include features related to recurrence times, and features originating in the complex network theory. To the best of our knowledge, no author has used them all for sleep apnea previously. The best performing feature subset is entered into a Linear Discriminant classifier. The best results in the “Apnea-ECG Physionet database” and the “HuGCDN2014 database” are, according to the area under the receiver operating characteristic curve, 0.93 (Accuracy: 86.33%) and 0.86 (Accuracy: 84.18%), respectively. Our system outperforms, using a relatively small set of features, previously existing studies in the context of sleep apnea. We conclude that working with dimensions around 7–8 and delays about 4–5, and using for the threshold distance the Fixed Amount of Nearest Neighbours (FAN) method with 5% of neighbours, yield the best results. Therefore, we would recommend these reference values for future work when applying RQA to the analysis of HRV in sleep apnea. We also conclude that, together with the commonly used vertical and diagonal RQA measures, there are newly used features that contribute valuable information for apnea minutes discrimination. Therefore, they are especially interesting for characterization purposes. Using two different databases supports that the conclusions reached are potentially generalizable, and are not limited by database variability.

In this paper, we thoroughly analyze the detection of sleep apnea events in the context of Obstructive Sleep Apnea (OSA), which is considered a public health problem because of its high prevalence and serious health implications. We especially evaluate patients who do not always show desaturations during apneic episodes (non-desaturating patients). For this purpose, we use a database (HuGCDN2014-OXI) that includes desaturating and non-desaturating patients, and we use the widely used Physionet Apnea Dataset for a meaningful comparison with prior work. Our system combines features extracted from the Heart-Rate Variability (HRV) and SpO₂, and it explores their potential to characterize desaturating and non-desaturating events. The HRV-based features include spectral, cepstral, and nonlinear information (Detrended Fluctuation Analysis (DFA) and Recurrence Quantification Analysis (RQA)). SpO₂-based features include temporal (variance) and spectral information. The features feed a Linear Discriminant Analysis (LDA) classifier. The goal is to evaluate the effect of using these features either individually or in combination, especially in non-desaturating patients. The main results for the detection of apneic events are: (a) Physionet success rate of 96.19%, sensitivity of 95.74% and specificity of 95.25% (Area Under Curve (AUC): 0.99); (b) HuGCDN2014-OXI of 87.32%, 83.81% and 88.55% (AUC: 0.934), respectively. The best results for the global diagnosis of OSA patients (HuGCDN2014-OXI) are: success rate of 95.74%, sensitivity of 100%, and specificity of 89.47%. We conclude that combining both features is the most accurate option, especially when there are non-desaturating patterns among the recordings under study.