- Open Access
Pitfalls of assessment of autonomic function by heart rate variability
Journal of Physiological Anthropologyvolume 38, Article number: 3 (2019)
Although analysis of heart rate variability is widely used for the assessment of autonomic function, its fundamental framework linking low-frequency and high-frequency components of heart rate variability with sympathetic and parasympathetic autonomic divisions has developed in the 1980s. This simplified framework is no longer able to deal with much evidence about heart rate variability accumulated over the past half-century. This review addresses the pitfalls caused by the old framework and discusses the points that need attention in autonomic assessment by heart rate variability.
Analysis of heart rate variability (HRV) is widely used as a standard method for assessing autonomic nervous functions. Particularly, low-frequency (LF, 0.04–0.15 Hz) and high-frequency (HF, 0.15–0.4 Hz) spectral components of HRV are used as the separate metrics of sympathetic and vagal (parasympathetic) functions [1, 2]. The association between HRV and autonomic function was first reported in the 1970s [3,4,5,6], and the currently prevailing frameworks of interpreting the HRV frequency components (LF and HF) in relation to the autonomic divisions were established in the 1980s [7,8,9]. Over the next half-century, more than 28,000 HRV studies were published, including a lot of evidence indicating important limitations of HRV frequency components as autonomic indices [10,11,12]. However, it does not seem that enough attention has been paid to such alarming evidence even by recent studies using HRV.
Analysis of HRV frequency components (LF and HF) provides a powerful tool for evaluating autonomic nervous function by simply recording an electrocardiogram (ECG). Furthermore, devices that automatically analyze HRV from ECGs, and also from pulse waves, are widely used in recent years. As a result, using the autonomic nervous index output there without criticism raises the risk of leading a wrong conclusion or a judgment that may harm health. Considering the recent knowledge of the mechanisms that generate HRV, the framework directly associating the HRV frequency components with the divisions of the autonomic nervous system (sympathetic and parasympathetic) has become too simplistic and it may be the time of reconstruction. This review will discuss the pitfalls where studies using HRV for autonomic assessment may fall based on the accumulated evidence of HRV.
General problems in the assessment of autonomic functions by HRV
Locality of autonomic neural regulations
The problem often seen in the studies of autonomic nerve function by HRV is interpretation as if HRV reflects the autonomic state of the whole body. The signals generating HRV are basically originated in the brain and mediated through the sympathetic and vagal nerves that innervate the sinoatrial node [13, 14]. Thus, the autonomic functions reflected in HRV are those regulating the pace-making function of the sinoatrial node. The autonomic nervous functions of other organ systems cannot be known from HRV. For example, food intake enhances parasympathetic (vagal) activity to stimulate digestive functions but it suppresses vagal activity to the heart to increase heart rate to a level that meets circulatory demand for digestion and absorption . Although sympathetic regulation may parallel among cardiac, renal, and muscle sympathetic nerves [16, 17], parasympathetic regulation is by function and therefore HRV cannot be used for estimating parasympathetic regulation other than that for cardiac pace-making function.
Short-term and long-term HRV analyses
Analysis of HRV is classified into two categories by the length of data recording . Short-term HRV is typically calculated over 5 min, and the magnitude of frequency components obtained by spectral analysis is used for autonomic function assessment. Long-term HRV is computed over nominal 24 h and time-domain  and frequency-domain  indices, nonlinear indices [20,21,22,23], and indices obtained by signal averaging [24,25,26] are used mainly for mortality risk prediction.
However, there is a more important difference between short-term and long-term HRV than the length of the data. For short-term HRV, data are usually recorded under monitoring while keeping subjects and measurement environment in specified conditions, while long-term HRV data are usually recorded by wearable sensors in subjects who are freely moving under daily activities. This difference is important for the interpretation of long-term HRV and the mechanisms by which it relates to clinical risks such as mortality.
For example, decreased 24-h standard deviation of normal-to-normal R-R interval (SDNN), representative long-term time-domain HRV, is an increased mortality risk after acute myocardial infarction [18, 20]. This relationship is often described as cardiac vagal dysfunction to increase the vulnerability of the heart to ventricular tachycardia and fibrillation, but the autonomic nervous function may not be a sole determinant of 24-h SDNN. In a recent study using the big data of simultaneously recorded 24-h ECG and tri-axial accelerograms, we found that the major determinant of 24-h SDNN is day/night difference in R-R interval that was associated with the level of daytime physical activity . This suggests that the association between reduced SDNN and adverse prognosis may be partly due to health conditions that restrict daily physical activity. Also, LF-to-HF power ratio (LF/HF) has been widely used as autonomic index, and it is well known that LF/HF increases with standing [7, 8, 28]. The major studies of the association between long-term HRV and disease prognosis reported that a decrease in LF/HF is an increased mortality risk [20, 22, 29]. In recent studies using big data of long-term HRV and accelerogram, we observed that the LF/HF reduces with the increasing lying period during 24 h (Fig. 1) [30, 31]. Decreased LF/HF in patients with adverse prognosis may reflect their daily life spending a long time in lying position due to health problems.
These indicate that long-term HRV is affected by physical activity and posture during the recording period and that the relationships with autonomic functions known for short-term HRV cannot be extended to the interpretation of long-term HRV. To use long-term HRV to evaluate autonomic nervous function during daily life, it is necessary to monitor at least physical activity and posture and analyze HRV in relation to them.
Separation of autonomic divisions by frequency
Why 0.15 Hz?
LF and HF component of HRV is divided at 0.15 Hz. Classification at this frequency was proposed by Sayers  and supported by animal experiments [32, 33]. Among those, a study by Berger et al.  performed the transfer function analysis of cardiac sinus rate response to the broad-frequency electric stimulation of the sympathetic and vagal nerves in dogs. They found that the sinus node responds as a low-pass filter to fluctuations in either sympathetic or vagal tone, while the filter for sympathetic stimulations has a lower corner frequency (0.15 Hz) than for vagal stimulations (> 0.5 Hz).
It is believed that the difference in frequency characteristics between sympathetic and vagal heart rate controls arises from the difference in the signaling mechanisms between beta-adrenergic and cholinergic receptors. Changes in sympathetic activity cause the effects through several phosphorylating enzymatic processes in the intracellular signal transduction mechanisms existing downstream of beta-adrenergic receptors. Consequently, the sympathetic modulation of heart rate cannot transfer fast fluctuations (> 0.15 Hz). On the other hand, changes in vagal activity cause the effects simply by the conformation change of the membrane potassium channels with Ach. Consequently, the vagal modulation of heart rate can transfer fluctuations up to a higher frequency range. The cutoff frequency of 0.15 Hz, however, is the value obtained from dogs, and there is no convincing evidence as to whether the value applies to humans of all ages uner all conditions.
Sinus arrhythmias constituting frequency component of HRV
The message brought about by studies showing the difference in frequency characteristics between sympathetic and vagal heart rate regulations is often over-interpreted. It is not saying that the magnitude of HF component of HRV reflects cardiac vagal function but is only saying that when heart rate fluctuations mediated by autonomic nerves are observed in HF band (> 0.15 Hz), it is mediated by the cardiac vagus. Indeed, even if HRV in HF band decreased or disappeared, it may not necessarily mean a withdrawal of cardiac vagal activity or cardiac vagal dysfunction. This could happen when the respiration frequency is out of the HF band, such as slow breathing < 0.15 Hz (9 breath/min) like during deep breathing, or fast breathing > 0.4 Hz (24 breath/min) during exercise and in children. Also, even when heart rate fluctuation is observed in HF band (> 0.15 Hz), if it is not mediated by autonomic nerves (as discussed in the next section), the framework of frequency separation of the autonomic heart rate regulation cannot be used.
Considering the knowledge of the mechanisms that generate HRV, the framework to directly associate the HRV frequency with the divisions of the autonomic nervous system (sympathetic and vagal) is too simplistic and already classical. HRV should be described at least dividing it into individual fluctuations from the corresponding physiological mechanisms. It is well known that there are two major components in short-term HRV; respiratory sinus arrhythmia (RSA) and a fluctuation at ~ 0.1 Hz called Mayer wave sinus arrhythmia (MWSA) [2, 7, 8]. The autonomic mechanisms of RSA are discussed in later sections. The mechanism mediating the heart rate fluctuation at ~ 0.1 Hz is believed to be the baroreceptor reflex cardiovascular regulation system, by which the fluctuation of arterial blood pressure at ~ 0.1 Hz that are known as Mayer wave  is reflected in the fluctuation of heart rate [12, 35, 36]. Therefore, we named this HRV component “Mayer wave sinus arrhythmia [15, 37].”
As mentioned above, when respiration frequency is < 0.15 Hz, RSA becomes a part of the LF component and the association between HF component and cardiac vagal function is lost. Additionally, when consciously breathing at ~ 0.1 Hz, it may cause resonance among respiration, heart rate, and blood pressure through the baroreceptor reflex mechanism, by which RSA and MWSA are merged into a large single oscillation at ~ 0.1 Hz [38, 39].
Non-vagal HRV component in HF band
Even if HRV is observed in the HF band, it may not necessarily be mediated by autonomic nerves. This phenomenon has been known for a long time and has been reported in different terms including complex HRV [40, 41], sinus node alternans , erratic heart rate [43, 44], and heart rate fragmentation (HRF; we will use this term in this review) [23, 45, 46]. This is a type of sinoatrial instability characterized by frequent (sometimes every beat) appearance of peak and valley in R-R interval time series despite that ECG shows sinus rhythm. By this phenomenon, switching between increases and decreases in heart rate occurs more frequently than those expected for that mediated by the vagus.
To quantify the level of HRF, Costa et al. [45, 46] developed a set of metrics including the percentage of the inflection point (PIP) and symbol dynamics . Using these measures, they demonstrated that the level of HRF increases with advancing age and in patients with coronary artery disease  and that its increase is an increased risk of cardiovascular events and death . Because the existence of this phenomenon appears as increases in the conventional HRV metrics such as the percentage of pairs of adjacent normal-to-normal R-R intervals differing by > 50 ms (pNN50), the square root of the mean of squared differences between adjacent normal-to-normal R-R intervals (RMSSD), and HF power, it may confound their associations with cardiac vagal function and with disease prognosis. In fact, this phenomenon is considered a potential cause of lesser prognostic relevance of the HF component than that of other HRV frequency components such as very-low-frequency (0.0033–0.04 Hz) and LF components . Although the exact mechanisms for this phenomenon are unclear, alternative pacemaker shift within the sinus node caused by sinoatrial degeneration or disorganization may be a potential mechanism.
RSA and cardiac vagal function
Non-vagal factors influencing RSA
Although the association between HF component of short-term HRV and cardiac vagal function depends on the neural mechanisms of RSA [47, 48], many factors are known to influence the magnitude of RSA besides the cardiac vagal function. An increase in respiratory frequency decreases the amplitude of RSA (the range of excursion of R-R interval with respiration) [49,50,51]. Also, an increase in tidal volume increases RSA amplitude [49,50,51,52,53]. These effects of respiratory parameters on the RSA operate independently of the level of cardiac vagal activity . Therefore, when using RSA amplitude to evaluate cardiac vagal function, it is desirable to control respiration parameters in experimental protocol or to adjust their effects statistically.
Also, RSA amplitude may be affected by sympathetic nervous activity. Even though the sympathetic nervous system cannot transfer HRV > 0.15 Hz , it may restrict the magnitude of cardiac vagal modulations of heart rate [55, 56]. Taylor et al.  observed that beta-adrenergic blockade increases RSA amplitude over a wide frequency range including > 0.15 Hz and concluded that cardiac sympathetic outflow reduces heart rate oscillations at all frequencies including RSA.
Neural mechanisms of RSA
Despite that RSA amplitude is widely used to estimate cardiac vagal activity that is controlling the level of heart rate, RSA amplitude and the level of heart rate are known to be controlled by different groups of vagal motor neurons . Preganglionic cardiac vagal motor neurons are known to locate separately in the dorsal motor nucleus of the vagus (DVN) and in the nucleus ambiguous (nA) . In cats and rats, C-fiber cardiac preganglionic neurons in the DVN shows regular ongoing activity that is unaffected by respiratory rhythm, while those in the nA fire with respiratory rhythm . During embryonic development, the cardiac preganglionic neurons forming nA migrate ventrolaterally from a dorsomedial position, possibly of the DVN . In axolotl, this neuron relocation occurs with the metamorphosis of the onset of air-breathing . These suggest that during the evolution of air-breathing vertebrate, the nA was developed to “generate” RSA in addition to the DVN to control the level of heart rate.
Biological function of RSA
Beside its relationship to cardiac vagal function, RSA itself is recognized as a physiological mechanism to ensure optimal ventilation-perfusion matching within the lungs of air-breathing mammals [11, 57, 61,62,63]. A similar matching of the counter-perfusions with water and blood of the gills is also known in aquatic animals without the nA, but this matching is accomplished simply by adjusting the level of heart rate with ventilatory gill movements . In air-breathing animals with reciprocal respiration, the presence of RSA may improve pulmonary gas exchange efficiency by matching the timing of alveolar air volume and alveolar perfusion within each respiratory cycle (Fig. 2) . Because this mechanism could reduce the intrapulmonary shunt ratio by reducing ineffective heartbeats when the alveoli are contracting and it could also suppress the functional dead space by reducing wasted inspiration when alveolar perfusion is small, RSA seems a biological function to save the energy for circulation and respiration at rest—the period when the energy saving of the pulmonary circulation is more important to survive against hunger than responding to the demand for systemic circulation .
Possible reason why RSA is associated with cardiac vagal activity
If RSA and heart rate are regulated by different vagal motor neurons, direct relationship may not be expected between the magnitude of RSA and cardiac vagal activity to control heart rate. In fact, Goldberger et al. [65, 66] demonstrated the evidence of dissociation between RSA magnitude and cardiac vagal activity. They observed that bradycardia caused by cardiac vagal activation through baroreceptor stimulation with phenylephrine that increases arterial blood pressure is accompanied by a paradoxical decrease in RSA amplitude. This indicates that the degree of respiratory modulation of cardiac vagal outflow is not always parallel to cardiac vagal tone (mean activity level).
RSA as an indicator of intrinsic resting function
Nevertheless, why is the RSA amplitude widely used as an indicator of cardiac vagal activity? This is probably because RSA is an intrinsic resting function that appears when an individual is about to rest. In such situations, the nA and DVN are expected to work cooperatively for the common purpose—energy saving (Fig. 3). This apparently indicates that RSA amplitude is a limited indicator of cardiac vagal activity, but at the same time, it suggests that RSA amplitude may be a useful indicator of resting function. The latter assumption seems supported by much evidence about the factors affecting RSA amplitude; i.e., RSA increases with supine posture , relaxation , and sleep , decreases with standing [7, 8, 28], physical exercise [69,70,71], mental stress , cardiovascular stress including baroreceptor stimulation [65, 66], advancing age , and compromised cardiac function [37, 74].
LF component and sympathetic function
Although the power of the LF component and its derived indices (normalized LF power and LF/HF) have been widely used as an indicator of sympathetic functions (activity, predominance, sympatho-vagal balance, etc.), the validity of this use has been denied by a lot of evidence accumulated by numerous studies . Unlike the studies specific to RSA in the HF component, this conclusion does not change, even if the discussion is specified to MWSA. Earlier studies proposed the LF as sympathetic component , but the power of LF spectrum is not eliminated completely by beta-adrenergic blockade and rather largely reduced by sinoatrial parasympathectomy . Although the increase in LF power with standing or head-up tilting is often cited as the evidence for its association with sympathetic activity, postural increase in LF power is observed only in one-third of healthy subjects, no significant change in other third, and rather decreases in the remaining third . While LF/HF increases with standing, this is mainly due to the postural decrease in HF power. Finally, LF power is not related to cardiac sympathetic innervation quantified by positron emission tomographic neuroimaging [36, 76]. There is no longer a convincing physiological basis justifying the assessment of cardiac sympathetic nervous functions whether by absolute or normalized LF power or LF/HF.
This paper described the major pitfalls of the assessment of autonomic function by HRV. Overall, HRV tends to be used beyond the limit of its capacity as a quantitative measure of autonomic function. Particularly, the framework to associate the HRV frequency components (LF and HF) with the divisions of the autonomic nervous system (sympathetic and parasympathetic) is already too simplistic, and it is the major cause of the pitfalls.
On the other hand, even if HRV is a limited measure of autonomic functions, it is still a fact that HRV is there. Except for some exception such as HRF, HRV is thought to originate from the brain and to transfer to the heart through the autonomic nervous system. Like the radio that extracts information from radio waves modulated by the broadcasting station, HRV analysis can intercept the information for circulatory regulation generated by the brain. We will be able to expect the continued development of HRV researches by aiming for a faithful and unbiased extraction of the information contained in HRV beyond the classic framework of LF and HF.
Dorsal motor nucleus of the vagus
Heart rate fragmentation
Heart rate variability
Mayer wave sinus arrhythmia
Percentage of inflection point
Percentage of pairs of adjacent normal-to-normal R-R intervals differing by > 50 ms
Square root of the mean of squared differences between adjacent normal-to-normal R-R intervals
Respiratory sinus arrhythmia
Standard deviation of normal-to-normal R-R interval
Appel ML, Berger RD, Saul JP, Smith JM, Cohen RJ. Beat to beat variability in cardiovascular variables: noise or music? J Am Coll Cardiol. 1989;14:1139–48.
Camm AJ, Malik M, Bigger JT Jr, Breithardt G, Cerutti S, Cohen RJ, Coumel P, Fallen EL, Kleiger RE, Lombardi F, Malliani A, Moss AJ, Rottman JN, Schmidt G, Schwartz PJ, Singer DH, 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(5):1043–65.
Katona PG, Poitras JW, Barnett GO, Terry BS. Cardiac vagal efferent activity and heart period in the carotid sinus reflex. Am J Phys. 1970;218(4):1030–7.
Hyndman BW, Kitney RI, Sayers BM. Spontaneous rhythms in physiological control systems. Nature. 1971;233:339–41.
Sayers BM. Analysis of heart rate variability. Ergonomics. 1973;16(1):17–32.
Hyndman BW, Gregory JR. Spectral analysis of sinus arrhythmia during mental loading. Ergonomics. 1975;18(3):255–70.
Pomeranz B, Macaulay RJ, Caudill MA, Kutz I, Adam D, Gordon D, Kilborn KM, Barger AC, Shannon DC, Cohen RJ, et al. Assessment of autonomic function in humans by heart rate spectral analysis. Am J Phys. 1985;248(1 Pt 2):H151–3.
Pagani M, Lombardi F, Guzzetti S, Rimoldi O, Furlan R, Pizzinelli P, Sandrone G, Malfatto G, Dell'Orto S, Piccaluga E, Turiel M, Baselli G, Cerutti S, Malliani A. Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ Res. 1986;59:178–93.
Billman GE. Heart rate variability - a historical perspective. Front Physiol. 2011;2:86.
Berntson GG, Bigger JT Jr, Eckberg DL, Grossman P, Kaufmann PG, Malik M, Nagaraja HN, Porges SW, Saul JP, Stone PH, Van der Molen MW. Heart rate variability: origins, methods, and interpretive caveats. Psychophysiology. 1997;34(6):623–48.
Hayano J, Yasuma F. Hypothesis: respiratory sinus arrhythmia is an intrinsic resting function of cardiopulmonary system. Cardiovasc Res. 2003;58(1):1–9.
Billman GE. The LF/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front Physiol. 2013;4:26.
Sands KEF, Appel ML, Lilly LS, Schoen FJ, Mudge GH Jr, Cohen RJ. Power spectrum analysis of heart rate variability in human cardiac transplant recipients. Circulation. 1989;79:76–82.
Hayano J, Sakakibara Y, Yamada A, Yamada M, Mukai S, Fujinami T, Yokoyama K, Watanabe Y, Takata K. Accuracy of assessment of cardiac vagal tone by heart rate variability in normal subjects. Am J Cardiol. 1991;67(2):199–204.
Hayano J, Sakakibara Y, Yamada M, Kamiya T, Fujinami T, Yokoyama K, Watanabe Y, Takata K. Diurnal variations in vagal and sympathetic cardiac control. Am J Phys. 1990;258(3 Pt 2):H642–6.
Kamiya A, Kawada T, Yamamoto K, Michikami D, Ariumi H, Miyamoto T, Shimizu S, Uemura K, Aiba T, Sunagawa K, Sugimachi M. Dynamic and static baroreflex control of muscle sympathetic nerve activity (SNA) parallels that of renal and cardiac SNA during physiological change in pressure. Am J Phys. 2005;289(6):H2641–8.
Kamiya A, Kawada T, Yamamoto K, Michikami D, Ariumi H, Miyamoto T, Uemura K, Sugimachi M, Sunagawa K. Muscle sympathetic nerve activity averaged over 1 minute parallels renal and cardiac sympathetic nerve activity in response to a forced baroreceptor pressure change. Circulation. 2005;112(3):384–6.
Kleiger RE, Miller JP, Bigger JT Jr, Moss AJ, the Multicenter Post-Infarction Research G. Decreased heart rate variability and its association with increased mortality after acute myocardial infarction. Am J Cardiol. 1987;59:256–62.
Bigger JT Jr, Fleiss JL, Steinman RC, Rolnitzky LM, Kleiger RE, Rottman JN. Frequency domain measures of heart period variability and mortality after myocardial infarction. Circulation. 1992;85:164–71.
Huikuri HV, Mäkikallio TH, Peng CK, Goldberger AL, Hintze U, Moller M, Grp DS. Fractal correlation properties of R-R interval dynamics and mortality in patients with depressed left ventricular function after an acute myocardial infarction. Circulation. 2000;101:47–53.
Kiyono K, Hayano J, Watanabe E, Struzik ZR, Yamamoto Y. Non-Gaussian heart rate as an independent predictor of mortality in patients with chronic heart failure. Heart Rhythm. 2008;5(2):261–8.
Suzuki M, Hiroshi T, Aoyama T, Tanaka M, Ishii H, Kisohara M, Iizuka N, Murohara T, Hayano J. Nonlinear measures of heart rate variability and mortality risk in hemodialysis patients. Clin J Am Soc Nephrol. 2012;7(9):1454–60.
Costa MD, Redline S, Davis RB, Heckbert SR, Soliman EZ, Goldberger AL. Heart rate fragmentation as a novel biomarker of adverse cardiovascular events: the multi-ethnic study of atherosclerosis. Front Physiol. 2018;9:1117.
Schmidt G, Malik M, Barthel P, Schneider R, Ulm K, Rolnitzky L, Camm AJ, Bigger JT Jr, Schomig A. Heart-rate turbulence after ventricular premature beats as a predictor of mortality after acute myocardial infarction. Lancet. 1999;353:1390–6.
Bauer A, Kantelhardt JW, Barthel P, Schneider R, Makikallio T, Ulm K, Hnatkova K, Schomig A, Huikuri H, Bunde A, Malik M, Schmidt G. Deceleration capacity of heart rate as a predictor of mortality after myocardial infarction: cohort study. Lancet. 2006;367(9523):1674–81.
Hayano J, Yasuma F, Watanabe E, Carney RM, Stein PK, Blumenthal JA, Arsenos P, Gatzoulis KA, Takahashi H, Ishii H, Kiyono K, Yamamoto Y, Yoshida Y, Yuda E, Kodama I. Blunted cyclic variation of heart rate predicts mortality risk in post-myocardial infarction, end-stage renal disease, and chronic heart failure patients. Europace. 2017;19:1392–400.
Hayano J, Yuda E, Furukawa Y, Yoshida Y. Association of 24-hour heart rate variability and daytime physical activity: ALLSTAR big data analysis. Int J Biosci Biochem Bioinformatics. 2018;8(1):61–7.
Hayano J, Mukai S, Fukuta H, Sakata S, Ohte N, Kimura G. Postural response of low-frequency component of heart rate variability is an increased risk for mortality in patients with coronary artery disease. Chest. 2001;120:1942–52.
La Rovere MT, Bigger JT Jr, Marcus FI, Mortara A, Schwartz PJ, Investigators A. Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. Lancet. 1998;351:478–84.
Yoshida Y, Ogasawara H, Yuda E, Hayano J. What does LF/HF of heart rate variability in ambulatory ECG mean? Effect of time in lying position during monitoring. Eur Heart J. 2016;37(suppl):96-97.
Yoshida Y, Furukawa Y, Ogasawara H, Yuda E, Hayano J, editors. Longer lying position causes lower LF/HF of heart rate variability during ambulatory monitoring. Kyoto: 2016 IEEE 5th Global Conference on Consumer Electronics (GCCE); 2016.
Akselrod S, Gordon D, Ubel FA, Shannon DC, Barger AC, Cohen RJ. Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science. 1981;213:220–2.
Berger RD, Saul JP, Cohen RJ. Transfer function analysis of autonomic regulation. I. Canine atrial rate response. Am J Phys. 1989;256(1 Pt 2):H142–52.
Penaz J. Mayer waves: history and methodology. Automedica. 1978;2:135-141.
Madwed JB, Albrecht P, Mark RG, Cohen RJ. Low-frequency oscillation in arterial pressure and heart rate: a simple computer model. Am J Phys. 1991;256:H1573–H9.
Rahman F, Pechnik S, Gross D, Sewell L, Goldstein DS. Low frequency power of heart rate variability reflects baroreflex function, not cardiac sympathetic innervation. Clin Auton Res. 2011;21(3):133–41.
Hayano J, Sakakibara Y, Yamada M, Ohte N, Fujinami T, Yokoyama K, Watanabe Y, Takata K. Decreased magnitude of heart rate spectral components in coronary artery disease. Its relation to angiographic severity. Circulation. 1990;81(4):1217–24.
Vaschillo EG, Vaschillo B, Lehrer PM. Characteristics of resonance in heart rate variability stimulated by biofeedback. Appl Psychophysiol Biofeedback. 2006;31(2):129–42.
Lehrer PM, Gevirtz R. Heart rate variability biofeedback: how and why does it work? Front Psychol. 2014;5:756.
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:565–9.
Woo MA, Stevenson WG, Moser DK, Trelease RB, Harper RM. Patterns of beat-to-beat heart rate variability in advanced heart failure. Am Heart J. 1992;123:704–10.
Binkley PF, Eaton GM, Nunziata E, Khot U, Cody RJ. Heart rate alternans. Ann Intern Med. 1995;122(2):115–7.
Stein PK, Domitrovich PP, Hui N, Rautaharju P, Gottdiener J. Sometimes higher heart rate variability is not better heart rate variability: results of graphical and nonlinear analyses. J Cardiovasc Electrophysiol. 2005;16(9):954–9.
Stein PK, Le Q, Domitrovich PP, Investigators C. Development of more erratic heart rate patterns is associated with mortality post-myocardial infarction. J Electrocardiol. 2008;41(2):110–5.
Costa MD, Davis RB, Goldberger AL. Heart rate fragmentation: a symbolic dynamical approach. Front Physiol. 2017;8:827.
Costa MD, Davis RB, Goldberger AL. Heart rate fragmentation: a new approach to the analysis of cardiac interbeat interval dynamics. Front Physiol. 2017;8:255.
Taylor EW. The evolution of efferent vagal control of the heart in vertebrates. Cardioscience. 1994;5(3):173–82.
Taylor EW, Hoar WS, Randall DJ, Farrell AP. Nervous control of the heart and cardiorespiratory interactions. Fish Physiology vol XII B. New York: Academic Press; 1992. p. 343–87.
Melcher A. Respiratory sinus arrhythmia in man: a study in heart rate regulating mechanisms. Acta Physiol Scand. 1976;Suppl 435:1–31.
Hirsch JA, Bishop B. Respiratory sinus arrhythmia in humans: how breathing pattern modulates heart rate. Am J Phys. 1981;241:H620–H9.
Brown TE, Beightol LA, Koh J, Eckberg DL. Important influence of respiration on human R-R interval power spectra is largely ignored. J Appl Physiol. 1993;75:2310–7.
Eckberg DL. Human sinus arrhythmia as an index of vagal cardiac outflow. J Appl Physiol. 1983;54:961–6.
Kollai M, Mizsei G. Respiratory sinus arrhythmia is a limited measure of cardiac parasympathetic control in man. J Physiol (Lond). 1990;424:329–42.
Hayano J, Mukai S, Sakakibara M, Okada A, Takata K, Fujinami T. Effects of respiratory interval on vagal modulation of heart rate. Am J Phys. 1994;267(1 Pt 2):H33–40.
Taylor JA, Myers CW, Halliwill JR, Seidel H, Eckberg DL. Sympathetic restraint of respiratory sinus arrhythmia: implications for vagal-cardiac tone assessment in humans. Am J Phys. 2001;280:H2804–14.
Cohen MA, Taylor JA. Short-term cardiovascular oscillations in man: measuring and modelling the physiologies. J Physiol. 2002;542(Pt 3):669–83.
Taylor EW, Jordan D, Coote JH. Central control of the cardiovascular and respiratory systems and their interactions in vertebrates. Physiol Rev. 1999;79:855–916.
Jones JFX, Wang Y, Jordan D. Activity of C-fiber cardiac vagal efferents in anaesthetized cats and rats. J Physiol (Lond). 1998;507:869–80.
Windle WF. Neurofibrillar development in the central nervous system of cat embryos between 8 and 12 mm long. J Comp Neurol. 1933;58:643–723.
Taylor EW, Al-Ghamdi MS, Ihmied IH, Wang T, Abe AS. The neuranatomical basis of central control of cardiorespiratory interactions in vertebrates. Exp Physiol. 2001;86(6):771–6.
Hayano J, Yasuma F, Okada A, Mukai S, Fujinami T. Respiratory sinus arrhythmia. A phenomenon improving pulmonary gas exchange and circulatory efficiency. Circulation. 1996;94(4):842–7.
Yasuma F, Hayano J. Respiratory sinus arrhythmia: why does the heartbeat synchronize with respiratory rhythm? Chest. 2004;125(2):683–90.
Ito S, Sasano H, Sasano N, Hayano J, Fisher JA, Katsuya H. Vagal nerve activity contributes to improve the efficiency of pulmonary gas exchange in hypoxic humans. Exp Physiol. 2006;91(5):935–41.
Satchell GH. The reflex co-ordination of the heart beats with respiration in dogfish. J Exp Biol. 1960;37:719–31.
Goldberger JJ, Ahmed MW, Parker MA, Kadish AH. Dissociation of heart rate variability from parasympathetic tone. Am J Phys. 1994;266:H2152–H7.
Goldberger JJ, Challapalli S, Tung R, Parker MA, Kadish AH. Relationship of heart rate variability to parasympathetic effect. Circulation. 2001;103(15):1977–83.
Sakakibara M, Takeuchi S, Hayano J. Effect of relaxation training on cardiac parasympathetic tone. Psychophysiology. 1994;31:223–8.
Bonnet MH, Arand DL. Heart rate variability: sleep stage, time of night, and arousal influences. Electroencephalogr Clin Neurophysiol. 1997;102:390–6.
Arai Y, Saul JP, Albrecht P, Hartley LH, Lilly LS, Cohen RJ, Colucci WS. Modulation of cardiac autonomic activity during and immediately after exercise. Am J Phys. 1989;256:H132–H41.
Yamamoto Y, Hughson RL, Peterson JC. Autonomic control of heart rate during exercise studied by heart rate variability spectral analysis. J Appl Physiol. 1991;71:1136–42.
Taylor JA, Hayano J, Seals DR. Lesser vagal withdrawal during isometric exercise with age. J Appl Physiol. 1995;79:805–11.
Sakakibara M, Kanematsu T, Yasuma F, Hayano J. Impact of real-world stress on cardiorespiratory resting function during sleep in daily life. Psychophysiology. 2008;45(4):667–70.
Shannon DC, Carley DW, Benson H. Aging of modulation of heart rate. Am J Phys. 1987;253:H874–H7.
Saul JP, Arai Y, Berger RD, Lilly LS, Colucci WS, Cohen RJ. Assessment of autonomic regulation in chronic congestive heart failure by heart rate spectral analysis. Am J Cardiol. 1988;61:1292–9.
Randall DC, Brown DR, Raisch RM, Yingling JD, Randall WC. SA nodal parasympathectomy delineates autonomic control of heart rate power spectrum. Am J Phys. 1991;260(3 Pt 2):H985–8.
Moak JP, Goldstein DS, Eldadah BA, Saleem A, Holmes C, Pechnik S, Sharabi Y. Supine low-frequency power of heart rate variability reflects baroreflex function, not cardiac sympathetic innervation. Heart Rhythm. 2007;4(12):1523–9.
There is no funding source for this paper.
Availability of data and materials
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.