numpy, xdf, matplotlib, re, resampy, mne, braindecode
- Implement quick data analysis pipeline:
- Load all runs and have a mechanism to keep them separate (e.g. keep them in a list)
- Create a downsampled (250 Hz) copy of the data using resampy. USe this copy for all next steps but the spectrum plots
- Make a comb filter to remove 50 Hz and 90 Hz and harmonics to make visual inspection possible
- Remove DC offset (subtract mean), keeping offset values somewhere
- Calls mne raw data viewer on the concatenated runs to perform visual inspection
- Calculate and topographically plot power spectrum for each electrode and each file
- Calculate and topographically plot variance of each electrode and each file
- Epoch the data based on user specified markers and time-windows
- Average the epochs for specified markers individually and plot all the marker specific lines with deviation shades for each electrode in a topographical arangement with axes and legend
- Calculate and topographically plot the average power spectrum for each electrode and each file with one line per marker
- Calculate time-frequency spectrograms for each epoch
- Calculate baseline (-500 ms - 0 ms before marker, median across time & trials in two steps, all markers pooled)
- Divide spectrograms by baseline and take natural log
- Plot median spectrogram topographically for each marker