Abstract
Background
Heart rate variability (HRV) is known to be impaired in patients with congestive heart failure (CHF). Timedomain analysis of ECG signals traditionally relies heavily on linear indices of an essentially nonlinear phenomenon. Poincaré plots are commonly used to study nonlinear behavior of physiologic signals. Lagged Poincaré plots incorporate autocovariance information and analysis of Poincaré plots for various lags can provide interesting insights into the autonomic control of the heart.
Methods
Using Poincaré plot analysis, we assessed whether the relation of the lag between heart beats and HRV is altered in CHF. We studied the influence of lag on estimates of Poincaré plot indices for various lengths of beat sequence in a public domain data set (PhysioNet) of 29 subjects with CHF and 54 subjects with normal sinus rhythm.
Results
A curvilinear association was observed between lag and Poincaré plot indices (SD1, SD2, SDLD and SD1/SD2 ratio) in normal subjects even for a small sequence of 50 beats (p value for quadratic term 3 × 10^{5}, 0.002, 3.5 × 10^{5 }and 0.0003, respectively). This curvilinearity was lost in patients with CHF even after exploring sequences up to 50,000 beats (p values for quadratic term > 0.5).
Conclusion
Since lagged Poincaré plots incorporate autocovariance information, these analyses provide insights into the autonomic control of heart rate that is influenced by the nonlinearity of the signal. The differences in lagresponse in CHF patients and normal subjects exist even in the face of the treatment received by the CHF patients.
Background
Poincaré plot is an intuitive and commonly used method to assess complex nonlinear behavior in the study of physiological signals [16]. In the assessment of heart rate variability (HRV), Poincaré plots are constructed by plotting duplets of successive RR intervals [1,2,46], with an implicit assumption that the next RR interval is significantly determined by the current one. This assumption lends itself to further generalization of Poincaré plots by plotting mlagged plots where m represents the distance (in number of beats) between the duplet beats, that is, the 'lag' of the second beat from the first [2]. It has been observed in the context of the short term variability that the current RR interval can influence up to approximately eight subsequent RR intervals [2]. Therefore, a series of lagged Poincaré plots can potentially provide more information about the behavior of Poincaré plot indices in health and disease than the conventional 1lagged plot does [2].
Heart rate variability analysis provides a noninvasive means to assess the autonomic status of the heart [68]. Under normal conditions, the feedback elements characterized by vagal and sympathetic activation of the heart combined with the cardiac automaticity determine the HRV [7]. In various clinical conditions in which the sympathovagal balance is disturbed, such as after an episode of myocardial infarction, in diabetic autonomic neuropathy and in congestive heart failure, HRV is usually reduced [710]. In the context of congestive heart failure (CHF), the decrease in HRV has also been observed to correlate with disease severity [1115]. However, two issues relating to the strategies employed for HRV analysis deserve closer scrutiny. First, the majority of methods of quantifying HRV (including conventional Poincaré plots) use successive RR interval duplets only, with the implicit assumption that the current beat is influenced by the immediately preceding beat. However, it has been reported that a heart beat influences not only the beat immediately following it, but also up to 6–10 beats downstream [2], possibly as a consequence of respiratory sinus arrhythmia. Thus, an analysis hinging on the use of only successive RR interval duplets will likely underestimate the role of the autocovariance function of RR intervals i.e., the ability of heart beats to influence a train of succeeding beats. Second, the Poincaré plot indices relating to shortterm and longterm variability in RR intervals do not capture the nonlinear disposition of HRV [16]. The autocovariance function of RR intervals captures the additional aspects of HRV (e.g. nonlinearity) that can be masked by the strong correlation between successive beats if 1lagged plots are used. Indeed, Brennan et al [16] argue that lagged Poincaré plots can fully describe the autocovariance as well as the power spectrum of HRV. Our proposed analysis uses lagged Poincaré plots to overcome the limitations of the present practice of timedomain analysis of HRV.
Therefore, we hypothesized that the lagresponse patterns of linear and nonlinear Poincaré plot indices would be different in a diseased heart as compared to a normal heart. To test our hypothesis, we compared longterm ECG recordings of a group of CHF subjects with those of normal subjects. We also explored whether linear versus nonlinear indices of HRV behaved differentially with respect to the lag. Finally, we compared the lagresponses of CHF patients and normal subjects in light of the fact that heart rate variability may have been potentially restored by the pharmacologic therapy in the CHF patients. Thus, we suggest an alternative analytic strategy that has the potential to provide unique insights into the pathophysiological basis of CHF.
Methods
Study subjects
We used the PhysioNet internet resource which is a large repository of various physiologic signals including electrocardiograms recorded by a 24hour Holter monitor [17,18]. From this repository, we chose the data sets which included information on interbeat RR intervals. The chf2db database included records on 29 individuals with congestive heart failure. These data come from two previous trials of longterm digoxin [19] and carvedilol [20] therapy. These subjects included eight men and 2 women (gender was unknown in the remaining 21 subjects) aged 34 – 79 years. There were 4 subjects belonging to the NYHA class I, 8 belonged to NYHA class II and 17 to NYHA class III. For comparison, we used a set of ECG recordings from 54 healthy subjects with normal sinus rhythm. These subjects included 30 men (aged 28.5 – 76 years) and 24 women (aged 58 – 73 years). Complete details of the study subjects are provided in the Supplementary Table 1 (see 1).
Additional file 1. This document contains supplementary tables 1–4. Supplementary Table 1 details characteristics of the subjects included in the present analysis. Supplementary Table 2 shows replication of analyses shown in Table 1 (Main Text) in the subset of subjects who had pNN50 ≤ 10 (for both CHF and NSR datasets) and mean heart rate <100 beats/min (for the NSR dataset). Supplementary Table 3 shows replication of analyses shown in Table 2 (Main Text) in the subset of subjects who had pNN50 ≥ 10 (for both CHF and NSR datasets) and mean heart rate <100 beats/min (for the NSR dataset). Supplementary Table 4 shows the distribution of traditional measures of HRV across the study groups.
Format: DOC Size: 205KB Download file
This file can be viewed with: Microsoft Word Viewer
Extraction of NN intervals
The databases in the PhysioNet repository contain 24hour Holter ECG recordings in two data formats; unaudited beatwise annotations sampled at a frequency of 128 samples per second and the corresponding header identification information. The WFDB software package provided by the PhysioNet resource was used to extract beatwise RR intervals. We used the "ann2rr" command to extract a sequence of up to 50,000 normal beats (approximately 12 hours of ECG recording) to estimate the normaltonormal (NN, by specifying the N option in the command) intervals in seconds.
Study variables
The standard measures of heart rate and its variability in the time domain are considered linear measures [16]. These include the average RR interval (Mean RR) and its standard deviation (SDRR) [21]. Another approach is to use the differences in successive RR intervals (successive deltas) and estimate the standard deviation (SDSD) of these deltas [21]. In a Poincaré plot, one of several quantitative approaches includes the fitting of an ellipse to the data points and estimating the short and long axes of the ellipse (Figure 1) [2,16]. The minor axis is proportional to the standard deviation of the successive differences (SD1) whereas the major axis is proportional to the standard deviation of the current beat (SD2) [2,16]. It has been recognized that the minor axis is a measure of the shortterm variability while the major axis is a measure of the longterm variability of the signal [2,16]. Thus, SD1 and SD2 together capture the total variability in a data set. Given the nonlinear disposition of the Poincaré plot, it has been a convention to treat SD1 and SD2 as the nonlinear indices of HRV. However, it can be shown mathematically that SD1 is a function of SDSD while SD2 is a function of SDRR and SDSD [16]. As stated earlier, both SDRR and SDSD are linear indices and, by extension, SD1 and SD2 also behave mathematically as measures of the linear component of variability even though they are based on Poincaré plots [16]. Nonlinear transformations of SD1 and SD2 have also been used to better capture the complete nature of HRV. For example, the ratio of SD1 and SD2 indicates the ratio of shortterm versus longterm variability and has been observed to be more informative than either of these indices alone [2].
Figure 1. A typical Poincaré plot. The abscissa represents the RR interval of the current normal beat and ordinate represents the RR interval of the succeeding normal beat. An ellipse is fitted to the data points and the Poincaré plot indices are calculated by estimating the short diameter (SD1), the long diameter (SD2) and the ratio of the short and long diameters (SD1/SD2 ratio) of the fitted ellipse.
For our analyses, we included only those indices that mathematically depend on the lag. As Mean RR and SDRR are theoretically independent of lag, we did not include these indices for studying the lagresponse. However, mean deltas and SDSD can be appropriately generalized to measure the lag response. For instance, we estimated deltas as the difference between the mlagged RR intervals. Thus, we were able to estimate mlagged deltas and their standard deviations (SDLD). In the special case of 1lagged beats, the mean deltas and SDLD are respectively equal to mean successive delta and SDSD. We, then, estimated the SD1 and SD2 for each of the mlagged beat sequence. Thus, we included the following four indices in our analysis: SDLD, SD1, SD2 and SD1/SD2 ratio.
Length of beat sequence used for analysis
Controversy exists with regard to the period of time for which ECG recordings should be monitored to best capture the HRV dynamics. While some researchers maintain that recordings less than 18 hours are insufficient [15], others have observed that shortterm RR interval recordings are as reliable and accurate as longterm recordings in analyzing HRV [22]. Therefore, we studied the lagresponse of HRV for different lengths of beat sequences. We used the following seven lengths of consecutive beats: 50, 100, 500, 1000, 5000, 10000, and 50000. These sequence lengths represent recordings ranging from ~1 minute to ~12 hours, depending on the overall heart rate. We, thus, used the following analytical strategy: for each subject included in the study, we used seven lengths of beat sequences. For each beat sequence, we used lag values from one to ten. For each value of lag we constructed a Poincaré plot and estimated SDLD, SD1, SD2, and SD1/SD2 ratio. For construction of Poincaré plots as well as for estimation of the Poincaré plot indices, we used normal beats only as annotated in the PhysioNet database resource.
Analysis of lagresponse
Considering the facts that mlagged Poincaré plots can describe the autocovariance function, the autocovariance function monotonically decreases with increasing lag for values of lags less than 10, and that the current beat influences only about six to eight successive beats in the context of a short range influence, we expected a pattern of lagresponse which is stronger at the lower values of lag and which attenuates with increasing lag. In normal subjects, therefore, we expected a curvilinear relationship that explains variations in Poincaré plot indices on the basis of lag. To test for curvilinearity we used a quadratic relationship model and plotted the estimates of SDLD, SD1, SD2, and SD1/SD2 ratio against lag. We then fitted a secondorder polynomial curve using the leastsquares method. We assessed the modelfit using R^{2 }values. We then compared the average coefficients of the quadratic term in the secondorder polynomial equations in subjects with and without CHF.
Statistical analyses
We compared the linear and nonlinear mlagged Poincaré plot indices in CHF and normal groups by using MannWhitney rank sum test. To study the potential contribution of varying beat sequence length to estimates of the Poincaré plot indices, we used the nonparametric method of Spearman's correlation coefficient (ρ). The curvilinearity of the lagresponse was assessed by using a second order polynomial regression and the modelfit was assessed using the R^{2 }value in order to quantify the amount of variation explained by the quadratic lagresponse of the Poincaré plot indices. We wrote dedicated routines in Visual Basic^{® }for estimation of mlagged Poincaré plot indices. Then, we used Stata 7.0^{® }for statistical analysis. For all statistical analysis, we assumed statistical significance as p < 0.05.
Results
Length of ECG recording and Poincaré plot indices
We first assessed the role of length of ECG recording in measurement of the Poincaré plot indices. Table 1 provides a summary of the estimates of Poincaré plot indices of the study population. We observed that the shortterm variability captured by SD1 consistently tended to be higher in subjects with CHF as compared to subjects with normal sinus rhythm across varying lengths of beat sequence. This difference reached statistical significance only when the length of the sequence was at least 5000 beats (p ≤ 0.0043). The longterm variability revealed by SD2, however, showed interesting relationship with the length of the beat sequence. For sequence lengths of 50 and 100, the longterm variability tended to be higher in subjects with CHF than that in normal subjects. For sequence length of 500 beats and higher, the longterm variability tended to be less in CHF patients as compared to that in normal subjects. This difference just reached statistical significance (p = 0.0449) at beat sequence length of 1000. However, for beat sequence lengths 5000 and higher, the difference was highly significant (p ≤ 0.0047). Thus, our results indicated a critical change in statistical significance at the sequence length of 5000 beats. Estimates of shortterm variability increased with sequence length both in CHF patients and in normal subjects and the increase was statistically comparable in these two groups (p = 0.6577). The longterm variability showed a stronger direct association with beat sequence length. Moreover, the Spearman correlation coefficient was significantly higher in normal subjects as compared to the CHF patients (p = 7 × 10^{5}). In spite of a nonsignificant difference between CHF patients and normal subjects for both SD1 and SD2 in sequences less than 5000 beats in length, the ratio of these two indices (SD1/SD2) in the same sequences showed a highly statistically significant difference with the ratio being consistently higher in CHF patients (p ≤ 2.5 × 10^{6}) (Table 1).
Table 1. Summary of Poincaré plot indices in study subjects.*
Lagged Poincaré plots in CHF and normal sinus rhythm
In general, we observed that the estimates of the traditional HRV measures in CHF patients were higher than those in the subjects with normal sinus rhythm. Figures 2A and 2B show representative Poincaré plots for two subjects – one with CHF (Figure 2A) and one normal subject (Figure 2B). It is apparent from the plots that the variability was more in the CHF patient as compared to the normal subject. To understand the reasons for the increased HRV in the CHF group, we first examined hourly Poincaré plots for all subjects using all the recorded beats. The results are highly interesting (supplementary material files CHF_Binder and NSR_Binder, see additional files 2 and 3). As per the classification proposed by Woo et al [4], most of the control subjects (51 out of 54) had the normal comet shaped Poincare plots during each hour in a 24hour period. However most (25 out of 29) of the CHF patients had complex Poincare plots for most of the 24hour recordings. As pointed out by Woo et al [5], CHF patients have either torpedo shaped or complex plots. The Poincaré plots of the CHF patients in this study best fall within the complex category.
Figure 2. Representative Poincaré plots from two subjects. Panels A, C, E, and G are from a patient with CHF and panels B, D, F, and H are from a normal subject. For each subject the Poincaré plots use varying lags of 1, 2, 4 and 8 beats. The spread of the data points can be observed to increase for increasing lag for each subject. All plots are based on normaltonormal RR intervals. It can also be noted that in this particular pair of subjects, the variability was more in the subject with CHF than that in the normal subject.
Additional file 2. This file has hourly Poincaré plots for the 29 CHF patients.
Format: PDF Size: 1.3MB Download file
This file can be viewed with: Adobe Acrobat Reader
Additional file 3. This file has hourly Poincaré plots for the 54 control subjects.
Format: PDF Size: 2.1MB Download file
This file can be viewed with: Adobe Acrobat Reader
To illustrate the putative influence of lag on estimates of Poincaré plot indices, we modified the conventional Poincaré plot by plotting (n+m)^{th }beat on the ordinate against the n^{th }beat on the abscissa. Figures 2C–2H depict this technique. It is evident from these Poincaré plots that the scatter of points increases with increasing lag. This indicates that increasing lag corresponds to increasingly unrelated beats. Our hypothesis attempted to characterize this behavior in normal subjects and CHF patients.
Curvilinearity of lagresponse
We analyzed the association between Poincaré plotrelated estimates of heart rate variability and lag for varying beat sequence lengths. Figure 3 shows these associations for two extremes of beat sequence lengths (50, Figures 3A–3D; 50000, Figures 3E–3H). We observed that SDLD, SD1, SD2, and SD1/SD2 ratio were laginvariant (R^{2 }range 0.0003–0.094) in CHF patients using the smallest sequence length of 50 beats. In striking contrast, smooth secondorder polynomial curves explained more than 99 percent of the variance in the estimates of these indices on the basis of lag in the normal subjects. We therefore compared the mean values of coefficients of the quadratic term in CHF patients versus that in normal subjects. It can be observed from Table 2 that the coefficient of the quadratic term was consistently insignificant in CHF patients for all the indices studied. In normal subjects, however, all the indices had significant curvilinearity as indicated by significant coefficients of the quadratic term (Table 2). These results indicated that the curvilinearity of the lagresponse is diminished in CHF. This seminal finding was replicated at the other extreme of sequence length of 50000 beats (Figures 3E–3H, Table 2). In fact, we observed the same finding regardless of sequence lengths.
Figure 3. Lag response of Poincaré plot indices in patients of CHF and normal subjects. Panels AD show the lag responses for sequences 50 beats long whereas panels EH show the lag responses for sequences 50000 beats long. All the panels on the left hand side relate to subjects with CHF whereas the panels on the right hand side relate to normal subjects. We first estimated the Poincaré plot indices for each subject in both groups. We, then, estimated the averages for each group and plotted the estimates against the lag. In the normal subjects, all the Poincaré plot indices showed a curvilinear association with lag. In patients of CHF the curvilinearity was lost.
Table 2. Coefficients and statistical significance of the quadratic term in equations regressing heart rate variability indices on lag.
We again examined if the increased HRV seen in the CHF patients (described in the previous section) confounded our interpretations of the lagged Poincaré plot analyses. We observed that (Supplementary Table 1 in 1) 10 of the 29 CHF patients and 8 of the 54 control subjects had pNN50 values exceeding 10%. In addition 7 of the control subjects had mean heart rates over 100. To exclude the possibility that these abnormal values for pNN50 and mean heart rates could be responsible for the higher HRV in CHF as against normal subjects, we excluded these 25 subjects and performed a subset analysis on 19 CHF patients and 39 control subjects. The results from these subset analyses are detailed in supplementary tables 2 and 3 (see 1). These results obtained are completely concordant with those obtained without excluding any subject from analysis (Tables 1 and 2). Thus, high pNN50 and mean heart rates in some subjects did not seem to influence our main inference about the loss of a curvilinear lag response in CHF.
Discussion
Our results proffer compelling evidence that normally there is a curvilinear relation between lag and Poincaré plot indices of heart rate variability and that this curvilinearity is impaired in CHF. The main implication of our work is that differences in heart rate variability of CHF patients and normal subjects can be detected when these are not apparent from traditional measures of HRV. In addition, a potential interpretation of these results may be that the short term influence of the lagged beats is importantly determined by the respiratory sinus arrhythmia. Although traditionally viewed as an autonomicmediated phenomenon, respiratory sinus arrhythmia may also involve nonautonomic mechanisms. The presence of complex Poincaré plots in the CHF group suggests that a combination of autonomic and nonautonomic control mechanisms may be responsible for the sinus arrhythmia in some of these patients. Further investigation is needed to better define the interpretation of these analyses in CHF patients.
Traditionally, HRV analysis has been used to assess the autonomic status of the heart and Poincaré plots have been widely used for HRV analysis. But the potential utility of lagged Poincaré plots has not been fully appreciated. Considering the fact that a heart beat is affected by its preceding beats, it follows logically that lag would affect the indices of HRV. We expect this effect to be stronger at the lower values of lag and to become weaker with increasing lag. Therefore, we expect a curvilinear relation between lag and indices of HRV normally. This curvilinearity is lost in CHF. This analysis, then, uncovers a potential means to further risk stratify patients with CHF.
Useful length of beat sequence
Our study also demonstrated that a sequence of 5000 beats and more is useful for analyzing both the shortterm and longterm HRV. Although our study was not designed to identify the useful minimum length of a beat sequence, we did observe the association of the estimates of Poincaré plot indices with varying lengths of beat sequences. We also found that the SD1/SD2 ratio was significantly different in CHF patients as compared to normal subjects for any length of the beat sequences. This observation is consistent with the recent notion that SDSD, SD1, and SD2 only capture the linear aspect of HRV whereas the SD1/SD2 ratio may better relate to the nonlinear component of HRV [2,16]. Therefore, in the timedomain analysis of ECG, it may be the most informative to use the SD1/SD2 ratio. Our method further capitalizes on this property of the SD1/SD2 ratio by using mlagged SD1/SD2 ratios. Such a modification, by virtue of an implicit property of the autocovariance function, can improve the use of SD1/SD2 ratio as a result of a closer approximation of frequencydomain analysis.
Differences in HRV in CHF and in normal sinus rhythm
The analysis conducted in this study makes use of the autocovariance function which allows for a stronger and more robust internal comparison for individual subjects as opposed to the more inaccurate method of group comparisons. We observed a loss of curvilinearity of lagresponse that prevailed in CHF patients with therapy in spite of an improvement in HRV attributable to the drug therapy. It is known that digoxin only improves the clinical profile of CHF patients but does little to improve survival [21,23]. Thus, even though the patients on therapy seem to have greater HRV than normal, they might still be at risk for early death. This can have great prognostic implications. We propose that patients in whom lagresponse curvilinearity is restored may have a better prognosis than those in whom it is not restored. Further prospective studies are needed to test this hypothesis. Nevertheless, these initial findings suggest that this method offers a unique advantage over the conventional method of 1lagged Poincaré plots.
Study limitations
Our study suffers from several limitations. First, we conducted a secondary data analysis. In addition to the inherent limitations of such analyses, we were restricted to the analysis of the ECG signals only. It would have been more informative to study other pathophysiological and sociodemographic correlates of CHF. Second, the small sample sizes of both the CHF and normal sinus rhythm datasets do not allow wider generalization of our results. Therefore, studies of larger magnitude that can provide more conclusive insights into the utility of these measures of HRV are needed. Third, rather surprisingly, we observed more heart rate variability in CHF patients as compared to the normal subjects, which contrasts the generally accepted notion that HRV is decreased in CHF as compared to normal sinus rhythm. This unexpected result can be partially explained by the facts that i) A large proportion of the CHF patients exhibited patterns of complex Poincaré plots; ii) It is possible that the differential distribution of age, gender, socioeconomic status and other unknown parameters in the comparison group may have confounded the observations on increased heart rate variability in CHF subjects; and iii) all the CHF patients in this study were receiving longterm digoxin therapy [19,20] which can potentially improve HRV in the CHF patients [23,24].
Furthermore, though SD1 and SD1/SD2 ratios were higher in CHF patients as compared to normal subjects, SD2 (an index of longterm variability) was higher in the normal subjects for beat sequences 500 and above. This is consistent with two sets of observations: first, CHF patients are more likely to have a high degree of sinus arrhythmia that is not of respiratory origin, that is, they are more likely to have sudden jumps in the NN interval; and second, a study by Huikuri et al reported that depression of long term variability is a predictor of mortality [25]. Thus, CHF patients in this study had lower SD2 and hence a lower longterm heart rate variability than normal subjects, which might be responsible for the poorer prognosis in the former. We also studied three other measures of HRV across the study groups: percentage of beats with change in successive NN intervals exceeding 50 ms (pNN50) and 20 ms (pNN20) as well as the standard deviation of the NN intervals (SDNN). The individual data is shown in Supplementary Table 1 (see 1). Using a MannWhitney test we observed that while in general the values for all the three parameters tended to be higher in the CHF group, these were not statistically significantly different (Supplementary Table 4, see 1). These observations further highlight the value of the proposed method to gain insights into HRV.
Conclusion
In summary, despite the presence of these limitations, our study provides strong evidence for the loss of curvilinearity in the lagresponse of Poincaré plot indices in CHF patients even after receiving longterm therapy and even in cases in which some traditional measures of HRV are similar to those seen in normal individuals. We propose that these approaches of analysis can be used to improve the timedomain analysis of ECG signals and can enhance the diagnostic and prognostic indicators of CHF.
Abbreviations
1. CHF – Congestive Heart Failure
2. ECG – Electrocardiogram
3. HRV – Heart Rate Variability
4. RR – RR interval (these are considered to be measured from normal beats, NN)
5. SD1 – Standard Deviation 1 (long axis of ellipse fitted to Poincaré plot)
6. SD2 – Standard Deviation 2 (short axis of ellipse fitted to Poincaré plot)
7. SDLD – Standard Deviation of Lagged Deltas
8. SDRR – Standard Deviation of RR intervals
9. SDSD – Standard Deviation of Successive Differences
10. pNN50 – percentage of beats with change in successive NN intervals exceeding 50 ms
11. pNN20 – percentage of beats with change in successive NN intervals exceeding 20 ms
12. SDNN – standard deviation of the NN intervals
Competing interests
The author(s) declare that they have no competing interests.
Authors' contributions
TPT conceptualized the study, performed the analysis of the data and wrote the initial and revised drafts of the manuscript. MLS contributed to the data analysis and the review and revision of the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors wish to thank Drs. Hemant Kulkarni and Manju Mamtani, Lata Medical Research Foundation, Nagpur, India for their constant help and discussion during conceptualization and conduct of the study. We also thank Dr. Phyllis Stein, Dr. Michal Javorka and Dr. Jacques Regnard for the insightful reviews.
References

