Next Article in Journal
Correlation Analysis of Noise, Vibration, and Harshness in a Vehicle Using Driving Data Based on Big Data Analysis Technique
Next Article in Special Issue
DFTrans: Dual Frequency Temporal Attention Mechanism-Based Transportation Mode Detection
Previous Article in Journal
An Ensemble Learning Model for COVID-19 Detection from Blood Test Samples
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Automatic Detection of the Cyclic Alternating Pattern of Sleep and Diagnosis of Sleep-Related Pathologies Based on Cardiopulmonary Resonance Indices

1
University of Chinese Academy of Sciences, Beijing 101408, China
2
CAS Institute of Healthcare Technologies, Nanjing 210046, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(6), 2225; https://doi.org/10.3390/s22062225
Submission received: 16 February 2022 / Revised: 6 March 2022 / Accepted: 7 March 2022 / Published: 14 March 2022
(This article belongs to the Special Issue Human-Centric Sensing Technology and Systems)

Abstract

:
The cyclic alternating pattern is the periodic electroencephalogram activity occurring during non-rapid eye movement sleep. It is a marker of sleep instability and is correlated with several sleep-related pathologies. Considering the connection between the human heart and brain, our study explores the feasibility of using cardiopulmonary features to automatically detect the cyclic alternating pattern of sleep and hence diagnose sleep-related pathologies. By statistically analyzing and comparing the cardiopulmonary characteristics of a healthy group and groups with sleep-related diseases, an automatic recognition scheme of the cyclic alternating pattern is proposed based on the cardiopulmonary resonance indices. Using the Hidden Markov and Random Forest, the scheme combines the variation and stability of measurements of the coupling state of the cardiopulmonary system during sleep. In this research, the F1 score of the sleep-wake classification reaches 92.0%. In terms of the cyclic alternating pattern, the average recognition rate of A-phase reaches 84.7% on the CAP Sleep Database of 108 cases of people. The F1 score of disease diagnosis is 87.8% for insomnia and 90.0% for narcolepsy.

1. Introduction

Sleep, which accounts for nearly a third of human life, is an important function that helps the body to recover. It has been proven that sleep could help to reduce stress, regulate hormone balance, stabilize appetite and cardiovascular function [1,2,3]. At the same time, sleep is essential for the recovery of the brain function, which is closely related to brain development, learning, memory and mental health of human beings [4]. A lack of sleep will cause different degrees of harm to the body and mind [5]. The monitoring of sleep and the detection of sleep-related diseases are of great significance in people’s daily life as well as in clinical treatment.
Sleep is known to rhythmically regulate autonomic nervous system activity. Quiet sleep is associated with increased parasympathetic arousal and activity, while rapid eye movement (REM) sleep is relevant to the increased sympathetic activity [6]. The sleep structure is based on the cyclical alternation of two main neurophysiological states: REM and NON-REM (NREM) sleep [7]. The alternations of non-REM and REM sleep constitute the sleep cycle, and its recurrence during the night determines the classical progressive sleep mode.
Many researchers suggest that both slow and fast electroencephalogram (EEG) activation complexes are involved in the structural organization of sleep [8]. Deep sleep is established and maintained through a process of periodic EEG instability accompanied by mild autonomic fluctuations in the wake state [9]. In contrast, the breakdown of slow-wave sleep and the introduction of REM sleep are mainly associated with EEG desynchronization and strong activation of muscle and autonomic nerve functions [10].
The Cyclic Alternating Pattern (CAP) is the periodic EEG activity occurring during NREM sleep. It is characterized by cyclic sequences of cerebral activation (phase A) followed by periods of deactivation (phase B) which separate two successive phase A periods with an interval <1 min [11]. Phase A period together with the following phase B period define a CAP cycle [12]. Detailed investigation has ascertained that the spontaneous EEG fluctuations of CAP are implicated in the subtle mechanisms that regulate the production and attenuation of slow-wave activities during sleep [13]. Different components of CAP have a sculpturing effect on the profile of the sleep cycle.
Phase A periods are subdivided into three subtypes. The abundance of A1 subtype in the descending branches and grooves may be an EEG expression of brain mechanisms involved in release activity, while the dominance of A2 and A3 subtypes in the ascending branches reflects the REM-on drive [14]. Therefore, in addition to a variety of EEG features, the activation complex also has non-random distribution characteristics at night, and has obvious periodicity in NREM sleep within the CAP framework, e.g., [15,16,17]. CAP is regarded as the main expression of sleep microstructure. CAP can be recognized in the sleep of both adults and children and it is a sensitive tool for studying sleep disorders throughout the life cycle [18,19].
Several efforts have been made to develop a reliable automatic CAP-scoring algorithm [20,21,22,23]. Most of these methods rely on the extraction of spectral features from the EEG and on the application of machine-learning algorithms, such as the k-nearest neighbor, support vector machine, artificial neural network, decision trees, and deep neural network [24,25]. However, EEG acquisition needs to be carried out under the guidance of experts, and the wearing of the equipment will also affect the sleep state of the subjects [26]. The preprocessing [27], recognition [28] and analysis [29] of EEG signal are complicated [30].
Some researchers and institutions have tried to classify sleep stages by physiological signals instead of EEG. As physiological signals that contain the physiological characteristics and autonomic nerve status information of the human body, Electrocardiograph (ECG), respiration and the three-axis acceleration signals on the chest have attracted much attention [31,32,33]. Wilhelm Daniel Scherz et al. used Heart Rate Variability (HRV) characteristics of three transformation domains to determine sleep stages [34]. Mourad Adnane et al. attempted to extract heart rate variability features from time domain, frequency domain, detrend fluctuation analysis and window detrend fluctuation analysis, and used a support vector machine to classify sleep and wake stage, with an average accuracy of 79.3% [35]. Martin Oswaldo Mendez et al. extracted HRV and body movement using a time-varying auto-regression model, and classified REM and NREM using hidden Markov model [36], and the time-varying relationship between EEG and HRV are studied [37]. Eline R.de Groot et al. researched on the value of cardiorespiratory parameters for sleep state classification in preterm infants [38].
Considering the connection between the heart and brain of people [39,40,41], the CAP, which derives from the transformation of EEG, affects the excitability of human nerves and then reflects not only in body movements and changes of heart rate, but also in respiration [42,43]. Cardiopulmonary characteristics of people show diverse manifestations in different stages of the sleep. Inspired by the causality analysis of the modulation of heart rate form respiration, cardiopulmonary resonance indices (CRI) are defined to measure the status of the cardiopulmonary resonance system by adopting the cardiopulmonary system circuit model based on the causality analysis in frequency domain [42].
In this research, we first expanded the index system of CRI, combined with the physiological significance of the cardiopulmonary coupling model. CRI in different sleep stages of different people were analyzed in detail in this paper. Inspired by the results of the statistical analysis of the cardiopulmonary characteristics of healthy people and patients with sleep-related diseases, we proposed an automatic recognition and disease diagnostic scheme of CAP based on CRI. The scheme combines the variation and stability of CRI using Hidden Markov [44] and Random Forest model [45]. In addition to the recognition of CAP, this scheme can also diagnose sleep-related diseases through the cardiopulmonary characteristics. Compared with the features of the directly acquired signals, such as ECG, three-axis acceleration signals, etc., the proposed method explores the deep rhythm in the cardiopulmonary system during human sleep, and is effective to conquer the weaknesses in signal acquisition and processing in clinical trials. At the same time, this study also proves the existence of the mind-brain connection in the human body.
In this paper, for materials and methods, Section 2 describes the database used in the study and the data preprocessing before calculating CRI. Section 3 describes the CRI and expands the index system, and statistical analysis of the cardiopulmonary characteristics of people with non-pathology and patients with sleep-related diseases, then gives the classification and recognition scheme. Section 3 and Section 4 give the details and effects of the methods. The last section is the conclusion and future work.

