Model-Based Evaluation of Methods for Respiratory Sinus Arrhythmia Estimation

Technical aspects of cardiorespiratory estimation using subspace projections and cross entropy

BACKGROUND

Respiratory sinus arrhythmia (RSA) is a form of cardiorespiratory coupling. Its quantification has been suggested as a biomarker to diagnose different diseases. Two state-of-the-art methods, based on subspace projections and entropy, are used to estimate the RSA strength and are evaluated in this paper. Their computation requires the selection of a model order, and their performance is strongly related to the temporal and spectral characteristics of the cardiorespiratory signals.

OBJECTIVE

To evaluate the robustness of the RSA estimates to the selection of model order, delays, changes of phase and irregular heartbeats as well as to give recommendations for their interpretation on each case.

APPROACH

Simulations were used to evaluate the model order selection when calculating the RSA estimates explained before, as well as 3 different scenarios that can occur in signals acquired in non-controlled environments and/or from patient populations: the presence of irregular heartbeats; the occurrence of delays between heart rate variability (HRV) and respiratory signals; and the changes over time of the phase between HRV and respiratory signals.

MAIN RESULTS

It was found that using a single model order for all the calculations suffices to characterize the RSA estimates correctly. In addition, the RSA estimation in signals containing more than 5 irregular heartbeats in a period of 5 minutes might be misleading. Regarding the delays between HRV and respiratory signals, both estimates are robust. For the last scenario, the two approaches tolerate phase changes up to 54°, as long as this lasts less than one fifth of the recording duration.

SIGNIFICANCE

Guidelines are given to compute the RSA estimates in non-controlled environments and patient populations.

Assessing stress variations in children during the strange situation procedure: comparison of three widely used respiratory sinus arrhythmia estimation methods

Objective.The respiratory sinus arrhythmia (RSA) is a well-known marker of vagal activity that can be exploited to measure stress changes. RSA is usually estimated from heart rate variability (HRV). This study aims to compare the RSA obtained with three widely adopted methods showing their strengths and potential pitfalls.Approach.The three methods are tested on 69 healthy preschoolers undergoing a stressful protocol, the strange situation procedure (SSP). We compare the RSA estimated by the Porges method, the univariate autoregressive (AR) spectral analysis of the HRV signal, and the bivariate AR spectral analysis of HRV and respirogram signals. We examine RSA differences detected across the SSP episodes and correlation between the estimates provided by each method.Main results.The Porges and the bivariate AR approaches both detected significant differences (i.e. stress variations) in the RSA measured across the SSP. However, the latter method showed higher sensitivity to stress changes induced by the procedure, with the mean RSA variation between baseline and first separation from the mother (the most stressful condition) being significantly different among methods: Porges, -17.5%; univariate AR, -18.3%; bivariate AR, -23.7%. Moreover, the performances of the Porges algorithm were found strictly dependent on the applied preprocessing.Significance.Our findings confirm the bivariate AR analysis of the HRV and respiratory signals as a robust stress assessment tool that does not require any population-specific preprocessing of the signals and warn about using RSA estimates that neglect breath information in more natural experiments, such as those involving children, in which respiratory frequency changes are extremely likely.

Cardiopulmonary coupling indices to assess weaning readiness from mechanical ventilation

The ideal moment to withdraw respiratory supply of patients under Mechanical Ventilation at Intensive Care Units (ICU), is not easy to be determined for clinicians. Although the Spontaneous Breathing Trial (SBT) provides a measure of the patients' readiness, there is still around 15-20% of predictive failure rate. This work is a proof of concept focused on adding new value to the prediction of the weaning outcome. Heart Rate Variability (HRV) and Cardiopulmonary Coupling (CPC) methods are evaluated as new complementary estimates to assess weaning readiness. The CPC is related to how the mechanisms regulating respiration and cardiac pumping are working simultaneously, and it is defined from HRV in combination with respiratory information. Three different techniques are used to estimate the CPC, including Time-Frequency Coherence, Dynamic Mutual Information and Orthogonal Subspace Projections. The cohort study includes 22 patients in pressure support ventilation, ready to undergo the SBT, analysed in the 24 h previous to the SBT. Of these, 13 had a successful weaning and 9 failed the SBT or needed reintubation -being both considered as failed weaning. Results illustrate that traditional variables such as heart rate, respiratory frequency, and the parameters derived from HRV do not differ in patients with successful or failed weaning. Results revealed that HRV parameters can vary considerably depending on the time at which they are measured. This fact could be attributed to circadian rhythms, having a strong influence on HRV values. On the contrary, significant statistical differences are found in the proposed CPC parameters when comparing the values of the two groups, and throughout the whole recordings. In addition, differences are greater at night, probably because patients with failed weaning might be experiencing more respiratory episodes, e.g. apneas during the night, which is directly related to a reduced respiratory sinus arrhythmia. Therefore, results suggest that the traditional measures could be used in combination with the proposed CPC biomarkers to improve weaning readiness.

Benchmarking Transfer Entropy Methods for the Study of Linear and Nonlinear Cardio-Respiratory Interactions

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.

Linear and Non-linear Quantification of the Respiratory Sinus Arrhythmia Using Support Vector Machines

Respiratory sinus arrhythmia (RSA) is a form of cardiorespiratory coupling. It is observed as changes in the heart rate in synchrony with the respiration. RSA has been hypothesized to be due to a combination of linear and nonlinear effects. The quantification of the latter, in turn, has been suggested as a biomarker to improve the assessment of several conditions and diseases. In this study, a framework to quantify RSA using support vector machines is presented. The methods are based on multivariate autoregressive models, in which the present samples of the heart rate variability are predicted as combinations of past samples of the respiration. The selection and tuning of a kernel in these models allows to solve the regression problem taking into account only the linear components, or both the linear and the nonlinear ones. The methods are tested in simulated data as well as in a dataset of polysomnographic studies taken from 110 obstructive sleep apnea patients. In the simulation, the methods were able to capture the nonlinear components when a weak cardiorespiratory coupling occurs. When the coupling increases, the nonlinear part of the coupling is not detected and the interaction is found to be of linear nature. The trends observed in the application in real data show that, in the studied dataset, the proposed methods captured a more prominent linear interaction than the nonlinear one.