Application of an automatic adaptive filter for Heart Rate Variability analysis

Unravelling COVID-19 waves in Rio de Janeiro city: Qualitative insights from nonlinear dynamic analysis

Since the COVID-19 pandemic was first reported in 2019, it has rapidly spread around the world. Many countries implemented several measures to try to control the virus spreading. The healthcare system and consequently the general quality of life population in the cities have all been significantly impacted by the Coronavirus pandemic. The different waves of contagious were responsible for the increase in the number of cases that, unfortunately, many times lead to death. In this paper, we aim to characterize the dynamics of the six waves of cases and deaths caused by COVID-19 in Rio de Janeiro city using techniques such as the Poincaré plot, approximate entropy, second-order difference plot, and central tendency measures. Our results reveal that by examining the structure and patterns of the time series, using a set of non-linear techniques we can gain a better understanding of the role of multiple waves of COVID-19, also, we can identify underlying dynamics of disease spreading and extract meaningful information about the dynamical behavior of epidemiological time series. Such findings can help to closely approximate the dynamics of virus spread and obtain a correlation between the different stages of the disease, allowing us to identify and categorize the stages due to different virus variants that are reflected in the time series.

Estimation of physiological exercise thresholds based on dynamical correlation properties of heart rate variability

Aerobic and anaerobic thresholds of the three-zone exercise model are often used to evaluate the exercise intensity and optimize the training load. Conventionally, these thresholds are derived from the respiratory gas exchange or blood lactate concentration measurements. Here, we introduce and validate a computational method based on the RR interval (RRI) dynamics of the heart rate (HR) measurement, which enables a simple, yet reasonably accurate estimation of both metabolic thresholds. The method utilizes a newly developed dynamical detrended fluctuation analysis (DDFA) to assess the real-time changes in the dynamical correlations of the RR intervals during exercise. The training intensity is shown to be in direct correspondence with the time- and scale-dependent changes in the DDFA scaling exponent. These changes are further used in the definition of an individual measure to estimate the aerobic and anaerobic threshold. The results for 15 volunteers who participated in a cyclo-ergometer test are compared to the benchmark lactate thresholds, as well as to the ventilatory threshods and alternative HR-based estimates based on the maximal HR and the conventional detrended fluctuation analysis (DFA). Our method provides the best overall agreement with the lactate thresholds and provides a promising, cost-effective alternative to conventional protocols, which could be easily integrated in wearable devices. However, detailed statistical analysis reveals the particular strengths and weaknessess of each method with respect to the agreement and consistency with the thresholds-thus underlining the need for further studies with more data.

A two-step pre-processing tool to remove Gaussian and ectopic noise for heart rate variability analysis

Artifacts in the Electrocardiogram (ECG) degrade the quality of the recorded signal and are not conducive to heart rate variability (HRV) analysis. The two types of noise most often found in ECG recordings are technical and physiological artifacts. Current preprocessing methods primarily attend to ectopic beats but do not consider technical issues that affect the ECG. A secondary aim of this study was to investigate the effect of increasing increments of artifacts on 24 of the most used HRV measures. A two-step preprocessing approach for denoising HRV is introduced which targets each type of noise separately. First, the technical artifacts in the ECG are eliminated by applying complete ensemble empirical mode decomposition with adaptive noise. The second step removes physiological artifacts from the HRV signal using a combination filter of single dependent rank order mean and an adaptive filtering algorithm. The performance of the two-step pre-processing tool showed a high correlation coefficient of 0.846 and RMSE value of 7.69 × 10^-5 for 6% of added ectopic beats and 6 dB Gaussian noise. All HRV measures studied except HF peak and LF peak are significantly affected by both types of noise. Frequency measures of Total power, HF power, and LF power and fragmentation measures; PAS, PIP, and PSS are the most sensitive to both types of noise.

Characterisation of neonatal cardiac dynamics using ordinal partition network