2. Materials and Methods

2.1. Data

This section describes the data set and the data processing in this research.

2.1.1. Data Set

The experimental data used in this study is obtained from the CAP Sleep Database [46] in MIT-BIH database. The CAP Sleep Database is a collection of 108 polysomnographic recordings registered at the Sleep Disorders Center of the Ospedale Maggiore of Parma, Italy. The waveforms include at least three EEG channels, Electro-Oculogram (EOG) (two channels), Electromyography (EMG) of the submentalis muscle, bilateral anterior tibial EMG, respiration signals, the acceleration and gyroscope signals in the chest and ECG [47].
The database includes 108 cases of people with non-pathology and several sleep-related pathologies. Expert neurologists who trained at the Sleep Center provided the scoring of the sleep macrostructure, according to the Rechtschaffen & Kales rules [48], including SLEEP−REM, SLEEP−S0, SLEEP−S1, SLEEP−S2, SLEEP−S3, SLEEP−S4, while the CAP was detected in agreement with Terzano’s reference atlas of rules [49], including MCAP−A1, MCAP−A2 and MCAP−A3 in the NREM period every 30 s. Sleep stages S0, S1 and S2 are collectively called light sleep stage, while S3 and S4 are referred to as deep sleep stage.
The Figure 1 below is an example of the EEG in sleep stage 2. Phases A1, A2 and A3 are framed in the diagram according to the labels from experts. As shown in Figure 1, the Phase A periods are subdivided into three subtypes [50]: Subtype A1: synchronized events with low impact on autonomic and somatomotor activities; Subtype A2: mixed synchronized–desynchronized EEG events with an intermediate influence on the autonomic and somatomotor activities; Subtype A3: predominantly desynchronized EEG events with heavy effects on the autonomic and somatomotor activities.

2.1.2. Data Selection and Preprocessing

The ECG and respiration signals of people with non-pathology, insomnia and narcolepsy are used in this study. For the signals in the database, stationary signals are selected for analysis.
Wavelet analysis is used to remove the baseline drift of the signal and a Pasteur band-pass filter is used to remove the noise. Then, for respiratory signals, we define the signal quality function as follows:
Q = 0.03 0.5 p s d ( r e s ) p s d ( r e s )
where res represents the respiratory signal, psd represents power spectral density function. Q value expresses the concentration degree of the power spectrum of respiratory signal, and the stability of respiratory frequency. Considering the respiratory signal is concentrated in the frequency domain from 0.03 Hz to 0.5 Hz [51], the signals which meet the requirement that Q > 0.85 are chosen [43,52]. This process ensures the quality of the respiratory signal to be representative of the breathing condition of the subject.
As for the ECG, to ensure RR intervals (intervals between adjacent R waves) in the normal range and without a mutation because of premature beat or other physiological phenomena [53], we use interpolation as a substitute for points that do not meet the following conditions based on the Gaussian model [54]:
| R R I i R R I ¯ | < 1.5 * S t d ( R R I ) 0.7 * R R I i 1 < R R I i < 1.3 * R R I i 1
where RRI is the RR intervals, RRIi−1 and RRIi are adjacent intervals.
After processing ECG and respiratory signals, we calculated CRI sequences taking the step size of 10 s with the window of 120 s.

2.2. Methods

2.2.1. Cardiopulmonary Resonance Indices (CRI)

Cardiopulmonary interaction is important in the circulation system to ensure efficient delivery of oxygen and nutrients, and that the efficiency is optimized at the state of cardiopulmonary resonance. Cardiopulmonary resonance indices (CRI) come from the bivariate autoregressive model of respiration series and RR intervals. The model could calculate respiratory along with the non-respiratory component effects on RR intervals in the frequency domain based on G-causality [42].
In the model, the change process of RRI was regarded as a Markov process, ignoring other factors affecting heart rate in the short term, the RRI and RSP (both of length T) were described by a bivariate auto-regressive model [43]:
R R I ( t ) = j = 1 p A 11 , j R R I ( t j ) + j = 1 p A 12 , j R S P ( t j ) + ε 1 ( t ) R S P ( t ) = j = 1 p A 21 , j R R I ( t j ) + j = 1 p A 22 , j R S P ( t j ) + ε 2 ( t )
where p is the maximum number of lagged observations included in the model (the model order, p < T). A contains the coefficients of the model, and ε1, ε2 are the residuals for each time series.
By the definition of G-causality, the magnitude of RSA can be measured by the log ratio of the prediction error variances for the restricted (R) and unrestricted (U) models:
G 2 1 = ln var ( ε 1 R ( 12 ) ) var ( ε 1 U )
where ε1R(12) is derived from the model omitting the A12,j (for all j) coefficients in the first equation and ε1U is derived from the full model [21]. In order to describe the effect of respiration on heart rate better, and to compare the model with Heart Rate Variability (HRV), the model conducted the operation above in the frequency domain. Thus, through the model, we obtained a curve G(f) which measures the cardiopulmonary interaction in the frequency domain.
Based on G(f) in the frequency domain, CRI are defined to measure the status of the cardiopulmonary resonance system by adopting the cardiopulmonary system circuit model [42,55,56] (Figure 2).
CRI include cardiopulmonary resonance amplitude (CRA), cardiopulmonary resonance bandwidth (CRB) and cardiopulmonary resonance quality factor (CRQ) shown in Figure 3. CRA represents the coupling depth of the cardiopulmonary system, CRB represents the width of the coupling frequency band, and CRQ represents the current quality of the cardiopulmonary system. CRI are the quantitative measurements for respiratory sinus arrhythmia (RSA) in the frequency domain [42]. Reflecting the modulation effect of breathing on heart rate changes, CRI well represent the degree of cardiopulmonary resonance and parasympathetic nerve activity level. They have been proved valid in the classification of NON-REM (NREM) and REM sleep [42] and in the antepartum autonomic nervous care in pregnant women [43].
CRQ is defined to measure the merit of the cardiopulmonary resonance system by adopting the quality factor measure for inductor, capacitor, and resistor LCR oscillator. Combined with the physiological significance of the cardiopulmonary coupling model, we expanded the index system with non-respiration factors, which is an analogy with the energy consumption element in the circuit model, as cardiopulmonary resonance resistance (CRR).
In the oscillating system, the quality factor of a circuit is the ratio between the reactance power and the average power consumed by the resistor, namely the ratio of reactive power and active power. In the non-radiative system, the capacitive reactance is equal to the inductive reactance at resonance, then we can deduce [57]:
CRQ = 1 2 π F A C R
where FA represents the cardiopulmonary resonance frequency. Then we define the cardiopulmonary resonance resistance (CRR) as follows [57]:
CRR = 2 π F A C R
Cardiopulmonary resonance resistance (CRR) represents the physiological phenomena of people that exist in the heart rate variability excluding the respiratory modulation, and is related to pathological state of the subject.
CRA, CRB, CRQ and CRR jointly constructed the evaluation index system of cardiopulmonary coupling. In the following analysis, we constructed the sleep analysis scheme by using the performance of the four indicators (including CRQ, CRR, CRB and CRA) in different groups. The significant difference of the CRI in A-phase and non-A phase (NA-phase) period were detected by the repeated one-way ANOVA test.

