Residual heart rate variability measures can better differentiate patients with acute myocardial infarction from patients with patent coronary artery

Abstract

Purpose

It has been shown that the power spectral density (PSD) of heart rate variability (HRV) can be decomposed into a power-law function and a residual PSD (rPSD) with a more prominent high-frequency component than that in traditional PSD. This study investigated whether the residual HRV (rHRV) measures can better discriminate patients with acute myocardial infarction (AMI) from patients with patent coronary artery (PCA) than traditional HRV measures.

Materials and methods

The rHRV and HRV measures of 48 patients with AMI and 69 patients with PCA were compared.

Results

The high-frequency power of rHRV spectrum was significantly enhanced while the low-frequency and very low-frequency powers of rHRV spectrum were significantly suppressed, as compared to their corresponding traditional HRV spectrum in both groups of patients. The normalized residual high-frequency power (nrHFP = residual high-frequency power/residual total power) was significantly greater than the corresponding normalized high-frequency power in both groups of patients. Between-groups comparison showed that the nrHFP in AMI patients was significantly smaller than that in PCA patients. Receiver operating characteristic curve analysis showed that the nrHFP or nrHFP + normalized residual very low-frequency power (residual very low-frequency power/rTP) had better discrimination capability than the corresponding HRV measures for predicting AMI.

Conclusions

Compared with traditional HRV measures, the rHRV measures can slightly better differentiate AMI patients from PCA patients, especially the nrHFP or nrHFP + normalized residual very low-frequency power.

Keywords:

• heart rate variability
• fractal
• residual power spectrum
• power-law function
• acute myocardial infarction

Introduction

Heart rate variability (HRV) refers to the continuous fluctuations of RR intervals (RRI) around its mean because of autonomic nervous modulation. Power spectral analysis of heart rate (HR) fluctuations provides a noninvasive and quantitative means to assess the sympathetic and vagal modulations of HR.^1,2 HRV analysis has been used in the assessment of autonomic dysfunction caused by many kinds of major clinical diseases including acute myocardial infarction (AMI), acute coronary syndrome, postmyocardial infarction status, orthotopic heart transplantation, etc.

Wolf et al³ were the first to describe the association between HRV reduction and increased postinfarction mortality in 1978. By analyzing the 1-minute echocardiogram (ECG) recording obtained in patients with AMI immediately upon admission to coronary unit, they found that patients with more pronounced sinus arrhythmia had a lower mortality rate than patients with less pronounced variability of sinus impulses.³ Later on, it was found that AMI is almost inevitably associated with considerable reduction in HRV,⁴ and that the HRV was significantly lower in patients with myocardial infarction even a year after the acute coronary event as compared to healthy age-matched subjects.⁵ It was also found that HRV measures were improved in patients with better ejection fraction and greater angiographic patency after thrombolysis, suggesting that early HRV assessment after myocardial infarction may be useful in noninvasive risk stratification.⁶ Though decreased HRV can independently predict poor prognosis after myocardial infarction, the cut-off points that should be used in clinical practice are still a matter for further investigation.⁷ Measurement of HR turbulence upon arrival at the emergency department may provide additional incremental value in the risk assessment for patients with non-ST elevation myocardial infarction or unstable angina.⁸ Abnormal HRV can predict both sudden and nonsudden cardiac death after AMI.⁹

It has been shown that the power spectral density (PSD) of HRV contains a power-law relation that can be obtained by plotting the logarithm of PSD against the logarithm of frequency.^10–21 It was further shown that the PSD of HRV can be decomposed mathematically into a power-law function representing the power-law relation part of HRV, and a residual HRV spectrum representing the residual HRV (rHRV).²² The residual part can give rise to a greatly enhanced high-frequency (HF) component as well as a greatly suppressed low-frequency (LF) and very low-frequency (VLF) components of HRV.²²

This study investigated whether the rHRV measures can better differentiate the patients with AMI from those with patent coronary artery (PCA) than traditional HRV measures.