Bergfeldt L, Haga Y: Power spectral and Poincare plot characteristics in sinus node dysfunction.
J Appl Physiol 2003, 94(6):22172224. PubMed Abstract  Publisher Full Text

Lerma C, Infante O, PerezGrovas H, Jose MV: Poincare plot indexes of heart rate variability capture dynamic adaptations after haemodialysis in chronic renal failure patients.
Clin Physiol Funct Imaging 2003, 23(2):7280. PubMed Abstract  Publisher Full Text

Tulppo MP, Makikallio TH, Seppanen T, Airaksinen JK, Huikuri HV: Heart rate dynamics during accentuated sympathovagal interaction.
Am J Physiol 1998, 274(3 Pt 2):H8106. PubMed Abstract  Publisher Full Text

Woo MA, Stevenson WG, Moser DK: Comparison of four methods of assessing heart rate variability in patients with heart failure.
Am J Crit Care 1996, 5(1):3441. PubMed Abstract

Woo MA, Stevenson WG, Moser DK, Middlekauff HR: Complex heart rate variability and serum norepinephrine levels in patients with advanced heart failure.
J Am Coll Cardiol 1994, 23(3):565569. PubMed Abstract

Woo MA, Stevenson WG, Moser DK, Trelease RB, Harper RM: Patterns of beattobeat heart rate variability in advanced heart failure.
Am Heart J 1992, 123(3):704710. PubMed Abstract  Publisher Full Text

