Next Article in Journal
Motor-Imagery EEG-Based BCIs in Wheelchair Movement and Control: A Systematic Literature Review
Next Article in Special Issue
A Novel Methodology for the Synchronous Collection and Multimodal Visualization of Continuous Neurocardiovascular and Neuromuscular Physiological Data in Adults with Long COVID
Previous Article in Journal
An Approach to Steady-State Power Transformer Modeling Considering Direct Current Resistance Test Measurements
Previous Article in Special Issue
Generation of Gait Events with a FSR Based Cane Handle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of ECG Data Length on Heart Rate Variability among Young Healthy Adults

1
Department of Biomedical Engineering, Henry Samueli School of Engineering, University of California at Irvine, Irvine, CA 92697, USA
2
School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ 85281, USA
3
Department of Physical Therapy, Crean College of Health and Behavioral Sciences, Chapman University, Irvine, CA 92618, USA
4
Department of Electrical and Computer Science Engineering, Fowler School of Engineering, Chapman University, Orange, CA 92866, USA
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(18), 6286; https://doi.org/10.3390/s21186286
Submission received: 1 August 2021 / Revised: 2 September 2021 / Accepted: 14 September 2021 / Published: 19 September 2021
(This article belongs to the Collection Sensors for Gait, Posture, and Health Monitoring)

Abstract

:
The relationship between the robustness of HRV derived by linear and nonlinear methods to the required minimum data lengths has yet to be well understood. The normal electrocardiography (ECG) data of 14 healthy volunteers were applied to 34 HRV measures using various data lengths, and compared with the most prolonged (2000 R peaks or 750 s) by using the Mann–Whitney U test, to determine the 0.05 level of significance. We found that SDNN, RMSSD, pNN50, normalized LF, the ratio of LF and HF, and SD1 of the Poincaré plot could be adequately computed by small data size (60–100 R peaks). In addition, parameters of RQA did not show any significant differences among 60 and 750 s. However, longer data length (1000 R peaks) is recommended to calculate most other measures. The DFA and Lyapunov exponent might require an even longer data length to show robust results. Conclusions: Our work suggests the optimal minimum data sizes for different HRV measures which can potentially improve the efficiency and save the time and effort for both patients and medical care providers.

1. Introduction

Heart rate variability (HRV) is a promising measure used to assess cardiovascular health by investigating heartbeat fluctuations over time. Electrocardiogram (ECG) is an autonomically controlled physiological vital signal that changes due to sympathetic or parasympathetic perturbations. As such, reduction in HRV is a well-established biomarker of diabetes [1,2], cardiovascular disease (CVD) [3,4], inflammation [5,6,7], obesity [8] and psychiatric disorders [9,10]. Over the past half-century, three groups of mathematical methods for determining HRV have been proposed (i) time domain, (ii) frequency domain, and (iii) nonlinear analyses [11].
Time- and frequency-domain HRV measures were standardized in 1996 by the task force of The European Society of Cardiology and the North American Society for Pacing and Electrophysiology [12]. Nonlinear variability measurements such as the Lyapunov exponent, fractal, entropy, and symbolic entropy are well-established [13]. However, the standardization in selecting a time series length for robustness in differentiating different populations and ailments is currently lacking. HRV measurements are influenced by data length, sampling frequency [14,15,16,17], noise [14,18] and, computation parameters used in specific methods such as the time delay and embedding dimensions [16,19,20,21,22,23]. The length of the heart-rate time series data are often considered a limitation for utilizing nonlinear analyses. The number of data points in time series is critical for nonlinear analysis since it is unknown whether fewer data sets can characterize the whole dynamics of the system.
The length of data is an essential factor considering shorter data acquisition time can improve patient throughput and efficiency of hospitals, healthcare providers, and home monitoring in general. Additionally, it improves patient’s adherence and overall experience. An essential part of the HRV analysis is knowing how many data points are needed to describe the system appropriately. An important rule of thumb suggested by researchers is to choose time series of at least 10d data points with ‘d’ as the system’s embedding dimension [24]. In such a scenario, if the embedding dimension is 6, then at least one million data points are required. However, sometimes obtaining such long time series data from human subjects in controlled clinical environments is practically impossible. For example, to collect one million heartbeats, one has to record ECG continuously for about 277 h. This makes human subject data collection practically impossible. Besides, it is not known if shorter time series can accurately characterize the system’s dynamics. An optimal data length should capture the essential dynamics of the system. Thus, it is imperative to understand the relationship between the robustness of the HRV to data lengths, such that the minimum data points necessary for accurate measurements can be determined.
Traditionally in HRV analysis, the data acquisition time is set to at least 5 min [12]. A limitation to longer data acquisition is that the hydrogel layer used in ECG electrodes can degrade and lower the signal-to-noise ratio [20,25,26,27]. Some studies evaluated the influence of shorter data acquisition time in different HRV measures [28]. For instance, Munoz and coworkers suggested that some time-domain HRV measures can be reliably obtained from less than a minute of recording [29]. Others investigated the influences of data length in both time- and frequency-domain HRV measures [20,30,31,32,33]. In the time-domain analysis, the standard deviation of normal-to-normal (NN) interval (SDNN), root mean square of standard deviation (RMSSD), and the percentage of successive NN intervals greater than 50 ms in all NN internals (pNN50) have been suggested as reliable measures for 5 min data lengths. However, frequency-domain HRV measures cannot produce consistent conclusions. For nonlinear dynamics, Entropy-based HRV measures [34,35,36] have been used to differentiate patients with cardiovascular disease from healthy controls using shorter time-series data. Sample entropy (SampEn) is reported to be less dependent on data length than approximate entropy (ApEn) [35]. In conclusion, the relationship between the robustness of HRV derived by linear and nonlinear methods to the required minimum data lengths has yet to be systematically evaluated and clearly understood.
In addressing such limitations, this study explores how the data length of the ECG signal affects the HRV measures (time and frequency domain variables and nonlinear variables). In this study, 14 healthy volunteers were monitored in resting-state. Various data lengths of ECG recordings were applied to eight (two time, two frequency, and four nonlinear) approaches, including 13 different methods (statistical and geometric methods in time domain, Welch’s and Lomb–Scargle in frequency domain, and Poincaré plot, recurrence quantification analysis (RQA), detrended fluctuation analysis (DFA), Wolf and Rosenstein’s Lyapunov exponents (LE), ApEn, SampEn, multiscale entropy (MSE), and composite multiscale entropy (CMSE) in nonlinear analyses), and 34 HRV measures (Figure 1), to determine the shortest data length that can keep the system dynamics intact. To understand the robustness of optimal minimum data length that can be utilized to quantify HRV, each data set size was compared with the most prolonged (2000 R peaks or 750 s) using the Mann–Whitney U test to determine the 0.05 level of significance.

2. Materials and Methods

2.1. Subjects

The data were collected from healthy young participants with normal ECG recordings. A total of 14 subjects (seven male) participated with the age of 23.8 ± 4.1 years (mean ± standard deviation) and body mass index (BMI) of 23.4 ± 5.1 kg/cm2. Participants were excluded if any neurological disorder or heart-related disease was reported. Only one male participant reported having a history of Kawasaki’s disease but did not present any symptoms during these data collection.

2.2. Experimental Protocol