Materials and methods

Study subjects

This study reanalyzed the RRI data obtained in our previous work²³ by using the new method,²² and compared the newly derived rHRV measures with the traditional HRV measures in the study and control groups, and compared the rHRV measures between both groups of patients. Forty-eight patients with PCA and 69 patients with AMI were included in this study. The general data and clinical characteristics are listed in our previous publication.²³

Patients admitted to the intensive care unit with documented AMI, and consecutive patients with PCA were studied. PCA was defined as without stenosis or with luminal narrowing <30%. Patients with diabetes mellitus, atrial fibrillation, coexisting valvular heart diseases, cardiac conduction abnormalities, and liver and/or kidney disease were excluded. Patients who had atrial fibrillation or those using class I antiarrhythmic medication were excluded from this study. All experimental protocols were approved by the Institutional Review Board Committee B of Changhua Christian Hospital, Taiwan (CCH IRB No 170809), which agreed for a waiver of documentation of informed consent of this study. All research was performed in accordance with relevant guidelines/regulations.

This study reanalyzed the RRI data obtained in our previous work²³ published in 2003 by using the new method published in 2016,²² and compared the newly derived rHRV measures with traditional HRV measures in patients with PCA and AMI, and compared the rHRV measures between PCA and AMI patients. Since the ECG signals and relevant clinical data were obtained more than 15 years ago and have been delinked with the patients, it is very difficult to obtain informed consent or consent to review their medical records from the patients themselves or their next of kin. So, the Institutional Review Board Committee B of Changhua Christian Hospital, Taiwan, gave us a waiver of documentation of informed consent for this study. The patient data confidentiality was ensured because the ECG signals and the clinical data were already delinked with the patients.

HRV analysis

The power spectra of the 512 RRI were obtained by means of fast Fourier transformation. The area under the curve (AUC) of the spectral peaks within the frequency range of 0.01–0.4, 0.01–0.04, 0.04–0.15, and 0.15–0.40 Hz were defined as the total power (TP), very low-frequency power (VLFP), low-frequency power (LFP), and high-frequency power (HFP), respectively. The normalized VLFP (nVLFP = VLFP/TP) was used as the index of vagal withdrawal, rennin–angiotensin modulation, and thermoregulation,^24–26 the normalized LFP (nLFP = LFP/TP) as the index of combined sympathetic and vagal modulation,²⁷ the normalized HFP (nHFP = HFP/TP) as the index of vagal modulation, and the low-/high-frequency power ratio (LHR = LFP/HFP) as the index of sympathovagal balance.²⁸

Power-law function and residual HRV

The power spectrum of HRV was decomposed into a power-law function and a residual part of HRV,²² $\text{PSD}={\mathrm{F}}_{\text{rg}}\cdot \text{rPSD}={10}^{\mathrm{Y}}\cdot {\text{Frq}}^{\mathrm{s}}\cdot \text{rPSD}$(1) where the PSD is the traditional power spectral density, F_rg is the function of linear regression between log(PSD) and log(Frq), the subscript “rg” stands for “regression,” the “r” denotes “residual”, rPSD is the residual PSD, and the “s” and “Y” are the slope and Y-intercept of linear regression between log(PSD) and log(Frq) within the frequency range from 0+ Hz to the Nyquist frequency, respectively.

Similar to the definition of traditional HRV measures, the AUC of the spectral peaks within the range of 0.01–0.4, 0.01–0.04, 0.04–0.15, and 0.15–0.40 Hz in the rPSD were defined as the residual total power (rTP), residual very low-frequency power (rVLFP), residual low-frequency power (rLFP), and residual high-frequency power (rHFP), respectively. The normalized rVLFP (nrVLFP = rVLFP/rTP), nor-malized rLFP (nrLFP = rLFP/rTP), normalized rHFP (nrHFP = rHFP/rTP), and residual low-/high-frequency power ratio (rLHR = rLFP/rHFP) were defined in similar ways to those of traditional HRV measures.

