EEGlab_time_frequency_Delorme.pdfVIP

Time-Frequency analysis of biophysical time series ? Arnaud Delorme CERCO, CNRS, France SCCN, UCSD, La Jolla, USA Frequency analysis 1 second 4-7 Hz Theta 9-11 Hz Alpha 18-21 Hz Beta 30-60 Hz Gamma 0.5-2 Hz Delta ? synchronicity of cell excitation determines amplitude and rhythm of the EEG signal Spectral phase and amplitude Real Imaginary Real I m a g . Fk(f,t) Real I m a g . Real Imaginary Spectral phase and amplitude Fk(f,t) function [a,b] = dft (y) % DFT - The Discrete Fourier Transform % [a, b] = DFT (y) % a, b are the cosine and sine components n = length (y); t = 2*pi*(0:n-1)/n; f = 2.0 / n; for j = 0:n2 cs = cos (j * t); ss = sin (j * t); a(j+1) = f * (cs * y); b(j+1) = f * (ss * y); end % boundaries n2 = floor (n / 2); a(1) = 0.5 * a(1); a(n2+1) = 0.5 * a(n2+1); b(1) = 0.0; b(n2+1) = 0.0; Loop on frequency Cosine component Sine component Multiply with signal From /pwessel/courses/gg691/dft.m Average of squared absolute values Spectral power 0 Hz 10 Hz 20 Hz 30 Hz 40 Hz 50 Hz P o w e r ( d B ) Frequency (Hz) Average of squared amplitude Average Overlap 50% padding Plot data spectrum using EEGLAB ‘winsize’, 256 (change FFT window length) ‘nfft’, 256 (change FFT padding) ‘overlap’, 128 (change window overlap) Average Overlap 50% Padding Exercice - Read dataset “eeglab_data_epoch_ica.set” - Plot spectral decomposition 1- default (15% of data) 2- with 100% of data 3- with 100% and 50% overlap (‘overlap’) 4- reduce window size (‘winsize’) 5- reduce fft length (‘nfft’) Spectrogram or ERSP 0 ms 10 ms 20 ms 30 ms 40 ms 50 ms 60 ms 5 Hz 10 Hz 20 Hz 30 Hz Spectrogram or ERSP 0 ms 10 ms 20 ms 30 ms 40 ms 50 ms 60 ms 5 Hz 10 Hz 20 Hz 30 Hz 5 Hz 10 Hz 20 Hz 30 Hz 0 ms 10 ms 20 ms 30 ms 40 ms 50 ms 60 ms Average of squared values ∑ = = n k k tfFn tfERSP 1 2),(1),( Complex number Scaled to dB 10Log10(ERSP) Power spectrum and event-related spectral perturbation From Delorme and Makeig, J Neuro Methods, 2004 Absolute versus relative power Absolute = ERS Relativ