The ECG signals were recorded with lead II placement at a sampling rate of 2 kHz from 14 subjects by BIOPAC MP36 System (BIOPAC Systems, Inc., Goleta, CA, USA). Subjects were asked to avoid taking caffeine two hours before the time of the study. Participants were provided enough rest on the chair (at least five minutes) after arrival to ensure the subjects were recorded in a relatively calm state during the studies. At least 10 min of ECG recordings were acquired while the subjects were seated on the chair. All subjects provided written informed consent for the study, which the Institutional Review Board of the University of California approved (IRB #2016-2924).

2.3. Data Preprocessing

The effect of data length on HRV measures was quantified by analyzing R-R intervals on the raw ECG recordings using time, frequency domain and nonlinear analyses methods. We segmented and standardized ECG data length using R peaks ranging from 60 to 2000, thus controlling heart rate (HR) differences [20]. In the study, three statistical (SDNN, RMSSD and pNN50) and one geometrical method (triangulation index) were included in time-domain measurements (Figure 1). Welch and Lomb–Scargle algorithms for each of the eight measures (total power, power and normalized power of very low-frequency, low-frequency, and high-frequency, and the ratio of low-frequency to high-frequency) were investigated in the frequency-domain method. Eight nonlinear methods were investigated to present different HRV measures. Besides the recurrence quantification analysis (RQA), all the measures were computed with consecutive discrete R-R intervals with R peaks ranging from 60 to 2000. RQA was estimated by raw ECG recordings ranging from 60 to 750 s. In the section below, we provided details for (i) R peak extraction, (ii) evaluation for HRV measures, and (iii) statistical analysis.

2.4. Extraction of R Peaks Using Wavelet Analysis

The consecutive discrete R-R intervals were processed initially from the raw ECG signal. Symlet 4 (Sym4) wavelet was chosen to enhance R peaks by using maximal overlap discrete wavelet transform (MODWT) due to its similarity with the QRS complex as shown in Figure 2. Potential artifacts as arrhythmic events were excluded. The identified R peaks were then extracted as R-R interval segments used for the HRV analysis.

2.5. Time-Domain Analysis

Three statistical measures (SDNN, RMSSD, pNN50) and a geometric measure (triangular index) were discussed in time-domain analysis. SDNN is the standard deviation of the normal-to-normal (NN) interval. RMSSD represents the square root of the mean of the sum of the squares of differences between adjacent NN intervals. pNN50 calculates the ratio of the counts of adjacent NN intervals that are more than 50 ms and the total number of NN intervals in the data set. The basic variable in the geometric method, HRV triangular index, is also computed (Equation (1)).
HRV   Triangular   Index   = total   number   of   NN   intervals number   of   NN   intervals   in   the   modal   bin
SDNN, RMSSD, pNN50 and HRV triangular index were computed as standards of measurement [12].

2.6. Frequency-Domain Analysis

According to the task force’s [12] guideline, frequency domain measures of HRV can be categorized into four bands, ultra-low-frequency (ULF, 0.003 Hz), very low-frequency (VLF, 0.003–0.04 Hz), low-frequency (LF, 0.04–0.15 Hz), and high-frequency (HF, 0.15–0.4 Hz). It is also suggested that VLF, LF, and HF rhythms are distinct components for a 2–5 min short-term ECG recording. Hence, there were eight parameters calculated in the study: (i) VLF power, (ii) LF power, (iii) HF power, (iv) total power, (v) normalized VLF (VLF norm), (vi) normalized LF (LF norm), (vii) normalized HF (HF norm), and (viii) the ratio of LF to HF (LF/HF). The total power is a sum of the VLF, LF, and HF absolute power. The normalized HRV power is defined as the proportion of one power range to the total power in absolute values and allows to compare individuals and various data lengths appropriately. To estimate HRV power spectral analysis, both Welch’s method [37] and Lomb–Scargle periodogram [38,39] were applied. Welch’s method utilizes fast Fourier transform (FFT) and is advantageous for low a computation workload. The original signal with N data points is split into K segments, X 1 ( j ) , , X K ( j ) . Each segment has a length of L and is apart from the previous segment in the distance of D such that N can be written as Equation (2).
N = ( K 1 ) D + L .
The overlapped segments are helpful to mitigate the loss. The windowed finite Fourier transform A k ( n ) is applied to each segment ( k = 1 ,   2 ,   ,   K ) with a data window W ( j ) where j = 0 ,   ,   L 1 as Equation (3):
A k ( n ) = 1 L j = 0 L 1 X k ( j ) W ( j ) e 2 k i j n L                 n = 0 ,   ,   L 2 ,
where X k ( j ) W ( j ) is the windowed segment sequences and i = 1 . Therefore, the modified periodogram, P k ( f n ) , for each segment, can be written as Equation (4).
P k ( f n ) = L 2 | A k ( n ) | 2 j = 0 L 1 W 2 ( j )             f n = n L             n = 0 ,   ,   L 2 .
And Welch’s method is estimated as the average of the periodogram values as Equation (5).
P ( f n ) = 1 K k = 1 K P k ( f n ) .
The Lomb–Scargle periodogram is inspired by the Fourier transform and a least-squares method known for identifying periodicity. It offers advantages in dealing with unevenly sampled data and allows ectopic or missing beats [40,41,42].

2.7. Nonlinear Methods

We investigated several nonlinear variability methods such as Poincaré plots, fractal dimension, Lyapunov exponent, and entropy in the study.

2.7.1. Poincaré Plot

The Poincaré plot is a nonlinear technique that can depict HRV in a two-dimensional graphic. It visualizes the beat-to-beat detail to dispersion and provides quantitative information on cardiac performances. Popular approaches to quantify Poincaré plot include ellipse fitting, histogram, and correlation coefficient. The study characterized the R-R interval as illustrated by the Poincaré plot by fitting an ellipse as in Figure 3 [43,44,45].
A new set of the coordinate plane, x 1 and x 2 , is formed at the intersection of the ellipse center (Equation (6)).
[ x 1 x 2 ] = [ cos θ sin θ sin θ cos θ ] [ R R n R R n + 1 ]
S D 1 and S D 2 represent the distribution of points perpendicular and parallel to the line-of-identity of the fitted ellipse, which indicates the level of short- and long-term variability [46,47]. Mathematically, S D 1 and S D 2 are the standard deviations around x 1 and x 2 , respectively. They were computed by using linear measures of HRV as shown in Equations (7) and (8) [47]:
S D 1 2 = V a r ( x 1 ) = V a r ( 1 2 R R n 1 2 R R n + 1 ) = 1 2 V a r ( R R n R R n + 1 ) = 1 2 S D S D 2
S D 2 2 = 2 S D N N 2 1 2 S D S D 2 ,
where the standard deviation of the differences between adjacent RR intervals is denoted by S D S D . The axis ratio S D 1 / S D 2 indicates the relationship of the instantaneous interval variation to the long-term variation.

2.7.2. Approximate Entropy

Approximate entropy (ApEn) is a method to quantify the complexity of time-series data. Its application is limited to data lengths greater than 100 data points [48,49]. Complexity can quantify variability to indicate the unpredictability of HR fluctuations. To assess complexity, ApEn was computed as the difference between the probability of the series of a vector with a fixed data length m and the probability of the series of another vector with a similar length ( m + 1 ) that both fall within a tolerance r as Equation (9):
A p E n ( m ,   r ,   N ) = Φ m + 1 ( r ) Φ m ( r ) ,
where Φ m ( r ) from the element in Eckmann–Ruelle (E–R) entropy [50] is defined in Equation (10):
Φ m ( r ) = 1 N m + 1 i = 1 N m + 1 l o g C i m ( r ) ,
where C i m ( r ) is the conditional probability of vector length m .
Pincus et. al. recommended that given 1000 data points, using parameters m = 2 and an r value between 0.1 and 0.25 of data standard deviation provide a robust result of ApEn [48,51]. Hence, ApEn was applied to different quantities of R-R interval time series with the values m = 2 and r = 0.2 in this study.

2.7.3. Sample Entropy

Sample entropy (SampEn) was initially introduced 9 years after ApEn [52]. With a similar goal, SampEn is a useful mathematical algorithm that measures the predictability in time series and can be viewed as a refinement of ApEn. SampEn addressed ApEn’s shortcoming as its ability to have independence for different data lengths and showed robustness for the change in data lengths [52,53,54]. Richman and Moorman defined the negative natural logarithm of the probability using vector length m , tolerance r , and signal data length N from ApEn (Equation (11)):
S a m p E n ( m ,   r ,   N ) = l n A B ,
where A and B are ( N m ) ( N m 1 ) 2 times the sum of all the conditional probabilities, C i m + 1 and C i m , divided by N m without considering the self-matches.

2.7.4. Multiscale Entropy

The traditional multiscale entropy (MSE) algorithm [55,56] is conducted in two parts: (i) a coarse-graining procedure for different scaled time series from an original signal; (ii) SampEn is used to calculate each time series scale. However, some disadvantages of traditional MSE had been issued [57]. Therefore, several newly developed MSE methods were reported [29,58,59,60]. However, most of the solutions suggested modifying the first step’s time-series scale with different algorithms or using other entropies for the second step instead [61]. Therefore, instead of the conventional MSE, a composite multiscale entropy (CMSE) [60] was used to evaluate cardiovascular complexity in the study. Traditional MSE computed SampEn with the first coarse-grained time series of each scale only (Equation (12)):
S a m p E n ( m ,   r ,   N ) = l n A B ,
where x is the original one-dimensional time series, τ is the scale factor, and the same parameters m and r from computing SampEn. Equation (13) is the first coarse-grained time series y 1 ( τ ) which is defined as:
y 1 , j ( τ ) = 1 τ i = ( j 1 ) τ + 1 j τ x i ,                           1 j N τ .
Nevertheless, the SampEn of the first coarse-grained time series derives poor reliability. The CMSE algorithm was then introduced to overcome the problem [60]. Therefore, instead of using the first coarse-grained time series, All the coarse-grained time series SampEn are considered in CMSE as shown in Equation (14):
C M S E ( x ,   τ ,   m ,   r ) = 1 τ k = 1 τ S a m p E n ( y k ( τ ) ,   m ,   r ) ,
where y k ( τ ) represents the coarse-graining procedure determined from x k , the k th data point from the original signal (Equation (15)).
y k , j ( τ ) = 1 τ i = ( j 1 ) τ + k j τ + k 1 x i ,                             1 j N τ ,   1 k τ .

2.7.5. Detrended Fluctuation Analysis

In the time series analysis, the detrended fluctuation analysis (DFA) is used to estimate self-similarity. It is calculated as the root-mean-square error of the least-squares line fitted in separate non-overlapping windows of the cumulative integral of the original time series. For example, given the interbeat intervals of an original ECG signal in time series x with length N , the cumulative integral of the signal y ( k ) is written as Equation (16):
y ( k ) = i = 1 k ( x ( i ) x ) ,
where x denotes the average of x , y ( k ) is then divided into segments of equal length n . The root-mean-square fluctuation F as a function of the window sizes n is computed by Equation (17):
F ( n ) = 1 N k = 1 N [ y ( k ) y n ( k ) ] 2 ,
where y n ( k ) is the y-coordinate of the linear fitted line in each of the segments.
The fluctuation is indicated as the slope of F ( n ) versus n in a logarithmic scale by a scaling exponent α . Peng’s work found that there is a significant difference in α over a wide range of window sizes ( 20 n 1000 ) among the interbeat interval time series of 15 severe heart failures and 12 controls [62].

2.7.6. Recurrence Quantification Analysis

The research of recurrences is commonly used to understand the complexity of a nonlinear dynamical system. The RQA is known for being able to handle short and nonstationary data. The recurrence plot (RP) and its quantification measures are the tools to visualize and quantify the recurrence behavior of the state space trajectory [63,64,65]. RP was first introduced in 1987, allowing the recurrences of higher dimensional phase space to be visualized by a two-dimensional representation [66]. It is an effective way to understand the behavior of a dynamic system. Mathematically, it shows the phase space vectors x i that recur at time i and another time j (Equation (18)):
R i , j = Θ ( ε i | | x i x j | | ) ,             x i   R m ,       i ,   j = 1 , ,   K ,
where K is considered size for the recurrence matrix, m is the highest dimension being investigated, ε i is a threshold distance, and Θ ( . ) is the Heaviside function.
To infer RQA, there are three parameters that need to be taken into consideration, a time delay τ , the embedding parameter D , and a threshold distance ε . The time-delay parameter of the Takens embedding theorem [67] brings the original one-dimensional time series into multiple dimensional manifolds. Here, the average mutual information (AMI) is used to estimate τ [68]. Essentially, AMI calculates the least dependent information in the time-delayed coordinates. Then, the embedding parameter comes in to reconstruct the phase space vector since a time delay is applied to the raw data. The choice of a deficient D may lead to unwanted results such as the false bifurcation points [68]. Hence, false nearest neighbors [69], which inspects whether there’s a significant change in the distance between two adjacent data points with embedding dimensions, was used to determine D here. In our work, the time delay and embedding parameters were chosen as the median values across 14 studies ( τ = 22 and D = 2 ) . Lastly, the threshold distance parameter defines who the RP neighbors are as the radius of a sphere. To select a sufficient threshold, distance can be critical, it may lose the key information of the recurrence structure or include a lot of artifacts if ε is chosen too small or too large [13,70]. It is suggested that a proper selection of ε should correspond to a specific range of the percent recurrence (%REC), a quantification measure is discussed in detail in the next paragraph [64,71].
In comparing various subsets of data, the recommended threshold parameter should be chosen so that a typical %REC is in the range of 5% to 10% and the minimum is at least 1% [71]. The threshold of 8 fulfilled the guideline in the study. The quantification measures are used to characterize the information in RPs. Four of them, %REC, percent determinism (%DET), the average diagonal line length (ADL), and the maximum diagonal line length (MDL), were reported in the paper. %REC quantifies the percentage of recurrent points in a RP (Equation (19)).
% REC   = sum   of   recurrent   points size   of   RP * 100 .
The minimum (0%) and maximum (100%) represent that no points and all the points are fallen into the ε -defined recurrent sphere, respectively. The second parameter, %DET, measures the percent of recurrent points occurring in connected trajectories, which is of the total counts of recurrent points (Equation (20)). The connected trajectories are formed by the continuous adjacent of two or more points that follow the diagonal lines.
% DET   = sum   of   diagonally   adjacent   recurrent   points sum   of   recurrent   points   in   RP * 100 .
ADL calculates the average length of connected trajectories. MDL simply counts the length of the longest connected trajectory in the RP.

2.7.7. Lyapunov Exponent

The Lyapunov exponent (LE) of a dynamical system is a quantity that measures how fast two infinitesimally close trajectories separate in phase space based on the initial condition. For instance, one point would exponentially diverge from another if the system is chaotic. Given two close trajectories, x ( t ) and y ( t ) are a function of time, Equation (21) shows the next iteration ε after time t separates their distance exponentially:
| x ( t + ε ) y ( t + ε ) | = | x ( t ) y ( t ) | e λ 1 ε ,  
where λ 1 represents as the first LE (Equation (22)). Hence, the first LE can be written as:
λ 1   = l i m ε l i m | x ( t ) y ( t ) | 0 1 ε l n ( | x ( t + ε ) y ( t + ε ) | | x ( t ) y ( t ) | ) .
The rate of separation, LE, may vary for different orientations of the initial two close trajectories. λ with a positive value ( λ > 0 ) means the trajectories diverge exponentially. On the contrary, two nearby points converge exponentially, which leads to a negative value of λ   ( λ < 0 ) . That is to say, the larger the value of LE, the lower the predictability for a dynamical system. The largest Lyapunov exponent (LLE) is commonly referred to as the indicator of chaos. There are various computational methods to quantify LLE. Two widely used algorithms, Wolf and Rosenstein methods, are included in the study. Both of them track the divergence of nearest neighbors over time. However, Wolf’s algorithm [72] takes one trajectory as the reference only (Equation (23)). Thus, each point on the reference trajectory iterates with its single nearest neighbor’s trajectory over time until their distance apart from each other grows beyond a threshold:
| | z ( t i ) x ( t i ) | | = L ( i ) ,
where i is the increment, L is the difference in two trajectories, x is the point on the reference trajectory, and z is the nearest neighbor of the corresponding point. The reference point re-evaluates the new single nearest neighbor once the previous trajectory’s separation is large. Following the nearest neighbor and replacing it with another trajectory completes at the end of the reference trajectory. The distance between the reference point and the beginning and last point of each new nearest neighbor trajectory is denoted as L and L , respectively. The LLE by Wolf’s algorithm, λ 1 , tracks all the L and L (Equation (24)):
λ 1   1 K i = 0 M 1 l n L i L i ,
where M represents the number of nearest neighbor trajectories and K is the total number of iterations in the reference trajectory.
Instead of focusing on a single nearest neighbor, Rosenstein’s method [73] finds the nearest neighbor overall points on the trajectory. The distance of a certain reference point, X j , to its nearest neighbor X j can be expressed as Equation (25):
d j ( 0 ) = m i n X j | | X j X j | | ,
where d j ( 0 ) means the initial distance between the j th point and the nearest neighbor. The LLE by Rosenstein’s work can be estimated as the average speed of nearest neighbor separation (Equation (26)).
d j ( i ) C j e λ 1 ( i Δ t ) ,
where Δ t is the period of each iteration, i is the number of iterations, and C i as the initial separation.
In our approach, the original one-dimensional signal was reconstructed using the methods mentioned above to define parameters, time delay τ and embedding dimension D (details in RQA section). To compare different data lengths, each τ and D in a particular data length group was obtained as the median values across all participants. The selected values of τ and D in each group are referenced in Table 1.

2.8. Statistical Analysis

A total of 34 HRV indices were computed in this study. The indices were divided into three groups: the effect of data length on (i) time domain, (ii) frequency domain, and (iii) nonlinear HRV measures. Descriptive data are presented as means and standard deviations (SDs) for continuous HRV variables. Normal distribution of all time/frequency domain, linear/nonlinear variables was tested using the Kolmogorov–Smirnov test and visually inspected histograms and Q-Q plots. HRV parameters from ECG recordings did not exhibit normal distributions and were analyzed as non-parametric. To understand the robust/optimal minimum data length that can be utilized to quantify HRV in the methods as mentioned above, each data set size was compared against the most extended (2000 R peaks or 750 s) using the Mann–Whitney U test in R. Moreover, HRV measures were then calculated with several randomly chosen segments and compared against each other using the Mann–Whitney U test to eliminate the possibility of biased results from short data length selection. The critical value was chosen at the 0.05 level of significance.

3. Results

3.1. Time-Domain HRV

Among the four time-domain HRV measures presented in the paper, we find SDNN, RMSSD, and pNN50 are consistent for a very short-term HRV analysis (Table 2). SDNN in 100 R peaks was not significantly different from SDNN at 2000 R peaks. Similarly, a minimum of 60 R peaks of RMSSD and pNN50 were found consistent with longer (2000 R peaks) HRV recordings. The HRV triangular index, however, is recommended to be used with at least 1000 R peaks. This means that if one were an adult with a normal resting HR, a 10 to 16-min recording would have no statistically significant difference compared with a 20 to 33-min recording.

3.2. Frequency-Domain HRV

Both Welch and Lomb–Scargle periodograms were included in the HRV power spectral analysis (Table 2). With Welch’s method, the appropriate minimum length for acquiring VLF power, LF power, total power, and VLF norm was 1000 R peaks. The HF norm analysis can use 750 R peaks. The HF power, LF norm, and LF/HF ratio had no change from 60 to 2000 R peaks. On the contrary, the lengths required to calculate HRV measures by the Lomb–Scargo periodogram were shorter than Welch’s in general. However, the recommended lengths to obtain VLF norm, LF norm, and LF/HF ratio remained same. Furthermore, using 200, 500, 500, and 750 R peaks for HF norm, LF power, total power, and HF power, respectively, were acceptable. VLF power with the Lomb–Scargle algorithm had no statistically significant difference using very short-term recordings.

3.3. Nonlinear HRV

The statistical results of nonlinear HRV indices were separated into two tables. The majority of nonlinear HRV indices are shown in Table 2. The measures of RQA are listed in Table 3. As a whole, the adequate lengths of nonlinearly-assessed HRV measures were longer than those of linearly-assessed methods.
Poincaré plots, when compared with different data lengths, did not show any differences with S D 1 . However, the minimum lengths of 1000 and 750 R peaks were not significantly different compared with the maximum (2000 R peaks) for S D 2 and S D 1 / S D 2 , respectively (Table 2). Furthermore, as entropy-based approaches, using 750 R peaks to quantify SampEn, ApEn, and CMSE showed no significant differences compared with maximum length (2000 R peaks); however, traditional MSE was significantly different until 750 R peaks and was robust for longer data lengths. Statistical results showed that DFA and Wolf and Rosenstein’s LE were affected due to the data length. The minimum data length for DFA and Wolf’s LE was 1500 R peaks such that no significant difference was found compared with 2000 R peaks (maximum data length). However, Rosenstein’s LE showed significantly different values in all the data lengths when compared with full data length (2000 R peaks). Four parameters of RQA (Recurrence (REC), determinism (DET), Maximum Diagonal Length (MDL) and Average Diagonal Length (ADL)) were computed for different data lengths to compare with the maximum (ECG data of 750 s or 12.5 min) data. The results did not show any significant differences among 1 min and 12.5 min of %REC, %DET, MDL and ADL.
Table 4 shows recommended minimum data length suggesting consistency and unbiased to maximum data length (2000 R peaks) using the Mann–Whitney U test. We reported time- and frequency-domain HRV measures and minimum data length with consistency. We found that most of nonlinear variables required a minimum data length of 1000 or 1500 R peaks. However, this was not the case for S D 1 of Poincaré plot and RQA measures.
We investigated the effects of data length HRV measures on 14 participants utilizing the Mann–Whitney U test. Box plots were used to summarize the distribution in each HRV measure per data length (Figure 4). Figure 4a shows how data length affects the complexity (ApEn) of the R-R interval data. With the data length increasing from 60 to 2000 R peaks, the overall value of ApEn showed a linear trend and reached a plateau around 750 R peaks. The median values and the 25–75 percentile range among 750 to 2000 R peaks remained consistent, contrasting with 60 to 500 R peaks (showed increasing linear trend). Unlike the pattern in ApEn, Rosenstein’s LE showed inconsistency while the data lengths increased. In Figure 4b, the box plots of Rosenstein’s LE show scattered values with the observation of the outlier quantity. Both MSE (Figure 4c) and CMSE (Figure 4d) had more extensive percentile ranges in shorter data lengths and a similar variation in median values after 200 R peaks.

4. Discussion

In this study, we investigated the effects of ECG data series length on the consistency of HRV parameters. Conventionally, short assessments of five minutes of ECG data are used for analysis (approximately 360 R-R intervals) [12], and in this study, we evaluated the statistical differences of HRV measures at different data lengths to maximum data length (2000 R peaks). We investigated 34 HRV measures in this study, including (a) time domain, (b) frequency domain, and (c) nonlinear analysis variables. We found that a length of 1000 R peaks or more can precisely estimate HRV for time and frequency domain variability features. Moreover, we found all variables were affected by ECG data length. For example, in the general frequency domain variables are more unstable for up to 750 R peaks of data length.

4.1. Use of HRV Measures in Pathology Differentiation

Previous studies suggested that HRV measures can be an indicator to differentiate the pathologies from the healthy controls or to predict the severity of disease. Patients with common neurodegenerative diseases such as Alzheimer’s disease and Parkinson’s disease were found to have significantly reduced time- and frequency-domain HRV measures regardless of the data length [74,75,76,77]. Valappil and coworkers discussed patients with REM sleep behavior disorder; considering premotor Parkinson’s disease, they used five-minute ECG recording and found significantly lower HRV in SDNN, pNN50, LF, HF, SD1 and SD2 in the Poincaré plot compared with controls [78]. The interest in studying the use of HRV measures in other pathologies is also widespread. For instance, linear and nonlinear variability such as complexity, measured by SampEn, reduced in patients with type 1 diabetes mellitus using 3200 R peaks [79]. Mussalo and coworkers demonstrated that both patients with mild and severe essential hypertension had lower time- and frequency- domain HRV measures computed with at least 250 R peaks in the 10-min recordings [80]. In addition, there was the study presented that a five-minute ECG recording can show differences in nonalcoholic fatty liver disease with diabetes and controls [81]. Most of the studies had focused on whether linear HRV measures can stratify pathologies. To further understand how to assess not only linear but nonlinear HRV in a short time, it is crucial to investigate the effect of data lengths.

4.2. Importance of Short Data Sets and R-R Intervals

Various studies have demonstrated very short-term (less than 5 min) or short-term (5 min) data lengths for HRV analysis [20,29,30,31,32,34,75,77,78,79,82,83,84,85]; others have vouched for using longer data lengths for HRV analysis [74,76,86]. However, the knowledge of how data length affects different HRV linear and nonlinear measures is presently unknown. This research is essential since variation in ECG recording length may result in differences in outcomes of HRV analysis in all temporal, frequency-based and linear/nonlinear analyses. When recordings with different duration are compared, it should be considered to use the most prolonged duration as a standard of comparison for stable HRV values. This allows us to identify the minimum length of ECG data that can capture system dynamics without significantly differing results from long data sets. A quick HRV analysis may serve as a promising diagnostic tool in healthcare. An effort to shorten ECG data recording is critical since HRV features add essential information for cardiac functioning. Our study highlights that most of the HRV measures are sensitive to changes in the data length. We found the sensitivity of each HRV measure was affected by the change in data lengths. It is important to note that a faster HR leads to a smaller HRV [20]. Hence, unlike most others using time as the length reference, the quantity of R peaks ranging from 60 to 2000 was used instead to avoid HR variation [87].

4.3. Linear ECG Variability Measures

Chen and coworkers conducted a study with 3387 adult participants with ECG recordings of longer than two minutes, but reported that such long ECG recording may not be required since the valid results of RMSSD and SDNN can be attained from 10 and 30 s of ECG recordings, respectively [29]. The robustness of RMSSD from 10 s recordings was also corroborated in Thong’s work [30]. However, 10 s based SDNN assessment was found inconsistent by both studies concluding that linear variability measure such as SDNN was more sensitive to data length than RMSSD. Similarly, some other studies reported similar results of RMSSD and SDNN when comparing data lengths of 3 and 5 min [20], 50-s [31] and 5-min [32,83,84] measurements as the reference. On the other hand, pNN50 evaluated from three and five minutes showed similar values [20]. Short ECG data sets of 20 s of pNN50 can reliably estimate similar to 150 s [31]. Thus, our results indicate that RMSSD and pNN50 were the least sensitive to data lengths. A shorter ECG data length of 30 s for SDNN evaluation can potentially replace existing guidelines by the task force [12]. The HRV triangular index was affected by data length, which was consistent with previous findings [83].

4.4. Frequency-Domain Analysis

The frequency-domain HRV measures were evaluated as per standard techniques defined by the task force [12]. We investigated two frequency domain HRV analysis methods using the Welch periodogram and Lomb–Scargle. Previously, Thong investigated 10-s HRV data for HF band variables and reported results unreliable for accuracy [30]. McNames and Aboy concluded that the performances of HF variable ranging from 10 s to 10 min compared with the 5-min estimation were comparable with the results in mean HR [83]. In addition, the study showed that 40-s HF and 50-s LF/HF, LF norm and HF norm were reliable to monitor mental stress under a mobile setting [31]. Similar to Salahuddin’s work, we found 60 R peaks (36- to 60-s) LF/HF and LF norm had no significant differences with 2000 R peaks using either Welch or Lomb–Scargle algorithm. In addition, the task force manual suggested that it can be inappropriate to assess VLF in short-term recordings (≤5 min). However, our findings show that the optimum data length to estimate VLF depends on the methods of power spectral density (for example, 1000 and 60 R peaks in using the Welch and Lomb–Scargle algorithm, respectively).

4.5. Nonlinear Variability Analysis

Most of the nonlinear methods were proposed 30 years ago. However, there are not literatures investigating the sensitivity on data length as the linear HRV analysis. A shorter data length of Poincaré plot, MSE, and CMSE were reported and discussed for different purposes without providing any suggestions to minimum data length for reliable measurements. For instance, the short-term assessment of the Poincaré plot was applied in different stress levels, yet was concluded as promising results [88,89]. For an entropy-based HRV analysis, SampEn and ApEn were compared and discussed together in most instances since SampEn was introduced to improve the unreliable outcome of ApEn due to data length. It was suggested that a minimum data length of 100 and 250 RR intervals of SampEn and ApEn can distinguish healthy from congestive heart failure patients [34]. Another group studied in the range of 2 min to the 20 min data length. The authors concluded that SampEn was a lot less sensitive to data length compared with ApEn [35]. McNames and Aboy considered several time domain and frequency domain HRV variables with ApEn and indicated that ApEn was the most unreliable one [83]. Although, SampEn was reported independent of data length compared with ApEn. However, we found 1000 R peaks to be optimum for estimating ApEn or SampEn.
The results of DFA ( α ) represent the relationship of F ( n ) and the window sizes n . Obviously, n cannot be larger than the length of the data set. Therefore, the data length is an essential factor of the accuracy of α . Moreover, a crossover phenomenon was observed when DFA was proposed [62]. Therefore, Peng suggested that at least a 24-h recording was required for diagnostic purposes since the crossover phenomena can play an essential role consistent with our results.
To calculate LLE, the parameters τ and D need to be defined first. In this study, both τ and D increases with the data lengths in HRV. Gao’s works suggested that D = 2 should be used when analyzing a finite HRV data set [90,91,92]. However, others reported larger parameters for longer data sets and shorter parameters for shorter data sets. For instance, Signorini and Cerutti calculated long-term HRV ( N = 20 , 000 ) with τ = 7 and D = 10 [86]. Li’s group used τ = 1 and D = 3 for less than 5 min data sets ( N = 200 355 ) [85]. Moreover, regarding the method differences in LE, Rosenstein’s LE is fast and easy to implement and applicable to small data sets [73]. A minimum of 200 consecutive R-R intervals was suggested as the optimum data length for calculating HRV [85]. However, Li and coworkers did not report any statistical evidence supporting their conclusion. On the contrary, we found that data length is significantly different from LLE of 2000 R peaks.
The varying data lengths required for nonlinear and chaos HR analyses can be partially explained by different autonomic heart control system aspects that complexity (ApEn, SampEn and MSE) and LLE can measure. For instance, MSE can measure deviations or differences at different time scales thus offering insightful information on temporal dynamical variations in the autonomic HR control system. In our laboratory pilot studies we found strenuous exercises significantly decreased the complexity of HR R-R interval data. However, LLE detects the presence of chaos in heart dynamical control systems by quantifying LEs (exponential divergence of initially close state space trajectories). The computation of LLE utilizing multiple LEs derived from divergence curves requires larger data sets for stable values. Additionally, DFA is an indicator of statistical persistence and antipersistence of HR time series. Persistence indicates the deviation in HR time series is statistically more likely followed by subsequent deviation in the same direction (increase in HR is followed by a subsequent increase in HR or decrease is followed by another decline of HR). On the other hand, antipersistence implies that the deviation is followed by subsequent deviation in the opposite direction (increase in HR is followed by a decrease in HR and vice-versa). Since supraspinal mechanisms involuntarily control HR, it is likely the HR control mechanism will naturally produce long-range correlated HR time series, thus requiring at least 1500 data points for robust and stable DFA values.

4.6. Limitations

The findings of this study must be seen together with the limitations. Firstly, our study is limited with sample size. Additionally, the subjects were within a limited range of age thus the study is limited in external validity to other age ranges. We will conduct a study with a larger number of participants and broader age ranges in the future. Secondly, this study investigated optimum data length for HRV measures only limited to a healthy group. Hence, a future study with pathological groups would be interesting to embolden our findings.

5. Conclusions

An ECG-based analysis of cardiac rhythm is critical for the diagnosis of a heart condition and disease management. In addition, novel clinical decision support systems require quick ECG analysis to assist clinicians. Our effort to shorten the ECG recording duration is vital to improve efficiency and save time and effort for patients and clinical care providers. This can be more critical for patients with frequent artifacts and HRV physiological features extracted from short ECG recordings with high confidence. Our study suggests that ECG data length collected from wearable devices must be optimized and selected such that more consistent and reliable results can be attained with existing laboratory-grade measurements. In conclusion, this study suggests ultra-short data sequences can be collected and analyzed for quick HRV assessments retaining the rich information from linear/ non-linear variability structure, but with caution since HRV variables are affected differentially to the data length. Chaotic HR analyses such as LLE (Rosenstein and Wolf) and long-range correlation through DFA required longer data set lengths than other nonlinear variability measures such as ApEn, SampEn, and MSE.

Author Contributions

T.L., M.K., R.S. conceived the idea and M.K. and E.-F.C. designed the experiments; E.-F.C. performed the experiments; R.S. along with E.-F.C. analyzed the data wrote the manuscript with support from T.L. and M.K. All authors discussed the results and contributed to the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of University of California Irvine (IRB Number 2016-2924 approved on 2 October 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Malpas, S.C.; Whiteside, E.A.; Maling, T.J. Heart Rate Variability and Cardiac Autonomic Function in Men with Chronic Alcohol Dependence. Heart 1991, 65, 84–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kudat, H.; Akkaya, V.; Sozen, A.; Salman, S.; Demirel, S.; Ozcan, M.; Atilgan, D.; Yilmaz, M.; Guven, O. Heart Rate Variability in Diabetes Patients. J. Int. Med. Res. 2006, 34, 291–296. [Google Scholar] [CrossRef] [PubMed]
  3. Stein, P.K.; Reddy, A. Non-Linear Heart Rate Variability and Risk Stratification in Cardiovascular Disease. Indian Pacing Electrophysiol. J. 2005, 5, 210–220. [Google Scholar]
  4. Thayer, J.F.; Yamamoto, S.S.; Brosschot, J.F. The Relationship of Autonomic Imbalance, Heart Rate Variability and Cardiovascular Disease Risk Factors. Int. J. Cardiol. 2010, 141, 122–131. [Google Scholar] [CrossRef]
  5. Sajadieh, A. Increased Heart Rate and Reduced Heart-Rate Variability Are Associated with Subclinical Inflammation in Middle-Aged and Elderly Subjects with No Apparent Heart Disease. Eur. Heart J. 2004, 25, 363–370. [Google Scholar] [CrossRef]
  6. Lampert, R.; Bremner, J.D.; Su, S.; Miller, A.; Lee, F.; Cheema, F.; Goldberg, J.; Vaccarino, V. Decreased Heart Rate Variability Is Associated with Higher Levels of Inflammation in Middle-Aged Men. Am. Heart J. 2008, 156, 759.e1–759.e7. [Google Scholar] [CrossRef] [Green Version]
  7. Williams, D.P.; Koenig, J.; Carnevali, L.; Sgoifo, A.; Jarczok, M.N.; Sternberg, E.M.; Thayer, J.F. Heart Rate Variability and Inflammation: A Meta-Analysis of Human Studies. Brain Behav. Immun. 2019, 80, 219–226. [Google Scholar] [CrossRef]
  8. Zahorska-Markiewicz, B.; Kuagowska, E.; Kucio, C.; Klin, M. Heart Rate Variability in Obesity. Int. J. Obes. Relat. Metab. Disord. 1993, 17, 21–23. [Google Scholar]
  9. Gorman, J.M.; Sloan, R.P. Heart Rate Variability in Depressive and Anxiety Disorders. Am. Heart J. 2000, 140, S77–S83. [Google Scholar] [CrossRef]
  10. Chalmers, J.A.; Quintana, D.S.; Abbott, M.J.-A.; Kemp, A.H. Anxiety Disorders Are Associated with Reduced Heart Rate Variability: A Meta-Analysis. Front. Psychiatry 2014, 5, 80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Shaffer, F.; Ginsberg, J.P. An Overview of Heart Rate Variability Metrics and Norms. Front. Public Health 2017, 5, 258. [Google Scholar] [CrossRef] [Green Version]
  12. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart Rate Variability: Standards of Measurement, Physiological Interpretation, and Clinical Use. Circulation 1996, 93, 1043–1065. [Google Scholar] [CrossRef] [Green Version]
  13. Henriques, T.; Ribeiro, M.; Teixeira, A.; Castro, L.; Antunes, L.; Costa-Santos, C. Nonlinear Methods Most Applied to Heart-Rate Time Series: A Review. Entropy 2020, 22, 309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rhea, C.K.; Silver, T.A.; Hong, S.L.; Ryu, J.H.; Studenka, B.E.; Hughes, C.M.L.; Haddad, J.M. Noise and Complexity in Human Postural Control: Interpreting the Different Estimations of Entropy. PLoS ONE 2011, 6, e17696. [Google Scholar] [CrossRef] [PubMed]
  15. Singh, B.; Singh, M.; Banga, V.K. Sample Entropy Based HRV: Effect of ECG Sampling Frequency. Biomed. Sci. Eng. 2014, 2, 68–72. [Google Scholar] [CrossRef]
  16. McCamley, J.; Denton, W.; Arnold, A.; Raffalt, P.; Yentes, J. On the Calculation of Sample Entropy Using Continuous and Discrete Human Gait Data. Entropy 2018, 20, 764. [Google Scholar] [CrossRef] [Green Version]
  17. Raffalt, P.C.; McCamley, J.; Denton, W.; Yentes, J.M. Sampling Frequency Influences Sample Entropy of Kinematics during Walking. Med. Biol. Eng. Comput. 2019, 57, 759–764. [Google Scholar] [CrossRef] [PubMed]
  18. Ramdani, S.; Bouchara, F.; Lagarde, J. Influence of Noise on the Sample Entropy Algorithm. Chaos 2009, 19, 013123. [Google Scholar] [CrossRef]
  19. Casaleggio, A.; Braiotta, S. Estimation of Lyapunov Exponents of ECG Time Series—The Influence of Parameters. Chaos Solitons Fractals 1997, 8, 1591–1599. [Google Scholar] [CrossRef]
  20. Xinnian, C.; Solomon, I.C.; Chon, K.H. Comparison of the Use of Approximate Entropy and Sample Entropy: Applications to Neural Respiratory Signal. In Proceedings of the 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference, Shanghai, China, 17–18 January 2006; pp. 4212–4215. [Google Scholar]
  21. Kaffashi, F.; Foglyano, R.; Wilson, C.G.; Loparo, K.A. The Effect of Time Delay on Approximate & Sample Entropy Calculations. Phys. D Nonlinear Phenom. 2008, 237, 3069–3074. [Google Scholar] [CrossRef]
  22. Singh, B.; Singh, D. Effect of Threshold Value r on Multiscale Entropy Based Heart Rate Variability. Cardiovasc. Eng. Tech. 2012, 3, 211–216. [Google Scholar] [CrossRef]
  23. Estrada, L.; Torres, A.; Sarlabous, L.; Jané, R. Influence of Parameter Selection in Fixed Sample Entropy of Surface Diaphragm Electromyography for Estimating Respiratory Activity. Entropy 2017, 19, 460. [Google Scholar] [CrossRef]
  24. Stergious, N. Nonlinear Analysis for Human Movement Variability, 1st ed.; CRC Press: Boca Raton, FL, USA, 2018; ISBN 978-1-315-36237-3. [Google Scholar]
  25. Yokus, M.A.; Jur, J.S. Fabric-Based Wearable Dry Electrodes for Body Surface Biopotential Recording. IEEE Trans. Biomed. Eng. 2016, 63, 423–430. [Google Scholar] [CrossRef] [PubMed]
  26. Arquilla, K.; Webb, A.; Anderson, A. Textile Electrocardiogram (ECG) Electrodes for Wearable Health Monitoring. Sensors 2020, 20, 1013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Crosby, J. Development of a Flexible Printed Paper-Based Battery; Western Michigan University: Kalamazoo, MI, USA, 2020. [Google Scholar]
  28. Smith, A.-L.; Owen, H.; Reynolds, K.J. Heart Rate Variability Indices for Very Short-Term (30 Beat) Analysis. Part 1: Survey and Toolbox. J. Clin. Monit. Comput. 2013, 27, 569–576. [Google Scholar] [CrossRef]
  29. Munoz, M.L.; van Roon, A.; Riese, H.; Thio, C.; Oostenbroek, E.; Westrik, I.; de Geus, E.J.C.; Gansevoort, R.; Lefrandt, J.; Nolte, I.M.; et al. Validity of (Ultra-)Short Recordings for Heart Rate Variability Measurements. PLoS ONE 2015, 10, e0138921. [Google Scholar] [CrossRef] [Green Version]
  30. Thong, T.; Li, K.; McNames, J.; Aboy, M.; Goldstein, B. Accuracy of Ultra-Short Heart Rate Variability Measures. In Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE Cat. No.03CH37439), Cancun, Mexico, 17–21 September 2003; pp. 2424–2427. [Google Scholar]
  31. Salahuddin, L.; Cho, J.; Jeong, M.G.; Kim, D. Ultra Short Term Analysis of Heart Rate Variability for Monitoring Mental Stress in Mobile Settings. In Proceedings of the 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Lyon, France, 22–26 August 2007; pp. 4656–4659. [Google Scholar]
  32. Nussinovitch, U.; Elishkevitz, K.P.; Katz, K.; Nussinovitch, M.; Segev, S.; Volovitz, B.; Nussinovitch, N. Reliability of Ultra-Short ECG Indices for Heart Rate Variability: Ultra-Short HRV Reliability. Ann. Noninvasive Electrocardiol. 2011, 16, 117–122. [Google Scholar] [CrossRef]
  33. Choi, W.-J.; Lee, B.-C.; Jeong, K.-S.; Lee, Y.-J. Minimum Measurement Time Affecting the Reliability of the Heart Rate Variability Analysis. Korean J. Health Promot. 2017, 17, 269. [Google Scholar] [CrossRef] [Green Version]
  34. Graff, B.; Graff, G.; Kaczkowska, A. Entropy Measures of Heart Rate Variability for Short ECG Datasets in Patients with Congestive Heart Failure. Acta Phys. Pol. B Proc. Suppl. 2012, 5, 153. [Google Scholar] [CrossRef]
  35. Singh, M.; Singh, B.; Singh, G. Optimal RR-Interval Data Length for Entropy Based Heart Rate Variability Analysis. IJCA 2015, 123, 39–42. [Google Scholar] [CrossRef]
  36. Lee, D.-Y.; Choi, Y.-S. Multiscale Distribution Entropy Analysis of Short-Term Heart Rate Variability. Entropy 2018, 20, 952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Welch, P. The Use of Fast Fourier Transform for the Estimation of Power Spectra: A Method Based on Time Averaging over Short, Modified Periodograms. IEEE Trans. Audio Electroacoust. 1967, 15, 70–73. [Google Scholar] [CrossRef] [Green Version]
  38. Lomb, N.R. Least-Squares Frequency Analysis of Unequally Spaced Data. Astrophys. Space Sci. 1976, 39, 447–462. [Google Scholar] [CrossRef]
  39. Estévez, M.; Machado, C.; Leisman, G.; Estévez-Hernández, T.; Arias-Morales, A.; Machado, A.; Montes-Brown, J. Spectral Analysis of Heart Rate Variability. Int. J. Disabil. Hum. Dev. 2016, 15, 5–17. [Google Scholar] [CrossRef] [Green Version]
  40. Press, W.H.; Flannery, B.P.; Teukolsky, A.A.; Vetterling, W.T. Numerical Recipes in C: The Art of Scientific Computing, 2nd ed.; Cambridge University Press: New York, NY, USA, 1992; ISBN 978-0-521-43108-8. [Google Scholar]
  41. Moody, G.B. Spectral Analysis of Heart Rate without Resampling. In Proceedings of the Proceedings of Computers in Cardiology Conference, London, UK, 5–8 September 1993; pp. 715–718. [Google Scholar]
  42. Fonseca, D.S.; Netto, A.D.; Ferreira, R.B.; de Sa, A.M.F.L.M. Lomb-Scargle Periodogram Applied to Heart Rate Variability Study. In Proceedings of the 2013 ISSNIP Biosignals and Biorobotics Conference: Biosignals and Robotics for Better and Safer Living (BRC), Rio de Janerio, Brazil, 18–20 February 2013; pp. 1–4. [Google Scholar]
  43. Marciano, F.; Migaux, M.L.; Acanfora, D.; Furgi, G.; Rengo, F. Quantification of Poincare’ Maps for the Evaluation of Heart Rate Variability. In Proceedings of the Computers in Cardiology 1994, Bethesda, MD, USA, 25–28 September 1994; pp. 577–580. [Google Scholar]
  44. Tulppo, M.P.; Makikallio, T.H.; Takala, T.E.; Seppanen, T.; Huikuri, H.V. Quantitative Beat-to-Beat Analysis of Heart Rate Dynamics during Exercise. Am. J. Physiol. Heart Circ. Physiol. 1996, 271, H244–H252. [Google Scholar] [CrossRef] [PubMed]
  45. D’Addio, G.; Acanfora, D.; Pinna, G.; Maestri, R.; Furgi, G.; Picone, C.; Rengo, F. Reproducibility of Short- and Long-Term Poincare Plot Parameters Compared with Frequency-Domain HRV Indexes in Congestive Heart Failure. In Proceedings of the Computers in Cardiology (Cat. No.98CH36292), Cleveland, OH, USA, 13–16 September 1998; Volume 25, pp. 381–384. [Google Scholar]
  46. Kamen, P.W.; Krum, H.; Tonkin, A.M. Poincaré Plot of Heart Rate Variability Allows Quantitative Display of Parasympathetic Nervous Activity in Humans. Clin. Sci. 1996, 91, 201–208. [Google Scholar] [CrossRef] [Green Version]
  47. Brennan, M.; Palaniswami, M.; Kamen, P. New Insights into the Relationship between Poincare Plot Geometry and Linear Measures of Heart Rate Variability. In Proceedings of the 2001 Conference Proceedings of the 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Istanbul, Turkey, 25–28 October 2001; Volume 1, pp. 526–529. [Google Scholar]
  48. Pincus, S.M. Approximate Entropy as a Measure of System Complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [Green Version]
  49. Pincus, S. Approximate Entropy (ApEn) as a Complexity Measure. Chaos 1995, 5, 110–117. [Google Scholar] [CrossRef]
  50. Eckmann, J.-P.; Ruelle, D. Ergodic theory of chaos and strange attractors. In The Theory of Chaotic Attractors; Hunt, B.R., Li, T.-Y., Kennedy, J.A., Nusse, H.E., Eds.; Springer: New York, NY, USA, 1985; pp. 273–312. ISBN 978-1-4419-2330-1. [Google Scholar]
  51. Pincus, S.M.; Cummins, T.R.; Haddad, G.G. Heart Rate Control in Normal and Aborted-SIDS Infants. Am. J. Physiol. Regul. Integr. Comp. Physiol. 1993, 264, R638–R646. [Google Scholar] [CrossRef]
  52. Richman, J.S.; Moorman, J.R. Physiological Time-Series Analysis Using Approximate Entropy and Sample Entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Yentes, J.M.; Hunt, N.; Schmid, K.K.; Kaipust, J.P.; McGrath, D.; Stergiou, N. The Appropriate Use of Approximate Entropy and Sample Entropy with Short Data Sets. Ann. Biomed. Eng. 2013, 41, 349–365. [Google Scholar] [CrossRef]
  54. Montesinos, L.; Castaldo, R.; Pecchia, L. On the Use of Approximate Entropy and Sample Entropy with Centre of Pressure Time-Series. J. NeuroEng. Rehabil. 2018, 15, 116. [Google Scholar] [CrossRef] [Green Version]
  55. Costa, M.; Goldberger, A.L.; Peng, C.-K. Multiscale Entropy Analysis of Complex Physiologic Time Series. Phys. Rev. Lett. 2002, 89, 068102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Costa, M.; Goldberger, A.L.; Peng, C.-K. Multiscale Entropy Analysis of Biological Signals. Phys. Rev. E 2005, 71, 021906. [Google Scholar] [CrossRef] [Green Version]
  57. Faes, L.; Porta, A.; Javorka, M.; Nollo, G. Efficient Computation of Multiscale Entropy over Short Biomedical Time Series Based on Linear State-Space Models. Complexity 2017, 2017, 1–13. [Google Scholar] [CrossRef]
  58. Amoud, H.; Snoussi, H.; Hewson, D.; Doussot, M.; Duchene, J. Intrinsic Mode Entropy for Nonlinear Discriminant Analysis. IEEE Signal Process. Lett. 2007, 14, 297–300. [Google Scholar] [CrossRef]
  59. Valencia, J.F.; Porta, A.; Vallverdu, M.; Claria, F.; Baranowski, R.; Orlowska-Baranowska, E.; Caminal, P. Refined Multiscale Entropy: Application to 24-h Holter Recordings of Heart Period Variability in Healthy and Aortic Stenosis Subjects. IEEE Trans. Biomed. Eng. 2009, 56, 2202–2213. [Google Scholar] [CrossRef] [PubMed]
  60. Wu, S.-D.; Wu, C.-W.; Lin, S.-G.; Wang, C.-C.; Lee, K.-Y. Time Series Analysis Using Composite Multiscale Entropy. Entropy 2013, 15, 1069–1084. [Google Scholar] [CrossRef] [Green Version]
  61. Humeau-Heurtier, A. The Multiscale Entropy Algorithm and Its Variants: A Review. Entropy 2015, 17, 3110–3123. [Google Scholar] [CrossRef] [Green Version]
  62. Peng, C.-K.; Havlin, S.; Stanley, H.E.; Goldberger, A.L. Quantification of Scaling Exponents and Crossover Phenomena in Nonstationary Heartbeat Time Series. Chaos 1995, 5, 82–87. [Google Scholar] [CrossRef]
  63. Zbilut, J.P.; Webber, C.L. Embeddings and Delays as Derived from Quantification of Recurrence Plots. Phys. Lett. A 1992, 171, 199–203. [Google Scholar] [CrossRef]
  64. Webber, C.L.; Zbilut, J.P. Dynamical Assessment of Physiological Systems and States Using Recurrence Plot Strategies. J. Appl. Physiol. 1994, 76, 965–973. [Google Scholar] [CrossRef] [PubMed]
  65. Trulla, L.L.; Giuliani, A.; Zbilut, J.P.; Webber, C.L. Recurrence Quantification Analysis of the Logistic Equation with Transients. Phys. Lett. A 1996, 223, 255–260. [Google Scholar] [CrossRef]
  66. Eckmann, J.-P.; Oliffson, S.K.; David, R. Recurrence plots of dynamical systems. In World Scientific Series on Nonlinear Science Series A; World Scientific Publishing Company: Singapore, 1995; pp. 441–446. ISBN 1793-1010. [Google Scholar]
  67. Takens, F. Detecting strange attractors in turbulence. In Dynamical Systems and Turbulence, Warwick 1980; Rand, D., Young, L.-S., Eds.; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1981; Volume 898, pp. 366–381. ISBN 978-3-540-11171-9. [Google Scholar]
  68. Gao, J.; Cai, H. On the Structures and Quantification of Recurrence Plots. Phys. Lett. A 2000, 270, 75–87. [Google Scholar] [CrossRef]
  69. Kennel, M.B.; Brown, R.; Abarbanel, H.D.I. Determining Embedding Dimension for Phase-Space Reconstruction Using a Geometrical Construction. Phys. Rev. A 1992, 45, 3403–3411. [Google Scholar] [CrossRef] [Green Version]
  70. Marwan, N.; Carmenromano, M.; Thiel, M.; Kurths, J. Recurrence Plots for the Analysis of Complex Systems. Phys. Rep. 2007, 438, 237–329. [Google Scholar] [CrossRef]
  71. Wallot, S. Recurrence Quantification Analysis of Processes and Products of Discourse: A Tutorial in R. Discourse Process. 2017, 54, 382–405. [Google Scholar] [CrossRef] [Green Version]
  72. Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov Exponents from a Time Series. Phys. D Nonlinear Phenom. 1985, 16, 285–317. [Google Scholar] [CrossRef] [Green Version]
  73. Rosenstein, M.T.; Collins, J.J.; De Luca, C.J. A Practical Method for Calculating Largest Lyapunov Exponents from Small Data Sets. Phys. D Nonlinear Phenom. 1993, 65, 117–134. [Google Scholar] [CrossRef]
  74. Zulli, R.; Nicosia, F.; Borroni, B.; Agosti, C.; Prometti, P.; Donati, P.; Vecchi, M.; Romanelli, G.; Grassi, V.; Padovani, A. QT Dispersion and Heart Rate Variability Abnormalities in Alzheimer’s Disease and in Mild Cognitive Impairment: Cardiovascular abnormalities in alzheimer’s disease and MCI. J. Am. Geriatr. Soc. 2005, 53, 2135–2139. [Google Scholar] [CrossRef]
  75. Giubilei, F.; Strano, S.; Imbimbo, B.P.; Tisei, P.; Calcagnini, G.; Lino, S.; Frontoni, M.; Santini, M.; Fieschi, C. Cardiac Autonomic Dysfunction in Patients with Alzheimer Disease: Possible Pathogenetic Mechanisms. Alzheimer Dis. Assoc. Disord. 1998, 12, 356–361. [Google Scholar] [CrossRef] [PubMed]
  76. Ke, J.-Q.; Shao, S.-M.; Zheng, Y.-Y.; Fu, F.-W.; Zheng, G.-Q.; Liu, C.-F. Sympathetic Skin Response and Heart Rate Variability in Predicting Autonomic Disorders in Patients with Parkinson Disease. Medicine 2017, 96, e6523. [Google Scholar] [CrossRef]
  77. Kallio, M.; Suominen, K.; Bianchi, A.M.; Mäkikallio, T.; Haapaniemi, T.; Astafiev, S.; Sotaniemi, K.A.; Myllylä, V.V.; Tolonen, U. Comparison of Heart Rate Variability Analysis Methods in Patients with Parkinson’s Disease. Med. Biol. Eng. Comput. 2002, 40, 408–414. [Google Scholar] [CrossRef] [PubMed]
  78. Valappil, R.A.; Black, J.E.; Broderick, M.J.; Carrillo, O.; Frenette, E.; Sullivan, S.S.; Goldman, S.M.; Tanner, C.M.; Langston, J.W. Exploring the Electrocardiogram as a Potential Tool to Screen for Premotor Parkinson’s Disease. Mov. Disord. 2010, 25, 2296–2303. [Google Scholar] [CrossRef]
  79. Javorka, M.; Trunkvalterova, Z.; Tonhajzerova, I.; Javorkova, J.; Javorka, K.; Baumert, M. Short-Term Heart Rate Complexity Is Reduced in Patients with Type 1 Diabetes Mellitus. Clin. Neurophysiol. 2008, 119, 1071–1081. [Google Scholar] [CrossRef] [PubMed]
  80. Mussalo, H.; Vanninen, E.; Ikäheimo, R.; Laitinen, T.; Laakso, M.; Länsimies, E.; Hartikainen, J. Heart Rate Variability and Its Determinants in Patients with Severe or Mild Essential Hypertension: HRV and Its Determinants in Severe and Mild Hypertension. Clin. Physiol. 2001, 21, 594–604. [Google Scholar] [CrossRef]
  81. Kumar, M.S.; Singh, A.; Jaryal, A.K.; Ranjan, P.; Deepak, K.K.; Sharma, S.; Lakshmy, R.; Pandey, R.M.; Vikram, N.K. Cardiovascular Autonomic Dysfunction in Patients of Nonalcoholic Fatty Liver Disease. Int. J. Hepatol. 2016, 2016, 1–8. [Google Scholar] [CrossRef]
  82. Nguyen Phuc Thu, T.; Hernández, A.I.; Costet, N.; Patural, H.; Pichot, V.; Carrault, G.; Beuchée, A. Improving Methodology in Heart Rate Variability Analysis for the Premature Infants: Impact of the Time Length. PLoS ONE 2019, 14, e0220692. [Google Scholar] [CrossRef] [Green Version]
  83. McNames, J.; Aboy, M. Reliability and Accuracy of Heart Rate Variability Metrics versus ECG Segment Duration. Med. Biol. Eng. Comput. 2006, 44, 747–756. [Google Scholar] [CrossRef]
  84. Baek, H.J.; Cho, C.-H.; Cho, J.; Woo, J.-M. 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]
  85. Li, L.; Liu, C.; Liu, C.; Zhang, Q.; Li, B. Physiological Signal Variability Analysis Based on the Largest Lyapunov Exponent. In Proceedings of the 2009 2nd International Conference on Biomedical Engineering and Informatics, Tianjin, China, 17–19 October 2009; pp. 1–5. [Google Scholar]
  86. Signorini, M.G.; Cerutti, S. Lyapunov Exponents Calculated from Heart Rate Variability Time Series. In Proceedings of the 16th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Baltimore, MD, USA, 3–6 November 1994; pp. 119–120. [Google Scholar]
  87. Lin, G.-H.; Chang, Y.-H.; Lin, K.-P. Comparison of Heart Rate Variability Measured by ECG in Different Signal Lengths. J. Med. Biol. Eng. 2005, 25, 67–71. [Google Scholar]
  88. Mayya, S.; Jilla, V.; Tiwari, V.N.; Nayak, M.M.; Narayanan, R. Continuous Monitoring of Stress on Smartphone Using Heart Rate Variability. In Proceedings of the 2015 IEEE 15th International Conference on Bioinformatics and Bioengineering (BIBE), Belgrade, Serbia, 2–4 November 2015; pp. 1–5. [Google Scholar]
  89. Pereira, T.; Almeida, P.R.; Cunha, J.P.S.; Aguiar, A. Heart Rate Variability Metrics for Fine-Grained Stress Level Assessment. Comput. Methods Programs Biomed. 2017, 148, 71–80. [Google Scholar] [CrossRef]
  90. Gao, J.B.; Hu, J.; Tung, W.W.; Cao, Y.H. Distinguishing Chaos from Noise by Scale-Dependent Lyapunov Exponent. Phys. Rev. E 2006, 74, 066204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Gao, J.; Cao, Y.; Tung, W.-W.; Hu, J. Multiscale Analysis of Complex Time Series: Integration of Chaos and Random Fractal Theory, and Beyong; John Wiley & Sons: Hoboken, NJ, USA, 2007; ISBN 978-0-471-65470-4. [Google Scholar]
  92. Hu, J.; Gao, J.; Tung, W. Characterizing Heart Rate Variability by Scale-Dependent Lyapunov Exponent. Chaos 2009, 19, 028506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Approaches, methods, and outputted measures used to calculate HRV in the study.