2.2.2. CAP Recognition and Disease Diagnostic Scheme

Inspired by the results of the statistical analysis of the cardiopulmonary characteristics of healthy people and patients with sleep-related diseases, we propose an automatic recognition and disease diagnostic scheme of CAP based on CRI. Using the Hidden Markov and Random Forest model, the scheme combines the values and stability of characteristics which measure the status of cardiopulmonary system of people during sleep. Precision and recall rates will be calculated in the next section.
The flow of the classification and diagnosis scheme is shown in the Figure 4 above. In our work, movement characteristics and cardiopulmonary characteristics are extracted every 30 s, along with improved Hidden Markov and Random Forest model, to classify the sleep stages. Support Vector Machine (SVM) is then used to diagnose diseases.
The classification algorithm in this paper uses a total of 21 features, consisting of three body movement (BM) features, 12 heart rate variability (HRV) features and six CRI features. Each kind of data is preprocessed.
Body movement (BM) features [58] include the variance, approximate entropy and sample entropy of three cycle acceleration signals. HRV features used in this paper can be split into three categories following existing research [59]. Three time-domain features: SDNN, RMSSD, pNN50. Four frequency-domain features: VLF power, LF power, HF power, LF/HF ratio. Five nonlinear features: SD1, SD2 and SD1/SD2 from Poincare plot, as well as sample entropy and approximate entropy. CRI features include follows: CRR, CRB, CRQ, CRA, the variance of CRA, the entropy of CRA.
During the training process, to improve the effectiveness of the extracted features, for the sleep stage WEAK, LIGHT SLEEP (stages S1 and S2) and DEEP SLEEP (stages S3 and S4), 20 consecutive tagged data is selected as the training data, that is, the data remaining unchanged in sleep status within 10 min is selected as the initial training data. For the identification of A1-phase, A2-phase and A3-phase, all the data are used. When the value of the calculated feature exceeds the normal range, the data fragment is preprocessed again, and the outliers are compensated until the eigenvalue is normal.
The scheme adopts the combination of Hidden Markov and Random Forest model to fully utilize the variation and stability of features.
The Hidden Markov Model is used to characterize both the stability and the variation of CRA between different sleep stages, and the final recognition scheme is realized by the Random Forest Model.
The Hidden Markov Model is a probabilistic model of timing sequences. It can learn the implicit characteristics of states by observing phenomena. The transformation process between sleep stage is a Markov process with implicit temporal information. The physiological signals and sleep states can be observed, while the numerical characteristics of physiological significance of the state transitions are unclear. Therefore, we choose the hidden Markov model as the initial model. In this study, for each subject, the transformation between the sleep stages and phases is a Markov chain.
We take the CRA sequence as the observation sequence and the label of the sleep stages as the state sequence. There are 18 kinds of sleep stages in this study, including WEAK, S1 (A1, A2, A3 and NA), S2 (A1, A2, A3 and NA), S3 (A1, A2, A3 and NA), S4 (A1, A2, A3 and NA) and REM. I is the state sequence of 18 dimension, and O is the corresponding observation sequence of CRA.
I = ( i 1 , i 2 , , i 18 )
A is the probability matrix of the observation transition:
A = [ a i j ] 18 × 18
a i j = P ( O t + 1 = q j | O t = q i ) , i = 1 , 2 , , 18 ; j = 1 , 2 , , N
where aij is the probability of the transformation to qj at time t + 1, given qi at time t.
In the CAP classification scheme, firstly, we train the HMM [44] for CRA and corresponding sleep stage tags. Then, the values along the diagonal of A and the sum of each column are calculated as the new eigenvalues of the states of the corresponding columns to represent the stability of CRA. Finally, we put the two new features and the 21 features mentioned earlier as the inputs of the Random Forest Model.
The Random Forest Model has natural advantages for classification tasks with many features and incomplete sample size, with high accuracy and fast training speed. It is suitable as our final classifier. The training and classification work are carried out in two layers, the first layer is carried out for identification stage WEAK, stages REM, S1, S2, S3 and stage S4, while the second layer is carried out for identifying A1-phase, A2-phase and A3-phase in S1/S2/S3/S4, respectively. The clarification of the two layers both used Random Forest Model.
SVM [60] is used in the scheme to diagnose insomnia and narcolepsy. The output of the Random Forest Model and the 23 characteristics above are used for training of the SVM. For subjects diagnosed with sleep-related disorders, their eigenvalues will be recalculated and retained in the training library.

3. Results

3.1. Results of the Statistical Analysis of CRI in People with Non-Pathology, Insomnia and Narcolepsy

This section describes the statistical analysis of the cardiopulmonary characteristics of people with non-pathology and several sleep-related pathologies.
The deep sleep stage period (including stages S3 and S4) during the whole sleep, in which the human body is least affected by external noise and voluntary activities, shows more remarkable physiological significance in analyzing statistical results [61]. Previous studies have shown that CRI characteristics are stable and have statistical significance during deep sleep stage [43].
Therefore, CRI characteristics (including CRQ, CRR, CRB and CRA) in people with non-pathology, insomnia and narcolepsy during the DEEP SLEEP period were analyzed. The results are shown in the Figure 5 below. The values of CRI are normalized in the figure for convenience of comparison. The blue bar shows the results of group with non-pathology. The red bar shows the results of group with insomnia. The green represents the results in subjects with narcolepsy.
In the bar chart Figure 5, we can see that, compared to people with non-pathology, insomniacs have higher CRQ and lower CRR, while narcolepsy patients show the opposite manifestations. The CRB of the healthy group is significantly higher than the CRB of insomniacs and narcolepsy. At the same time, the CRA of insomniac patients is higher, while the CRA of narcolepsy patients is significantly lower than the CRA of the healthy group. The overall coupling of depth and width of narcolepsy patients are low.
For the analysis of CAP mode, the repeated one-way ANOVA, followed by Dunnett’s post hoc test [62] is used to represent the significant difference of CRA in the A-phase and non-A phase (NA-phase) period during the sleep stages S3 and S4 in Table 1 below as an example.
It is a method of comparing means in analysis of variance to judge whether the influence of CRA on the sleep stage is significant [62]. For repeated one-way ANOVA, if p < 0.05, there is significant difference between the groups. We calculated the p value of CRA in S3 and S4: p = 0.004. For Dunnett’s post hoc test, if the value of difference of the mean > LSR, there is a significant difference between the groups being compared.
In the Table 1, it can be seen that CRA, as a measure of coupling depth, has significant differences between A-phase and NA-phase. For the different groups, measurements of people with non-pathology, insomnia and narcolepsy were analyzed and compared to each other below.
In the groups with non-pathology, insomnia and narcolepsy, by comparing the performance of CRA in period A1-phase, A2-phase and A3-phase in sleep stages S1, S2, S3 and S4 shown in Figure 6, the following results are obtained.
In the non-pathology group, in stage S4, there is a significant drop in CRA from A2 to A3; in stage S3, there is a huge drop between A1 and A2, whereas S1 and S2 are inherently less coupled. CRA of A3 in stages S3 and S4 is lower than in stages S1 and S2 in the NA-phase. In insomniacs, CRA of A3-phase in stage S3 is very low, while it is very high in the NA-stage. We can see that CRA fluctuates greatly between phase A and non-Phase A. This may be a compensatory mechanism of the autonomic nervous system in insomniacs. In narcolepsy, the CRA in S3 and S4 is small regardless of whether they are in an A-phase or NA-phase period.
To facilitate the analysis and comparison between three groups, Figure 7 is drawn. From Figure 7, compared with the healthy group, the cardiopulmonary performance of insomniacs fluctuates greatly with the generation of A-phase. The value of CRA is very large in the non-A phase period and very small in the A-phase period. However, the value of CRA of narcolepsy patients is generally small. Combined with the overall analysis of CRB and CRR, it can be concluded that the cardiopulmonary coupling damping of narcolepsy patients is relatively high.
Overall, from the detailed analysis of the cardiopulmonary characteristics of different groups during the sleep period, we could see that CRI characteristics (including CRQ, CRR, CRB and CRA) in people with non-pathology, insomnia and narcolepsy are different during the whole sleep. The CRI could help diagnose insomnia and narcolepsy. In the CAP mode of different groups, both the variation and stability of CRA characteristics are significantly different. Constructing a classifier combining variation and stability can help us to recognize CAP more accurately.

