Spectral power of the electrodermal signal in the sympathetic band (0.045-0.25 Hz by default), following Posada-Quintero et al. (2016). The signal is optionally downsampled, its Welch power spectral density is estimated, and the band power is reported both in absolute units and normalised by the total spectral power.
Usage
edaSympatheticIndex(
x,
band = c(0.045, 0.25),
assay_name = NULL,
resample_hz = 2,
seg_sec = NULL
)Arguments
- x
A PhysioExperiment object with an electrodermal assay.
- band
Numeric
c(low, high)sympathetic band in Hz (defaultc(0.045, 0.25)).- assay_name
Input assay (default:
defaultAssay(x)).- resample_hz
Rate (Hz) the signal is downsampled to before spectral estimation (default 2;
NULLkeeps the native rate). Ignored if it is not below the native sampling rate.- seg_sec
Welch segment length in seconds (
NULL= one whole-signal Hann window).
Value
A data frame with one row per channel: channel, edasymp
(absolute band power), edasymp_norm (band power / total power, in
[0, 1], or NA for a constant / zero-power channel), total_power, and
the band_low/band_high used.
References
Posada-Quintero, H.F. et al. (2016). Power spectral density analysis of electrodermal activity for sympathetic function assessment. Ann Biomed Eng 44(10):3124-3135.
Examples
pe <- edaSimulate(n_time = 1200, sr = 4, scr_count = 6, seed = 1)
edaSympatheticIndex(pe)
#> channel edasymp edasymp_norm total_power band_low band_high
#> 1 1 0.006002253 0.1195644 0.05020099 0.045 0.25