Statistical analysis

Kruskal–Wallis rank sum test was used to compare the HRV and rHRV measures between PCA and AMI patients. Wilcoxon signed rank test was employed to compare the traditional HRV measures with the corresponding rHRV measures in both groups of patients. All data are presented as median (25%–75%). Receiver operating characteristic (ROC) curves with AUC analysis was conducted to compare the predictive value of HRV and rHRV measures for predicting AMI. The optimal cutpoint for the classification of high-risk patient for each HRV or rHRV measure is determined based on the Youden’s index in ROC analysis to produce an analysis with sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy. A P<0.05 was considered statistically significant.

Results

Figure 1 shows the tachogram, traditional HRV spectrum, log–log plots of traditional HRV spectra, power-law function of traditional HRV spectra, and the residual HRV spectra of a representative patient with PCA and a representative patient with AMI. The mean RRI and amplitude of RRI oscillation in the RRI tachogram and the powers in the HRV spectrum of the AMI patient are smaller than those of the PCA patient. The negative linear correlation between log(PSD) and log(Frq) in the HRV spectrum indicates that a power-law relation exists between PSD and Frq in the HRV spectra in both PCA and AMI patients. After removal of the power-law function, the rHRV spectra have greatly enhanced HF components and greatly reduced VLF and LF components in both PCA and AMI patients, as compared with the traditional HRV spectra. The VLF and HF components of the HRV spectrum in the AMI patient are smaller than those of the PCA patient.

Figure 1 The tachogram (A and F), traditional HRV spectrum (B and G), log–log plots of traditional HRV spectra (C and H), power-law function inside the traditional HRV spectra (D and I), and the residual HRV spectra (E and J) of a representative patient with PCA and a representative patient with AMI.

Notes: Both mean RRI and amplitude of RRI oscillation in the AMI patient (F) are smaller than those of the PCA patient (A). The powers in the HRV spectrum of the AMI patients (G) are also smaller than those of the PCA patient (I). The negative linear correlation between log (PSD) and log(Frq) in the HRV spectrum indicates that there is a power-law relation between PSD and Frq in the HRV spectra in both PCA and AMI patients. After removal of the power-law relation between PSD and Frq (D and I), the resultant rPSD (E and J) has relatively prominent HF component compared with the PSD in the traditional HRV spectra (B and G) in both PCA and AMI patients. The VLF and HF components in the rHRV spectrum of the AMI patient (J) are smaller than those of the PCA patient (E).
Abbreviations: AMI, acute myocardial infarction; Frq, frequency; HF, high frequency; HRV, heart rate variability; PCA, patent coronary artery; PSD, power spectral density; RRI, RR interval; VLF, very low frequency.

Table 1 compares the HRV and rHRV measures between the PCA and AMI groups of patients. The AMI patients had significantly greater HR, nVLFP, and nrLFP/nLFP than the PCA patients, and significantly smaller mean, root mean squared successive difference, HFP, nHFP, and nrHFP than the PCA patients. Table 1 also compares the rHRV measures with their corresponding HRV measures in the same group of PCA or AMI patients. The rTP, rVLFP, rLFP, rHFP, nrVLFP, nrLFP, and rLHR were all significantly smaller than their corresponding traditional TP, VLFP, LFP, HFP, nVLFP, nLFP, and LHR, while only the nrHFP was significantly greater than its corresponding nHFP in both groups of patients. The finding that the nrHFP/nHFP was >1 while the nrVLFP/nVLFP, nrLFP/nLFP, and rLHR/LHR were all smaller than 1 indicated that the removal of the power-law function from the PSD in the HRV spectrum could result in an enhanced HF part and suppressed low- and very low-frequency parts in the rHRV spectrum. The significantly greater nrLFP/nLFP of AMI patients than that of PCA patients indicates that the reduction in the LF components in AMI patients was smaller than that in PCA patients after the removal of the fractal part of HRV spectra.