Stys A, Stys T: Current clinical applications of heart rate variability.
Clin Cardiol 1998, 21(10):719724. PubMed Abstract

Casolo G, Balli E, Taddei T, Amuhasi J, Gori C: Decreased spontaneous heart rate variability in congestive heart failure.
Am J Cardiol 1989, 64(18):11621167. PubMed Abstract  Publisher Full Text

Casolo G, Balli E, Fazi A, Gori C, Freni A, Gensini G: Twentyfourhour spectral analysis of heart rate variability in congestive heart failure secondary to coronary artery disease.
Am J Cardiol 1991, 67(13):11541158. PubMed Abstract  Publisher Full Text

Casolo GC, Stroder P, Sulla A, Chelucci A, Freni A, Zerauschek M: Heart rate variability and functional severity of congestive heart failure secondary to coronary artery disease.
Eur Heart J 1995, 16(3):360367. PubMed Abstract

Guzzetti S, Cogliati C, Mezzetti S, Turiel M, Lombardi F, Malliani A: Heart rate and systolic arterial blood pressure variabilities in the progression of chronic heart failure.
Clin Sci (Lond) 1996, 91 Suppl:3739. PubMed Abstract

Guzzetti S, Magatelli R, Borroni E, Mezzetti S: Heart rate variability in chronic heart failure.
Auton Neurosci 2001, 90(12):102105. PubMed Abstract  Publisher Full Text

