NeuroAnalyzer tutorials: Plot band powers

using NeuroAnalyzer
using Plots
eeg = load("files/eeg.hdf");

Plot band power in the δ, θ, α and β ranges

Collect frequency bounds:

bands = [:delta, :theta, :alpha, :beta]
frq = Vector{Tuple{Float64, Float64}}()
for idx in bands
    push!(frq, band_frq(eeg, band=idx))
end

Calculate power for channel 2, epoch 1:

bp = Vector{Float64}()
ch = 2
ep = 1
for idx in frq
    push!(bp, band_power(eeg, ch="all", frq_lim=idx)[ch, ep])
end

Plot PSD and add band ranges:

plot_psd(eeg,
         ch="Fp2",
         ep=1,
         db=false,
         frq_lim=(0, 30),
         mono=true)
vline!(frq,
       lc=:black,
       alpha=0.5,
       lw=0.5)

Plot band powers:

plot_bar(bp,
         xlabels=string.(bands),
         xlabel="",
         ylabel="Power [μV²/Hz]",
         title="Band powers\n[channel: $ch, epoch: $ep]")

Plot PSD relative to alpha band power

plot_psd(eeg,
         ch="Fp2",
         ep=ep,
         ref=:alpha,
         frq_lim=(0, 30),
         mono=true)
vline!(frq,
       lc=:black,
       alpha=0.5,
       lw=0.5)