The maturation of the autonomic nervous system (ANS) starts in the gestation period and it is completed after birth in a variable time, reaching its peak in adulthood. However, the development of ANS maturation is not entirely understood in newborns. Clinically, the ANS condition is evaluated with monitoring of gestational age, Apgar score, heart rate, and by quantification of heart rate variability using linear methods. Few researchers have addressed this problem from the perspective nonlinear data analysis. This paper proposes a new data-driven methodology using nonlinear time series analysis, based on complex networks, to classify ANS conditions in newborns. We map 74 time series given by RR intervals from premature and full-term newborns to ordinal partition networks and use complexity quantifiers to discriminate the dynamical process present in both conditions. We obtain three complexity quantifiers (permutation, conditional, and global node entropies) using network mappings from forward and reverse directions, and considering different time lags and embedding dimensions. The results indicate that time asymmetry is present in the data of both groups and the complexity quantifiers can differentiate the groups analysed. We show that the conditional and global node entropies are sensitive for detecting subtle differences between the neonates, particularly for small embedding dimensions (m < 7). This study reinforces the assessment of nonlinear techniques for RR interval time series analysis. Graphical Abstract.

Comparing stress prediction models using smartwatch physiological signals and participant self-reports

Recent advances in wearable technology have facilitated the non-obtrusive monitoring of physiological signals, creating opportunities to monitor and predict stress. Researchers have utilized machine learning methods using these physiological signals to develop stress prediction models. Many of these prediction models have utilized objective stressor tasks (e.g., a public speaking task or solving math problems). Alternatively, the subjective user responses with self-reports have also been used for measuring stress. In this paper, we describe a methodological approach (a) to compare the prediction performance of models developed using objective markers of stress using participant-reported subjective markers of stress from self-reports; and (b) to develop personalized stress models by accounting for inter-individual differences. Towards this end, we conducted a laboratory-based study with 32 healthy volunteers. Participants completed a series of stressor tasks-social, cognitive and physical-wearing an instrumented commercial smartwatch that collected physiological signals and participant responses using timed self-reports. After extensive data preprocessing using a combination of signal processing techniques, we developed two types of models: objective stress models using the stressor tasks as labels; and subjective stress models using participant responses to each task as the label for that stress task. We trained and tested several machine learning algorithms-support vector machine (SVM), random forest (RF), gradient boosted trees (GBT), AdaBoost, and Logistic Regression (LR)-and evaluated their performance. SVM had the best performance for the models using the objective stressor (i.e., stressor tasks) with an AUROC of 0.790 and an F-1 score of 0.623. SVM also had the highest performance for the models using the subjective stress (i.e., participant self-reports) with an AUROC of 0.719 and an F-1 score of 0.520. Model performance improved with a personalized threshold model to an AUROC of 0.751 and an F-1 score of 0.599. The performance of the stress models using an instrumented commercial smartwatch was comparable to similar models from other state-of-the-art laboratory-based studies. However, the subjective stress models had a lower performance, indicating the need for further research on the use of self-reports for stress-related studies. The improvement in performance with the personalized threshold-based models provide new directions for building stress prediction models.

Heart Rate Variability Assessment and Its Application for Autonomic Function Evaluation in Patients with Hyperhidrosis

BACKGROUND

Pathophysiology mechanism of primary focal hyperhidrosis (PFHH) is controversial. Heart rate variability (HRV) could explain if there is a systemic component present. We aimed to investigate the functions of the autonomic nervous system in patients diagnosed with PFHH compared to controls using the analysis of HRV in the domains of time, frequency, and nonlinearity, as well as analysis of the recurrence plots (RPs).

METHODS

We selected 34 patients with PFHH (29.4 ± 10.2 years) and 34 controls (29.2 ± 9.6 years) for HRV analysis. Heart beats were recorded with Polar RS800CX monitor (20 min, at rest, in supine position), and RR intervals were analyzed with Kubios Premium HRV software. RPs were constructed with Visual Recurrence Analysis software. Statistical analysis included unpaired t test (p < 0.05).

RESULTS

Our results showed that HRV parameters in the 3 domains evaluated did not show any differences between the groups. The same was observed with RPs.

CONCLUSIONS

