Plot HRV Power Spectral Density with Frequency Bands

Plots the RR-interval power spectrum with the VLF, LF and HF bands shaded and their integrated powers annotated. The PSD and band integration reuse the same internals as ecgHRVfreq, so the shaded band areas equal the ecgHRVfreq band powers. The band powers are attached as the "band_power" attribute.

Usage

plotHRVpsd(
  rr,
  method = c("welch", "ar", "lomb"),
  detrend = FALSE,
  channel = NULL,
  vlf_band = c(0.003, 0.04),
  lf_band = c(0.04, 0.15),
  hf_band = c(0.15, 0.4)
)

Arguments

rr: An rr data.frame or an object from ecgProcess.
method: Spectral method ("welch" default, "ar" or "lomb"), matching ecgHRVfreq.
detrend: Smoothness-priors detrending, matching ecgHRVfreq.
channel: Channel to plot (default: first).
vlf_band, lf_band, hf_band: Band edges (Hz).

Value

A ggplot object with a "band_power" attribute.

References

Task Force (1996). Heart rate variability. Circulation, 93(5), 1043-1065.

Examples

pe <- make_ecg(40000, sr = 500)
rr <- ecgRRintervals(pe, ecgDetectRpeaks(pe))
plotHRVpsd(rr)