Table 1 Comparison of HRV and rHRV measures between patients with PCA and patients with AMI

Parameters	PCA (N=48)	AMI (N=69)	P-value
General data
Age (years)	68 (60–72)	68 (53–74)	0.346
Gender (M/F)	38/10	59/10	0.370
HR (bpm)	65 (54–69)	71 (64–80)	<0.001
Traditional HRV measures
Mn (ms)	921 (864–1113)	845 (752–933)	<0.001
SDRR (ms)	26 (19–45)	27 (16–37)	0.368
CVRR	0.03 (0.02–0.04)	0.03 (0.02–0.04)	0.846
RMSSD (ms)	21 (16–36)	19 (10–29)	0.037
TP (ms^Citation2)	199 (144–807)	254 (91–516)	0.316
VLFP (ms^Citation2)	71 (42–220)	106 (45–24)	0.702
LFP (ms^Citation2)	56 (28–151)	52 (18–101)	0.292
HFP (ms^Citation2)	103 (51–240)	74 (22–171)	0.036
nVLFP (nu)	33.8 (24.7–45.0)	45.9 (29.8–58.8)	0.007
nLFP (nu)	19.7 (12.7–30.7)	19.7 (12.6–26.9)	0.587
nHFP (nu)	39.7 (24.8–61.1)	29.3 (18.9–46.5)	0.045
LHR	0.59 (0.21–1.18)	0.74 (0.36–1.03)	0.308
Residual HRV measures
rTP (ms^Citation2)	1.17 (1.02–1.45)a	1.16 (0.95–1.35)a	0.192
rVLFP (ms^Citation2)	0.08 (0.05–0.13)a	0.09 (0.07–0.15)a	0.161
rLFP (ms^Citation2)	0.16 (0.12–0.20)a	0.16 (0.13–0.20)a	0.929
rHFP (ms^Citation2)	0.91 (0.77–1.10)a	0.85 (0.70–0.99)a	0.115
nrVLFP (nu)	6.43 (4.84–10.02)a	8.45 (5.51–11.03)a	0.053
nrLFP (nu)	13.5 (9.2–18.0)a	13.7 (10.3–19.3)a	0.463
nrHFP (nu)	77.6 (72.6–82.9)a	75.3 (71.1–77.9)a	0.039
rLHR	0.17 (0.12–0.23)a	0.18 (0.14–0.26)a	0.321
Traditional/residual ratio
nrVLFP/nVLFP	0.19 (0.14–0.37)	0.19 (0.14–0.29)	0.820
nrLFP/nLFP	0.65 (0.54–0.79)	0.73 (0.64–0.92)	0.042
nrHFP/nHFP	1.98 (1.43–3.15)	2.45 (1.64–3.77)	0.085
rLHR/LHR	0.28 (0.19–0.56)	0.25 (0.19–0.44)	0.475

Notes:

a Significant difference between rHRV and the corresponding HRV measure in the same group of patients. Numeric values are expressed as the median (interquartile range).

Abbreviations: AMI, acute myocardial infarction; CV_RR, coefficient of variation of RR intervals; HFP, high-frequency power; HR, heart rate; HRV, heart rate variability; LFP, low-frequency power; LHR, low-/high-frequency power ratio; Mn, mean RR interval; PCA, patent coronary artery; rHRV, residual HRV; RMSSD, root mean squared successive differences; SD_RR, standard deviation of RR intervals; TP, total power; VLFP, very low-frequency power; bpm, beats per minute; ms, millisecond; nHFP, normalized HFP; nLFP, normalized LFP; nVLFP, normalized VLFP; nrHFP, normalized rHFP; nrLFP, normalized rLFP; nrVLFP, normalized rVLFP; nu, normalized unit; rHFP, residual HFP; rLFP, residual LFP; rLHR, residual LHR; rTP, residual TP; rVLFP, residual VLFP.

