Load EEG object:
= import_edf("test/eeg-test-edf.edf")
eeg = 1
channel = 1 epoch
Calculate PSD data using Welch periodogram:
= psd(eeg, norm=true) p, f
Draw scatter plot of PSD data:
scatter(f, p[channel, :], ms=1, legend=false)
Calculate and plot PSD slope:
= psd_slope(eeg, norm=true)
f, psd_slope, frq plot!(frq, lf[channel, :, epoch], title="slope = $(round(psd_slope[channel, epoch], digits=2))")
Calculate PSD slope of the alpha band:
= psd_slope(eeg, f=(8, 14), norm=true) f, psd_slope, frq