NeuroAnalyzer tutorials: Analyze EEG (5)

Load data:

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

[ Info: Loaded: EEG (24 × 308480 × 1; 1204.996 s)

Lateralization index is calculated as log(A / B), where A is average power at given frequency (default is 10 Hz, alpha) for the right hemisphere and B is average power at that frequency for the left hemisphere.

Lateralization index (default for alpha 10 Hz frequency):

lidx = lat_idx(eeg)

1.517379806560375

Lateralization index for the range 8.0-12.5 Hz:

lidx = lat_idx(eeg, frq=(8.0, 12.5))

1.4688999711396038

Calculate frontal alpha asymmetry index:

# left
ch_left = channel_cluster(eeg, cluster=:f1)
# or
ch_left = channel_pick(eeg, p=[:left, :frontal])
# right
ch_right = channel_cluster(eeg, cluster=:f2)
# or
ch_right = channel_pick(eeg, p=[:right, :frontal])

# alpha frequency
alpha = band_frq(eeg, band=:alpha)

# asymmetry index
band_asymmetry(eeg, frq_lim=alpha, ch1=labels(eeg)[ch_left], ch2=labels(eeg)[ch_right])

(ba = -2.6469560812055706,
 ba_norm = -0.8676462643674967,)