Forum: help

RE: Use CalculateEnergyInPSDBands or CalculatePowerBand for LF and HF in 60s record [ Reply ]
By: Abraham Otero on 2023-03-12 18:26

[forum:49625]

We added this function after the publication of the book to facilitate studies where you do not want to analyze the temporal evolution of frequency bands. You can find the code here:

https://rdrr.io/github/constantino-garcia/RHRV_tests/src/R/StationaryFrequencyAnalysis.R

In principle, it should not be discarding 10% of the samples at the beginning or end of the record (unless they have been filtered by FilterNIHR).

RE: Use CalculateEnergyInPSDBands or CalculatePowerBand for LF and HF in 60s record [ Reply ] By: Robert Zieringer on 2023-03-10 15:23	[forum:49624]
Thanks Abraham! Regarding the taper parameter in the CalculatePSD function: I can't find much information about it in your book nor in R help functions. I noticed, that the results change, when I play around with it. Does it mean, that by default 10% at the beginning and 10% at the end of the recording are are discarded?

RE: Use CalculateEnergyInPSDBands or CalculatePowerBand for LF and HF in 60s record [ Reply ]
By: Abraham Otero on 2023-03-08 06:50

[forum:49621]

If you want to calculate the spectral power of a whole recording (and being recordings as short as 60 seconds it would not make too much sense to do otherwise) indeed CalculateEnergyInPSDBands is the best option. CalculatePowerBand calculates the energy per window, returning a time series with the energy value for each of the frequency bands within each window. This is useful to see changes in frequency along the record.

Regarding the units, yes, they would be ms^2.

RE: Use CalculateEnergyInPSDBands or CalculatePowerBand for LF and HF in 60s record [ Reply ] By: Robert Zieringer on 2023-03-06 16:01	[forum:49620]
Thanks Abraham for your quick reply. So I understand, that CalculateEnergyInPSDBands would be the correct command to use. Which unit does the output have though? is it ms^2?

RE: Use CalculateEnergyInPSDBands or CalculatePowerBand for LF and HF in 60s record [ Reply ] By: Abraham Otero on 2023-03-06 14:43	[forum:49619]
After filtering and interpolating the RR series do: HRVData=CreateFreqAnalysis(HRVData) HRVData=CalculatePSD(HRVData, "pgram") # get Energy in the default ULF, VLF and LF frequency bands. CalculateEnergyInPSDBands(HRVData)

Use CalculateEnergyInPSDBands or CalculatePowerBand for LF and HF in 60s record [ Reply ]
By: Robert Zieringer on 2023-03-04 10:51

[forum:49618]

I have several recordings of a length of 60s. I want to get the LF and the HF components of those recordings. However, I am not clear which command I should use: CalculateEnergyInPSDBands or CalculatePowerBand?

Also: if CalculatePowerBand is the way to go, what would make sense for the parameters shift and size, given that it is such a short recording?
The results appear to be quite different depending on which value I choose for size.

Grateful for you help