Synthesis of HRV signals characterized by predetermined time-frequency structure by means of time-varying ARMA models

Model-Based Evaluation of Methods for Respiratory Sinus Arrhythmia Estimation

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 real-life applications, 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.

Quantifying the lagged Poincaré plot geometry of ultrashort heart rate variability series: automatic recognition of odor hedonic tone

The application of Poincaré plot analysis to characterize inter-beat interval dynamics has been successfully proposed in the scientific literature for the assessment of humans' physiological states and related aberrations. In this study, we proposed novel descriptors to trace the evolution of Poincaré plot shape over the lags. Their reliability in ultra-short cardiovascular series analysis was validated on synthetic inter-beat series generated through a physiologically plausible integral pulse frequency modulation model. Furthermore, we used the proposed approach for the investigation of the direct relationship between autonomic nervous system (ANS) dynamics and hedonic olfactory elicitation, in a group of 30 healthy subjects. Participants with a similar olfactory threshold were selected, and were asked to score 5-s stimuli in terms of arousal and valence levels according to the Russell's circumflex model of affect. Their ANS response was investigated in 35-s windows after the elicitation. Experimental results showed a gender-specific, high discriminant power of the proposed approach, discerning between pleasant and unpleasant odorants with an accuracy of 83.33% and 73.33% for men and for women, respectively. Graphical Abstract Olfaction plays a crucial role in our life and is strictly related to the Autonomic Nervous System (ANS) activity, which can be monitored studying Heart Rate Variability. We used the Lagged Poincare Plot approach to recognize gender-specific ANS response in 35-second windows after the elicitation through pleasant/unpleasant odorants.

The Effect of Emotional Valence on Ventricular Repolarization Dynamics Is Mediated by Heart Rate Variability: A Study of QT Variability and Music-Induced Emotions

Background

Emotions can affect cardiac activity, but their impact on ventricular repolarization variability, an important parameter providing information about cardiac risk and autonomic nervous system activity, is unknown. The beat-to-beat variability of the QT interval (QTV) from the body surface ECG is a non-invasive marker of repolarization variability, which can be decomposed into QTV related to RR variability (QTVrRRV) and QTV unrelated to RRV (QTVuRRV), with the latter thought to be a marker of intrinsic repolarization variability.

Aim

To determine the effect of emotional valence (pleasant and unpleasant) on repolarization variability in healthy volunteers by means of QTV analysis.

Methods

75 individuals (24.5 ± 3.2 years, 36 females) without a history of cardiovascular disease listened to music-excerpts that were either felt as pleasant (n = 6) or unpleasant (n = 6). Excerpts lasted about 90 s and were presented in a random order along with silent intervals (n = 6). QTV and RRV were derived from the ECG and the time-frequency spectrum of RRV, QTV, QTVuRRV and QTVrRRV as well as time-frequency coherence between QTV and RRV were estimated. Analysis was performed in low-frequency (LF), high frequency (HF) and total spectral bands.

Results

The heart rate-corrected QTV showed a small but significant increase from silence (median 347/interquartile range 31 ms) to listening to music felt as unpleasant (351/30 ms) and pleasant (355/32 ms). The dynamic response of QTV to emotional valence showed a transient phase lasting about 20 s after the onset of each musical excerpt. QTV and RRV were highly correlated in both HF and LF (mean coherence ranging 0.76–0.85). QTV and QTVrRRV decreased during listening to music felt as pleasant and unpleasant with respect to silence and further decreased during listening to music felt as pleasant. QTVuRRV was small and not affected by emotional valence.

Conclusion

Emotional valence, as evoked by music, has a small but significant effect on QTV and QTVrRRV, but not on QTVuRRV. This suggests that the interaction between emotional valence and ventricular repolarization variability is mediated by cycle length dynamics and not due to intrinsic repolarization variability.

Novel characterization method of impedance cardiography signals using time-frequency distributions

The purpose of this document is to describe a methodology to select the most adequate time-frequency distribution (TFD) kernel for the characterization of impedance cardiography signals (ICG). The predominant ICG beat was extracted from a patient and was synthetized using time-frequency variant Fourier approximations. These synthetized signals were used to optimize several TFD kernels according to a performance maximization. The optimized kernels were tested for noise resistance on a clinical database. The resulting optimized TFD kernels are presented with their performance calculated using newly proposed methods. The procedure explained in this work showcases a new method to select an appropriate kernel for ICG signals and compares the performance of different time-frequency kernels found in the literature for the case of ICG signals. We conclude that, for ICG signals, the performance (P) of the spectrogram with either Hanning or Hamming windows (P = 0.780) and the extended modified beta distribution (P = 0.765) provided similar results, higher than the rest of analyzed kernels. Graphical abstract Flowchart for the optimization of time-frequency distribution kernels for impedance cardiography signals.