Table 2 shows the HRV and rHRV measures the AUC of which is >0.5 in predicting AMI. The Youden’s index suggests that the high-risk patients with AMI are those whose TP <95.40, rTP <0.96, nVLFP >42.14, nrVLFP >7.18, nHFP <40.43, and nrHFP <78. Figure 2 shows the ROC curves comparing the discrimination capability of various cutpoints of HRV and rHRV measures in predicting AMI. The patients with nrHFP <78 have an AUC of 0.627 with an accuracy rate of 64.96% in predicting AMI, while the patients with nHFP <40.43 have an AUC of 0.612 with an accuracy rate of 63.25% in predicting AMI. The nrHFP was better than nHFP in the differentiation between AMI and PCA. The patients with nrVLFP >7.18 and nrHFP <78 had an AUC of 0.665 with an accuracy rate of 64.96% in predicting AMI, while the patients with nVLFP <40.43 and nHFP >42.14 had an AUC of 0.649 with an accuracy rate of 61.54% in predicting AMI. The combination of nrHFP and nrVLFP had the highest accuracy rate compared with the other rHRV measures in the differentiation between AMI and PCA patients, and was found to be better than the corresponding combination of nHFP and nVLFP. Thus, the rHRV measures have better differential capability than the HRV measures in the differentiation between AMI and PCA patients, especially the nrHFP or nrHFP + nrVLFP.

Figure 2 The ROC curves with optimal cutpoints for those HRV and rHRV measures the AUC of which is >0.5 in predicting AMI.

Notes: (A) The high-risk patients with AMI are those whose TP <95.40 and rTP <0.96, with an AUC of 0.617 and 0.578, respectively. (B) The high-risk patients with AMI are those whose nHFP <40.43 and nrHFP <78, with an AUC of 0.612 and 0.627, respectively. (C) The high-risk patients with AMI are those whose nVLFP >42.14 and nrVLFP >7.18, with an AUC of 0.654 and 0.618, respectively. (D) The high-risk patients with AMI are those whose nHFP <40.43 + nVLFP >42.14 and nrHFP <78 + nrVLFP >7.18, with an AUC of 0.649 and 0.665, respectively. The nrHFP <78 + nrVLFP >7.18 is the most accurate diagnostic criteria for AMI, with an AUC of 0.665 (95% CI, 0.57–0.77; P=0.002) and accurate rate of 64.96%.
Abbreviations: AMI, acute myocardial infarction; AUC, area under the curve; HRV, heart rate variability; nHFP, normalized high-frequency power; nVLFP, normalized very low-frequency power; nrHFP, normalized residual high-frequency power; nrVFLP, normalized very low-frequency power; ROC, receiver operating characteristic; rTP, residual total power; TP, total power.

Table 2 The HRV and rHRV measures that have Youden’s cutpoints with an AUC >0.5 in predicting AMI

Predictor	Cutpoint	AUC	Sensitivity	Specificity	PPV	NPV	Accurate rate
TP	TP<95.40	0.617a	27.54%	95.83%	90.48%	47.92%	55.56%
rTP	TP<0.96	0.578	26.09%	89.58%	78.26%	45.74%	52.14%
nVLFP	nVLFP<42.14	0.654**	57.97%	72.92%	75.47%	54.69%	64.10%
nrVLFP	nrVLFP<7.18	0.618a	65.22%	58.33%	69.23%	53.85%	62.39%
nHFP	nHFP<40.43	0.612a	72.46%	50.00%	67.57%	55.81%	63.25%
nrHFP	nrHFP<78	0.627a	75.36%	50.00%	68.42%	58.54%	64.96%
nHFP + nVLFP	nHFP<40.43+ nVLFP<42.14	0.649**	46.38%	83.33%	80.00%	51.95%	61.54%
nrHFP + nrVLFP	nrHFP<78+ nrVLFP<7.18	0.665**	57.97%	75.00%	76.92%	55.38%	64.96%

Notes:

a P<0.05;

** P<0.01.