Figure 1. Approaches, methods, and outputted measures used to calculate HRV in the study.
Sensors 21 06286 g001
Figure 2. R peak detection technique. (a) Sym4 resembles the QRS complex that can be used for the wavelet transform. (b) A representative raw ECG signal with extracted R peaks in red circles.
Figure 2. R peak detection technique. (a) Sym4 resembles the QRS complex that can be used for the wavelet transform. (b) A representative raw ECG signal with extracted R peaks in red circles.
Sensors 21 06286 g002
Figure 3. A representative Poincaré plot with a new set of a coordinate plane. The axes of the plane are x 1 and x 2 . S D 1 and S D 2 represent the radii of a fitted ellipse on x 1 - and x 2 -axis.
Figure 3. A representative Poincaré plot with a new set of a coordinate plane. The axes of the plane are x 1 and x 2 . S D 1 and S D 2 represent the radii of a fitted ellipse on x 1 - and x 2 -axis.
Sensors 21 06286 g003
Figure 4. Bar plots of (a) ApEn, (b) Rosenstein’s LE, (c) MSE, and (d) CMSE with different numbers of R peaks. Different R-peaks are represented with different colors for four nonlinear methods. The error bars represent standard deviation (SD) of the values among 14 participants.
Figure 4. Bar plots of (a) ApEn, (b) Rosenstein’s LE, (c) MSE, and (d) CMSE with different numbers of R peaks. Different R-peaks are represented with different colors for four nonlinear methods. The error bars represent standard deviation (SD) of the values among 14 participants.
Sensors 21 06286 g004
Table 1. The values of time delay and embedding dimension used in different data lengths.
Table 1. The values of time delay and embedding dimension used in different data lengths.
Length (R Peaks)Time DelayEmbedding Dimension
6022
10032
15023
20033
30023
40034
50034
75034
100044
150034
200055
Table 2. Mann–Whitney U test results for comparing HRV measures at 2000 R peaks with shorter data lengths. Data length is in R peaks. Statistical significant differences (p < 0.05) and statistical highly significant differences (p < 0.001) are color labeled in lighter and darker gray with bold font, respectively.
Table 2. Mann–Whitney U test results for comparing HRV measures at 2000 R peaks with shorter data lengths. Data length is in R peaks. Statistical significant differences (p < 0.05) and statistical highly significant differences (p < 0.001) are color labeled in lighter and darker gray with bold font, respectively.
Length6010015020030040050075010001500
Time-domain HRV
  Geometric measure
