Phase-amplitude coupling (cross-frequency coupling) for EEG

Measures how the phase of a slow oscillation modulates the amplitude of a fast oscillation — within a single channel (amp_channel = NULL) or across two channels. Delegates the estimator to PhysioCrossModal::phaseAmplitudeCoupling().

Usage

eegPAC(
  pe,
  phase_channel = 1L,
  amp_channel = NULL,
  phase_band = c(4, 8),
  amp_band = c(30, 80),
  method = c("tort", "canolty", "ozkurt", "plv"),
  n_bins = 18L,
  assay_name = NULL
)

Arguments

pe: A PhysioExperiment (continuous data; epoched data is concatenated).
phase_channel: Channel (index or label) supplying the modulating phase.
amp_channel: Channel supplying the modulated amplitude; NULL (default) uses phase_channel (within-channel PAC).
phase_band: Phase frequency band in Hz (default theta c(4, 8)).
amp_band: Amplitude frequency band in Hz (default gamma c(30, 80)).
method: "tort" (modulation index), "canolty" (mean vector length), "ozkurt", or "plv".
n_bins: Number of phase bins for the Tort modulation index (default 18).
assay_name: Assay to use (default: the object's default assay).

Value

The PAC result from PhysioCrossModal::phaseAmplitudeCoupling() (a list with pac and, for "tort", the phase-amplitude distribution).

References

Tort et al. 2010; Canolty et al. 2006.

Examples

if (FALSE) { # \dontrun{
pe <- make_eeg(n_time = 5000, n_channels = 4, sr = 250)
eegPAC(pe, phase_channel = 1, phase_band = c(4, 8), amp_band = c(30, 80))
} # }