Compute power spectral density of bands listed in the bands argument.

Source: R/eeg.R

`bands_psd` calculates power spectral densities estimates on bands. Bands are computed from spectrogram bands equal or greater than lower limit and inferior to the upper limit.

bands_psd(signal, sRate, bands, normalize = FALSE, method = "pwelch")

Arguments

signal

Numerical vector of the signal.

sRate

Signal sample rate in Hertz.

bands

A list of bands to compute with lower and upper limits in the form `list(c(0,4),c(4,8))``

normalize

A band to normalize (divide) by. Defaults to `c(0.5,40)`. Can be set up to FALSE for raw results. Defaults to FALSE.

method

Character. Method to use to compute power spectral density. "pwelch" or "psm". Defaults to "pwelch".

Value

A list of bands powers.

Examples

signal <- sin(seq(0,100,0.01))
bands_psd(bands = list(c(0,4),c(4,8)), signal = signal, sRate = 200)
#> [[1]]
#> [1] -169.1249
#> 
#> [[2]]
#> [1] -237.7298
#>