Investigation of Lagged Poincaré Plot reliability in ultra-short synthetic and experimental Heart Rate Variability series

This study reports on the reliability of Lagged Poincaré Plot (LPP) parameters calculated from ultra-short cardiovascular time series (from 30 to 180 seconds). ity (HRV) signals, whereas a few studies have studied nonlinear approaches. Particularly, methods derived from the phase-space theory, especially the ones employing multi-lag analyses, are usually considered to be inaccurate with a low number of samples. Here we propose a comprehensive study about LPP, using both synthetic and real RR series. Specifically, we considered 109 5-minutes HRV series: 60 synthetic series generated through the Integral Pulse Frequency Modulation (IPFM) model and 49 experimental series acquired from healthy subjects during resting-state. Three parameters have been extracted through the ellipse-fitting method, SD1, SD2 and S, using ten values of lag. All LPP parameters were estimated by averaging estimates gathered from segments of 30, 120 and 180 seconds, and compared with the once from 5-minute series. Results showed Spearman's correlation coefficients higher than 0.9 in both synthetic and real series. In conclusion, SD1 gave promising results in terms of percentage absolute error, when it was extracted from series with a duration less than three minutes.

Assessment of Quadratic Nonlinear Cardiorespiratory Couplings During Tilt-Table Test by Means of Real Wavelet Biphase

OBJECTIVE

In this paper, a method for assessment of quadratic phase coupling (QPC) between respiration and heart rate variability (HRV) is presented.

METHODS

First, a method for QPC detection is proposed named real wavelet biphase (RWB). Then, a method for QPC quantification is proposed based on the normalized wavelet biamplitude (NWB). A simulation study has been conducted to test the reliability of RWB to identify QPC, even in the presence of constant delays between interacting oscillations, and to discriminate it from quadratic phase uncoupling. Significant QPC was assessed based on surrogate data analysis. Then, quadratic cardiorespiratory couplings were studied during a tilt-table test protocol of 17 young healthy subjects.

RESULTS

Simulation study showed that RWB is able to detect even weak QPC with delays in the range of [Formula: see text] s, which are usual in the autonomic nervous system (ANS) control of heart rate. Results from the database revealed a significant reduction ([Formula: see text]) of NWB between respiration and both low and high frequencies of HRV in head-up tilt position compared to early supine.

CONCLUSION

The proposed technique detects and quantifies robustly QPC and is able to track the coupling between respiration and various HRV components during ANS changes.

SIGNIFICANCE

The proposed method can help to assess alternations of nonlinear cardiorespiratory interactions related to ANS dysfunction and physiological regulation of HRV in cardiovascular diseases.

A Time-Varying Nonparametric Methodology for Assessing Changes in QT Variability Unrelated to Heart Rate Variability

OBJECTIVE

To propose and test a novel methodology to measure changes in QT interval variability (QTV) unrelated to RR interval variability (RRV) in nonstationary conditions.

METHODS

Time-frequency coherent and residual spectra representing QTV related (QTVrRRV) and unrelated (QTVuRRV) to RRV, respectively, are estimated using time-frequency Cohen's class distributions. The proposed approach decomposes the nonstationary output spectrum of any two-input one-output model with uncorrelated inputs into two spectra representing the information related and unrelated to one of the two inputs, respectively. An algorithm to correct for the bias of the time-frequency coherence function between QTV and RRV is proposed to provide accurate estimates of both QTVuRRV and QTVrRRV. Two simulation studies were conducted to assess the methodology in challenging nonstationary conditions and data recorded during head-up tilt in 16 healthy volunteers were analyzed.

RESULTS

In the simulation studies, QTVuRRV changes were tracked with only a minor delay due to the filtering necessary to estimate the nonstationary spectra. The correlation coefficient between theoretical and estimated patterns was even for extremely noisy recordings (signal to noise ratio (SNR) in QTV dB). During head-up tilt, QTVrRRV explained the largest proportion of QTV, whereas QTVuRRV showed higher relative increase than QTV or QTVrRRV in all spectral bands ( for most pairwise comparisons).