The findings suggest that PFHH, from the pathophysiological point of view, may be caused by peripheral involvement of the sympathetic nervous system (glandular level or nerve terminals), as there was no difference between the groups studied. More specific studies should help elucidate this issue.

Adaptive Filtering Improved Apnea Detection Performance Using Tracheal Sounds in Noisy Environment: A Simulation Study

Objective

Tracheal sounds were used to detect apnea on various occasions. However, ambient noises can contaminate tracheal sounds which result in poor performance of apnea detection. The objective of this paper was to apply the adaptive filtering (AF) algorithm to improve the quality of tracheal sounds and examine the accuracy of the apnea detection algorithm using tracheal sounds after AF.

Method

Tracheal sounds were acquired using a primary microphone encased in a plastic bell, and the ambient noises were collected using a reference microphone resting outside the plastic bell in quiet and noisy environments, respectively. Simultaneously, the flow pressure signals and thoracic and abdominal movement were obtained as the standard signals to determine apnea events. Then, the normalized least mean square (NLMS) AF algorithm was applied to the tracheal sounds mixed with noises. Finally, the algorithm of apnea detection was used to the tracheal sounds with AF and the tracheal sounds without AF. Sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), accuracy, and Cohen's kappa coefficient of apnea detection were calculated.

Results

Forty-six healthy subjects, aged 18-35 years and with BMI < 21.4, were included in the study. The apnea detection performance using tracheal sounds was as follows: in the quiet environment, the tracheal sounds without AF detected apnea with 97.2% sensitivity, 99.9% specificity, 99.8% PPV, 99.4% NPV, 99.5% accuracy, and 0.982 kappa coefficient. The tracheal sounds with AF detected apnea with 98.2% sensitivity, 99.9% specificity, 99.4% PPV, 99.6% NPV, 99.6% accuracy, and 0.985 kappa coefficient. While in the noisy environment, the tracheal sounds without AF detected apnea with 81.1% sensitivity, 96.9% specificity, 85.1% PPV, 96% NPV, 94.2% accuracy, and 0.795 kappa coefficient and the tracheal sounds with AF detected apnea with 91.5% sensitivity, 97.4% specificity, 88.4% PPV, 98.2% NPV, 96.4% accuracy, and 0.877 kappa coefficient.

Conclusion

The performance of apnea detection using tracheal sounds with the NLMS AF algorithm in the noisy environment proved to be accurate and reliable. The AF technology could be applied to the respiratory monitoring using tracheal sounds.

Extended Central Tendency Measure and difference plot for heart rate variability analysis

Heart rate variability (HRV) is a non-invasive alternative to analyze the role of the autonomic nervous system (ANS) on heart functioning. Many tools have been developed to analyze collected cardiac data. Among them, the Central Tendency Measure (CTM) is a quantitative method for variability analysis of RR intervals. The values of the CTM must be between 0 and 1 (inclusive) for different radius, which follows the intrinsic characteristics of each time series. Using the conventional CTM, the successive differences of the time series may be calculated, and it can classify and differentiate the disturbances in the ANS involving HRV. This method was extended (e-CTM) to analyze the differences between RR interval time series. In this extension, a new parameter is added, which allows analysis of long time intervals, instead of successive and adjacent RR intervals. The ability of the e-CTM to differentiate the groups of the RR interval time series was verified with 145 RR interval time series divided into three groups: subjects with congestive heart failure, healthy subjects, and nurses during one hour of their workday. Results evidence that the new parameter added differentiates the group with pathology (and subsequent impairment of ANS) and group under stress at work (temporary impairment of ANS). These results suggest that the e-CTM is capable of detection long-term variations in the HRV according to the ANS impairment.

A preprocessing tool for removing artifact from cardiac RR interval recordings using three-dimensional spatial distribution mapping