Guzzetti S, Mezzetti S, Magatelli R, Porta A, De Angelis G, Rovelli G, Malliani A: Linear and nonlinear 24 h heart rate variability in chronic heart failure.
Auton Neurosci 2000, 86(12):114119. PubMed Abstract  Publisher Full Text

Stefenelli T, BerglerKlein J, Globits S, Pacher R, Glogar D: Heart rate behaviour at different stages of congestive heart failure.
Eur Heart J 1992, 13(7):902907. PubMed Abstract

Wijbenga JA, Balk AH, Meij SH, Simoons ML, Malik M: Heart rate variability index in congestive heart failure: relation to clinical variables and prognosis.
Eur Heart J 1998, 19(11):17191724. PubMed Abstract  Publisher Full Text

Brennan M, Palaniswami M, Kamen P: Do existing measures of Poincare plot geometry reflect nonlinear features of heart rate variability?
IEEE Trans Biomed Eng 2001, 48(11):13421347. PubMed Abstract  Publisher Full Text

Goldberger AL, Amaral LA, Glass L, Hausdorff JM, Ivanov PC, Mark RG, Mietus JE, Moody GB, Peng CK, Stanley HE: PhysioBank, PhysioToolkit, and PhysioNet: components of a new research resource for complex physiologic signals.
Circulation 2000, 101(23):E21520. PubMed Abstract  Publisher Full Text