3.2. Results of the Recognition and Disease Diagnostic Scheme

To illustrate the classification performance and generalization ability of the scheme, we conducted the 7-fold cross-validation [63]. We proved the validity of the sleep stage classification method by comparing our estimates with the expert diagnosis results.
The sleep stage classification results are shown in the confusion matrix Table 2. Each value in the table represents the number of samples. Each column represents the prediction category. Each row represents the actual category the data belong to. The results include 18 categories: Wake, REM, S1, S2 and S3, S4. In the S1, S2, S3 and S4, there are A1, A2, A3 and NA in every stage. According to Table 2, the accuracy and recall rates of the different classes could be calculated.
In order to express classification accuracy and generalization performance form Table 2, the F1 score is defined as the average of accuracy and recall rates of all the eighteen classes [64].
F 1 i = 2 * p r e c i s i o n i * r e c a l l i p r e c i s i o n i + r e c a l l i
In this formula, i stands for class i, including WEAK, REM and A1, A2, A3 and NA from stage 1 to stage 4. The precision refers to the specific gravity of the positive sample in the positive example determined by the classifier; recall refers to the proportion of the total positive cases that are predicted to be positive. Then the final F1 score in CAP experiments is [64]:
F 1 = i = 1 18 F 1 i 18
According to the Table 2, we could calculate that F1 score of sleep-wake classification is 92.0%. F1 score in the sleep stage classification (WEAK, REM, S1, S2, S3 and S4) is 83.8%, while the F1 score in the CAP experiments is 80.4%.
In detail, about the identification of A-phase, the accuracy in stages S1, S2, S3 and S4 are 84.4%, 90.1%, 84.2% and 79.9% respectively. The average recognition rate of A-phase reaches 84.7%. Similarly, in terms of disease diagnosis, the F1 score of diagnosis and recognition is 87.8% for insomnia and 90.0% for narcolepsy. For real cases, this table could help to analysis differences between the results of sleep classification and the actual sleep structure easily.

4. Discussion

CRI characteristics (including CRQ, CRR, CRB and CRA) in people with non-pathology, insomnia and narcolepsy are different in all stages during the whole sleep.
Compared to people with non-pathology, insomniacs show higher CRQ and lower CRR, while narcolepsy patients have the opposite manifestations. The CRB of the healthy group is significantly higher than the CRB of insomniacs and narcolepsy. At the same time, the CRA of insomniac patients is higher than the healthy group, while the CRA of narcolepsy patients is significantly lower than that of the healthy group. This phenomenon can be considered as a compensation of autonomic nerve of insomniac patients. Meanwhile, the coupling depth CRA and coupling width CRB of narcolepsy patients are always low during the whole sleep.
In the CAP, in stage S4, there is a significant drop in CRA from A2 to A3; in stage S3, there is a sharp drop between A1 and A2, whereas S1 and S2 are inherently less coupled. CRA of A3 in stages S3 and S4 is lower than in stages S1 and S2 in the NA-phase. These performances support the conclusions including the presence of heart-brain connection [39] and the EEG performance of A2 and A3 type [50]. Compared to people with non-pathology, the cardiopulmonary performance of insomniacs fluctuates greatly with the generation of A-phase. The CRA is very large in the non-A phase period and very small in the A-phase period. However, the CRA of narcolepsy patients is generally low. Combined with the analysis of CRB and CRR, it can be concluded that the cardiopulmonary coupling damping of narcolepsy patients is relatively high.
In the CAP mode of different groups, both the variation and stability of CRA are significantly different. Our scheme using improved HMM and RF combines these two points well. The average recognition rate of A-phase reaches 84.7% and the F1 score in the CAP experiments reaches 80.4%.
In order to illustrate the role of CRI features in the scheme, as the control group, the features related to CRI were removed while the remaining features were used. After the training of the Random Forest Model and 7-fold cross-validation test, the F1 score of CAP classification is 73.8%. The F1 score of disease diagnosis is 71.6%. Besides, as comparison, features in references [34,35,36,37,38] are used in our scheme instead of CRI. Meanwhile, research on the feature extraction methods of ECG in recent ten years [65,66,67,68] are also considered. The results are shown in Table 3.
It can be seen that the performance of the CRI surpasses that of the features in other studies in sleep stage classification, especially in CAP pattern recognition. In the diagnosis of sleep-related disorders, CRI significantly outperforms other features.
The performances of ‘removing CRI’ and ‘CRI’ are better than the other benchmark papers, such as [34,35,36,37,38]. These papers start with the heart rate and body movement signals, try to find features that can characterize sleep stages. The HRV features in our classifying framework, include: three time-domain features: SDNN, RMSSD, pNN50; four frequency-domain features: VLF power, LF power, HF power, LF/HF ratio; five nonlinear features: SD1, SD2, SD1/SD2, sample entropy and approximate entropy. Body movement features include the variance, approximate entropy and sample entropy of three cycle acceleration signals. These are more comprehensive than those in these papers. ‘CRI’ performances better than [67,68], which uses wavelet filter bank to analysis the heart rate characteristics. The wavelet analysis could dig out heart rate characteristics in depth, which are closely related to the activities of cardiopulmonary system, to some extent. Compared with the features of the directly acquired signals, such as ECG, three-axis acceleration signals, etc., CRI deeply explores the degree of cardiopulmonary coupling in different stages of sleep, shows the physiological significance of the autonomic nerve regulation.
CRI characteristics quantify the coupling status in the cardiopulmonary system, and could express the activity characteristics of autonomic nerve regulation. These autonomic nerve regulation characteristics are determined by EEG, and ultimately manifest as tiny periodic changes in sleep stages. At the same time, using Hidden Markov and Random Forest model, the scheme proposed combines the values and stability of these characteristics during sleep. This makes the results of classification more accurate. For the groups with non-pathology, insomnia and narcolepsy, CRI show significant differences in various stages of sleep and express the pathological features of the disease. Thus CRI make the framework more effective in diagnosing diseases.
The results of the scheme, on the one hand, prove the existence of a heart-brain connection, and on the other hand, verify the validity of CRI in the expression of the cardiopulmonary coupling and autonomic nervous activity of human being. The scheme is also helpful to identify insomniac patients and narcolepsy patients.