Artifact is common in cardiac RR interval data that is recorded for heart rate variability (HRV) analysis. A novel algorithm for artifact detection and interpolation in RR interval data is described. It is based on spatial distribution mapping of RR interval magnitude and relationships to adjacent values in three dimensions. The characteristics of normal physiological RR intervals and artifact intervals were established using 24-h recordings from 20 technician-assessed human cardiac recordings. The algorithm was incorporated into a preprocessing tool and validated using 30 artificial RR (ARR) interval data files, to which known quantities of artifact (0.5%, 1%, 2%, 3%, 5%, 7%, 10%) were added. The impact of preprocessing ARR files with 1% added artifact was also assessed using 10 time domain and frequency domain HRV metrics. The preprocessing tool was also used to preprocess 69 24-h human cardiac recordings. The tool was able to remove artifact from technician-assessed human cardiac recordings (sensitivity 0.84, SD = 0.09, specificity of 1.00, SD = 0.01) and artificial data files. The removal of artifact had a low impact on time domain and frequency domain HRV metrics (ranging from 0% to 2.5% change in values). This novel preprocessing tool can be used with human 24-h cardiac recordings to remove artifact while minimally affecting physiological data and therefore having a low impact on HRV measures of that data.

Cardiovascular and nervous system changes during meditation

Background: A number of benefits have been described for the long-term practice of meditation, yet little is known regarding the immediate neurological and cardiovascular responses to meditation. Wireless sensor technology allows, for the first time, multi-parameter and quantitative monitoring of an individual's responses during meditation. The present study examined inter-individual variations to meditation through continuous monitoring of EEG, blood pressure, heart rate and its variability (HRV) in novice and experienced meditators.

Methods: Participants were 20 experienced and 20 novice meditators involved in a week-long wellness retreat. Monitoring took place during meditation sessions on the first and last full days of the retreat. All participants wore a patch that continuously streamed ECG data, while half of them also wore a wireless EEG headset plus a non-invasive continuous blood pressure monitor.

Results: Meditation produced variable but characteristic EEG changes, significantly different from baseline, even among novice meditators on the first day. In addition, although participants were predominately normotensive, the mean arterial blood pressure fell a small (2–3 mmHg) but significant (p < 0.0001) amount during meditation. The effect of meditation on HRV was less clear and influenced by calculation technique and respiration. No clear relationship between EEG changes, HRV alterations, or mean blood pressure during meditation was found.

Conclusion: This is the first study to investigate neurological and cardiovascular responses during meditation in both novice and experienced meditators using novel, wearable, wireless devices. Meditation produced varied inter-individual physiologic responses. These results support the need for further investigation of the short- and long-term cardiovascular effects of mental calm and individualized ways to achieve it.

A fixed mass method for the Kramers-Moyal expansion--application to time series with outliers

Extraction of stochastic and deterministic components from empirical data-necessary for the reconstruction of the dynamics of the system-is discussed. We determine both components using the Kramers-Moyal expansion. In our earlier papers, we obtained large fluctuations in the magnitude of both terms for rare or extreme valued events in the data. Calculations for such events are burdened by an unsatisfactory quality of the statistics. In general, the method is sensitive to the binning procedure applied for the construction of histograms. Instead of the commonly used constant width of bins, we use here a constant number of counts for each bin. This approach-the fixed mass method-allows to include in the calculation events, which do not yield satisfactory statistics in the fixed bin width method. The method developed is general. To demonstrate its properties, here, we present the modified Kramers-Moyal expansion method and discuss its properties by the application of the fixed mass method to four representative heart rate variability recordings with different numbers of ectopic beats. These beats may be rare events as well as outlying, i.e., very small or very large heart cycle lengths. The properties of ectopic beats are important not only for medical diagnostic purposes but the occurrence of ectopic beats is a general example of the kind of variability that occurs in a signal with outliers. To show that the method is general, we also present results for two examples of data from very different areas of science: daily temperatures at a large European city and recordings of traffics on a highway. Using the fixed mass method, to assess the dynamics leading to the outlying events we studied the occurrence of higher order terms of the Kramers-Moyal expansion in the recordings. We found that the higher order terms of the Kramers-Moyal expansion are negligible for heart rate variability. This finding opens the possibility of the application of the Langevin equation to the whole range of empirical signals containing rare or outlying events. Note, however, that the higher order terms are non-negligible for the other data studied here and for it the Langevin equation is not applicable as a model.