Response: Commentary: Decomposition of Heart Rate Variability Spectrum into a Power-Law Function and a Residual Spectrum.

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.

Author Contributions

JK: drafting of the manuscript. C-DK: conception, design and finalization of the work.

Conflict of Interest Statement

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.