Many of the existing electrocardiogram (ECG) toolboxes focus on the derivation of heart rate variability features from RR-intervals. By doing so, they assume correct detection of the QRS-complexes. However, it is highly likely that not all detections are correct. Therefore, it is recommended to visualize the actual R-peak positions in the ECG signal and allow manual adaptations. In this paper we present R-DECO, an easy-to-use graphical user interface (GUI) for the detection and correction of R-peaks. Within R-DECO, the R-peaks are detected by using a detection algorithm which uses an envelope-based procedure. This procedure flattens the ECG and enhances the QRS-complexes. The algorithm obtained an overall sensitivity of 99.60% and positive predictive value of 99.69% on the MIT/BIH arrhythmia database. Additionally, R-DECO includes support for several input data formats for ECG signals, three basic filters, the possibility to load other R-peak locations and intuitive methods to correct ectopic, wrong, or missed heartbeats. All functionalities can be accessed via the GUI and the analysis results can be exported as Matlab or Excel files. The software is publicly available. Through its easy-to-use GUI, R-DECO allows both clinicians and researchers to use all functionalities, without previous knowledge.

The presence of noise sources could reduce the diagnostic capability of the ECG signal and result in inappropriate treatment decisions. To mitigate this problem, automated algorithms to detect artefacts and quantify the quality of the recorded signal are needed. In this study we present an automated method for the detection of artefacts and quantification of the signal quality. The suggested methodology extracts descriptive features from the autocorrelation function and feeds these to a RUSBoost classifier. The posterior probability of the clean class is used to create a continuous signal quality assessment index. Firstly, the robustness of the proposed algorithm is investigated and secondly, the novel signal quality assessment index is evaluated. Data were used from three different studies: a Sleep study, the PhysioNet 2017 Challenge and a Stress study. Binary labels, clean or contaminated, were available from different annotators with experience in ECG analysis. Two types of realistic ECG noise from the MIT-BIH Noise Stress Test Database (NSTDB) were added to the Sleep study to test the quality index. Firstly, the model was trained on the Sleep dataset and subsequently tested on a subset of the other two datasets. Secondly, all recording conditions were taken into account by training the model on a subset derived from the three datasets. Lastly, the posterior probabilities of the model for the different levels of agreement between the annotators were compared. AUC values between 0.988 and 1.000 were obtained when training the model on the Sleep dataset. These results were further improved when training on the three datasets and thus taking all recording conditions into account. A Pearson correlation coefficient of 0.8131 was observed between the score of the clean class and the level of agreement. Additionally, significant quality decreases per noise level for both types of added noise were observed. The main novelty of this study is the new approach to ECG signal quality assessment based on the posterior clean class probability of the classifier.

Objective: Respiratory sinus arrhythmia (RSA) refers to heart rate oscillations synchronous with respiration, and it is one of the major representations of cardiorespiratory coupling. Its strength has been suggested as a biomarker to monitor different conditions and diseases. Some approaches have been proposed to quantify the RSA, but it is unclear which one performs best in specific scenarios. The main objective of this study is to compare seven state-of-the-art methods for RSA quantification using data generated with a model proposed to simulate and control the RSA. These methods are also compared and evaluated on a real-life application, for their ability to capture changes in cardiorespiratory coupling during sleep. Methods: A simulation model is used to create a dataset of heart rate variability and respiratory signals with controlled RSA, which is used to compare the RSA estimation approaches. To compare the methods objectively in a real-life application, regression models trained on the simulated data are used to map the estimates to the same measurement scale. Results and conclusion: RSA estimates based on cross entropy, time-frequency coherence and subspace projections showed the best performance on simulated data. In addition, these estimates captured the expected trends in the changes in cardiorespiratory coupling during sleep similarly. Significance: An objective comparison of methods for RSA quantification is presented to guide future analyses. Also, the proposed simulation model can be used to compare existing and newly proposed RSA estimates. It is freely accessible online.

This study aims at investigating the development of premature infants' autonomic nervous system (ANS) based on a quantitative analysis of the heart-rate variability (HRV) with a variety of novel features. Additionally, the role of heart-rate drops, known as bradycardias, has been studied in relation to both clinical and novel sympathovagal indices. ECG data were measured for at least 3 h in 25 preterm infants (gestational age ≤32 weeks) for a total number of 74 recordings. The post-menstrual age (PMA) of each patient was estimated from the RR interval time-series by means of multivariate linear-mixed effects regression. The tachograms were segmented based on bradycardias in periods after, between and during bradycardias. For each of those epochs, a set of temporal, spectral and fractal indices were included in the regression model. The best performing model has R² = 0.75 and mean absolute error MAE = 1.56 weeks. Three main novelties can be reported. First, the obtained maturation models based on HRV have comparable performance to other development models. Second, the selected features for age estimation show a predominance of power and fractal features in the very-low- and low-frequency bands in explaining the infants' sympathovagal development from 27 PMA weeks until 40 PMA weeks. Third, bradycardias might disrupt the relationship between common temporal indices of the tachogram and the age of the infant and the interpretation of sympathovagal indices. This approach might provide a novel overview of post-natal autonomic maturation and an alternative development index to other electrophysiological data analysis.

Impedance pneumography has been suggested as an ambulatory technique for the monitoring of respiratory diseases. However, its ambulatory nature makes the recordings more prone to noise sources. It is important that such noisy segments are identified and removed, since they could have a huge impact on the performance of data-driven decision support tools. In this study, we investigated the added value of machine learning algorithms to separate clean from noisy bio-impedance signals. We compared three approaches: a heuristic algorithm, a feature-based classification model (SVM) and a convolutional neural network (CNN). The dataset consists of 47 chronic obstructive pulmonary disease patients who performed an inspiratory threshold loading protocol. During this protocol, their respiration was recorded with a bio-impedance device and a spirometer, which served as a gold standard. Four annotators scored the signals for the presence of artefacts, based on the reference signal. We have shown that the accuracy of both machine learning approaches (SVM: 87.77 ± 2.64% and CNN: 87.20 ± 2.78%) is significantly higher, compared to the heuristic approach (84.69 ± 2.32%). Moreover, no significant differences could be observed between the two machine learning approaches. The feature-based and neural network model obtained a respective AUC of $92.77\pm 2.95$% and $92.51\pm 1.74$%. These findings show that a data-driven approach could be beneficial for the task of artefact detection in respiratory thoracic bio-impedance signals.

Transfer entropy ($TE$) has been used to identify and quantify interactions between physiological systems. Different methods exist to estimate $TE$, but there is no consensus about which one performs best in specific applications. In this study, five methods (linear, k-nearest neighbors, fixed-binning with ranking, kernel density estimation and adaptive partitioning) were compared. The comparison was made on three simulation models (linear, nonlinear and linear + nonlinear dynamics). From the simulations, it was found that the best method to quantify the different interactions was adaptive partitioning. This method was then applied on data from a polysomnography study, specifically on the ECG and the respiratory signals (nasal airflow and respiratory effort around the thorax). The hypothesis that the linear and nonlinear components of cardio-respiratory interactions during light and deep sleep change with the sleep stage, was tested. Significant differences, after performing surrogate analysis, indicate an increased $TE$ during deep sleep. However, these differences were found to be dependent on the type of respiratory signal and sampling frequency. These results highlight the importance of selecting the appropriate signals, estimation method and surrogate analysis for the study of linear and nonlinear cardio-respiratory interactions.