The PhysioNet Database [http://www.physionet.org] webcite

Krum H, Bigger JTJ, Goldsmith RL, Packer M: Effect of longterm digoxin therapy on autonomic function in patients with chronic heart failure.
J Am Coll Cardiol 1995, 25(2):289294. PubMed Abstract  Publisher Full Text

Goldsmith RL, Bigger JT, Bloomfield DM, Krum H, Steinman RC, SacknerBernstein J, Packer M: Longterm carvedilol therapy increases parasympathetic nervous system activity in chronic congestive heart failure.
Am J Cardiol 1997, 80(8):11011104. PubMed Abstract  Publisher Full Text

Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology
Circulation 1996, 93(5):10431065. PubMed Abstract  Publisher Full Text

Jiang W, Hathaway WR, McNulty S, Larsen RL, Hansley KL, Zhang Y, O'Connor CM: Ability of heart rate variability to predict prognosis in patients with advanced congestive heart failure.
Am J Cardiol 1997, 80(6):808811. PubMed Abstract  Publisher Full Text

Flapan AD, Goodfield NE, Wright RA, Francis CM, Neilson JM: Effects of digoxin on time domain measures of heart rate variability in patients with stable chronic cardiac failure: withdrawal and comparison group studies.
Int J Cardiol 1997, 59(1):2936. PubMed Abstract  Publisher Full Text

Binkley PF, Haas GJ, Starling RC, Nunziata E, Hatton PA, Leier CV, Cody RJ: Sustained augmentation of parasympathetic tone with angiotensinconverting enzyme inhibition in patients with congestive heart failure.
J Am Coll Cardiol 1993, 21(3):655661. PubMed Abstract

Huikuri HV, Makikallio TH, Airaksinen KE, Seppanen T, Puukka P, Raiha IJ, Sourander LB: Powerlaw relationship of heart rate variability as a predictor of mortality in the elderly.
Circulation 1998, 97(20):20312036. PubMed Abstract  Publisher Full Text
Prepublication history
The prepublication history for this paper can be accessed here: