SCM

Forum: open-discussion

Monitor Forum | Start New Thread Start New Thread
RE: Power calculation limits and plotting [ Reply ]
By: Abraham Otero on 2023-06-16 06:51
[forum:49688]
By default in the function clips the frequency calculation in 0.5. Add the to parameter to modify this:

hrv.data=CalculatePSD(hrv.data,indexFreqAnalysis = 1, ULFmin=0, ULFmax=0.003, VLFmin = 0.0033, VLFmax = 0.04,
LFmin = 0.04, LFmax = 0.15, HFmin = 0.15, HFmax = 0.8,
main="Lomb Periodogram Pac 30 14? day",
method = "lomb",doPlot=T,alpha=0.05, to =0.9)

Power calculation limits and plotting [ Reply ]
By: Tomas FariƱa on 2023-06-15 17:22
[forum:49687]

IMG_5180.HEIC (5) downloads
I want to calculate the HF upper limit up to 0.6 Hz.

I use the CalculatePSD function but when I plot it, I get the following warning message:
In PlotPSD (hav.data, indexFreqAnalysis=1, ULFmin=NULL,ULFmax=NULL, : Plotting larger frequency range than computed. Increase the "to" parameter in CalculatePSD().
I get a plot. The x-axis label has 0.6 Hz but there is nothing in the graphic beyond 0.4 Hz.

I have the same values of HF using CalculateEnergyInPSDBands function. So I suppose RHRV only calculate HF up to 0.4 Hz.

I try to find the "to" parameter in help sections but can't find a solution.

Any clue to change this setting??

My code is:
hrv.data=CalculatePSD(hrv.data,indexFreqAnalysis = 1, ULFmin=0, ULFmax=0.003, VLFmin = 0.0033, VLFmax = 0.04,
LFmin = 0.04, LFmax = 0.15, HFmin = 0.15, HFmax = 0.6,
main="Lomb Periodogram Pac 30 14? day",method = "lomb",normalize="press",doPlot=T,alpha=0.05)

PlotPSD(hrv.data,indexFreqAnalysis=1,
ULFmin=NULL,ULFmax=NULL,VLFmin = 0.0033, VLFmax = 0.04,
LFmin = 0.04, LFmax = 0.15, HFmin = 0.15, HFmax = 0.6,
main="PSD estimated using Lomb Periodogram,log="")

Thanks to:
Vienna University of Economics and Business Powered By FusionForge