5. Conclusions

Cyclic Alternating Pattern (CAP) is a sensitive tool for studying sleep disorders throughout the life cycle. Several efforts have been made to develop a reliable automatic CAP-scoring algorithm. Most of these methods rely on the extraction of spectral features from the EEG.
Considering the connection between the heart and brain of people, we analyze the cardiopulmonary characteristics of cardiopulmonary resonance indices (CRI) in people with non-pathology, insomnia and narcolepsy under the CAP mode. We analyze and identify key fields that contribute to insomnia and narcolepsy. The results show that CRI have different manifestations in different groups. The conclusion is drawn that CRI are capable of representing the cardiopulmonary coupling degree, autonomic nerve state and sensitivity of the subjects, and thus are able to measure the health status of human body.
Inspired by the results of the statistical analysis of the cardiopulmonary characteristics, our CAP Recognition and Disease Diagnostic scheme uses CRI in the frequency domain as the feature to recognize the A-phase stage during the whole sleep.
There are many trials in pattern recognition using Hidden Markov and Random Forest model. Different methods have shown superiority in different applications. The scheme in this article takes both variation and stability of measurements of coupling state of the cardiopulmonary system during sleep into account. The precision and recall show the good performance of the scheme.
The proposed method is effective to conquer the weakness about the signal acquisition and processing in clinical trials, particularly when we consider the stress of the examinee and the cost of measurements using the conventional Rechtschaffen & Kales rules [50]. Results show that the scheme could automatically recognize the Cyclic Alternating Pattern accurately, and hence help to diagnose insomnia and narcolepsy. Further, we would conduct more clinical trials to validate and improve our model.
To develop a classification and diagnosis scheme for a specific hospital, the clinical data and expert opinions should be added into the scheme. It could be interesting to see how these contribute to the prediction accuracy.
Furthermore, the meaning and application of CRI are also worth exploring, such as the psychological stress, inflammation, etc.

Author Contributions

J.C., J.W. and Z.H. designed the study and performed the research. J.C. analyzed the data and wrote the paper. Z.H. and J.W. reviewed edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by National Natural Science Foundation of China, contract 61431017.

Institutional Review Board Statement

This study was conducted in an open database and did not require ethical approval.

Informed Consent Statement

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

Data Availability Statement

The data for this article is publicly available on the website https://www.physionet.org/content/capslpdb (accessed on 1 January 2022).

Acknowledgments