Abbreviations: AUC, area under the curve; NPV, negative predictive value; PPV, positive predictive value; nHFP, normalized high-frequency power; nVLFP, normalized very low-frequency power; nrHFP, normalized residual high-frequency power; nrVLFP, normalized residual very low-frequency power; rTP, residual TP; TP, total power.

Discussion

Mathematically, the power spectrum can be decomposed into a power-law function and a residual part. The PSD underneath the power-law function in the LF range is higher than that in the HF range. The clinically concerned vagal and sympathetic activity as shown in the HF, LF, and VLF ranges are often obscured in the traditional HRV measures. After removal of the power-law function from the PSD, the autonomic nervous activities of the subject could be better reflected by the rHRV measures than by the traditional HRV measures.

The guidelines on HRV state that the regression line should be estimated within the frequency range below 0.04 Hz to exclude interference from the oscillations in the LF and HF bands.² However, the power-law relation between the PSD and the frequency does not exist below 0.04 Hz only, it exists within the whole frequency range below the Nyquist frequency. Therefore, in this study, the whole frequency range was included for regression analysis so that the powers within the whole frequency range could be analyzed. Because all powers are included in the linear regression analysis, the slope and intercept of the linear regression line depend on every individual power in the HRV spectrum. If the harmonic oscillations of various physiological systems above 0.04 Hz are removed before the application of linear regression analysis, then the decomposition of the PSD into a power-law function and a residual PSD will be inapplicable because the power-law function in that case will be valid only below the frequency of 0.04 Hz while the HRV spectrum is valid within the whole frequency range below the Nyquist frequency. Thus, the linear regression analysis between log(PSD) and log(Frq) was performed within the whole frequency range below the Nyquist frequency in this study.

In this study the HFP was normalized by the true TP rather than by TP–VLFP as suggested in the guidelines on HRV.² If the HRV or rHRV measures used in the analysis contained only HF and LF parts, it may be acceptable to normalize the HFP and LFP by TP–VLFP. If, however, the VLF parts and even ultralow-frequency parts are used in the analysis, normalization by the true TP is needed. Furthermore, if the powers in the whole range of frequency from 0+ to the Nyquist frequency were included in the regression analysis, it may not be logical to normalize the individual powers by using TP–VLF, as suggested by the guidelines on HRV.

It has been reported that the HRV measures are decreased in survivors of cardiac arrest.²⁹ Some studies found that the SD of normal-to-normal QRS intervals is a significant parameter for long-term prognosis in patients with AMI.^30,31 In the frequency domain, it has been found that the nHFP was the lowest in patients with AMI, followed in increasing order by patients with coronary artery disease and patients with PCA in three recumbent positions; whereas the nLFP and LHR were the largest in patients with AMI, followed in decreasing order by patients with coronary artery disease and patients with PCA in three recumbent positions.³² In this study, we found similarly that that the AMI patients had significantly greater HR and nVLFP than the PCA patients, and significantly smaller mean, root mean squared successive difference, HFP, nHFP, and nrHFP than the PCA patients. Our findings are compatible with previous reports that the patients with AMI have overactive sympathetic modulation and suppressed vagal modulation.^33–36

HRV is a powerful predictor of arrhythmia-associated complications in patients surviving AMI.^37,38 The predictive value of depressed HRV after AMI is independent of other risk factors, such as left ventricular ejection fraction, frequency of ventricular premature complexes on the 24-hour Holter recording,³⁷ and the presence of late potentials on signal-averaged ECGs.³⁹ In previous studies, the HR assessed from standard short-term ECGs has been shown to predict mortality in postmyocardial infarction patients.^40–42 In accordance with those studies, we found that the AMI patients had significantly greater HR and nVLFP, and significantly smaller mean, root mean squared successive difference, HFP, nHFP, and nrHFP than the PCA patients. In the ROC analysis, we found that the nrHFP or nrHFP + nrVLFP had better differential capability than the corresponding HRV measures for the differentiation between AMI and PCA patients. Our findings suggested that the rHRV measures can be used in the clinical diagnosis and monitoring of AMI in the future.