CONCLUSION

The proposed approach accurately tracks changes in QTVuRRV. Head-up tilt induced a slightly greater increase in QTVuRRV than in QTVrRRV.

SIGNIFICANCE

The proposed index QTVuRRV may represent an indirect measure of intrinsic ventricular repolarization variability, a marker of cardiac instability associated with sympathetic ventricular modulation and sudden cardiac death.

Reliability of Lagged Poincaré Plot Parameters in Ultrashort Heart Rate Variability Series: Application on Affective Sounds

The number of studies about ultrashort cardiovascular time series is increasing because of the demand for mobile applications in telemedicine and e-health monitoring. However, the current literature still needs a proper validation of heartbeat nonlinear dynamics assessment from ultrashort time series. This paper reports on the reliability of the Lagged Poincaré Plot (LPP) parameters-calculated from ultrashort cardiovascular time series. Reliability is studied on simulated as well as on real RR series. Simulated RR series are generated and LPP parameters estimated for ultrashort time series (from 15 to 60 s) are compared to those estimated from 1 h. All LPP parameters estimated from time series longer than 35 s presented a Spearman's correlation coefficient higher than 0.99. RR series acquired from 32 healthy subjects during 5-min resting state sessions are used to test the LPP approach in experimental data. The usefulness of ultrashort term parameters in real data is accomplished also studying their ability to discriminate positive and negative valence of auditory stimuli taken from the International Affective Digitized Sound System (IADS) dataset. The achieved accuracies in the recognition of elicitation along the valence dimension, using only the LPP parameters, were of 77.78% for 1 min 28 s series, and of 79.17% for 35 s series.

A multivariate time-frequency approach for tracking QT variability changes unrelated to heart rate variability

The beat-to-beat variability of the QT interval (QTV) is a marker of ventricular repolarization (VR) dynamics and it has been suggested as an index of sympathetic ventricular outflow and cardiac instability. However, QTV is also affected by RR (or heart rate) variability (RRV), and QTV due to RRV may reduce QTV specificity as a VR marker. Therefore, it would be desirable to separate QTV due to VR dynamics from QTV due to RRV. To do that, previous work has mainly focused on heart rate corrections or time-invariant autoregressive models. This paper describes a novel framework that extends classical multiple inputs/single output theory to the time-frequency (TF) domain to quantify QTV and RRV interactions. Quadratic TF distributions and TF coherence function are utilized to separate QTV into two partial (conditioned) spectra representing QTV related and unrelated to RRV, and to provide an estimates of intrinsic VR dynamics. In a simulation study, a time-varying ARMA model was used to generate signals representing realistic RRV and VR dynamics with controlled instantaneous frequencies and powers. The results demonstrated that the proposed methodology is able to accurately track changes in VR dynamics, with a correlation between theoretical and estimated patterns higher than 0.88. Data from healthy volunteers undergoing a tilt table test were analyzed and representative examples are discussed. Results show that the QTV unrelated to RRV dynamics quickly increased during orthostatic challenge.

Influence of Heart Rate in Non-linear HRV Indices as a Sampling Rate Effect Evaluated on Supine and Standing

The purpose of this study is to characterize and attenuate the influence of mean heart rate (HR) on nonlinear heart rate variability (HRV) indices (correlation dimension, sample, and approximate entropy) as a consequence of being the HR the intrinsic sampling rate of HRV signal. This influence can notably alter nonlinear HRV indices and lead to biased information regarding autonomic nervous system (ANS) modulation. First, a simulation study was carried out to characterize the dependence of nonlinear HRV indices on HR assuming similar ANS modulation. Second, two HR-correction approaches were proposed: one based on regression formulas and another one based on interpolating RR time series. Finally, standard and HR-corrected HRV indices were studied in a body position change database. The simulation study showed the HR-dependence of non-linear indices as a sampling rate effect, as well as the ability of the proposed HR-corrections to attenuate mean HR influence. Analysis in a body position changes database shows that correlation dimension was reduced around 21% in median values in standing with respect to supine position (p < 0.05), concomitant with a 28% increase in mean HR (p < 0.05). After HR-correction, correlation dimension decreased around 18% in standing with respect to supine position, being the decrease still significant. Sample and approximate entropy showed similar trends. HR-corrected nonlinear HRV indices could represent an improvement in their applicability as markers of ANS modulation when mean HR changes.

