Forum: open-discussion


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]![]() |
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="") |