Next Article in Journal
Applying Convolutional Neural Networks to Predict the ICD-9 Codes of Medical Records
Next Article in Special Issue
Development of a Self-Powered Piezo-Resistive Smart Insole Equipped with Low-Power BLE Connectivity for Remote Gait Monitoring
Previous Article in Journal
Intelligent Perception System of Robot Visual Servo for Complex Industrial Environment
Previous Article in Special Issue
Leveraging Wearable Sensors for Human Daily Activity Recognition with Stacked Denoising Autoencoders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Error Estimation of Ultra-Short Heart Rate Variability Parameters: Effect of Missing Data Caused by Motion Artifacts

1
Department of Computer Science, University of Pisa, 56126 Pisa, Italy
2
Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford OX1 2JD, UK
3
Huma Therapeutics Limited, London SW1P 4QP, UK
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(24), 7122; https://doi.org/10.3390/s20247122
Submission received: 19 November 2020 / Revised: 9 December 2020 / Accepted: 10 December 2020 / Published: 11 December 2020

Abstract

:
Application of ultra–short Heart Rate Variability (HRV) is desirable in order to increase the applicability of HRV features to wrist-worn wearable devices equipped with heart rate sensors that are nowadays becoming more and more popular in people’s daily life. This study is focused in particular on the the two most used HRV parameters, i.e., the standard deviation of inter-beat intervals (SDNN) and the root Mean Squared error of successive inter-beat intervals differences (rMSSD). The huge problem of extracting these HRV parameters from wrist-worn devices is that their data are affected by the motion artifacts. For this reason, estimating the error caused by this huge quantity of missing values is fundamental to obtain reliable HRV parameters from these devices. To this aim, we simulate missing values induced by motion artifacts (from 0 to 70%) in an ultra-short time window (i.e., from 4 min to 30 s) by the random walk Gilbert burst model in 22 young healthy subjects. In addition, 30 s and 2 min ultra-short time windows are required to estimate rMSSD and SDNN, respectively. Moreover, due to the fact that ultra-short time window does not permit assessing very low frequencies, and the SDNN is highly affected by these frequencies, the bias for estimating SDNN continues to increase as the time window length decreases. On the contrary, a small error is detected in rMSSD up to 30 s due to the fact that it is highly affected by high frequencies which are possible to be evaluated even if the time window length decreases. Finally, the missing values have a small effect on rMSSD and SDNN estimation. As a matter of fact, the HRV parameter errors increase slightly as the percentage of missing values increase.

1. Introduction

Heart Rate Variability (HRV) is a physiological phenomenon of fluctuation between adjacent Inter-Beat Intervals (IBI) [1,2]. The heart beats (HR) average and HRV are regulated by two main mechanisms. First, they are controlled by Autonomic Nervous Systems (ANS), which affect the activity of the sinoatrial node (SAN) [3,4,5]. In particular, ANS controls SAN by two signaling pathways that permit increasing a decrease in HR by sympathetic and parasympathetic nervous systems stimulations, respectively [2]. The second mechanism that regulates HR and HRV is the coupled-clock system which drives the automaticity of human SAN pacemaker cells even without neural input [3,6]. Aging and clinical condition (e.g., cardiovascular disease) affect both ANS and SAN responses inducing an alteration in HR and HRV [2,3]. For this reason, monitoring HR and HRV response permits assessing possible deterioration of cardiovascular responses. In particular, abnormal hourly heart rate patterns associated with an increase in HRV values have been found to be strongly correlated to an increased risk of mortality (particularly among the elderly) causing, for example, atrial fibrillation [7]. Moreover, HRV parameters are useful parameters to provide insight about sympathetic and parasympathetic nervous system responses of cardiac vagal tone that was found to be linked to cognitive, emotional, social, and health status [8].
SDNN, i.e., standard deviation of the normal-to-normal interval between QRS complexes, rMSSD, i.e., the root mean squared of successive inter-beat intervals, are the two widely used HRV parameters to predict many health issues [2]. SDNN is considered the “gold standard” for cardiovascular risks when recorded over 24 h (SDNN24) [9]. As a matter of fact, it has been found that SDNN24 values lower than 50 milliseconds (ms), between 50 and 100 ms and higher than 100 ms, indicate people with unhealthy, compromised health and healthy status, respectively [10]. In particular, Keiger et al. [10] demonstrated that people with SDNN24 higher than 100 ms have a 5.3 times lower risk of mortality at follow-up compared to those with SDNN24 lower than 50 ms. Moreover, it was found that SDNN24 values lower than 80 ms are predictive of cardiac events [11,12]. On the contrary, RMSSD reflects the IBI variance that is sensitive to high-frequency (HF) heart period fluctuations such as in the respiratory frequency range and is used as an index of vagal cardiac control [13]. It is considered the most useful marker for the parasympathetic neural system activation [9,14]. RMSSD is widely used by professional athletes in order to monitor cardiac activity and consequently their health status. In particular, it has been found that acute decreases in HRV have been linked to intense endurance training [15], resistance training [16], and competition [17]. Therefore, low HRV reflects acute fatigue from high intensity physical effort. Moreover, reduced resting HF is also linked to altered activation of prefrontal cortex and consequently to low cognitive functions that were found to be associated with stress, panic, anxiety, or worry [18]. Actually, in the last few decades, researchers investigated the association between HRV and depression reporting a strong association between reduced rMSSD and depression [19,20].
HRV has been traditionally obtained from 5 min or 24 h IBI time series [2]. Thanks to the recent technological advancements of wrist-worn wearable devices equipped with a heart rate sensor, it is now possible to semi-continuously obtain passive cardiac response measurements during people’s daily lives. For this reason and because of the low cost of these devices, they might have a great impact on the medical field due to the fact that it is possible to assess HRV features 24 h per day, 7 days per week [21]. The main problem of wrist-worn wearable devices is the artifacts induced by external stimuli that produce inconsistent IBI values. Usually, when IBI are recorded with gold standard technology (i.e., electrocardiography), the number of abnormal values is close to 1% [22], while, when it is recorded by these low cost instruments, is more than 10% [23]. Hence, estimating the error caused by this huge quantity of missing values is fundamental to obtain reliable values from wrist-worn wearable devices equipped with heart rate sensors avoiding misleading results [24]. To obtain reliable HRV data from these devices, shorter recordings (lower than 5 min) of IBI data during real-life scenarios are required. The ultra-short HRV parameters permit reducing the influence of motion artifacts on data recording due to the low time window length and consequently low percentage of missing values. Several studies estimate the reliability of HRV analysis from IBI time series shorter than 5 min for the clinical purposes demonstrating that 10 and 30 s recordings are the minimum time required to obtain accurate measures of RMSSD and SDNN, respectively [25,26,27,28]. However, to the best of our knowledge, no study are conducted in order to estimate the effect of missing values induced by wrist-worn devices equipped with HR sensors on ultra-short HRV parameters.
Hence, the aim of this study is to quantify the error induced by ultra-short IBI time series (i.e., from 4 min to 30 s) and missing values (from 0 to 70%) in order to detect if it is possible to accurately estimate SDNN and rMSSD by using wrist-worn wearable devices. To this aim, we simulate missing values induced by motion artifacts by a Gilbert burst model which creates a burst-error with a two-state Markov chain in 22 young healthy subjects.