Triangular index0.0000.0020.0040.0040.0050.0080.0170.0480.1480.800
  Statistical measure
SDNN0.0160.0560.0560.0690.0560.0620.1040.0940.2650.946
RMSSD0.9820.8040.6670.7340.7340.7010.7010.6350.7691.000
pNN500.9460.9090.7300.8720.8360.7650.9090.8360.9090.982
Frequency-domain HRV
  Welch’s periodogram
VLF0.0000.0020.0020.0020.0010.0030.0060.0090.1040.734
LF0.0350.0500.1040.0560.0440.0480.0440.0440.2100.839
HF0.6030.8040.9100.8390.8740.9100.8740.9460.9820.982
Total power0.0160.0390.0500.0500.0440.0350.0570.0500.2100.910
VLF norm0.0000.0000.0000.0000.0000.0010.0050.0090.1030.646
LF norm0.9460.7010.7340.9270.8040.7690.7480.9460.7340.890
HF norm0.0080.0030.0080.0140.0190.0240.0440.0850.1370.734
LF/HF0.1370.0620.0850.0620.0690.0850.1130.1830.3060.839
  Lomb–Scargle’s periodogram
VLF0.9450.1210.1880.1050.2560.3060.4181.0000.0690.728
LF0.0000.0000.0060.0040.0130.0500.1880.1880.7980.694
HF0.0000.0000.0010.0010.0030.0110.0300.1120.9820.963
Total power0.0000.0000.0010.0000.0010.0130.0730.1400.6450.890
VLF norm0.0000.0000.0000.0020.0020.0100.0560.0400.3580.804
LF norm0.1370.4350.8390.8071.0000.9820.9100.7690.8900.982
HF norm0.0270.0770.0480.0620.0810.0940.1780.1980.5110.818
LF/HF0.8390.5410.4820.5110.5200.4010.5200.8040.7340.982
Nonlinear HRV
  Poincaré plot
