Depression affects autonomic system of the body? Yes, it does!
Sheena Singla, Shikha Jhamb, Kamal Dev Singh, Avnish Kumar
Department of Physiology, Government Medical College, Patiala, Punjab, India
|Date of Submission||19-Oct-2019|
|Date of Acceptance||18-Mar-2020|
|Date of Web Publication||31-Aug-2020|
Dr. Shikha Jhamb
House No. E-8, Government Medical College, Patiala - 147 001, Punjab
Source of Support: None, Conflict of Interest: None
BACKGROUND: Depression is a disorder of the brain, a state of low mood and aversion to activity that can affect a person's thoughts, behavior, feelings, and sense of well-being and can cause alteration in sympathetic activity of the body, thus affecting heart rate variability (HRV).
AIM: The study was conducted to determine the effects of depression on HRV parameters in clinically known cases of depression.
MATERIALS AND METHODS: A cross-sectional (observational) study was conducted on fifty known patients of depression, aged 18–65 years, and fifty healthy, age-matched, normal controls. HRV test was conducted, and the results were analyzed statistically using Student's “t”-test for equal variance, for various parameters.
RESULTS: The study showed that the differences in mean R-R interval(s), mean HR (beats/min), very low frequency (VLF) power (milliseconds squared [ms2]), and low frequency/high frequency (LF/HF) power (%) (LF/HF ratio) among controls and cases were statistically highly significant, while the differences in LF power (ms2) in controls and cases were statistically significant. However, the differences in basic anthropometric parameters, STD (standard deviation of the normal to normal)(s), root mean square successive difference (ms), NN50 (count), pNN50 (%), VLF peak (Hz), LF peak (Hz), HF peak (Hz), HF power (ms2), VLF power (%), LF power (%), HF power (%), LF power (nu), and HF power (nu) were found to be statistically insignificant.
CONCLUSION: Depression leads to changes in autonomic control of the body and changes the autonomic balance in favor of an increased sympathetic tone, which can be detected with fair accuracy with HRV analysis.
Keywords: Cases, controls, depression, heart rate variability, patients
|How to cite this article:|
Singla S, Jhamb S, Singh KD, Kumar A. Depression affects autonomic system of the body? Yes, it does!. J Edu Health Promot 2020;9:217
|How to cite this URL:|
Singla S, Jhamb S, Singh KD, Kumar A. Depression affects autonomic system of the body? Yes, it does!. J Edu Health Promot [serial online] 2020 [cited 2020 Oct 28];9:217. Available from: https://www.jehp.net/text.asp?2020/9/1/217/293948
| Introduction|| |
Depression (major depressive disorder or clinical depression) is a common but serious mood disorder. It causes severe symptoms that affect how you feel, think, and handle daily activities, such as sleeping, eating, or working.
Depression is a disorder of the brain. There are a variety of causes, including genetic, biological, environmental, and psychological factors. Depression can happen at any age, but it often begins in teens and young adults. It is much more common in women. Women can also get postpartum depression after the birth of a baby. Some people get seasonal affective disorder in the winter. Depression is one part of bipolar disorder.
Anxiety disorders increase the risk of future cardiovascular disease (CVD) and mortality, even after controlling possible confounding factors such as smoking, lifestyle, and socioeconomic status, irrespective of a history of medical disorders. While impaired vagal function, indicated by reductions in heart rate variability (HRV), may be one mechanism linking anxiety disorders to CVD, prior studies have reported inconsistent findings highlighting the need for meta-analysis.
To be diagnosed with depression, the symptoms must be present for at least 2 weeks.
Symptoms can include:
- Feeling sad or “empty”
- Loss of interest in favorite activities
- Overeating, or not wanting to eat at all
- Not being able to sleep, or sleeping too much
- Feeling very tired
- Feeling hopeless, irritable, anxious, or guilty
- Aches or pains, headaches, cramps, or digestive problems
- Thoughts of death or suicide.
Depression is most common in ages 18 to 25 (10.9%) and in individuals belonging to two or more races (10.5%). Women are twice as likely as men to have had a depressive episode, according to the NIMH and the World Health Organization (WHO). The WHO estimates that more than 300 million people worldwide suffer from depression. It is also the world's leading cause of disability.
Depression is also associated with elevated HR and reduced HRV, which are known risk factors for cardiac morbidity and mortality that may explain the increased risk associated with depression.
Heart rate variability
The human heart beat in a healthy individual is neither absolutely regular nor completely random. It varies based on the interplay between many factors including physical and mental stress, exercise, respiration, thermoregulation, blood pressure regulation, renin–angiotensin system, circadian rhythm, and other unknown complex mechanisms. This subtle fluctuation in sinus rhythm is known as HRV. HRV has proved to be a more sensitive tool for the detection of autonomic balance than mean HR. It is seen as reflective and predictive of general health and overall physiological wellness. HRV is one of the promising markers of autonomic activity. Variations in the HR are dependent on the activities of the sinoatrial (SA) node, which is the natural pacemaker of the heart innervated by the sympathetic and parasympathetic branches of the autonomic nervous system. The degree of balance between sympathetic and vagal nerve activity determines HRV.
Currently, HRV is recognized as an effective and noninvasive tool for evaluating autonomic nervous system regulations of the heart. HRV describes the variation between consecutive heart beats. The rhythm of the heart is controlled by SA node, which is modulated by both sympathetic and parasympathetic branches of the autonomic nervous system. Sympathetic activity tends to increase HR and its response is slow (few seconds). Parasympathetic activity, on the other hand, tends to decrease HR and mediates faster (0.2–0.6 s).
The parasympathetic activity's influence on HR is modulated by acetylcholine released by the vagus nerve on the SA node and the sympathetic influence by the release of epinephrine and norepinephrine.
The continuous modulation of the sympathetic and parasympathetic innervations results in variation in HRs. The most conspicuous periodic component of HRV is the so-called respiratory sinus arrhythmia which is considered to range from 0.15 to 0.4 Hz. In addition to the physiological influence of breathing on HRV, this high-frequency (HF) component is generally believed to be of parasympathetic origin. Another widely studied component of HRV is the low-frequency (LF) component usually ranging from 0.04 to 0.15 Hz including the component referred to as the 10-s rhythm or the Mayer wave. The rhythms within the LF band have been thought to be of both sympathetic and parasympathetic origin, even though some researchers have suggested them to be mainly of sympathetic origin. The very low frequency (VLF) component accounts for all other HR changes including those associated with thermoregulation and hormonal and local factors.
The HRV is evaluated by two ways: time domain analysis and frequency domain analysis.
Frequency domain method, which is also known as power spectral density, provides the basic information on how power as an expression of variance distributes in the function of frequency. VLF, LF, and HF powers are usually measured in absolute values of power (milliseconds squared [ms2]). LF and HF can be also measured in normalized units (nu) to emphasize the controlled and balanced behavior of the two branches of the autonomic nervous system, as well as baroreflex responsiveness to beat-to-beat variation in arterial blood pressure. Normalization of LF and HF power tends to minimize the effect of the changes in the total power on the values of these two components. Normalized units and absolute values of LF and HF power should both be calculated to provide a better measurement of the degree of autonomic modulation rather than just the level of autonomic tone.
The following parameters are calculated in time domain: mean R-R interval (sec), the standard deviation of intervals between normal beats (SDNN), mean HR (per minute), the root mean square differences of successive NN intervals (RMSSD), the number of interval differences of successive NN intervals >50 ms (NN50), a measure of the number of adjacent NN intervals which differ by >50 ms (pNN50). On a standard electrocardiogram (ECG), the maximum upward deflection of a normal QRS complex is at the peak of the R-wave, and the duration between two adjacent R-wave peaks is termed as the R-R interval.
According to a research conducted by the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, in a time-domain analysis, higher values of the SDNN intervals, RMSSD, and percentage of adjacent normal R–R intervals 50 ms (pNN50) indicate stronger parasympathetic dominance. Among the frequency-domain measurement indicators, a HF denotes parasympathetic activity, and a LF is related to the activity of both the sympathetic nervous system and the parasympathetic nervous system.
The term HRV refers, in general, to changes in heart beat interval which is a reciprocal of the HR. The starting point for HRV analysis is the ECG recording from which the HRV time series can be extracted. In the formulation of HRV time series, a fundamental issue is the determination of heart beat period.
A study conducted in 2002 shows that in physically healthy depressed adults, HRV does not differ from healthy controls.
A study conducted in 2008 shows that depression is associated with significantly lowered HRV. The present study was necessitated because of conflicting reports. The aim of the present study was to assess the effect of depression on HRV, by conducting study on clinically known patients of depression in the age group of 18–65 years.
| Materials and Methods|| |
A cross-sectional (observational) study was conducted in the department of physiology in the institute during November 2017–November 2018. The participants were divided into two groups: (i) fifty known patients of depression (cases) and (ii) fifty age-matched, normal healthy subjects (controls). The study protocol was approved by the ethical committee of the institute. Written consent was obtained from all the enrolled individuals after explaining them the details of the study in their own language. The participants for the study (cases) were enrolled from patients attending the outdoor clinic of the Psychiatry Department of Government Medical College and Rajindra Hospital, Patiala, whereas the controls were volunteers from among the paramedical and lower staff of the institute.
The subject selection is based on exclusion–inclusion criteria.
Inclusion criteria (Patients were on regular treatment with antidepressants like escitalopram, sertraline).
- Individuals in the age range of 18–65 years
- Diagnosed and known cases of depression. Patients were on antidepressants like escitalopram, sertraline, fluoxetine.
Exclusion criteria (Hypertension, renal, endocrine, cardiovascular and respiratory disorders are known to affect autonomic system of the body in many ways; that's why these were chosen as exclusion criteria).
- Patients with hypertension
- Patients with renal and endocrine disorders
- Patients with cardiovascular disorders
- Patients with respiratory disorders
In every case selected, thorough and detailed history of depression was taken. The participants were diagnosed in accordance with the WHO's depression diagnostic Criteria ICD-10 and Hamilton's Rating Scale by the consultant psychiatrist.
Heart rate variability analysis
The HR of each participant was recorded by ECG monitoring, in RR mode (beat to beat), for 5 min at rest, in supine position, using “physiopac hardware” by Medicaid (manufacturers at SAS Nagar, Mohali, Punjab, India); HRV analysis was done as per the procedure given in the manual supplied by RMS Recorders and Medicare Systems (P) Ltd. (Chaunki, Haryana, India). The participant was allowed to relax on a comfortable chair with his/her back toward the recording machine. The following anthropometric parameters were measured according to standardized techniques: height (cm), weight (kg), and body mass index (BMI): BMI was calculated using Quetlet's Index: BMI = weight (kg)/height2 (in m); body surface area (BSA): BSA was calculated using Dubois and Dubois formula: BSA in m2 = 0.007184 × weight0.425 (kg) × height (m)0.725
The results were compiled in the form of charts and graphs and statistically analyzed by using Student's “t”-test. Evaluation of the data was carried out by Microsoft Excel and SPSS/PC+(a statistical program to carry out analysis). “P” < 0.05 and 0.01 was considered statistically significant and highly significant, respectively.
| Results|| |
The comparison of age (in years), height (in centimeters), and weight (in kg) of controls and cases was done as shown in [Table 1]. There were no statistically significant differences between the basic data: age, height, and weight of both controls and cases.
|Table 1: Comparison of mean and standard deviation of anthropometric parameters in controls and cases|
Click here to view
The comparison of baseline variables, i.e., BMI (kg/m2) and BSA (m2) of controls and cases, was done as shown in [Table 2]. The differences in BSA and BMI of both controls and cases were found to be statistically insignificant.
The comparison of time domain analysis, i.e., mean RR(s), STD (SDNN)(s), mean HR (beats/min), RMSSD (ms), NN50 (count), and pNN50(%) of controls and cases, was done as shown in [Table 3]. The differences in mean RR and mean HR of controls and cases were statistically highly significant, with the values of mean RR higher among controls as compared to cases, and the values of mean HR higher among cases as compared to controls. The differences in STD (SDNN), RMSSD (ms), NN50 (count), and pNN50 (%) were found to be statistically insignificant.
The comparison of frequency domain analysis, i.e., VLF peak (Hz), LF peak (Hz), HF peak (Hz), VLF power (ms2), LF power (ms2), HF power (ms2), VLF power (%), LF power (%), HF power (%), LF/HF ratio, LF power (nu), and HF power (nu), of controls and cases was done as shown in [Table 4]. The differences in VLF (power) (ms2) and LF/HF (power %) (LF/HF ratio) of controls and cases were statistically highly significant, with values being much higher in cases as compared to controls, while the differences in LF (power) (ms2) were statistically significant, with the values higher among cases. The differences in VLF (peak) (Hz), LF (peak) (Hz), HF (peak) (Hz), HF (power) (ms2), VLF power (%), LF power (%), HF power (%), LF (power) (nu), and HF (power) (nu) were statistically insignificant.
| Discussion|| |
The study was conducted on fifty patients of depression and fifty normal healthy controls to observe the effect of depression on HR, in the age group of 18–65 years.
We studied mean RR, STD (SDNN), mean HR, RMSSD, NN50, and pNN50 in time domain analysis and VLF peak, LF peak, HF peak, VLF power, LF power, HF power, VLF power%, LF power%, HF power%, LF/HF ratio, LF power nu, and HF power nu in frequency domain analysis.
The LF component of HRV analysis reflects the sympathetic parameters, whereas the HF component of HRV analysis shows the parasympathetic activity of the body. If the values of LF component increase, it shows that the balance has shifted in favor of sympathetic system, and hence the sympathetic system in the body is overactive.
Higher values of LF/HF ratio also show increased sympathetic activity in the body.
A detailed comparison between our study and those derived from various other well-documented studies was done.
Comparison of anthropometric parameters
[Table 1] and [Table 2] show the comparison of anthropometric parameters of controls and cases. There was no statistically significant difference in the mean age (years), mean height (cm), mean weight (kg), mean BMI (kg/m2), and mean BSA (m2) of controls and cases (P > 0.05). This may be due to methodological limitations and smaller sample size.
Comparison of time domain analysis
[Table 3] shows the comparison of time domain analysis of controls and cases. Our study observed statistically highly significant differences in the mean RR of controls and cases (P < 0.01). It showed increase in the mean RR in controls as compared to cases. The values of the mean RR in the present study support the results obtained by Agelink et al. (2002), Jonathan et al. (2005), Kemp et al. (2012), Brunoni et al. (2013), Regan et al. (2014), Sgoifo et al. (2016), Paniccia et al. (2017), and Chen et al. (2017). The values of the mean RR in the present study are not in agreement with Sayar et al. (2002) and Gehi et al (2005). These contradictory findings may be due to variation in recording procedure, sample size, manual error, and sensitivity of machine.
No statistically significant differences were observed in the mean STD (SDNN) of controls and cases.
Statistically highly significant differences were observed in the mean HR of controls and cases (P > 0.05). Our study showed increase in mean HR of cases as compared to decrease in controls. The values of the mean HR in the present study are in agreement with those of Carney et al. (2001), Agelink et al. (2002), Jonathan et al. (2005), Kemp et al. (2012), Brunoni et al. (2013), Regan et al. (2014), Sgoifo et al. (2016), Paniccia et al. (2017), and Chen et al. (2017). Our study is not in agreement with the findings of Sayar et al. (2002) and Gehi et al. (2005). These contradictory findings may be due to variation in recording procedure, sample size, manual error, and sensitivity of machine.
No statistically significant differences were found in RMSSD, NN50, and pNN50 of controls and cases (P > 0.05).
Comparison of frequency domain analysis
[Table 4] shows the comparison of frequency domain analysis parameters in controls and cases. Our study observed no statistically significant differences in mean VLF peak, mean LF peak, and mean HF peak of controls and cases (P > 0.05). Our study observed highly statistically significant differences in mean VLF power (ms2) of controls and cases (P < 0.01), with the values higher among cases as compared to controls. The values of the VLF power in the present study are in agreement of the findings of Carney et al. (2001), Ha et al. (2015), and Chen et al. (2017). Statistically significant differences were observed in the mean LF power (ms2) of controls and cases (P < 0.05), with values higher among cases as compared to controls. The values of the LF power in the present study are in agreement with those of Carney et al. (2001), Wang et al. (2013), and Chen et al. (2017). No statistically significant differences were observed in mean HF power (ms2) of controls and cases. No statistically significant differences were observed in mean HF power %, mean VLF power %, and mean LF power % of controls and cases. The study observed decrease in the mean LF/HF ratio of controls and increase in cases. Statistically highly significant differences were observed in the mean LF/HF ratio of controls and cases (P < 0.01). The values of the mean LF/HF ratio in the present study are in agreement with the findings of Carney et al. (2001), Wang et al. (2013), and Chen et al. (2017). The values of the mean LF/HF ratio in the present study are not in agreement with those of Royster et al. (2012). These contradictory findings may be due to variation in recording procedure, sample size, manual error, and sensitivity of machine. LF/HF ratio is considered to mirror the sympathovagal balance or reflect the sympathetic modulations. No statistically significant (P > 0.05) differences were observed in the mean LF (nu) and mean HF (nu) of controls and cases.
Our study observed that the values of mean LF power (ms2) were very high in known cases of depression as compared to controls. In addition, the LF/HF ratio was very high in cases of depression as compared to controls. Both these findings demonstrate that the autonomic nervous system shifts in favor of sympathetic overactivity.
Alterations of autonomic nervous system functioning that promote vagal withdrawal are reflected in the reductions of HRV indices. A high sympathetic and a low cardiovagal activity in patients with major depression (MD) may contribute to the higher cardiac morbidity and mortality of MD patients. Depressed patients showed a higher HR and a significantly lower modulation of cardiovagal activity compared to controls. Sympathetic activity increases in depressed patients, and this enhancement is more obvious which expresses higher complexity of HRV time series. Parasympathetic nervous system of depressed patients is influenced by mental state, and they become dysfunctional under long-term depression.
Keeping in mind the prevalence of depression and its effects on the autonomic nervous system, it is suggested that HRV analysis may also be done as a part of routine screening procedures which are of vital importance in patients with depression.
| Conclusion|| |
The results of the present study demonstrate that the depressed patients influence HRV and change the autonomic balance in favor of an increased sympathetic tone. A high sympathetic and/or a low cardiovagal activity in patients with MD may contribute to the higher cardiac morbidity and mortality of MD patients. Parasympathetic nervous system of MD patients is influenced by mental state, and they become dysfunctional under long-term depression. Alterations of autonomic nervous system functioning that promote vagal withdrawal are reflected in the reductions of HRV indices.
Awareness should be created among health professionals involved directly or indirectly in the care of patients of depression and to inform the general public about the association of depression with the effects on HRV and autonomic nervous system. Education and treatment strategies have to be planned accordingly. Special care must be given to patients with depression before a significant change in HRV occurs. In addition, a close follow-up is necessary in patients with depression.
- Small sample size of the study
- Methodological limitations
- Sensitivity of the machine
- Possible manual errors.
The authors are grateful to all the participants without whose co-ordination and active participation, this work would not have been fruitful. We are also thankful to the clinical labs, for their help and support.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest (None).
| References|| |
Chalmers JA, Quintana DS, Abbott MJ, Kemp AH. Anxiety disorders are associated with reduced heart rate variability: A meta-analysis. Front Psychiatry 2014;5:80.
Carney RM, Freedland KE, Stein PK, Skala JA, Hoffman P, Jaffe AS. Change in heart rate and heart rate variability during treatment for depression in patients with coronary heart disease. Psychosom Med 2000;62:639-47.
Sathish B, Bhat R, Damodara M. Autonomic modulation in different phases of menstrual cycle in young and older Indian women. Int J Med Res 2013;1:12-6.
Lee CL, Chang WD. The effect of cigarette smoking on aerobic and anaerobic capacity and heart rate variability among female university students. Int J Womens Health 2013;5:667-79.
Tarvainen MP, Kubois NJ. HRV Users Guide. Ver. 2.1. Biosignal Analysis and Medical Imaging Group (BSAMIG), Department of Applied Physics, University of Eastern Finland, Kuopio, FINLAND. MATLAB. Copyright 1984-2012 The MathWorks, Inc. MATLAB is a registered trademark of The MathWorks, Inc. 2012. p. 8-12.
Gritti I, Defendi S, Mauri C, Banfi G, Duca P, Roi GS. Heart rate variability, standard of measurement, physiological interpretation and clinical use in mountain marathon runners during sleep and after acclimatization at 3480 m. J Behav Brain Sci 2013;3:26-48.
Doss SS, Anandhalakshmi S, Rekha K, Akhil AK. Effect of Smoking on heart rate variability in normal healthy volunteers. Asian J Pharm Clin Res 2016;9:230-34.
Kaya H, Suner A, Koroglu S, Akcay A, Turkbeyler IH and Koleoglu M. Heart rate variability in familial Mediterranean fever patients. Europ J Rheumatol 2014;1:58-61.
Karim N, Hasan JA, Syed SA. Heart rate variability – A review. J Basic Appl Sci 2011;7:71-7.
(Authors not listed) 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. Eur Heart J 1996;17:354-81.
Sayar K, Güleç H, Gökçe M, Ismail AK. Heart rate variability in depressed patients. Bull Clin Psychopharmacol 2002;12:130-33.
Licht CM, de Geus EJ, Zitman FG, Hoogendijk WJ, van Dyck R, Penninx BW. Association between major depressive disorder and heart rate variability in the Netherlands study of depression and anxiety (NSEDA). Arch Gen Psychiatry 2008;65:1358-67.
Janca A, Hiller W. ICD-10 checklist – A tool for clinicians use of the ICD-10 classification of mental and behavioral disorder. Compr Psychiatry 1996;37:180-7.
BMI Classification, Global Database on Body Mass Index. World Health Organization; 2006.
Du Bois D, Du Bois EF. A formula to estimate the approximate surface area if height and weight be known. Arch Intern Med. 1916;17:863-71.
Agelink MW, Boz C, Ullrich H, Andrich J. Relationship between major depression and heart rate variability. Clinical consequences and implications for antidepressive treatment. Psychiatry Res 2002;113:139-49.
Davidson J, Watkins L, Owens M, Krulewicz S, Connor K, Carpenter D, et al
. Effects of paroxetine and venlafaxine XR on heart rate variability in depression. J Clin Psychopharmacol 2005;25:480-4.
Kemp AH, Quintana DS, Felmingham KL, Matthews S, Jelinek HF. Depression, comorbid anxiety disorders, and heart rate variability in physically healthy, unmedicated patients: Implications for cardiovascular risk. PLoS One 2012;7:e30777.
Brunoni AR, Kemp AH, Dantas EM, Goulart AC, Nunes MA, Boggio PS, et al
. Heart rate variability is a trait marker of major depressive disorder: Evidence from the sertraline vs. electric current therapy to treat depression clinical study. Int J Neuropsychopharmacol 2013;16:1937-49.
O'Regan C, Kenny RA, Cronin H, Finucane C, Kearney PM. Antidepressants strongly influence the relationship between depression and heart rate variability: Findings from the Irish longitudinal study on ageing (TILDA). Psychol Med 2015;45:623-36.
Sgoifo A, Carnevali L, Alfonso Mde L, Amore M. Autonomic dysfunction and heart rate variability in depression. Stress 2015;18:343-52.
Paniccia M, Paniccia D, Thomas S, Taha T, Reed N. Clinical and non-clinical depression and anxiety in young people: A scoping review on heart rate variability. Auton Neurosci 2017;208:1-14.
Chen X, Yang R, Kuang D, Zhang L, Lv R, Huang X, et al
. Heart rate variability in patients with major depression disorder during a clinical autonomic test. Psychiatry Res 2017;256:207-11.
Gehi A, Mangano D, Pipkin S, Browner W, Whooley M. Depression and heart rate variability in patients with stable coronary heart disease. Arch Gen Psychiatry 2005;62:661-6.
Carney RM, Blumenthal JA, Stein PK, Watkins L, Catellier D, Berkman LF, et al
. Depression, heart rate variability, and acute myocardial infarction. Circulation 2001;104:2024-8.
Ha JH, Park S, Yoon D, Kim B. Short-term heart rate variability in older patients with newly diagnosed depression. Psychiatry Res 2015;226:484-8.
Wang Y, Zhao X, O'Neil A, Turner A, Liu X, Berk M. Altered cardiac autonomic nervous function in depression. BMC Psychiatry 2013;13:187.
Royster EB, Trimble LM, Cotsonis G, Schmotzer B, Manatunga A, Rushing NN, et al
. Changes in heart rate variability of depressed patients after electroconvulsive therapy. Cardiovasc Psychiatry Neurol 2012;2012:794043.
Askelrod S, Gordon D, Ubel FA, Shannon DC, Berger AC, Cohen RJ. Power spectrum analysis of heart rate fluctuation: A quantitative probe of beat-to-beat cardiovascular control. Science 1981;213:220-2.
[Table 1], [Table 2], [Table 3], [Table 4]