In all three equations, QT interval is in milliseconds and the
RR interval is in seconds.
2.5 QTV indices
2.5.1 Time domain indices
Time domain indices include mean and standard deviation
(SD)ofQT,QTc,aswellasQTatimeseries,meanofT
amplitudes and QTVI, a metric that incorporates HRV in
QTV using the following formula [10]:
QTVI ¼ log
10
QT
v
=QT
2
m
HR
v
=HR
2
m
: ð4Þ
Here, QT
v
and QT
m
are the variance and mean of QT inter-
vals, respectively. HR
v
and HR
m
are the variance and mean of
heart rate series, respectively. QTVI was calculated only for
the raw QT time series.
2.5.2 Frequency-domain indices
Before the power spectral analysis, QT and QTa time series
were evenly resampled to 4 Hz by spline interpolation.
Power spectral density (PSD) was estimated by
autoregressive (AR) model based parametric method. In
this study, the Burg’s method was used to estimate the
model coefficients with an order of 16 [17].Basedonthe
frequency-domain characteristics of HRV, two frequency
bands were defined: (1) a low-frequency band from 0.04
to 0.15 Hz and (2) a high-frequency band from 0.15 to
0.40 Hz. The frequency-domain i ndices for QT and QTa
time series include the power of low frequency bands (LF)
and high frequency bands (HF).
2.5.3 Sample entropy
For the time series {u
i
= u(i), 1 ≤ i ≤ N}, form N − m + 1 vectors
X
m
i
¼ uiðÞ; uiþ 1ðÞ; …; uiþ m−1ðÞ; 1≤i≤N−m
.Define
the distance between all possible pair of vectors X
m
i
and X
m
j
by d
i, jm
=max{|u(i + k) − u(j + k)|, 0 ≤ k ≤ m − 1}. Denote the
average number of j that meets d
i, jm
≤ r for all 1 ≤ j ≤ N − m,
j ≠ i by B
m
i
rðÞ. Then increase the dimension m to m +1 and
repeat the above-described procedures to calculate B
mþ1
i
rðÞ.
The SampEn for u
i
can be defined by [16]:
SampEn m; r; NðÞ¼−ln ∑
N−m
t¼1
B
mþ1
i
rðÞ= ∑
N−m
t¼1
B
m
i
rðÞ
: ð5Þ
The threshold r value has been recommended to choose
between 0.1 × SD and 0.25 × SD [27]. In t his study, the
threshold was set at r =0.2×SDin order to minimize the pos-
sibility of getting invalid SampEn. The dimension was set at
m = 2 in accordance with the previous study [40].
2.5.4 Permutation entropy and dynamical patterns
For a time series {u
i
= u(i), 1 ≤ i ≤ N}, form N − m + 1 vectors
X
m
i
¼ uiðÞ; uiþ 1ðÞ; …; uiþ m−1ðÞ; 1≤i≤N−m þ 1
.A
symbolic counterpart is then obtained by coarse-graining the
time series using a given resolution Δ,i.e.,:
φ iðÞ¼floor
uiðÞ−min X
m
i
Δ
: ð6Þ
The N − m + 1 vectors can thus be represented by
φ
m
i
¼ φ iðÞ; φ i þ 1ðÞ; …; φ i þ m−1ðÞðÞ. The ordinal pattern
corresponding to φ
m
i
can be defined by the index vector when
rearranging φ
m
i
in ascending order. For instance, if
φ
m
i
¼ 10; 30; 20½, the index vector when rearranging it in
ascending order should be (132) which is defined as the ordi-
nal pattern for φ
m
i
.
Five dynamical patterns were defined by categorizing the
ordinal patterns into five different groups. The five dynamical
patterns for m = 3 were illustrated in Table 2.
The probability (frequency) of each dynamical pattern was
counted and denoted by {p(D
i
), 1 ≤ i ≤ 5}. PE is then calculat-
ed by:
PE m; ΔðÞ¼−∑
i
pD
i
lnpD
i
: ð7Þ
Table 1 Demographical and
clinical characteristics of all
subjects
Variables Healthy group CAD group CHF group
N 29 (16/13) 29 (22/7)** 20 (9/11)
Age (year) 56 ± 8 58 ± 8 59 ± 9
Body mass index (BMI) (kg/m
2
)24±826±325±4
Systolic blood pressure (SBP) (mmHg) 114 ± 13 117 ± 14 118 ± 14
Diastolic blood pressure (DBP) (mmHg) 72 ± 9 74 ± 9 75 ± 9
Left ventricular ejection fraction (LVEF) (%) 65 ± 4 60 ± 6 39 ± 7*
Data are expressed as number (male/female) or mean ± SD.
*p < 0.05 vs. healthy group as revealed by t test
**p < 0.01 vs. healthy group as revealed by chi-squared test
Med Biol Eng Comput (2019) 57:389–400 391