Dr. Paolo Castiglioni (1) comments on our paper “Decomposition of heart rate variability spectrum into a power-law function and a residual spectrum” (2). In that commentary, he pointed out that both the rg (= s · log(Frq) + Y) in Equation (1) in his commentary and the PSDrg have the units of PSD, or ms2/Hz, so that the spectral ratio rPSD in Equation (3) in his commentary is a dimensionless number. In the above argument, the “rg” denotes “regression,” “PSD” denotes “power spectral density,” and “rPSD” denotes “residual PSD.” We are grateful to Dr. Castiglioni for his careful examination of our equations. Unfortunately, the argument of Dr. Catiglioni was based on the misunderstanding caused by the omission of the units of the variables in our Equation (2). To avoid confusion and misunderstanding, we shall re-express our equations in the followings with the units of the variables placed in suitable positions in the equations.
To facilitate the decomposition of the PSD in the whole heart rate variability (HRV) spectrum, the power-law relation of HRV can be obtained by plotting against within the frequency range from >0 Hz to the Nyquist frequency. Both and are dimensionless. The 0 Hz point must be excluded because log(0) is not defined mathematically. The linear regression relation between and in the PSD of HRV can be expressed as
where the “log” denotes logarithm, the subscript “rg” denotes “regression,” and the “s” and “Y” are the “slope” and “Y-intercept” of linear regression between and , respectively.
Equation (1) can be re-written as
Equation (2) shows that the is a power-law function of . The difference between and is the logarithm of the residual power spectral density, , that cannot be accounted for by the :
Equations (2, 3) give immediately the following expression of rPSD:
Since the PSDrg has the unit of ms2/Hz, the rPSD also has the same unit of ms2/Hz as that of PSD, according to Equation (4). Thus, the rPSD is not dimensionless, as asserted by Dr. Castiglioni (1). Because the rPSD is obtained from the original PSD by removing the power law constituent in it, its frequency components are all significantly smaller than those of the original PSD. This is comprehensible.
Castiglioni also commented on our fitting of the regression line up to the Nyquist frequency (1). He asserted that the PSDrg cannot be considered the true “fractal component” of the spectrum, because the regression slope s is influenced by the oscillations in the low-frequency (LF) and high-frequency (HF) bands. Though the PSDrg can be influenced by the oscillations in the LF and HF bands, its overall behavior is self-similar with respect to frequency because the PSDrg is a power-law function of frequency, according to Equation (2). It is doubtful that the PSDrg cannot be interpreted as a fractal. If the PSDrg cannot be interpreted as a “fractal,” then call it a power-law function of frequency. It is not important whether the power-law function PSDrg obtained in our decomposition method can be regarded as a fractal or not.
The purpose of decomposing the PSD of HRV into a power-law function of frequency and a residual PSD is to examine more closely the heart rate oscillations in the LF and HF regions, rather than to examine the frequency region below 0.04 Hz, because LF and HF regions are mostly concerned in many clinical settings that are associated with autonomic dysfunction of the patients. If we adhere to the international guidelines on HRV to analyze the PSD within the frequency range of <0.04 Hz (3–7), then the heart rate oscillations in the LF and HF regions cannot be examined in more details by using the decomposition techniques introduced by us in our previous article (2).
The PSD of HRV obtained from either short-term or long-term (6, 7) recordings of heart periods can be decomposed into a power-law function and a residual spectrum by using the mathematical technique devised by our group. The point is that if we are going to examine the LF and HF components of HRV in more details by using the decomposition method, then fitting the regression line between and beyond the upper frequency limit (0.04 Hz) of the very low frequency (VLF) component is a necessity. Otherwise, the PSD within the frequency range of LF and HF will not be decomposed, and the LF and HF components of HRV cannot be examined in more details by using the decomposition method developed by our group. Thus, the decomposition of the PSD of HRV within the whole frequency range from >0 to the Nyquist frequency is the simplest way of performing the decomposition analysis.
Statements
Author contributions
JK: drafting of the manuscript. C-DK: conception, design and finalization of the work.
Funding
This work was supported by the grant 106-CCH-IRP-100 from the Changhua Christian Hospital, Changhua, Taiwan.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1.
CastiglioniP. Commentary: decomposition of heart rate variability spectrum into a power-law function and a residual spectrum. Front Cardiovasc Med. (2018) 5:94. 10.3389/fcvm.2018.00094
2.
KuoJKuoCD. Decomposition of heart rate variability spectrum into a power-law function and a residual spectrum. Front Cardiovasc Med. (2016) 3:16. 10.3389/fcvm.2016.00016
3.
KobayashiMMushaT. 1/f fluctuation of heartbeat period. IEEE Trans Biomed Eng. (1982) 29:456–7.
4.
OtsukaKNakajimaSYamanakaT. Vagal tone and its association with a new index of heart rate variability called 1/f fluctuations. J Ambul Monit. (1994) 7:213–18.
5.
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. Eur Heart J. (1996) 17:354–81.
6.
SassiRCeruttiSLombardiFMalikMHuikuriHVPengCKet al. Advances in heart rate variability signal analysis: joint position statement by the e-Cardiology ESC Working Group and the European Heart Rhythm Association co-endorsed by the Asia Pacific Heart Rhythm Society. Europace (2015) 17:1341–53. 10.1093/europace/euv015
7.
KazumaNNozakiMNakamuraEMatsuokaIOtaniT. Biological rhythm in 1/f fluctuations of heart rate in asthmatic children. Allergol Intern (2004) 53:265–9. 10.1111/j.1440-1592.2004.00343.x
Summary
Keywords
decomposition, heart rate variability, power-law function, residual spectrum, power spectral analysis, fractal
Citation
Kuo J and Kuo C-D (2019) Response: Commentary: Decomposition of Heart Rate Variability Spectrum into a Power-Law Function and a Residual Spectrum. Front. Cardiovasc. Med. 5:191. doi: 10.3389/fcvm.2018.00191
Received
31 July 2018
Accepted
17 December 2018
Published
17 January 2019
Volume
5 - 2018
Edited by
Shimon Rosenheck, Meir Medical Center, Israel
Reviewed by
Yael Yaniv, Technion Israel Institute of Technology, Israel
Updates
Copyright
© 2019 Kuo and Kuo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cheng-Deng Kuo cdkuo23@gmail.com
This article was submitted to Cardiac Rhythmology, a section of the journal Frontiers in Cardiovascular Medicine
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.