CalculatePowerBand {RHRV} | R Documentation |
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
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)
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 |
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
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)
M. Lado, A. Mendez, D. Olivieri, L. Rodriguez, X. Vila
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)
CalculateSpectrogram
for spectrogram calculation