SD10.6671.0000.8900.9630.9080.8740.9100.7690.8540.910
SD20.0090.0310.0210.0270.0140.0140.0310.0440.1980.910
SD1/SD20.0040.0020.0080.0110.0090.0140.0390.0850.1500.667
  Entropy
SampEn0.0460.0090.0060.0070.0350.0210.0310.0770.3290.734
ApEn0.0000.0000.0000.0000.0000.0000.0110.5410.7010.511
MSE0.0050.1250.0120.0030.0020.0030.0070.0350.1640.839
CMSE0.0000.0020.0070.0020.0090.0040.0090.0620.1640.839
  Fractal
DFA0.0000.0000.0000.0000.0000.0000.0020.0000.0020.667
Lyapunov exponent
Wolf0.0000.0000.0000.0350.0000.0010.0690.0850.0020.701
Rosenstein0.0000.0000.0000.0000.0000.0000.0000.0000.0020.000
Table 3. Mann–Whitney U test results for comparing measures of RQA at 750 s with shorter data lengths. No significant differences show in any length considered.
Table 3. Mann–Whitney U test results for comparing measures of RQA at 750 s with shorter data lengths. No significant differences show in any length considered.
Length (s)60100150200300400500600
%REC0.9100.9101.0000.9461.0000.8740.9460.982
%DET1.0000.9820.9461.0001.0000.9461.0000.946
MDL0.9820.7690.8390.9820.8041.0001.0000.982
ADL0.7010.9100.9820.9820.9460.9460.9820.946
Table 4. The recommended minimum data length of each HRV measure.
Table 4. The recommended minimum data length of each HRV measure.
HRV ParametersRecommended Minimum Data Length (R Peaks)
Time-domain HRV
   Geometric measure
    Triangular index1000
   Statistical measure
    SDNN100
    RMSSD60
    pNN5060