Though PCA patients were found to have patent coronary arteries by using coronary angiography, they were not truly healthy people because they had cardiac symptoms and were willing to receive invasive coronary angiography examination. This may explain why the rHRV measures can give only slightly better discrimination capability between PCA and AMI patients than traditional HRV measures. If truly normal subjects were used as the control subjects, the rHRV measures might have better discrimination capability than traditional HRV measures in the differentiation between AMI patients and truly normal subjects.

The main purpose of this study was to investigate whether the isolation of the rHRV spectrum from traditional HRV spectrum could give rise to better parameters that can identify high-risk patients with AMI. Though the discrimination capability of the rHRV measures is not tremendously greater than that of traditional HRV measures, some improvements were observed. This small improvement in the discrimination between PCA and AMI suggests that with more refinements in the technique of decomposition, we may be able to obtain better rHRV measures to identify high-risk AMI patients.

The use of 24-hour ECG recordings in spectral HRV analysis has its advantage over using 10-minute ECG recordings; however, it also has some disadvantages. First, the physician caring for the patient with cardiac symptoms, the patient, and his/her family cannot wait for 24 hours to obtain the HRV data for diagnosis. Second, interferences from diurnal changes, therapeutic activities of medical personnel, and body motions of the patient will be recorded in the ECG tracing such that a lot of artifacts and interferences will be present in the HRV spectrum. The HRV and rHRV measures obtained from 24-hour ECG recordings may not reflect the true cardiac autonomic activities of the patient. Therefore, the reliability of HRV measures obtained from 10-minute ECG recording may not be less than that obtained from 24-hour ECG recording.

Though most of the bibliography cited, as well as the HRV guidelines, recommended the use of frequency range from 10⁻² to 10⁻⁴ Hz, some studies used the frequency ranges from 10⁻³ to 1 Hz,¹³ 10⁻² to 0.7 Hz,¹⁹ and 10⁻⁴ to 1 Hz.²⁰ The purpose of decomposing the power spectrum of HRV into a power-law function of frequency and a residual power spectrum is to examine more closely the HR oscillations in the VLF, LF, and HF regions, because the VLF, LF, and HF regions are associated with autonomic modulation of the patients and are the parts most concerned in many clinical settings. If we adhere to the guidelines on HRV to analyze the PSD within the frequency range from 10⁻⁴ to 10⁻² Hz, then the HR oscillations in the VLF, LF, and HF regions cannot be analyzed by using the decomposition technique. Furthermore, in Figure 1C and H, we can see that a negative linear correlation between log(PSD) and log(Frq) exists within the whole frequency range. There is no reason to do linear regression analysis between log(PSD) and log(Frq) within the frequency range from 10⁻⁴ to 10⁻² Hz only. Thus, linear regression analysis between log(PSD) and log(Frq) within the whole frequency range from 0+ to the Nyquist frequency was performed in this study.

One limitation of this study was that the relationship between rHRV measures and ventricular arrhythmia risk in AMI patients was not assessed, including late potentials on signal-averaged ECG and ECG markers of sudden cardiac death (QT interval, T_peak–T_end interval, QRS duration, and fragmented QRS, ST segment changes).⁴³ Further studies are needed to delineate the relationship between rHRV measures and the abovementioned ECG risk factors in AMI patients.

Conclusion

The HF part of the rHRV spectrum is augmented while the VLF and LF parts of the rHRV spectrum are suppressed compared to traditional HRV spectrum. The combination of nrHFP <78 nu and nrVLFP >7.18 nu can slightly better differentiate AMI patients from PCA patients than the combination of corresponding traditional HRV measures.

Acknowledgments

This study was supported by a grant (V98C1-008) from the Taipei Veterans General Hospital, and a grant (NSC100-2314- B-075-026) from the Ministry of Science and Technology, Taiwan.

Disclosure

The authors report no conflicts of interest in this work.