Skip to contents

Automatic Sleep Staging

Source: R/eeg-sleep.R
eegSleepStage.Rd

Classifies EEG epochs into sleep stages (Wake, N1, N2, N3, REM) using spectral power analysis based on simplified AASM criteria. Each epoch is scored by computing spectral band powers (delta, theta, alpha, sigma, beta) and applying rule-based classification.

Usage

eegSleepStage(
  x,
  method = c("spectral", "rule_based"),
  epoch_sec = 30,
  assay_name = NULL
)

Arguments

x

A PhysioExperiment object with EEG data.

method

Staging method: "spectral" (FFT power-based) or "rule_based" (synonym, same algorithm).

epoch_sec

Epoch duration in seconds (default: 30, AASM standard).

assay_name

Input assay name. If NULL, uses the default assay.

Value

A data.frame with columns:

epoch

Integer epoch number.

stage

Character sleep stage: "W", "N1", "N2", "N3", or "REM".

start_sample

Integer start sample of the epoch.

end_sample

Integer end sample of the epoch.

delta_power

Numeric delta band power (0.5-4 Hz).

theta_power

Numeric theta band power (4-8 Hz).

alpha_power

Numeric alpha band power (8-13 Hz).

sigma_power

Numeric sigma band power (12-16 Hz).

beta_power

Numeric beta band power (16-30 Hz).

The result is also stored in metadata(x)$sleep_stages.

References

Berry, R. B., et al. (2017). AASM Scoring Manual Updates for 2017. Journal of Clinical Sleep Medicine, 13(5), 665-666.

See also

eegSpindleDetect(), eegKcomplexDetect(), eegSlowWaveDetect(), eegSleepMetrics(), eegPlotHypnogram()

Examples

pe <- make_eeg_sleep(n_time = 90000, n_channels = 2, sr = 500)
stages <- eegSleepStage(pe, epoch_sec = 30)
head(stages)
#>   epoch stage start_sample end_sample  delta_power theta_power alpha_power
#> 1     1     W            1      15000 4.635160e+02    484.5215 842732.6567
#> 2     2   REM        15001      30000 4.186072e+02 239693.1601    704.3590
#> 3     3    N3        30001      45000 3.530500e+04   3501.0902   7650.3457
#> 4     4    N3        45001      60000 1.274289e+07    500.4129    594.3589
#> 5     5    N3        60001      75000 3.420374e+04   3550.2539   7842.9891
#> 6     6   REM        75001      90000 4.190141e+02  94328.4560    665.6874
#>   sigma_power beta_power
#> 1    503.6035   1668.842
#> 2    540.5154   1593.042
#> 3  13495.6811   1773.220
#> 4    430.2320   1684.233
#> 5  13784.6707   1694.592
#> 6    516.8163  35357.397