2. Materials and Methods

2.1. Data

In this paper, we used the Multilevel Monitoring of Activity and Sleep in Healthy people (MMASH) dataset [29,30,31] providing 24 h of continuous Inter-Beat Intervals data (IBI), triaxial accelerometer data, sleep quality, physical activity, and psychological characteristics (i.e., anxiety status, stress events, and emotions) for 22 healthy young males (age = 27.29 ± 4.21 yrs; height = 179.91 ± 8.22 cm; weight = 75.05 ± 12.79 kg). Participants’ anthropomorphic characteristics (i.e., age, height, and weight) were recorded at the start of the data recording. Moreover, participants filled in questionnaires that provide information about participants’ psychological status, i.e., chronotype, anxiety status, and sleep quality. During the 24 h of the data recording, participants wore two devices that continuously recorded heart response (Polar H7 heart rate monitor—Polar Electro Inc., Bethpage, NY, USA) and the actigraphy data (ActiGraph wGT3X-BT—ActiGraph LLC, Pensacola, FL, USA). Moreover, the perceived moods were recorded at different times of the day and the daily stress was provided before sleeping in order to summarize the individual’s stressful events of the day. Finally, twice a day (i.e., before going to bed and when they woke up), the subjects collected saliva samples at home in appropriate vials in order to assess Melatonin and Corisol saliva concentration. More details about the experimental setup of the MMASH dataset are provided in the data descriptor paper written by Rossi et al. and published on MDPI Data journal in 2020 [29].
In this study, we used only IBI data continuously recorded over 24 h on 22 healthy young males by a Polar H7 chest strap (Polar Electro Inc., Bethpage, NY, USA). This device is a Bluetooth low energy chest strap with an ECG sensor that provides valid information of inter-beat intervals (IBI) [32]. Participants wore the Polar H7 chest strap for 24 h (between 9:00 a.m. and 9:00 p.m. on the next day), and were instructed to wear the chest straps during both day (during physical activities too) and night. For the aim of this study, the IBI timeseries of each participants were split into 5 min time windows in order to compute gold standard HRV parameters. Moreover, to assess the effect of length of the time window on HRV estimation, ultra-short time windows (from 4 min to 30 s) were selected in each 5 min time window, and the HRV parameters were computed. Finally, we introduce missing values to assess their effect on HRV features. More details about the dataset preprocessing and the HRV estimation were provided in Section 2.2 and Section 2.3, respectively. Finally, the data analysis details are provided in Section 2.4.
In accordance with the Helsinki Declaration as revised in 2013, the study was approved by the Ethical Committee of the University of Pisa ( # 0077455 / 2018 ).

2.2. Data Preprocessing

First of all, by using the Python hrv-analysis library (https://pypi.org/project/hrv-analysis), 2.18 ± 1.29% of RR-intervals in the dataset were detected as ectopic beats (i.e., disturbance of the cardiac rhythm frequently related to the electrical conduction system of the heart) or missing values induced by motion artifacts. These missing values were reconstructed via quadratic interpolation applied on the time domain, i.e., the heartbeats timestamp, instead of the duration domain, i.e., the duration of the heartbeats, as suggested by Morelli et al. [23].
In order to compute gold standard SDNN and rMSSD HRV features, the participants’ 24 h RR-intervals time series were split into 5 min windows (6130 5 min time windows were extracted). Moreover, to detect the influence of ultra-short HRV features in each 5 min window, ultra-short windows (4 min, 3 min, 2 min, 1 min and 30 s) were randomly selected.
Finally, in order to simulate motion artifacts observable during IBI recording by wrist-worn devices equipped with PPG sensors, artificial missing values (i.e., 5%, 10%, 15%, 30%, 50%, and 70%) were created for each 5 min and ultra-short time windows in accordance with a random walk Gilbert burst model which simulates burst-error with a two-state Markov chain (i.e., good as 0 and bed as 1) [33] as showed in Morelli at al.’s paper [23].

2.3. HRV Parameters

HRV is usually calculated from 5 min IBI time series [34]. For the aim of our study, we focused on the two most popular time-domain HRV measurements:
  • SDNN refers to the standard deviation of IBI. It estimates overall power spectrum of IBI timeseries. The SDNN is defined as the “gold standard” to assess both morbidity and mortality in the population [2].
  • rMSSD is the root mean square of the successive IBI differences estimates short-term components of HRV [35].
These two HRV metrics are computed in each 5 min window and in each ultra-short window with and without missing values.
Moreover, to better investigate the effect of missing values in the ultra-short HRV parameters, the power spectrum of the IBI time series was computed. The power spectrum S x x ( f ) of a time series x ( t ) describes the distribution of power into frequency components ( ω ) composing that signal. The minimal frequency that is possible to observe for each ultra-short time series was defined as a ratio between 1 and its length expressed in seconds (1/s), while the maximum is set at 0.4 in accordance with the literature. In order to derive the Power Spectrum Density (PSD) of the time series x ( t ) , the Lomb–Scargle Periodogram was used instead of Fourier transformation because the inter-beats intervals (IBI) are not uniformly distributed [36,37]. VLF (power in very low-frequency ranges, i.e., ≤0.04 Hz), LF (power in low-frequency ranges, i.e., 0.04–0.15 Hz), HF (Power in high-frequency ranges, i.e., 0.15–0.4 Hz), and total power (Power in all the frequency ranges, i.e., ≤0.4) were obtained by the integration of the power in the relevant frequency range in the spectrum as showed in Equation (1). Finally, the ratio between LF and HF (LF/HF) was also computed:
P S D f r e q = m i n f r e q m a x f r e q S x x ( ω ) d ω

2.4. Data Analysis

The relationship between “gold standard” HRV parameters computed in a 5 min window and ultra-short HRV parameters with and without missing values are assessed by Pearson’s correlation coefficient. Values higher than 0.8 reflect a strong correlation. The estimation error between HRV parameters compute on 5 min and in ultra-short time window with and without missing values were computed by the Root Mean Squared Error (RMSE). The magnitude of the difference was assessed by Cohen’s d effect size (ES). It is computed as the ratio between the mean difference between HRV gold standard parameters and ultra-short ones, and the pooled standard deviation as showed in Equation (2). ES < 0.2, (0.2–0.5], (0.5–0.8], >0.8 are considered trivial, small, medium, and large effect size. For example, if the ES is equal to 0.04, the two groups’ means differ by 0.04 standard deviations, and this difference should be considered trivial:
E S = H R V 5 m i n ¯ H R V u l t r a s h o r t ¯ ( σ H R V 5 m i n + σ H R V u l t r a s h o r t ) / 2
Bland and Altman plot (BA)was assessed in order to quantify the agreement between two quantitatives. Due to the fact that a huge number of Bland–Altman plots are required to describe all the comparisons, we decided to not insert all of these plots. We have resumed the results of BA plots by Bias, i.e., mean ± standard deviation of the differences, and systematic error, i.e., Pearson’s correlation coefficient between difference and means, of ultra-short HRV parameters with and without missing values compared to “gold standard” HRV parameters.
Finally, Spearman’s rank correlation coefficient was computed in order to assess the relationship between HRV features.
All the data preprocessing, feature extraction and data analyses are conducted by using the Python 3 programming language.

3. Results

Descriptive statistics of all the participants’ psycho-physiological features provided in MMASH dataset are provided in Rossi et al.’s data descriptor paper [29].
Table 1 and Table 2 report the descriptive statistics and statistical results of SDNN and rMSSD for all of the windows and missing values percentage, respectively. Small ES and moderate correlation (i.e., ES < 0.38, r > 0.69; see Table 1) are detected for ultra-short SDNN values with a time window longer than 2 min showing a bias lower than 10.70 ms (RMSE < 26.84 ms) without a statistically significant systematic error. On the contrary, moderate ES (i.e., ES > 0.45) is detected for time windows shorter than 1 min. On the contrary, rMSSD shows a trivial ES and moderate/strong correlation (i.e., ES < 0.20, r > 0.75; see Table 2) for almost all of the time windows showing a bias lower than 4.83 ms (RMSE < 21.47 ms) with a small systematic error (not statistically significant). In particular, only 70% of the missing values in 1 min and 30 s windows show small ES and moderate correlation (i.e., ES > 0.20, r < 0.62) with a bias higher than 5.37 ms (RMSE > 24.01 ms), while other time series lengths with a different percentage of missing values show trivial ES.
Figure 1 compares the power spectrum density (PSD) of the 5 min time series and the ultra-short ones with missing values. The lower the time series length, the lower the very low frequencies (VLF) that are possible to evaluate. In particular, frequencies lower than 0.004 (1/240 s), 0.006 (1/180 s), 0.008 (1/120 s), 0.017 (1/60 s), and 0.03 (1/30 s) Hz are not possible to evaluate for 4 min, 3 min, 2 min, 1 min, and 30 s, respectively. However, even if the time series decrease, it is possible to assess both LF and HF. Due to the fact that rMSSD reflects the HF (a strong correlation was detected between rMSSD and HF in 5 min time series—Table 3), it is possible to accurately estimate rMSSD even if the length of the time series decreases due to the fact that all HF could be assessed. On the contrary, missing very low frequencies negatively affect SDNN estimation as the length of the time series decreases. In particular, the moderate correlation detected between SDNN and Total power (Table 3) demonstrates that short time series do not permit assessing all the frequencies in the power spectrum negatively affecting the SDNN estimation.
In line with the results found in Table 1 and Table 2 for SDNN and rMSSD, respectively, Table 4 shows trivial and small ES for LF and HF for all of the ultra-short time series with almost all of the missing values (except for 30 s time series with 70% of missing value). Moreover, trivial and small ES for time series longer than 1 min and of 4 min were detected for total power and VLF, respectively. In particular, only total power of 4 min time series with all of the missing values percentage and 3 min time series with 0% of missing values show a trivial ES. Additionally, Table 4 and Figure 1 show that VLF and LF are underestimated as the percentage of missing values increases, while HF is overestimated. For these reasons, the length of the time series and the missing values strongly affect LF/HF HRV parameters. As a matter of fact, it is possible to estimate LF/HF values only in 4 min time windows without missing values where LF and HF are accurately detected.

4. Discussion

HRV is a widely used tool in many research areas, e.g., cardiovascular, preventive medicine, and sport [38,39,40,41,42,43,44], due to its non-invasive and reliable characteristics that permit assessing ANS and SAN activation. Nowadays, long- and short-term recordings (i.e., 24 h and 5 min, respectively) are considered the more reasonable options for HRV analysis. However, in real-world applications, shorter IBI recordings than 5-min are required which permit the mobile devices to instantaneously display results avoiding any possible misleading insight induced by motion artifacts. Therefore, representative studies are required in literature to provide reference values for short-term HRV analysis with missing values. The main finding of this study is that our results are in line with the ones reported for other previous studies [25,26,27,28]. As a matter of fact, we found that rMSSD requires at least 30 s to obtain reliable values (trivial effect size; Table 2). On the contrary, SDNN requires IBI timeseries longer than 2 min to have reliable results (small effect size; Table 1). Moreover, the novel results provided in this study are that rMSSD is more highly affected by missing values compared to SDNN (i.e., 0.12 ± 0.04 and 0.04 ± 0.03, respectively) showing, however, a small increment in effect size from all of the time series lengths as their percentage of missing values increase.
To truly estimate SDNN, we need precise information from the whole power spectrum, i.e., from VLF to HF [45]. As a matter of fact, we have found a moderate/strong relationship (r = 0.73) between total power and SDNN obtained from 5 min time series. In particular, SDNN is more correlated with VLF than with LF and HF (Table 3), and the power of VLF (49.92%) is larger than the power of LF (31.88%) and HF (18.20%) by several orders of magnitude (Table 4) resulting in a smaller contribution of HF and LF to SDNN than VLF [45]. Hence, the small ES and RMSE higher than 22.88 ms between SDNN obtained from 5 min time series and ultra-short ones are induced by the fact that it is not possible to evaluate all the VLF required to estimate SDNN (Figure 1). For example, in 2 min IBI time series, it is possible to evaluate the power of frequencies higher than 0.008 losing the information of frequencies between 0.0033 and 0.008. For this reason, SDNN shows a bias of about −9.51 ± 23.28 ( RMSE = 25.15) with a small effect size (ES = 0.34). For these reasons, SAN is lower reflected in ultra-short SDNN estimation compared to 5-min time window ones. Moreover, a very small change in ES (0.04 ± 0.03) between 0 to 70% of missing values in all the time windows demonstrated that SDNN is accurate even if motion artifacts affect the IBI time series. Hence, SDNN is more affected by the length of the time window than the missing values that could be detected throughout the IBI time series.
On the contrary, rMSSD is more sensitive to high-frequency (HF) fluctuations [13,45] compared to LF, VLF, and total power. As a matter of fact, strong correlation is detected between rMSSD and HF obtained in 5 min IBI time series (r = 0.80, Table 3). Hence, even if the length of the IBI time series decreases, it still possible to accurately assess the HF. As a matter of fact, Beak et al. [26] have found that rMSSD and in particular HF calculated using ultra-short IBI time series could be a good surrogate to those calculated using standard 5-min intervals. For this reason and as shown in Table 2, rMSSD is found to be accurate for all the IBI timeseries without missing values (ES < 0.11, RMSE < 15.13 ms, r < 0.79 and bias <−2.34 ms). For all of the time windows, ES shows a small increment with the missing values up to 50% (0.05 ± 0.02), while an increase of about 0.15 ± 0.03 for rMSSD estimated with 70% of missing values. However, trivial ES was observed for almost all the time windows with 70% of missing values. As a matter of fact, Table 4 and Figure 1 show trivial ES in HF for up to 50% of missing values for all the time windows, while a small ES was detected for time windows with 70% of missing values shorter than 1 min. On the contrary, missing values and the length of the time windows strongly affect the ratio between HF and LF. In particular, only a 4 min time window without missing values permits accurately estimating LF/HF HRV parameters.
These results highlight the fact that approaches to interpolating missing values are not required for rMSSD, but they are needed for SDNN estimation as also suggested by Morelli et al. [23]. As a matter of fact, rMSSD captures fast changes in heart activity (HF), while SDNN is an index of the entire power spectrum and in particular of VLF (Table 3 and Table 4). For this reason, interpolation methods that act as low pass filters affecting the power spectrum density of IBI time series have a negative impact on rMSSD, while a small effect on SDNN [23]. On the contrary, no interpolation approach introduced white errors, thus minimizing the impact on IBI successive differences that were the first computation step of rMSSD. Moreover, due to the fact that SDNN can be estimated from the entire power spectrum [45], the estimation of missing frequencies is required to improve accuracy of SDNN. Future works are scheduled in order to assess the effect of interpolation approaches and missing frequencies’ estimation from known ones. Moreover, due to the fact that low cost wearable devices equipped with heart rate sensors are becoming more and more accurate in IBI recording, the reduction of noisy data could be the right choice to improve the quality of HRV parameters’ estimation. However, this paper suggests that, even if the IBI time series contain a large amount of noise, a small error is expected in the estimation of both SDNN and rMSSD. In conclusion, wrist-worn low-cost devices could be a good surrogate of a good quality simple ECG even if they could have a huge quantity of missing values.
A limitation of this study is that the participants are all young males. Future works could be performed in order to detect if individual characteristics such as gender, age, and individual’s habits could affect the estimation of ultra-short SDNN and rMSSD with and without missing values. Moreover, due to the fact that SDNN24 is an important HRV parameter to detect people’s health status, an investigation is required to assess the effect of missing values on this HRV parameter.

5. Conclusions

As showed in previous studies, it is possible to accurately estimate rMSSD and SDNN by using 30 s and 2 min time windows, respectively. In particular, in this paper, we highlight the fact that ultra-short IBI time series do not permit assessing all the VLF increasing the bias on SDNN. Moreover, to the best of our knowledge, this is the first study that evaluates the effect of missing values on ultra-short HRV data. We found that the missing values have a low effect on both SDNN estimations, while having a moderate effect on rMSSD. In particular, the higher the missing values, the higher the bias in both rMSSD and SDNN showing trivial and small magnitudes of the differences, respectively.

Author Contributions

Conceptualization, A.R and D.M.; methodology, A.R. and D.M.; validation, A.R.; formal analysis, A.R.; investigation, A.R.; resources, D.M.; data curation, A.R.; writing—original draft preparation, A.R.; writing—review and editing, D.M., D.P., and D.A.C.; visualization, A.R.; supervision, D.M. and D.A.C.; project administration, D.M. and D.A.C.; funding acquisition, D.M. and D.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the European Community’s H2020 Program under the funding scheme H2020-INFRAIA-2019-1 Research Infrastructures grant agreement 871042, www.sobigdata.eu, SoBigData++: European Integrated Infrastructure for Social Mining and Big Data Analytics. Moreover, the research was also funded supported by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC). The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McCraty, R.; Shaffer, F. Heart rate variability: New perspectives on physiological mechanisms, assessment of self-regulatory capacity, and health risk. Glob. Adv. Health Med. 2015, 4, 46–61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Shaffer, F.; Ginsberg, J. An Overview of Heart Rate Variability Metrics and Norms. Front. Public Health 2017, 5, 258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rosenberg, A.; Weiser-Bitoun, I.; Billman, G.; Yaniv, Y. Signatures of the autonomic nervous system and the heart’s pacemaker cells in canine electrocardiograms and their applications to humans. Sci. Rep. 2020, 10, 9971. [Google Scholar] [CrossRef] [PubMed]
  4. Petretta, M.; Bonaduce, B.; Scalfi, L.; De Filippo, E.; Marciano, F.; Migaux, M.; Themistoclakis, S.; Ianniciello, A.; Contaldo, F. Heart rate variability as a measure of autonomic nervous system function in anorexia nervosa. Clin. Cardiol. 1997, 20, 219–224. [Google Scholar] [CrossRef] [PubMed]
  5. Karemaker, J. An introduction into autonomic nervous function. Physiol. Meas. 2017, 38, 89–118. [Google Scholar] [CrossRef]
  6. Tsutsui, K.; Monfredi, O.; Sirenko-Tagirova, S.; Maltseva, L.; Bychkov, R.; Kim, M.S.; Ziman, B.D.; Tarasov, K.V.; Tarasova, Y.S.; Wang, M.; et al. A coupled-clock system drives the automaticity of human sinoatrial nodal pacemaker cells. Sci. Signal. 2018, 11, eaap7608. [Google Scholar] [CrossRef] [Green Version]
  7. Stein, P.; Domitrovich, P.; Hui, N.; Rautaharju, P.; Gottdiener, J. Sometimes higher heart rate variability is not better heart rate variability: Results of graphical and nonlinear analyses. J. Cardiovasc. Electrophysiol. 2005, 16, 954–959. [Google Scholar] [CrossRef]
  8. Laborde, S.; Mosley, E.; Thayer, J. Heart Rate Variability and Cardiac Vagal Tone in Psychophysiological Research—Recommendations for Experiment Planning, Data Analysis and Data Reporting. Front. Psychol. 2017, 8, 213. [Google Scholar] [CrossRef] [Green Version]
  9. Camm, A.J.; Malik, M.; Bigger, J.T.; Breithardt, G.; Cerutti, S.; Cohen, R.J.; Coumel, E.L.; Fallen, H.L.; Kennedy, R.E.; Lombardi, F.; et al. Heart rate variability: Standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Circulation 1996, 93, 1043–1065. [Google Scholar]
  10. Kleiger, R.; Miller, J.; Bigger, J.; Moss, A.J. Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. Am. J. Cardiol. 1987, 59, 256–262. [Google Scholar] [CrossRef]
  11. Karcz, M.; Chojnowska, L.W.Z.R. Prognostic significance of heart rate variability in dilated cardiomyopathy. EP Eur. 2003, 87, 75–81. [Google Scholar] [CrossRef]
  12. Gao, L.; Chen, Y.; Shi, Y.; Xue, H.; Wang, J. Value of DC and DRs in prediction of cardiovascular events in acute myocardial infarction patients. Zhonghua Yi Xue Za Zhi 2016, 96. [Google Scholar] [CrossRef]
  13. Shaffer, F.; McCraty, R.; Zerr, C. A healthy heart is not a metronome: An integrative review of the heart’s anatomy and heart rate variability. Front. Psychol. 2014, 5, 1040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zygmunt, A.; Stanczyk, J. Methods of evaluation of autonomic nervous system function. Arch. Med. Sci. 2010, 6, 11–18. [Google Scholar] [CrossRef]
  15. Plews, D.; Laursen, P.; Kilding, A.; Buchheit, M. Heart rate variability and training intensity distribution in elite rowers. Int. J. Sport. Physiol. Perform. 2014, 9, 1026–1032. [Google Scholar] [CrossRef]
  16. Chen, J.; Yeh, D.; Lee, J.; Chen, C.; Huang, C.; Lee, S.D.; Chen, C.C.; Kuo, T.B.J.; Kao, C.L.; Kuo, C.H. Parasympathetic nervous activity mirrors recovery status in weightlifting performance after training. J. Strength Cond. Res. 2011, 25, 1546–1552. [Google Scholar] [CrossRef] [Green Version]
  17. Edmonds, R.; Sinclair, W.; Leicht, A. Effect of a training week on heart rate variability in elite youth rugby league players. Int. J. Sport. Med. 2013, 34, 1087–1092. [Google Scholar] [CrossRef]
  18. Jennings, J.; Sheu, L.; Kuan, D.; Manuck, S.; Gianaros, P. Resting state connectivity of the medial prefrontal cortex covaries with individual differences in high-frequency heart rate variability. Psychophysiology 2016, 53, 444–454. [Google Scholar] [CrossRef] [Green Version]
  19. Kemp, A.; Quintana, D.; Gray, M.; Felmingham, K.; Brown, K.; Gatt, J. Impact of depression and antidepressant treatment on heart rate variability: A review and meta-analysis. Biol. Psychiatry 2010, 67, 1067–1074. [Google Scholar] [CrossRef]
  20. Agelink, M.; Majewski, T.; Wurthmann, C.; Postert, T.; Linka, T.; Rotterdam, S. Autonomic neurocardiac function in patients with major depression and effects of antidepressive treatment with nefazodone. J. Affect. Disord. 2001, 62, 187–198. [Google Scholar] [CrossRef]
  21. Haghi, M.; Thurow, K.; Stoll, R. Wearable Devices in Medical Internet of Things: Scientific Research and Commercially Available Devices. Healthc. Inf. Res 2017, 23, 4–15. [Google Scholar] [CrossRef] [PubMed]
  22. Karlsson, M.; Hörnsten, R.; Rydberg, A.; Wiklund, U. Automatic filtering of outliers in RR intervals before analysis of heart rate variability in Holter recordings: A comparison with carefully edited data. Biomed. Eng. 2012, 11, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Morelli, D.; Rossi, A.; Cairo, M.; Clifton, D. Analysis of the Impact of Interpolation Methods of Missing RR-Intervals Caused by Motion Artifacts on HRV Features Estimations. Sensors 2019, 19, 3163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kim, K.; Lim, Y.; Kim, J.; Park, K. Effect of missing RR-interval data on heart rate variability analysis in the time domain. Physiol. Meas. 2007, 28, 1485–1494. [Google Scholar] [CrossRef]
  25. Munoz, M.; Van Roon, A.; Riese, H.; Thio, C.; Oostenbroek, E.; Westrik, I.; De Geus, E.J.C.; Gansevoort, R.; Lefrandt, J.; Nolte, I.; et al. Validity of (Ultra-)Short Recordings for Heart Rate Variability Measurements. PLoS ONE 2015, 10, e0138921. [Google Scholar] [CrossRef] [Green Version]
  26. Beak, H.; Cho, C.; Cho, J.; Woo, J. Reliability of Ultra-Short-Term Analysis as a Surrogate of Standard 5-Min Analysis of Heart Rate Variability. Telemed. e-Health 2015, 21, 404–414. [Google Scholar] [CrossRef]
  27. Castaldo, R.; Montesinos, P.; Melillo, C.; Pecchia, L. Ultra-short term HRV features as surrogates of short term HRV: A case study on mental stress detection in real life. BMC Med. Inform. Decis. Mak. 2019, 19, 12. [Google Scholar] [CrossRef] [Green Version]
  28. Pecchia, L.; Castaldo, R.; Montesinos, L.; Melillo, P. Are ultra-short heart rate variability features good surrogates of short-term ones? State-of-the-art review and recommendations. Healthc. Technol. Lett. 2018, 5, 94–100. [Google Scholar] [CrossRef]
  29. Rossi, A.; Da Pozzo, E.; Menicagli, D.; Tremolanti, C.; Priami, C.; Sirbu, A.; Clifton, D.; Martini, C.; Morelli, D. A Public Dataset of 24-h Multi-Levels Psycho-Physiological Responses in Young Healthy Adults. Data 2020, 5, 91. [Google Scholar] [CrossRef]
  30. Rossi, A.; Da Pozzo, E.; Menicagli, D.; Tremolanti, C.; Priami, C.; Sirbu, A.; Clifton, D.; Martini, C.; Morelli, D. Multilevel Monitoring of Activity and Sleep in Healthy People (version 1.0.0). PhysioNet 2020. [Google Scholar] [CrossRef]
  31. Goldberger, A.; Amaral, L.; Glass, L.; Hausdorff, J.; Ivanov, P.; Mark, R.; Stanley, H. PhysioBank, PhysioToolkit, and PhysioNet: Components of a new research resource for complex physiologic signals. Circulation 2000, 101, e215–e220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Plews, D.; Scott, B.; Altini, M.; Wood, M.; Kilding, A.; Laursen, P. Comparison of Heart-Rate-Variability Recording With Smartphone Photoplethysmography, Polar H7 Chest Strap, and Electrocardiography. Int. J. Sport. Physiol. Perform. 2017, 12, 1324–1328. [Google Scholar] [CrossRef] [PubMed]
  33. Hideki, I. Essentials of Error-Control Coding Techniques; Academic Press: Cambridge, MA, USA, 1990. [Google Scholar]
  34. Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology. Heart Rate Variability: Standards of Measurment, Physiological Interpretation, and Clinical Use. Circulation 1996, 93, 1043–1065. [Google Scholar] [CrossRef] [Green Version]
  35. Liao, D.; Cai, J.; Rosamond, W.; Barnes, R.; Hutchinson, R.; Whitsel, E. Cardiac autonomic function and incident coronary heart disease: A population-based case-cohort study. The ARIC study. Am. J. Epidemiol. 1997, 145, 696–706. [Google Scholar] [CrossRef] [Green Version]
  36. Lomb, N. Least-squares frequency analysis of unequally spaced data. Astrophys. Space Sci. 1976, 39, 447–462. [Google Scholar] [CrossRef]
  37. Scargle, J.D. Studies in Astronomical Time Series Analysis. V. Bayesian Blocks, a New Method to Analyze Structure in Photon Counting Data. Am. Astron. Soc. 1998, 504, 405. [Google Scholar] [CrossRef] [Green Version]
  38. Kleiger, R.; Stein, P.; Bigger, J. Heart rate variability: Measurement and clinical utility. Ann. Noninvasive Electrocardiol. 2005, 10, 88–101. [Google Scholar] [CrossRef]
  39. Baek, H.; Lee, H.; Kim, J.; Choi, J.; Kim, K.; Park, K. Non intrusive biological signal monitoring in a car to evaluate a driver’s stress and health state. Telemed. J. e-Health 2009, 15, 182–189. [Google Scholar] [CrossRef] [PubMed]
  40. Gouin, J.; Wenzel, K.; Deschenes, S.; Dang-Vu, T. Heart rate variability predicts sleep efficiency. Sleep Med. 2013, 14, e142. [Google Scholar] [CrossRef]
  41. Patel, M.; Lal, S.; Kavanagh, D.; Rossiter, P. Applying neural network analysis on heart rate variability data to assess deriver fatigue. Expert Syst. Appl. 2011, 38, 7235–7242. [Google Scholar] [CrossRef]
  42. Kaikkonen, P.; Hynynen, E.; Mann, T.; Rusko, H.; Nummela, A. Heart rate variability is related to training load variables in interval running exercises. Eur. J. Appl. Physiol. 2012, 112, 829–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Stein, P.; Pu, Y. Heart rate variability, sleep and sleep disorders. Sleep Med. Rev. 2012, 16, 47–66. [Google Scholar] [CrossRef] [PubMed]
  44. Thayer, J.; Ahs, F.; Fredrikson, M.; Sollers, J.; Wager, T. A meta-analysis of heartrate variability and neuromonitoring studies: Implications for heart rate variability as a marker of stress and health. Neurosci. Biobehav. Rev. 2012, 36, 747–756. [Google Scholar] [CrossRef] [PubMed]
  45. Malik, M.; Bigger, J.; Camm, A.; Kleiger, R.; Malliani, A.; Moss, A.; Schwartz, P. Heart rate variability: Standards of measurement, physiological interpretation and clinical use. Circulation 1996, 93, 927–934. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Power Spectrum Density plot for timeseries with different time series length and missing values percentage. The blue dotted line in each plot reflects the power spectrum of the 5 min time window.
Figure 1. Power Spectrum Density plot for timeseries with different time series length and missing values percentage. The blue dotted line in each plot reflects the power spectrum of the 5 min time window.
Sensors 20 07122 g001
Table 1. Descriptive statistics and statistical analysis of SDNN for all of the time window and missing values. The statistical analysis is focused to compare 5 min HRV parameters to ultra-short ones. * refers to small ES of ultra-short HRV features compared to 5 min ones, respectively. r, RMSE, and ES refer to the correlation, the Root Mean Square Error, and the effect size between the ‘gold standard’ and the ultra-short with and without missing values on SDNN, respectively. Bias refers to mean ± standard deviation of the difference of ultra-short HRV parameters compared to ‘gold standard’ HRV parameters (i.e., 5 min time window), while systematic error is the Pearson’s correlation coefficient between means and differences between 5 min and ultra-short HRV estimation.
Table 1. Descriptive statistics and statistical analysis of SDNN for all of the time window and missing values. The statistical analysis is focused to compare 5 min HRV parameters to ultra-short ones. * refers to small ES of ultra-short HRV features compared to 5 min ones, respectively. r, RMSE, and ES refer to the correlation, the Root Mean Square Error, and the effect size between the ‘gold standard’ and the ultra-short with and without missing values on SDNN, respectively. Bias refers to mean ± standard deviation of the difference of ultra-short HRV parameters compared to ‘gold standard’ HRV parameters (i.e., 5 min time window), while systematic error is the Pearson’s correlation coefficient between means and differences between 5 min and ultra-short HRV estimation.
Time WindowMissing ValuesSDNNrRMSEBiasSystematic ErrorES
5 min72.70 ± 30.23
4 min0%69.32 ± 35.830.7923.17−6.92 ± 22.110.080.20 *
5%69.15 ± 36.300.7923.20−6.94 ± 22.170.080.20 *
10%69.29 ± 36.420.7923.25−7.02 ± 22.310.080.20 *
15%69.18 ± 36.360.7823.31−7.08 ± 22.450.080.21 *
30%69.02 ± 36.160.7823.87−7.20 ± 22.750.080.21 *
50%68.62 ± 36.260.7724.71−7.48 ± 23.550.090.21 *
70%67.64 ± 36.660.7426.02−8.27 ± 24.670.100.23 *
3 min0%65.96 ± 31.990.7622.88−7.64 ± 21.570.040.26 *
5%65.92 ± 31.970.7622.95−7.66 ± 21.590.040.26 *
10%65.9 ± 31.790.7622.98−7.69 ± 21.610.040.26 *
15%65.79 ± 31.950.7623.06−7.73 ± 21.650.040.26 *
30%65.72 ± 32.180.7623.19−7.84 ± 21.820.050.26 *
50%65.31 ± 32.330.7523.71−8.19 ± 22.250.050.27 *
70%64.54 ± 32.980.7324.65−8.76 ± 23.040.080.29 *
2 min0%63.57 ± 31.390.7125.15−9.51 ± 23.280.030.34 *
5%63.52 ± 31.420.7125.19−9.56 ± 23.510.030.34 *
10%63.5 ± 31.370.7125.25−9.59 ± 23.570.030.35 *
15%63.4 ± 31.440.7125.31−9.60 ± 23.550.040.35 *
30%63.22 ± 31.560.7025.66−9.76 ± 23.730.040.35 *
50%62.82 ± 31.850.7025.98−10.60 ± 23.950.050.37 *
70%61.94 ± 32.350.6926.84−10.70 ± 24.620.070.38 *
1 min0%58.95 ± 30.350.6329.03−13.50 ± 25.700.010.45 *
5%58.9 ± 31.090.6329.16−13.59 ± 25.780.010.45 *
10%58.83 ± 30.930.6329.21−13.64 ± 25.810.010.45 *
15%58.69 ± 31.150.6329.28−13.70 ± 25.840.010.45 *
30%58.54 ± 30.530.6329.46−13.84 ± 26.000.020.45 *
50%58.04 ± 30.780.6229.98−14.18 ± 26.420.020.47 *
70%56.61 ± 31.310.6030.98−15.28 ± 26.940.040.52
30 s0%51.41 ± 28.620.5334.18−19.97 ± 27.74−0.010.67
5%51.36 ± 28.70.5334.21−20.01 ± 27.88−0.010.67
10%51.14 ± 28.540.5334.36−20.06 ± 27.91−0.010.69
15%51.23 ± 28.860.5234.48−20.08 ± 27.93−0.010.69
30%50.77 ± 28.800.5234.73−20.49 ± 28.04−0.010.70
50%49.95 ± 29.050.5135.42−21.16 ± 28.410.000.71
70%48.11 ± 29.630.4836.85−22.63 ± 29.080.020.77
Table 2. Descriptive statistics (mean ± standard deviation) and statistical analysis of rMSSD for all of the time window and missing values. The statistical analysis is focused to compare 5 min HRV parameters to ultra-short ones. ** and * refer to the trivial and small ES of ultra-short HRV features compared to 5 min ones, respectively. r, RMSE, and ES refer to the correlation, the Root Mean Square Error and the effect size between the ‘gold standard’ and the ultra-short with and without missing values on rMSSD, respectively. Bias refers to mean ± standard deviation of the difference of ultra-short HRV parameters compared to ‘gold standard’ HRV parameters (i.e., 5 min time window), while systematic error is the Pearson’s correlation coefficient between means and differences between 5 min and ultra-short HRV estimation.
Table 2. Descriptive statistics (mean ± standard deviation) and statistical analysis of rMSSD for all of the time window and missing values. The statistical analysis is focused to compare 5 min HRV parameters to ultra-short ones. ** and * refer to the trivial and small ES of ultra-short HRV features compared to 5 min ones, respectively. r, RMSE, and ES refer to the correlation, the Root Mean Square Error and the effect size between the ‘gold standard’ and the ultra-short with and without missing values on rMSSD, respectively. Bias refers to mean ± standard deviation of the difference of ultra-short HRV parameters compared to ‘gold standard’ HRV parameters (i.e., 5 min time window), while systematic error is the Pearson’s correlation coefficient between means and differences between 5 min and ultra-short HRV estimation.
Time WindowMissing ValuesrMSSDrRMSEBiasSystematic ErrorES
5 min42.49 ± 20.75
4 min0%42.34 ± 22.120.937.87−0.56 ± 7.850.160.03 **
5%42.18 ± 22.080.928.12−0.58 ± 7.910.160.03 **
10%42.11 ± 22.060.928.48−0.62 ± 8.180.180.03 **
15%42.07 ± 22.460.918.55−0.67 ± 8.540.190.04 **
30%42.02 ± 22.960.909.97−0.87 ± 9.930.210.05 **
50%41.35 ± 24.060.8512.78−1.50 ± 12.690.260.07 **
70%39.04 ± 26.090.7218.34−3.69 ± 17.960.320.16 **
3 min0%42.15 ± 22.360.928.82−0.57 ± 8.810.190.04 **
5%42.05 ± 22.540.909.01−0.66 ± 8.980.220.04 **
10%41.94 ± 23.010.899.56−0.71 ± 9.060.230.04 **
15%41.78 ± 23.070.899.88−0.85 ± 9.810.250.04 **
30%41.73 ± 23.210.8810.97−0.96 ± 10.930.240.05 **
50%41.08 ± 24.450.8313.86−1.56 ± 13.760.290.08 **
70%39.00 ± 26.890.7119.21−3.48 ± 18.890.360.16 **
2 min0%41.83 ± 22.700.909.85−0.82 ± 9.820.190.04 **
5%41.79 ± 22.780.8810.51−1.16 ± 10.050.250.05 **
10%41.81 ± 22.180.8810.36−1.12 ± 9.980.250.05 **
15%41.79 ± 22.810.8611.76−1.16 ± 10.800.270.06 **
30%41.29 ± 23.600.8612.05−1.31 ± 11.980.250.07 **
50%40.40 ± 24.750.8014.80−2.17 ± 14.640.280.11 **
70%38.38 ± 27.600.6820.63−4.00 ± 20.240.380.17 **
1 min0%41.26 ± 23.470.8612.02−1.22 ± 11.960.240.06 **
5%40.97 ± 23.550.8512.48−1.50 ± 12.870.240.07 **
10%40.90 ± 23.610.8412.69−1.53 ± 13.130.260.07 **
15%40.69 ± 23.770.8413.17−1.71 ± 13.240.260.09 **
30%40.54 ± 24.570.8114.42−1.90 ± 14.290.290.09 **
50%39.52 ± 26.230.7517.72−2.84 ± 17.500.350.13 **
70%36.73 ± 29.700.6224.01−5.37 ± 23.400.440.22 *
30 s0%38.73 ± 24.130.7915.13−2.34 ± 14.950.300.11 **
5%38.71 ± 24.180.7915.38−2.55 ± 15.860.310.13 **
10%38.52 ± 24.670.7815.99−2.57 ± 15.970.310.13 **
15%37.74 ± 24.980.7617.05−3.13 ± 17.110.330.13 **
30%37.66 ± 25.340.7317.83−3.37 ± 17.510.330.15 **
50%36.07 ± 27.200.6521.41−4.83 ± 20.860.390.19 **
70%32.63 ± 30.470.5327.14−7.80 ± 26.000.470.31 *
Table 3. Spearman’s rank correlation coefficient between HRV features in 5 min time series.
Table 3. Spearman’s rank correlation coefficient between HRV features in 5 min time series.
SDNNrMSSD
VLF0.620.19
LF0.510.43
HF0.450.80
LF/HF−0.04−0.47
Total Power0.730.45
Table 4. Descriptive statistics (mean ± standard error of the mean) of VLF, LF, HF, LF/HF, and Total power of the power spectrum density for all of the time window and missing values. The statistical analysis is focused to compare 5 min HRV parameters to ultra-short ones. ** and * refer to the trivial and small ES of ultra-short HRV features compared to 5 min ones, respectively.
Table 4. Descriptive statistics (mean ± standard error of the mean) of VLF, LF, HF, LF/HF, and Total power of the power spectrum density for all of the time window and missing values. The statistical analysis is focused to compare 5 min HRV parameters to ultra-short ones. ** and * refer to the trivial and small ES of ultra-short HRV features compared to 5 min ones, respectively.
Time WindowMissing ValuesVLFLFHFLF/HFTotal Power
5 min1891.44 ± 46.101208.11 ± 30.84689.49 ± 24.843.31 ± 0.043907.50 ± 86.95
4 min0%1533.88 ± 37.96 **1208.11 ± 30.84 **689.48 ± 24.84 **3.31 ± 0.04 **3549.94 ± 80.76 **
5%1462.12 ± 35.36 **1183.27 ± 31.12 **736.11 ± 25.53 **2.72 ± 0.033496.31 ± 80.57 **
10%1397.1 ± 34.58 **1147.97 ± 30.16 **782.07 ± 25.4 **2.31 ± 0.023438.69 ± 79.8 **
15%1326.12 ± 33.39 **1112.15 ± 24.97 **828 ± 24.78 **2.02 ± 0.023374.25 ± 71.88 **
30%1128.62 ± 27.88 *1027.24 ± 17.10 **954.92 ± 22.97 *1.44 ± 0.013208.93 ± 59.54 **
50%853.10 ± 22.20 *917.61 ± 19.01 **1123.03 ± 30.88 *0.99 ± 0.012976.95 ± 66.27 **
70%567.70 ± 15.31819.50 ± 27.86 **1313.26 ± 32.97 *0.69 ± 0.0042778.61 ± 73.01 **
3 min0%1127.11 ± 25.10 *1110.41 ± 19.26 **547.88 ± 25.18 **3.41 ± 0.042897.76 ± 58.20 **
5%1079.4 ± 24.31 *1083.32 ± 18.99 **591.89 ± 26.46 **2.81 ± 0.032861.37 ± 59.38 *
10%1023.22 ± 22.57 *1047.71 ± 15.22 **625.92 ± 17.50 **2.4 ± 0.022799.93 ± 46.28 *
15%975.01 ± 21.49 *1020.62 ± 16.45 **666.64 ± 18.70 **2.09 ± 0.022761.14 ± 48.79 *
30%827.45 ± 18.66 *937.27 ± 13.08 **785.34 ± 13.43 **1.48 ± 0.012640.77 ± 40.29 *
50%625.49 ± 13.81822.96 ± 12.18 *960.21 ± 17.47 *1.01 ± 0.012484.46 ± 40.24 *
70%429.55 ± 9.17719.21 ± 13.36 *1137.89 ± 20.49 *0.70 ± 0.0042350.27 ± 41.60 *
2 min0%1035.59 ± 24.53 *1082.92 ± 18.27 **518.80 ± 16.84 **3.42 ± 0.042743.70 ± 48.45 *
5%989.72 ± 23.62 *1054.75 ± 17.79 **554.49 ± 18.46 **2.86 ± 0.032702.02 ± 49.17 *
10%946.56 ± 22.56 *1028.3 ± 17.47 **592.71 ± 16.89 **2.47 ± 0.032667.16 ± 47.16 *
15%893.91 ± 21.12 *997.58 ± 16.83 **639.02 ± 19.08 **2.14 ± 0.022626.85 ± 48.55 *
30%753.04 ± 17.71904.07 ± 14.64 **743.62 ± 14.54 **1.51 ± 0.012486.35 ± 41.52 *
50%572.83 ± 13.52786.10 ± 12.41 *895.60 ± 15.57 *1.02 ± 0.012326.40 ± 38.32 *
70%394.79 ± 9.66682.95 ± 14.42 *1067.63 ± 28.46 *0.71 ± 0.0052204.56 ± 50.53 *
1 min0%616.19 ± 16.251060.77 ± 17.50 **491.83 ± 10.75 **3.46 ± 0.042272.02 ± 38.02 *
5%588.58 ± 15.481033.71 ± 16.93 **527.44 ± 11.94 **2.93 ± 0.032249.53 ± 38.05 *
10%559.61 ± 14.65999.09 ± 16.35 **557.29 ± 9.69 **2.51 ± 0.032212.12 ± 35.94 *
15%534.03 ± 13.78964.84 ± 15.65 **589.81 ± 10.8 **2.21 ± 0.022181.41 ± 35.56 *
30%453.45 ± 11.82880.46 ± 14.18 *693.12 ± 12.15 **1.57 ± 0.022110.12 ± 34.48 *
50%337.15 ± 8.56741.32 ± 12.09 *831.16 ± 14.37 **1.05 ± 0.011976.68 ± 32.74 *
70%226.16 ± 5.78623.91 ± 11.57 *964.02 ± 17.28 *0.71 ± 0.0051867.88 ± 33.62
30 s0%1043.22 ± 19.46 **497.03 ± 10.04 **3.47 ± 0.051634.33 ± 27.33
5%1009.19 ± 19.01 **518.67 ± 10.16 **3.01 ± 0.041618.65 ± 27.22
10%973.86 ± 18.21 **536.36 ± 9.59 **2.64 ± 0.041597.54 ± 26.45
15%932.54 ± 17.52 **571.69 ± 10.78 **2.29 ± 0.031588.48 ± 26.87
30%820.17 ± 15.69 *628.70 ± 11.05 **1.65 ± 0.021521.87 ± 26.04
50%671.04 ± 13.42 *720.85 ± 13.05 **1.08 ± 0.011451.23 ± 25.83
70%539.98 ± 11.27839.28 ± 16.60 *0.72 ± 0.011425.51 ± 27.72
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rossi, A.; Pedreschi, D.; Clifton, D.A.; Morelli, D. Error Estimation of Ultra-Short Heart Rate Variability Parameters: Effect of Missing Data Caused by Motion Artifacts. Sensors 2020, 20, 7122. https://doi.org/10.3390/s20247122

AMA Style

Rossi A, Pedreschi D, Clifton DA, Morelli D. Error Estimation of Ultra-Short Heart Rate Variability Parameters: Effect of Missing Data Caused by Motion Artifacts. Sensors. 2020; 20(24):7122. https://doi.org/10.3390/s20247122

Chicago/Turabian Style

Rossi, Alessio, Dino Pedreschi, David A. Clifton, and Davide Morelli. 2020. "Error Estimation of Ultra-Short Heart Rate Variability Parameters: Effect of Missing Data Caused by Motion Artifacts" Sensors 20, no. 24: 7122. https://doi.org/10.3390/s20247122

APA Style

Rossi, A., Pedreschi, D., Clifton, D. A., & Morelli, D. (2020). Error Estimation of Ultra-Short Heart Rate Variability Parameters: Effect of Missing Data Caused by Motion Artifacts. Sensors, 20(24), 7122. https://doi.org/10.3390/s20247122

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