Inclusion of Respiratory Frequency Information in Heart Rate Variability Analysis for Stress Assessment

Respiratory rate and heart rate variability (HRV) are studied as stress markers in a database of young healthy volunteers subjected to acute emotional stress, induced by a modification of the Trier Social Stress Test. First, instantaneous frequency domain HRV parameters are computed using time-frequency analysis in the classical bands. Then, the respiratory rate is estimated and this information is included in HRV analysis in two ways: 1) redefining the high-frequency (HF) band to be centered at respiratory frequency; 2) excluding from the analysis those instants where respiratory frequency falls within the low-frequency (LF) band. Classical frequency domain HRV indices scarcely show statistical differences during stress. However, when including respiratory frequency information in HRV analysis, the normalized LF power as well as the LF/HF ratio significantly increase during stress ( p-value 0.05 according to the Wilcoxon test), revealing higher sympathetic dominance. The LF power increases during stress, only being significantly different in a stress anticipation stage, while the HF power decreases during stress, only being significantly different during the stress task demanding attention. Our results support that joint analysis of respiration and HRV obtains a more reliable characterization of autonomic nervous response to stress. In addition, the respiratory rate is observed to be higher and less stable during stress than during relax ( p-value 0.05 according to the Wilcoxon test) being the most discriminative index for stress stratification (AUC = 88.2 % ).

A reweighted ℓ1-minimization based compressed sensing for the spectral estimation of heart rate variability using the unevenly sampled data

In this paper, a reweighted ℓ₁-minimization based Compressed Sensing (CS) algorithm incorporating the Integral Pulse Frequency Modulation (IPFM) model for spectral estimation of HRV is introduced. Knowing as a novel sensing/sampling paradigm, the theory of CS asserts certain signals that are considered sparse or compressible can be possibly reconstructed from substantially fewer measurements than those required by traditional methods. Our study aims to employ a novel reweighted ℓ₁-minimization CS method for deriving the spectrum of the modulating signal of IPFM model from incomplete RR measurements for HRV assessments. To evaluate the performance of HRV spectral estimation, a quantitative measure, referred to as the Percent Error Power (PEP) that measures the percentage of difference between the true spectrum and the spectrum derived from the incomplete RR dataset, was used. We studied the performance of spectral reconstruction from incomplete simulated and real HRV signals by experimentally truncating a number of RR data accordingly in the top portion, in the bottom portion, and in a random order from the original RR column vector. As a result, for up to 20% data truncation/loss the proposed reweighted ℓ₁-minimization CS method produced, on average, 2.34%, 2.27%, and 4.55% PEP in the top, bottom, and random data-truncation cases, respectively, on Autoregressive (AR) model derived simulated HRV signals. Similarly, for up to 20% data loss the proposed method produced 5.15%, 4.33%, and 0.39% PEP in the top, bottom, and random data-truncation cases, respectively, on a real HRV database drawn from PhysioNet. Moreover, results generated by a number of intensive numerical experiments all indicated that the reweighted ℓ₁-minimization CS method always achieved the most accurate and high-fidelity HRV spectral estimates in every aspect, compared with the ℓ₁-minimization based method and Lomb's method used for estimating the spectrum of HRV from unevenly sampled RR data.

Cardiorespiratory dynamic response to mental stress: a multivariate time-frequency analysis

Mental stress is a growing problem in our society. In order to deal with this, it is important to understand the underlying stress mechanisms. In this study, we aim to determine how the cardiorespiratory interactions are affected by mental arithmetic stress and attention. We conduct cross time-frequency (TF) analyses to assess the cardiorespiratory coupling. In addition, we introduce partial TF spectra to separate variations in the RR interval series that are linearly related to respiration from RR interval variations (RRV) that are not related to respiration. The performance of partial spectra is evaluated in two simulation studies. Time-varying parameters, such as instantaneous powers and frequencies, are derived from the computed spectra. Statistical analysis is carried out continuously in time to evaluate the dynamic response to mental stress and attention. The results show an increased heart and respiratory rate during stress and attention, compared to a resting condition. Also a fast reduction in vagal activity is noted. The partial TF analysis reveals a faster reduction of RRV power related to (3 s) than unrelated to (30 s) respiration, demonstrating that the autonomic response to mental stress is driven by mechanisms characterized by different temporal scales.