The authors would like to thank the doctors of China-Japan Friendship Hospital for their carefully interpretation.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Knutson, K.L. Sleep duration and cardiometabolic risk: A review of the epidemiologic evidence. Best Pr. Res. Clin. Endocrinol. Metab. 2010, 24, 731–743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Campbell, R.; Vansteenkiste, M.; Delesie, L.; Tobback, E.; Mariman, A.; Vogelaers, D.; Mouratidis, A. Reciprocal associations between daily need-based experiences, energy, and sleep in chronic fatigue syndrome. Health Psychol. 2018, 37, 1168. [Google Scholar] [CrossRef] [PubMed]
  3. Oliveira, M.M.; Akerman, S.; Hoffmann, J.; Goadsby, P. EHMTI-0287. Peptides involved in sleep and appetite homeostatic regulation and its effects in the modulation of trigeminovascular nociceptive activation. J. Headache Pain 2014, 15, 1. [Google Scholar] [CrossRef] [Green Version]
  4. Becker, S.P.; Jarrett, M.A.; Luebbe, A.M.; Garner, A.A.; Burns, G.L.; Kofler, M.J. Sleep in a large, multi-university sample of college students: Sleep problem prevalence, sex differences, and mental health correlates. Sleep Health 2018, 4, 174–181. [Google Scholar] [CrossRef] [PubMed]
  5. Kaniecki, R.G. Sleep and tension-type headache. Curr. Neurol. Neuroen. Rep. 2019, 6, 556. [Google Scholar]
  6. Lu, J.; Sherman, D.M.; Devor, M.; Saper, C.B. A putative flip–flop switch for control of REM sleep. Nature 2006, 441, 589–594. [Google Scholar] [CrossRef] [PubMed]
  7. McCarley, R.W. Neurobiology of REM and NREM sleep. Sleep Med. 2007, 8, 302–330. [Google Scholar] [CrossRef]
  8. Merica, H.; Blois, R.; Gaillard, J.-M. Spectral characteristics of sleep EEG in chronic insomnia. Eur. J. Neurosci. 1998, 10, 1826–1834. [Google Scholar] [CrossRef]
  9. Putilov, A.A. Spectral EEG indicator of pressure to enter into deep sleep: Its responsiveness to closing the eyes for just a few minutes exhibits a pure exponential buildup during sleep deprivation. Biol. Rhythm Res. 2016, 47, 773–786. [Google Scholar] [CrossRef]
  10. Sharon, O.; Nir, Y. Attenuated Fast Steady-State Visual Evoked Potentials During Human Sleep. Cereb. Cortex 2018, 28, 1297–1311. [Google Scholar] [CrossRef]
  11. Terzano, M.G.; Parrino, L.; Sherieri, A.; Chervin, R.; Chokroverty, S.; Guilleminault, C.; Hirshkowitz, M.; Mahowald, M.; Moldofsky, H.; Rosa, A.; et al. Atlas, rules, and recording techniques for the scoring of cyclic alternating pattern (CAP) in human sleep. Sleep Med. 2001, 2, 537–553. [Google Scholar] [CrossRef]
  12. Parrino, L.; Boselli, M.; Spaggiari, M.C.; Smerieri, A.; Terzano, M.G. Cyclic alternating pattern (CAP) in normal sleep: Polysomnographic parameters in different age groups. Electroencephalogr. Clin. Neurophysiol. 1998, 107, 439–450. [Google Scholar] [CrossRef]
  13. Feinberg, I.; Floyd, T.C. Systematic Trends across the Night in Human Sleep Cycles. Psychophysiology 1979, 16, 283–291. [Google Scholar] [CrossRef] [PubMed]
  14. Elliott, R.; Mckinley, S.; Cistulli, P. The quality and duration of sleep in the intensive care setting: An integrative review. Int. J. Nurs. Stud. 2011, 48, 384–400. [Google Scholar] [CrossRef]
  15. Priano, L.; Bigoni, M.; Albani, G.; Sellitti, L.; Giacomotti, E.; Picconi, R.; Cremascoli, R.; Zibetti, M.; Lopiano, L.; Mauro, A. Sleep microstructure in Parkinson’s disease: Cycling alternating pattern (CAP) as a sensitive marker of early NREM sleep instability. Sleep Med. 2019, 61, 57–62. [Google Scholar] [CrossRef]
  16. Angriman, M.; Caravale, B.; Novelli, L.; Ferri, R.; Bruni, O. Sleep in Children with Neurodevelopmental Disabilities. Neuropediatrics 2015, 46, 199–210. [Google Scholar] [CrossRef] [Green Version]
  17. Khachatryan, S.G.; Prosperetti, C.; Rossinelli, A.; Pedrazzi, P.; Agazzi, P.; Ratti, P.L.; Manconi, M. Sleep-onset central apneas as triggers of severe nocturnal seizures. Sleep Med. 2015, 16, 1017–1019. [Google Scholar] [CrossRef]
  18. Bruni, O.; Novelli, L.; Miano, S.; Parrino, L.; Terzano, M.G.; Ferri, R. Cyclic alternating pattern: A window into pediatric sleep. Sleep Med. 2010, 11, 628–636. [Google Scholar] [CrossRef]
  19. Shusterman, V.; Troy, W.C.; Abdelmessih, M.; Hoffman, S.; Nemec, J.; Strollo, P.J.; London, B.; Lampert, R. Nighttime instabilities of neurophysiological, cardiovascular, and respiratory activity: Integrative modeling and preliminary results. J. Electrocardiol. 2015, 48, 1010–1016. [Google Scholar] [CrossRef] [Green Version]
  20. Medina, D.I.; Mendez, M.O.; Murguia, J.S.; Chouvarda, I. Wavelet Singularity Analysis for CAP Sleep Delineation. In Proceedings of the 2019 IEEE 19th International Conference on Bioinformatics and Bioengineering (BIBE), Athens, Greece, 28–30 October 2019; pp. 764–769. [Google Scholar]
  21. Tripathy, R.K.; Acharya, U.R. Use of features from RR-time series and EEG signals for automated classification of sleep stages in deep neural network framework. Biocybern. Biomed. Eng. 2018, 38, 890–902. [Google Scholar] [CrossRef]
  22. Hartmann, S.; Baumert, M. Automatic A-Phase Detection of Cyclic Alternating Patterns in Sleep Using Dynamic Temporal Information. IEEE Trans. Neural Syst. Rehabil. Eng. 2019, 27, 1695–1703. [Google Scholar] [CrossRef] [PubMed]
  23. Largo, R.; Munteanu, C.; Rosa, A. CAP event detection by wavelets and GA tuning. In Proceedings of the IEEE International Workshop on Intelligent Signal Processing, Faro, Portugal, 1–3 September; 2005; pp. 44–48. [Google Scholar]
  24. Zhou, X.; Yao, D.; Zhu, M.; Zhang, X.; Qi, L.; Pan, H.; Zhu, X.; Wang, Y.; Zhang, Z. Vigilance detection method for high-speed rail using wireless wearable EEG collection technology based on low-rank matrix decomposition. IET Intell. Transp. Syst. 2018, 12, 819–825. [Google Scholar] [CrossRef]
  25. Robbins, K.; Su, K.-M.; Hairston, W.D. An 18-subject EEG data collection using a visual-oddball task, designed for benchmarking algorithms and headset performance comparisons. Data Brief 2018, 16, 227–230. [Google Scholar] [CrossRef] [PubMed]
  26. Mingai, L.; Shuoda, G.; Guoyu, Z.; Yanjun, S.; Jinfu, Y. Removing ocular artifacts from mixed EEG signals with fastkica and DWT. J. Intell. Fuzzy Syst. 2015, 28, 2851–2861. [Google Scholar] [CrossRef]
  27. Gao, Y.; Li, H.; Himed, B. Adaptive Subspace Tests for Multichannel Signal Detection in Auto-Regressive Disturbance. IEEE Trans. Signal Process. 2018, 66, 5577–5587. [Google Scholar] [CrossRef]
  28. Hua, X.; Ono, Y.; Peng, L.; Cheng, Y.; Wang, H. Target detection within nonhomogeneous clutter via total bregman divergence-based matrix information geometry detectors. IEEE Trans. Signal Process. 2021, 69, 4326–4340. [Google Scholar] [CrossRef]
  29. Nie, D.; Xie, K.; Zhou, F.; Qiao, G. A Correlation Detection Method of Low SNR Based on Multi-Channelization. IEEE Signal Process. Lett. 2020, 27, 1375–1379. [Google Scholar] [CrossRef]
  30. Dora, C.; Biswal, P.K. Efficient detection and correction of variable strength ECG artifact from single channel EEG. Biomed. Signal Process. Control 2019, 50, 168–177. [Google Scholar] [CrossRef]
  31. Olsen, M.; Mignot, E.; Jennum, P.J.; Sorensen, H.B.D. Robust, ECG-based detection of Sleep-disordered breathing in large population-based cohorts. Sleep 2020, 43, 276. [Google Scholar] [CrossRef]
  32. Meng, C.X.; Di, L.; Xiao, Z.; Wang, S.Q.; Ji, P. Study on eeg changes after occlusal trauma during sleep in rats. Shanghai Kou Qiang Yi Xue Shanghai J. Stomatol. 2014, 23, 26–29. [Google Scholar]
  33. Sadr, N.; de Chazal, P. A comparison of three ECG-derived respiration methods for sleep apnoea detection. Biomed. Phys. Eng. Express 2019, 5, 025027. [Google Scholar] [CrossRef]
  34. Scherz, W.D.; Fritz, D.; Velicu, O.R.; Seepold, R.; Madrid, N.M. Heart rate spectrum analysis for sleep quality detection. EURASIP J. Embed. Syst. 2017, 2017, 26. [Google Scholar] [CrossRef] [Green Version]
  35. Adnane, M.; Jiang, Z.; Yan, Z. Sleep–wake stages classification and sleep efficiency estimation using single-lead electrocardiogram. Expert Syst. Appl. 2012, 39, 1401–1413. [Google Scholar] [CrossRef]
  36. Mendez, M.O.; Matteucci, M.; Castronovo, V.; Ferini-Strambi, L.; Cerutti, S.; Bianchi, A. Sleep staging from heart rate variability: Time-varying spectral features and hidden Markov models. Int. J. Biomed. Eng. Technol. 2010, 3, 246–263. [Google Scholar] [CrossRef]
  37. Rothenberger, S.D.; Krafty, R.; Taylor, B.J.; Cribbet, M.R.; Thayer, J.F.; Buysse, D.J.; Kravitz, H.M.; Buysse, E.D.; Hall, M.H. Time-varying correlations between delta EEG power and heart rate variability in midlife women: The SWAN Sleep Study. Psychophysiology 2015, 52, 572–584. [Google Scholar] [CrossRef] [Green Version]
  38. de Groot, E.R.; Knoop, M.S.; Hoogen, A.V.D.; Wang, X.; Long, X.; Pillen, S.; Benders, M.; Dudink, J. The value of cardiorespiratory parameters for sleep state classification in preterm infants: A systematic review. Sleep Med. Rev. 2021, 58, 101462. [Google Scholar] [CrossRef]
  39. Wolters, F.J.; Ikram, M.A. Response by Wolters and Ikram to Letter Regarding Article, “Cerebral Perfusion and the Risk of Dementia: A Population-Based Study”. Circulation 2018, 137, 1416–1417. [Google Scholar] [CrossRef]
  40. Roger, V.L. Heart Failure as a Risk Factor for Stroke: Another Facet of the Heart–Brain Connection. J. Card. Fail. 2019, 25, 448–449. [Google Scholar] [CrossRef]
  41. Hillman, C.H.; Erickson, K.I.; Kramer, A.F. Be smart, exercise your heart: Exercise effects on brain and cognition. Nat. Rev. Neurosci. 2008, 9, 58–65. [Google Scholar] [CrossRef]
  42. Cui, J.; Huang, Z.; Wu, J.; Jiang, H. Cardiopulmonary Resonance Function and Indices—A Quantitative Measurement for Respiratory Sinus Arrhythmia. Front. Physiol. 2020, 11, 867. [Google Scholar] [CrossRef]
  43. Cui, J.; Huang, Z.; Jiaerken, D.; Fan, Y.; Zhao, S.; Zhang, L.; Wu, J. A wearable system for cardiopulmonary assessment and personalized respiratory training. Futur. Gener. Comput. Syst. 2020, 112, 1131–1140. [Google Scholar] [CrossRef]
  44. Ding, J.; Ye, M. Radar emitter identification using hidden markov model. Comput. Model. Eng. Ences 2019, 123, 401–426. [Google Scholar]
  45. El Haouij, N.; Poggi, J.M.; Ghozi, R.; Sevestre-Ghalila, S.; Jaïdane, M. Random forest-based approach for physiological functional variable selection for driver’s stress level classification. Stat. Methods Appl. 2019, 28, 157–185. [Google Scholar] [CrossRef]
  46. Goldberger, A.L.; Amaral, L.A.; Glass, L.; Hausdorff, J.M.; Ivanov, P.C.; Mark, R.G.; Mietus, J.E.; Moody, G.B.; Peng, C.K.; Stanley, H.E. PhysioBank, PhysioToolkit, and PhysioNet: Components of a new research resource for complex physiologic signals. Circulation 2000, 101, 215–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Parrino, L.; Ferri, R.; Bruni, O.; Terzano, M.G. Cyclic alternating pattern (CAP): The marker of sleep instability. Sleep Med. Rev. 2012, 16, 27–45. [Google Scholar] [CrossRef] [PubMed]
  48. Bounhoure, J.P.; Galinier, M.; Didier, A.; Leophonte, P. Sleep apnea syndromes and cardiovascular disease. Bull. l’Acad. Natl. Méd. 2005, 189, 445–459; discussion 460–464. [Google Scholar]
  49. Poryazova, R.; Werth, E.; Parrino, L.; Terzano, M.G.; Bassetti, C.L. Cyclic alternating pattern in narcolepsy patients and healthy controls after partial and total sleep deprivation. Clin. Neurophysiol. 2011, 122, 1788–1793. [Google Scholar] [CrossRef]
  50. Rechtscahffen, A.; Kales, A. A manual of standardized terminology, techniques and scoring system for sleep stages in human subjects. Arch. Gen. Psychiatry 1968, 20, 246–247. [Google Scholar]
  51. Negro, C.A.D.; Morgado-Valle, C.; Feldman, J.L. Respiratory rhythm: An emergent network property? Neuron 2002, 34, 821–830. [Google Scholar] [CrossRef] [Green Version]
  52. Dong, F.; Langford, W. Models of Cheyne-Stokes respiration with cardiovascular pathologies. J. Math. Biol. 2008, 57, 497–519. [Google Scholar] [CrossRef]
  53. Janse, M.J. The premature beat. Cardiovasc. Res. 1992, 26, 89–100. [Google Scholar] [CrossRef] [PubMed]
  54. Epskamp, S.; Waldorp, L.J.; Mõttus, R.; Borsboom, D. The Gaussian graphical model in cross-sectional and time-series data. Multivar. Behav. Res. 2018, 53, 453–480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lee, R.M.; Tsai, N.C. Dynamic model of integrated cardiovascular and respiratory systems. Math. Methods Appl. Sci. 2013, 36, 2224–2236. [Google Scholar] [CrossRef]
  56. Barbini, P.; Brighenti, C.; Cevenini, G.; Gnudi, G. A dynamic morphometric model of the normal lung for studying expiratory flow limitation in mechanical ventilation. Ann. Biomed. Eng. 2005, 33, 518–530. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, H.; Xie, X.; Zhang, C.; Li, Y.; Liu, H.; Hu, Y. Quantitative SSR analysis of series-compensated DFIG-based wind farms using aggregated RLC circuit model. IEEE Trans. Power Syst. 2016, 32, 474–483. [Google Scholar] [CrossRef]
  58. Lan, L.; Qian, X.L.; Lian, Z.W.; Lin, Y.B. Local body cooling to improve sleep quality and thermal comfort in a hot environment. Indoor Air 2018, 28, 135–145. [Google Scholar] [CrossRef] [PubMed]
  59. Zhuang, Z.; Gao, X.; Gao, S. The relationship of HRV to sleep EEG and sleep rhythm. Int. J. Neurosci. 2005, 115, 315–327. [Google Scholar] [CrossRef] [PubMed]
  60. Saunders, M.O.; Stitson, J.; Weston, R.H.; Bottou, L.; Scholkopf, B. Support vector machine. Comput. Sci. 2002, 1, 1–28. [Google Scholar]
  61. Fattinger, S.; de Beukelaar, T.T.; Ruddy, K.L.; Volk, C.; Heyse, N.C.; Herbst, J.A.; Hahnloser, R.H.; Wenderoth, N.; Huber, R. Deep sleep maintains learning efficiency of the human brain. Nat. Commun. 2017, 8, 15405. [Google Scholar] [CrossRef] [Green Version]
  62. Brown, A.M. A new software for carrying out one-way ANOVA post hoc tests. Comput. Methods Programs Biomed. 2015, 79, 89–95. [Google Scholar] [CrossRef]
  63. Jung, Y. Multiple predicting k-fold cross-validation for model selection. J. Nonparametric Stat. 2018, 30, 197–215. [Google Scholar] [CrossRef]
  64. Crestani, F.; Lalmas, M.; van Rijsbergen, C.J. (Eds.) Information Retrieval: Uncertainty and Logics: Uncertainty and Logics: Advanced Models for the Representation and Retrieval of Information; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  65. Kirk, P.A.; Garfinkel, S.; Robinson, O.J. RapidHRV: An open-source toolbox for extracting heart rate and heart rate variability. PsyArXiv 2021. [Google Scholar] [CrossRef]
  66. Kondo, H.; Ozone, M.; Ohki, N.; Sagawa, Y.; Yamamichi, K.; Fukuju, M.; Yoshida, T.; Nishi, C.; Kawasaki, A.; Mori, K.; et al. Association between heart rate variability, blood pressure and autonomic activity in cyclic alternating pattern during sleep. Sleep 2014, 37, 187–194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Sharma, M.; Dhiman, H.S.; Acharya, U.R. Automatic identification of insomnia using optimal antisymmetric biorthogonal wavelet filter bank with ECG signals. Comput. Biol. Med. 2021, 131, 104246. [Google Scholar] [CrossRef] [PubMed]
  68. Tzabazis, A.; Eisenried, A.; Yeomans, D.C.; Hyatt, M.I. Wavelet analysis of heart rate variability: Impact of wavelet selection. Biomed. Signal Proc. Control 2018, 40, 220–225. [Google Scholar] [CrossRef]