Frequency-domain HRV
   Welch’s periodogram
    VLF1000
    LF1000
    HF60
    Total power1000
    VLF norm1000
    LF norm60
    HF norm750
    LF/HF60
   Lomb–Scargle’s periodogram
    VLF60
    LF500
    HF1000
    Total power500
    VLF norm1000
    LF norm60
    HF norm500
    LF/HF60
Nonlinear HRV
   Poincaré plot
    SD160
    SD21000
    SD1/SD21000
   Entropy
    SampEn1000
    ApEn1000
    MSE1000
    CMSE1000
   Fractal
    DFA1500
   Lyapunov exponent
    Wolf1500
    Rosenstein-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chou, E.-F.; Khine, M.; Lockhart, T.; Soangra, R. Effects of ECG Data Length on Heart Rate Variability among Young Healthy Adults. Sensors 2021, 21, 6286. https://doi.org/10.3390/s21186286

AMA Style

Chou E-F, Khine M, Lockhart T, Soangra R. Effects of ECG Data Length on Heart Rate Variability among Young Healthy Adults. Sensors. 2021; 21(18):6286. https://doi.org/10.3390/s21186286

Chicago/Turabian Style

Chou, En-Fan, Michelle Khine, Thurmon Lockhart, and Rahul Soangra. 2021. "Effects of ECG Data Length on Heart Rate Variability among Young Healthy Adults" Sensors 21, no. 18: 6286. https://doi.org/10.3390/s21186286

APA Style

Chou, E. -F., Khine, M., Lockhart, T., & Soangra, R. (2021). Effects of ECG Data Length on Heart Rate Variability among Young Healthy Adults. Sensors, 21(18), 6286. https://doi.org/10.3390/s21186286

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