getNormSpectralUnits: Normalized Spectral Units

View source: R/CalculatePowerBand.R

getNormSpectralUnitsR Documentation

Normalized Spectral Units

Description

Calculates the spectrogram bands in normalized units

Usage

getNormSpectralUnits(HRVData,
  indexFreqAnalysis = length(HRVData$FreqAnalysis), VLFnormalization = T)

Arguments

HRVData

Data structure that stores the beats register and information related to it

indexFreqAnalysis

Reference to the data structure that contains the spectrogram analysis

VLFnormalization

Logical value. If TRUE (default), the function normalizes LF and HF power series by its sum. If FALSE, the function computes VLF, LF and HF power series by its sum.

Details

The default behaviour of this function computes the normalized power time series in the LF and HF bands following the Task Force recommendations:

normalized_LF = LF_power / (total_power - VLF_power - ULF_power)

normalized_HF = HF_power / (total_power - VLF_power -ULF-power)

If VLFnormalization is set to FALSE, the functions computes:

normalized_VLF = VLF_power / (total_power - ULF_power)

normalized_LF = LF_power / (total_power - ULF_power)

normalized_HF = HF_power / (total_power - ULF_power)

The resulting time series are returned in a list. Note that before using this function, the spectrogram should be computed with the CalculatePowerBand function.

Value

The getNormSpectralUnits returns a list storing the resulting normalized power-band series. Note that this list is not stored in the HRVData structure.

References

Camm, A. J., et al. "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 93.5 (1996): 1043-1065.

Examples

## Not run: 
# load some data...
data(HRVProcessedData)
hd = HRVProcessedData
# Perform some spectral analysis and normalize the results
hd = CreateFreqAnalysis(hd)
hd = CalculatePowerBand(hd,indexFreqAnalysis = 1,shift=30,size=60)
normUnits = getNormSpectralUnits(hd)
# plot the normalized time series
par(mfrow=c(2,1))
plot(normUnits$Time, normUnits$LF, xlab="Time", ylab="normalized LF",
     main="normalized LF",type="l")
plot(normUnits$Time, normUnits$HF, xlab="Time", ylab="normalized HF",
     main="normalized HF",type="l")
par(mfrow=c(1,1))


## End(Not run)

RHRV documentation built on Nov. 1, 2022, 1:05 a.m.