Figure 1. An example of the EEG in cyclic alternating pattern (CAP) in sleep stage 2. The horizontal axis represents the sampling point, and the vertical axis represents the amplitude of the signal. The shapes of the EEG signals in red, blue, and green boxes correspond to the CAP−A1, CAP−A2, and CAP−A3, respectively.
Figure 1. An example of the EEG in cyclic alternating pattern (CAP) in sleep stage 2. The horizontal axis represents the sampling point, and the vertical axis represents the amplitude of the signal. The shapes of the EEG signals in red, blue, and green boxes correspond to the CAP−A1, CAP−A2, and CAP−A3, respectively.
Sensors 22 02225 g001
Figure 2. Schematic diagram of cardiopulmonary system and its circuit model. Interaction between lung and heart resembles the energy flow between inductor and capacitor. The non-respiration factors are equivalent to resistor, damping the resonance.
Figure 2. Schematic diagram of cardiopulmonary system and its circuit model. Interaction between lung and heart resembles the energy flow between inductor and capacitor. The non-respiration factors are equivalent to resistor, damping the resonance.
Sensors 22 02225 g002
Figure 3. Cardiopulmonary Resonance Indices (CRI). CRA is the maximum amplitude of the curve, and CRB is the bandwidth of the curve. FA is the cardiopulmonary resonance frequency.
Figure 3. Cardiopulmonary Resonance Indices (CRI). CRA is the maximum amplitude of the curve, and CRB is the bandwidth of the curve. FA is the cardiopulmonary resonance frequency.
Sensors 22 02225 g003
Figure 4. The classification and diagnosis scheme.
Figure 4. The classification and diagnosis scheme.
Sensors 22 02225 g004
Figure 5. Cardiopulmonary characteristics CRQ, CRR, CRB and CRA during deep sleep of people with non-pathology, insomnia and narcolepsy.
Figure 5. Cardiopulmonary characteristics CRQ, CRR, CRB and CRA during deep sleep of people with non-pathology, insomnia and narcolepsy.
Sensors 22 02225 g005
Figure 6. The change line diagram of CRA of different people in CAP. CRA in period A1-phase, A2-phase, A3-phase and NA-phase in sleep stages S1, S2, S3 and S4 are shown for every group.
Figure 6. The change line diagram of CRA of different people in CAP. CRA in period A1-phase, A2-phase, A3-phase and NA-phase in sleep stages S1, S2, S3 and S4 are shown for every group.
Sensors 22 02225 g006
Figure 7. CRA in NA, A1, A2 and A3 period of different groups in the sleep stages S1, S2, S3, and S4 during the whole sleep.
Figure 7. CRA in NA, A1, A2 and A3 period of different groups in the sleep stages S1, S2, S3, and S4 during the whole sleep.
Sensors 22 02225 g007
Table 1. Repeated one-way ANOVA and Dunnett’s post hoc test. If the value of difference of the mean > LSR, there is a significant difference between the groups being compared (p = 0.004 < 0.05).
Table 1. Repeated one-way ANOVA and Dunnett’s post hoc test. If the value of difference of the mean > LSR, there is a significant difference between the groups being compared (p = 0.004 < 0.05).
CRA
in Deep Sleep
S3S4
Difference of the MeanLSR
(p < 0.05)
Difference of the MeanLSR
(p < 0.05)
A and NA0.1340.0470.2160.096
Table 2. The confusion matrix of the sleep stage classification.
Table 2. The confusion matrix of the sleep stage classification.
Pre S1S2S3S4
WRNAA1A2A3NAA1A2A3NAA1A2A3NAA1A2A3
ActW36221025748723161127058492332231420
R708242028454032231926133122124
S1NA52162789138164119293322111418111330
A17328538381483472111131113
A279301868403392242401130014
A38630613173091142700140005
S2NA5526250931289487645621310424814
A152232742020243421602321122
A249201710316322002611321133
A364331415218172126811130115
S3NA3081006200161148119110
A12231001110446527100
A21021001110243136100
A31411101200123336111
S4NA217000151001920015111105
A11130003000611024431
A21121001100231124312
A31311001100232211437
The value of row i and column j represents the number of the sections whose truth value is class i and is predicted to be class j. Each column represents the prediction category; Each row represents the actual category of data belonging to. ‘W’ represents Wake, ‘R’ represents REM. ‘Act’ represents the number of actual sleep stage; ‘Pre’ represents the number of the predicted sleep stage.
Table 3. The F1 scores of CRI and features of other studies in experiments of sleep stage classification and disease diagnosis.
Table 3. The F1 scores of CRI and features of other studies in experiments of sleep stage classification and disease diagnosis.
MethodSleep-Wake ClassificationS1, S2, S3, S4 and Wake Stage ClassificationCAP RecognitionDisease Diagnosis
Heart rate spectrum analysis [34,65]77.6%72.6%66.7%70.5%
detrended fluctuation analysis [35]78.6%71.4%66.3%64.7%
time-varying spectral features [36,37]82.0%76.6%70.3%72.5%
Heart rate fluctuations [38,66]79.9%73.1%66.7%70.5%
wavelet filter bank [67,68]90.1%82.6%76.7%80.9%
Removing CRI85.9%77.7%73.8%71.6%
CRI92.0%83.8%80.4%88.9%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cui, J.; Huang, Z.; Wu, J. Automatic Detection of the Cyclic Alternating Pattern of Sleep and Diagnosis of Sleep-Related Pathologies Based on Cardiopulmonary Resonance Indices. Sensors 2022, 22, 2225. https://doi.org/10.3390/s22062225

AMA Style

Cui J, Huang Z, Wu J. Automatic Detection of the Cyclic Alternating Pattern of Sleep and Diagnosis of Sleep-Related Pathologies Based on Cardiopulmonary Resonance Indices. Sensors. 2022; 22(6):2225. https://doi.org/10.3390/s22062225

Chicago/Turabian Style

Cui, Jiajia, Zhipei Huang, and Jiankang Wu. 2022. "Automatic Detection of the Cyclic Alternating Pattern of Sleep and Diagnosis of Sleep-Related Pathologies Based on Cardiopulmonary Resonance Indices" Sensors 22, no. 6: 2225. https://doi.org/10.3390/s22062225

APA Style

Cui, J., Huang, Z., & Wu, J. (2022). Automatic Detection of the Cyclic Alternating Pattern of Sleep and Diagnosis of Sleep-Related Pathologies Based on Cardiopulmonary Resonance Indices. Sensors, 22(6), 2225. https://doi.org/10.3390/s22062225

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