CalculatePowerBand {RHRV}R Documentation

Calculates power per band

Description

Calculates the spectrogram of the heart rate signal and the power of the spectrogram of the heart rate signal at the HF, LF, VLF bands

Usage

CalculatePowerBand(HRVData, indexFreqAnalysis = -1, size, shift, sizesp = 1024, scale = "linear", ULFmin = 0, ULFmax = 0.03, VLFmin = 0.03, VLFmax = 0.05, LFmin = 0.05, LFmax = 0.15, HFmin = 0.15, HFmax = 0.4, verbose = FALSE)

Arguments

HRVData Data structure that stores the beats register and information related to it
indexFreqAnalysis Reference to the data structure that will contain the variability analysis
size Size of window for calculating spectrogram (seconds)
shift Displacement of window for calculating spectrogram (seconds)
sizesp Seconds for calculating spectrogram (zero padding)
scale Scale used to plot spectrogram, linear or logarithmic
ULFmin Lower limit ULF band
ULFmax Upper limit ULF band
VLFmin Lower limit VLF band
VLFmax Upper limit VLF band
LFmin Lower limit LF band
LFmax Upper limit LF band
HFmin Lower limit HF band
HFmax Upper limit HF band
verbose Boolean argument that allows to specify if the function returns additional information

Value

Returns HRVData, the structure that contains beat positions register, associated heart rate instantaneous values, filtered heart rate signal equally spaced, and the analysis structure including the power of the spectrogram at different bands of the heart rate signal

Note

An example including all the necessary steps to obtain the power bands of a wfdb register is giving below:

##Reading a wfdb register and storing into a data structure:
md = CreateHRVData(verbose = TRUE)
md = LoadBeatWFDB(md, RecordName = "register_name", RecordPath = "register_path", verbose = TRUE)
##Calculating heart rate signal:
md = BuildNIHR(md, verbose = TRUE)
##Filtering heart rate signal:
md = FilterNIHR(md, verbose = TRUE)
##Interpolating heart rate signal:
md = InterpolateNIHR(md, verbose = TRUE)
##Calculating spectrogram and power per band:
md = CreateFreqAnalysis(md, verbose = TRUE)
md = CalculatePowerBand(md, indexFreqAnalysis = 1, size = 120, shift = 10, sizesp = 1024, verbose = TRUE)

Author(s)

M. Lado, A. Mendez, D. Olivieri, L. Rodriguez, X. Vila

References

L. Rodriguez-Linares, X. Vila, A. Mendez, M. Lado, D. Olivieri, "RHRV: An R-based software package for heart rate variability analysis of ECG recordings," 3rd Iberian Conference in Systems and Information Technologies (CISTI 2008), Proceedings I, 565-573, ISBN: 978-84-612-4476-8 (2008)

See Also

CalculateSpectrogram for spectrogram calculation


[Package RHRV version 2.1 Index]