EEGLABで脳波解析に挑戦④

のろのろですが、分析をすすめています。

本を読んでもわからないところは、ハズバンドにも相談しています。

ありがとう。

これまでの分析はこちら↓を見てね。

１．イベントファイルの作成

itoy82.hatenablog.com

２．脳波ファイルの整理

itoy82.hatenablog.com

３．EEGでの分析①②③

itoy82.hatenablog.com

前回までにエポックを区切るところまでやっているので、今日はエポックリジェクションから始めます！

脳波解析入門を読みながら分析を進めているのですが、なぜか私のEEGLABでは、エポックリジェクションがうまくいかず、仕方がないのでコードを書くことにしました…。

[商品価格に関しましては、リンクが作成された時点と現時点で情報が変更されている場合がございます。]

脳波解析入門　Windows10対応版 EEGLABとSPMを使いこなす [ 開　一夫 ]
価格：4180円（税込、送料無料) (2022/11/8時点)

１．エポックリジェクション

ハズバンドに見てもらいながら作ったコードがこちら。

i = 0;
rej_filename = ['sub10',num2str(i,'%02d'),'_rej.set'];
EEG = pop_eegthresh(EEG,1,[1:3] ,-100,100,-0.2,0.8,2,1);
EEG = eeg_rejsuperpose( EEG, 1, 1, 1, 1, 1, 1, 1, 1);
EEG = pop_rejepoch( EEG, EEG.reject.rejglobal,0);
EEG = pop_saveset( EEG, 'filename',rej_filename,'filepath','C:\Users\itoyu_000\OneDrive\Documents\MATLAB');
EEG = eeg_checkset( EEG );
eeglab redraw

これを、MATLABのコマンドウィンドウに貼り付けてエンターするだけです。

ここでようやく、最初からコードが書ければ、めちゃくちゃ分析がらくちんなんじゃないか？と、コードを最初から書いちゃう人の気持ちを想像することができました笑

最初の i = 0;　のところは、参加者の番号を入れるようにしました。

今回は、予備実験のデータを使ってやってみているので、参加者番号00としました。

2行目は、ファイルの名前の付け方です。sub1000‗rej.setというファイル名で保存されるように設定しています。1行目のi=0をi=1に変更すると、ファイル名sub1001_rej.setとなります。

3行目は、どういうルールでリジェクトするかを書いています。EEGの中に、今EEGLABで選択しているファイルが入っていて、次の1はローデータということ(コンポーネントデータであれば0になります）、その次はリジェクトする周波数なので、-100～100として、範囲は、私はｰ200ms～800msなので秒に直して-0.2，0.8としています。その次の2は、データをそのまま上書きして、最後の1ですぐにリジェクトするという感じです。

4行目で、これらのデータを上書きして、

6行目で、ファイル名を付けて、保存場所を設定しておしまいという感じです。

7行目のeeg_checksetっていうのは、よくわからないのだけど、定期的に出てくるので、EEGの構造に問題がないかチェックしてるみたいです！

eeglab redrawは、eeglabのGUIに現在の状態を反映させるために再読み込みするコードなようです。

MATLABでなんかやっていて、わからなくなったら、わからないコードのところにカーソルを持って行ってF1を押すと、ヘルプが出てきてくれるので、安心です🐍

２．刺激の選択

ここから、それぞれの刺激を選んで、刺激ごとの加算平均波形を出力することが出来ます！