DPSVOX.NKのブログ

dipfit -> could not locate MEX for ...

Fri, 13 Aug 2010 15:40:21 +0900

Matlab versionは2007a（7.4.0）

eeglab7_1_6_18b
eeglab8_0_3_4b
で検証

OSはWinXP 64bit

タイトルにあるようにMEXファイルがみつからない的なことを言われます．
この問題は極単純ですがいろんなところの問題がこんがらがっているのでややこしかったです．

詳しく書いておいていつか何かの足しになれば．

まずEEGLABでdipfitをやろうとすると出ました．
以前に
Mex File のコンパイル
と題してWindows 64bit環境でdipfitをやろうとするときの注意点を書きました．

これに付随してまた新たに問題が出たのが今回の報告．
EEGLAB bugzillaでも問題になっているので，
以下の記述中の言葉で適当に検索すれば出てきます．

EEGLAB error in function ptriproj()@ line 29
ERROR using ==> ptriproj at 29 could not locate MEX for ptriproj

ここで問題になっているMEXファイルはfieldtripの関数なので，
ここのページ
が紹介されたりして，以前の記事のようにコンパイルのやり方が示されています．
つまりコンパイルされた適切なファイルがないと，
あらかじめ用意されたmファイルを読みに行く形になっていて
そこに記述されているエラーメッセージが返されるということです．

しかし今回はコンパイル作業を完璧にこなしてもエラーが出続けました．
前回もコンパイルだけではダメで，
コンパイルした.mexw64ファイル群をより階層が上のフォルダに入れて，
.mexw32との競合を避ける必要がありました．

今回もそれだろうと，
・workフォルダ
・マイドキュメントのMATLABフォルダ
・そのときのカレントフォルダ
とさまざま上位階層と思われるところに入れてみましたがダメ．

何か別のパスと競合しているんじゃないかと思って，
MATLAB上のパスを一回全部消して，設定しなおしてみてもダメ．

古いバージョンで使っていたfieldtrip-liteを別から持ってきて，
そっちを読みにいくようにしてもダメ．

とお手上げな感じでした．

バタバタ色々やっていてやっと動いたのが，非常に基本的な問題でした．eeglabのパス設定です．普段私はeeglabを立ち上げていない状態でも
eeglabの関数を利用したりするので基本「サブフォルダも追加」でeeglabのフォルダ以下全てのフォルダをMATLABに登録しています．それがまずいみたいですね．
eeglabは「eeglab」とタイプして立ち上げを行ったとき，スタートアップで自分の使うフォルダを自分で設定しなおします．つまりeeglabを使うには，eeglabの最上位のフォルダのみをMATLABにパス登録しておけばきちんと動くようになっているのです．逆に言うとそうやって作られているので「サブフォルダも追加」で全部登録されてしまうと困る仕様になっているということです．今回のMEXが見つからないバグも，この二重登録が混乱を生じて起こったようで，パスを最上位のみにしたら直りました．

直接どこが競合したか突き止めていませんが，とにかくパス全部登録はダメらしい．特にMEXを使って処理するもので，パス階層の順位で読むファイルが変わるような時には注意する必要がありそうです．関数を使いたいときはその都度パスを通すか，逆にdipfitをやるときなどに限ってパスを外すかなど，で対処しなくてはならない，ということみたいです．

差分トポなど③応用編

Tue, 01 Jun 2010 14:12:28 +0900

差分トポなど②改造編の続きです

std_chantopo(...

はmファイルですから，

std_chantopo

とコマンドウィンドウに打ってやれば，
処理をしてくれるはずです．

しかしそれを適切に動かすには，

std_chantopo(erpdata, 'condnames', STUDY.condition, 'plottopo', fastif(length(allinds)==1, 'off', 'on'), ...
'datatype', 'spec', 'plotmode', opt.plotmode, 'groupnames',STUDY.group, 'unitx', '\muV', ...
'groupstats', pgroup, 'condstats', pcond, 'interstats', pinter, ...
'chanlocs', locs, 'plotsubjects', opt.plotsubjects, 'topovals',titlestr, plotcurveopt{:});

の括弧内にあるオプションを全部理解して，
適切なデータを引っ張ってこないといけません．

これはめんどい

そこで②で作ったmatファイルです．
先ほどの状態でstd_chantopoを動かすのに
必要な変数を全部保存しておいたんです．

これをMATLABワークスペース上にオープンしてやれば，
まさに

std_chantopo(erpdata, 'condnames', condition, 'plottopo', fastif(length(allinds)==1, 'off', 'on'), ...
'datatype', 'spec', 'plotmode', opt.plotmode, 'groupnames', group, 'unitx', '\muV', ...
'groupstats', pgroup, 'condstats', pcond, 'interstats', pinter, ...
'chanlocs', locs, 'plotsubjects', opt.plotsubjects, 'topovals',titlestr, plotcurveopt{:});

をそのままベタ貼りするだけで，
トポの表示ができてしまいます．

（ここでは元のプログラムを少し改変していますが，
conditionとgroupの前のSTUDY.を抜き取っただけです．
状況によってはSTUDY構造体が非常に重くなっている可能性があるので，
上記の対処を取りました．）

でもそれでは一度行ったトポ表示を繰り返すだけ．
今回は差分を行いたいわけですから，
表示するトポを各条件のモノではなく，
条件差のデータに変えたいということになります．

一目瞭然でわかると思いますが，

std_chantopo(...

のあとの一つ目の引数に

erpdata

というものがあります．

まさにこれがデータです．
中身をみるとセル扱いですが，
条件分のデータが格納されています．
3条件であれば，

< 3 x 1 cell >

です．

そのままでは使えないので

a=erpdata{1};

のように抽出してみると，

1 x 64 x 13

のように，
1，電極，人数
の3次元データになっています．

あとは色々いじって差分にするなり，
なんなり（といっても差分意外にないか）できます．

3条件なら以下のスクリプトを通すと，
条件1-条件2，条件2-条件3，条件3-条件1
の三つを出すようにできます．

for i=1:length(condition)
a(i).data=erpdata{i};
end
erpdata2=erpdata;
for i=1:length(condition)
if i==length(condition), erpdata{i}=a(i).data-a(1).data;
else erpdata{i}=a(i).data-a(i+1).data;
end
end

erpdata2として元データを残しているので，
両方出して確認してみましょう．

上記を応用すればERP波形を書く，

std_plotcurve(alltimes, erpdata, 'condnames', STUDY.condition, 'plottopo', fastif(length(allinds)==1, 'off', 'on'), ...
'datatype', 'spec', 'plotmode', opt.plotmode, 'groupnames', STUDY.group, 'unitx', 'ms', ...
'groupstats', pgroup, 'condstats', pcond, 'interstats', pinter, ...
'chanlocs', locs, 'plotsubjects', opt.plotsubjects, plotcurveopt{:});

も同じように変数書き出し法で，色々できます．

差分トポなど②改造編

Tue, 01 Jun 2010 13:24:12 +0900

差分トポなど①準備編の続きです．

改造編と称していますが，
これはEEGLAB既存のmファイルを書き換えるからです．
基本的に悪さをしない書き換えをしていますが，
嫌な人はSTUDYファイルを作る際に，
差分したデータセットを用意して，
それをつかって①をやる，という感じにすればよいと思います．

では本題．

「①準備編」で紹介したことができている前提です．
先ほどのParam設定の中に

[Plot scalp map at latency [ms]]

という設定がありました．

さきほどはスルーしていましたが，
ここに数値を入れるとトポグラフィー描画モードになります．

画像は400-600msの平均電位のトポを出したい場合．
単一の数値でも機能します．

ただ①で紹介した[Sel. all]をやってからでないと，
Plot ERPsを押してもエラーが出ます．
プログラムで電極が5個以上選択されていないと
出してくれないようになっているからです．

しかしここで出せるのは各条件ごとのトポ．
差分など好きなデータでトポを出したい場合は
いくつかのmファイルを書き換えます．

ポイントになるのは

std_erpplot
std_chantopo
std_readdata

の三つです．

このうち今回はできるだけ書き込みがバグにならないように，
std_erpplotのみを書き換えて，
後は手打ちで算出するという方法にします．

which std_erpplot

とかでどこにあるかを探してeditorで見てみます．

長いので細かくは見ませんが，

約330-340行目くらいにある以下の記述を変更します．

% compute statistics and plot
% ---------------------------
[pcond pgroup pinter] = std_stat(erpdata, 'groupstats', opt.groupstats, 'condstats', opt.condstats, ...
'statistics', opt.statistics, 'naccu', opt.naccu, 'threshold', opt.threshold);
locs = eeg_mergelocs(ALLEEG.chanlocs);
locs = locs(std_chaninds(STUDY, opt.channels));

if ~isempty(opt.topotime) & ~isnan(opt.topotime)
std_chantopo(erpdata, 'condnames', STUDY.condition, 'plottopo', fastif(length(allinds)==1, 'off', 'on'), ...
'datatype', 'spec', 'plotmode', opt.plotmode, 'groupnames',STUDY.group, 'unitx', '\muV', ...
'groupstats', pgroup, 'condstats', pcond, 'interstats', pinter, ...
'chanlocs', locs, 'plotsubjects', opt.plotsubjects, 'topovals',titlestr, plotcurveopt{:});
else
std_plotcurve(alltimes, erpdata, 'condnames', STUDY.condition, 'plottopo', fastif(length(allinds)==1, 'off', 'on'), ...
'datatype', 'spec', 'plotmode', opt.plotmode, 'groupnames', STUDY.group, 'unitx', 'ms', ...
'groupstats', pgroup, 'condstats', pcond, 'interstats', pinter, ...
'chanlocs', locs, 'plotsubjects', opt.plotsubjects, plotcurveopt{:});
end;

ここはデータに基づいて統計解析とデータのプロットをするところです．
注目するのは

std_chantopo(erpdata, 'condnames', STUDY.condition, 'plottopo', fastif(length(allinds)==1, 'off', 'on'), ...
'datatype', 'spec', 'plotmode', opt.plotmode, 'groupnames',STUDY.group, 'unitx', '\muV', ...
'groupstats', pgroup, 'condstats', pcond, 'interstats', pinter, ...
'chanlocs', locs, 'plotsubjects', opt.plotsubjects, 'topovals',titlestr, plotcurveopt{:});

です．

ここで

std_chantopo

はトポを表示するmファイルなので，
括弧内で指定したパラメータを引数として持って，
表示処理に向かう感じです．
そこでここで必要になるパラメータをmatファイルに保存してしまいます．

そのためのスクリプトが以下，

%% added by Noriaki Kanayama
condition=STUDY.condition;group=STUDY.group;
save('std_topo_set.mat','erpdata','condition', 'group', 'allinds', 'opt', 'pgroup', 'pcond', 'pinter', 'locs', 'titlestr', 'plotcurveopt','-mat');
%% added by Noriaki Kanayama

これを

std_chantopo(...

の直前に貼り付けます．
（もちろん%%がついている行はコメントなので任意で．
でもmファイル改変なので，後で検索しやすいようにしておくとよいです）

すると
std_topo_set.mat
というマットファイルを「トポを表示する」たびに
カレントフォルダに算出してくれます．

これで改造は終了です．

プログラムの途中で変数をファイルに保存しているだけなので，
基本的に他のどんな作業にも影響を与えないはずです．

このmatファイルを使うと，
後々楽に色々できます．

では続きは別で

差分トポなど①準備編

Tue, 01 Jun 2010 13:12:55 +0900

今回はSTUDYファイルを作った状態が前提です．

EEGLAB versionは7.1.3.11b.

STUDYファイルは何かと便利なので作っておきましょう．
これまではICAクラスタリングで使う印象でしたが，
Channel Baseの解析にも使えます．

今回は差分トポを出す方法．

差分波形はEEGLABでも実装されていますが，
差分波形をトポ表示したり，色々自由にやるのは難しいです．

多分開発者は元々のデータの情報量を削減するような
ERPの差分表示の方法があまり好きじゃないので（興味がない？）
きちんと実装していないのかなぁと思います．

個人的にもERPレベルで差分してしまうと，
（+10μV）-（+8μV）
と
（+1μV）-（-1μV）
の結果が同じになってしまうところとか，
ちょっと怖いなぁと思いますが．

取りあえず準備編

まずはメニューから

[Study]-[Precompute channel measures]

でERPを算出しておきます．

Grand Averageを出すだけだと思いますが，
カレントフォルダに

.daterp

という妙な拡張子のファイルがデータセット分できます．
これがデータです．

ここまでが下準備．

そして

[Study]-[Plot channel measures]

でEdit/Plot画面を出します．
（ちなみにこのウィンドウは特殊な構造で，
これが出ている間はデバック作業や，
コマンドの打ち込みがストップされます）

Plot ERPsの横のParamを押します．
すると設定画面が開きます．
（ここではERPしかPrecomputeしていないので，
その他のボタンはアクティブになっていません）

とかいう感じで適当に設定して，OKを押し，
先ほどの画面に戻ったら

Plot ERPsを押すとERPを表示してくれます．
この時真ん中上あたりの[Sel. all]ボタンを押すと，
電極が全て選択され，
その状態でPlotすれば全電極のデータを並べて表示してくれます．

ではここまでが準備編．

改造編に続く．

STUDYを作る際のエラーメモ

Tue, 23 Mar 2010 19:53:13 +0900

適当なメモなのであとで対処してちゃんと書いておきます

①pop_savestudy line93, 97でエラー？

②conditionが数値で表されてた場合,
std_specのline268で

boundaries = strmatch('boundary', lower({ EEG.event.type }));

赤字の部分が文字列じゃないと使えない関数みたいでエラーが起こる．

boundaries = strmatch('boundary', lower({ char(EEG.event.type) }));

のように無理矢理文字列にしてやれば通過できます．
マッチングして'boundary'って書いてあるか見るだけなので,
この変更をしても問題ないと思います．
（boundaryと一致するなら必ず文字列になっているはずだし,
一致しないならすべて数字なので,charにしたところで比較結果は変わらない．
単純にエラー回避）

コンタープロット修正版

Fri, 19 Mar 2010 12:46:41 +0900

完全内輪向け．

EEGLAB versionは7.1.3.11b.

前回のコンタープロットがまさに処理していたデータに限定的だったので，
少し汎用性を持たせて見ました．
ただどのデータをプロットするのか指定するだけですが…

以前のもの

contourf(ALLTFA.timec,ALLTFA.freqs,ALLTFA.ERSP.avdata, ...
[min(min(ALLTFA.ERSP.avdata)):0.1:max(max(ALLTFA.ERSP.avdata))])

では単極，単一条件で取りあえず出してみた時にしか機能しませんが，
普通多くの電極，いくつかの条件でデータを算出するはずですので，

ALLTFA.ERSP.avdata

は

ALLTFA.ERSP(電極数,条件数).avdata

になるはずです．

よって以下のようにしました．

%------------(Please set your request for plotting)--------------
ElecLabel='';
ElecNum=0;
ConditionNumber=0;

%------------(The program follows)--------------
if isempty(ElecLabel)
if ElecNum=0, fprintf('please set the Elec Label or number.\n');
elseif ~isempty(ElecNum), fprintf('please set the Elec Label or number.\n');
else elec=ElecNum;
end
else
for i=1:size(ALLTFA.channel,2)
if strcmp(ALLTFA.channel(1,i),ElecLabel)==1
ElecNum=i;break
end
if i==size(ALLTFA.channel,2)
error('please check whether ElecLabel is correct, or use the number')
end
end
end

cond=ConditionNumber;

contourf(ALLTFA.timec,ALLTFA.freqs,ALLTFA.ERSP(elec,cond).avdata, ...
[min(min(ALLTFA.ERSP.avdata)):0.1:max(max(ALLTFA.ERSP.avdata))])

で作ったmatファイルのデータをコンタープロット表示できます．

Conditionもラベルで指定できるようにしたら,
プロットに条件名も表示できていいですが,
それは元のTFA_EEGLAB_revisitedを少し改変する必要があるのでまた今度．
とりあえずどんな感じか見てみるだけってことで．

少しずつ汎用性を高めていきます．

ICAプロットのパワースペクトラム

Tue, 02 Mar 2010 12:19:12 +0900

完全内輪向け.
EEGLAB versionは7.1.3.11b.

ICAの

Reject components by map

でパワースペクトラムの表示周波数を変えたい場合.

pop_selectcomps.m

の

Line251にある

command = sprintf('pop_prop( %s, 0, %d, %3.15f, { ''freqrange'', [1 50] });', inputname(1), ri, button);

の青部分で決定されているのでここに数値を直接打ち込む.

東日本でノッチかけたり，
早めに高周波フィルタかけた場合など，
フィルタ部分に引っ張られて肝心の周波数帯域が
見れない時などに.

使ったらそのつど戻すか，
変数組み込んで外から制御できるように
しておくことをお勧めします.

STUDYファイル③

Fri, 26 Feb 2010 20:41:51 +0900

ひたすらプログラムの解説です

currdir='E:\FaceEmoCh\FECStudy';

はカレントディレクトリを打ち込んで設定しておきます.
間違えないように初めに打ち込んでおいて，後は触らないようにします.

[ALLEEG EEG CURRENTSET ALLCOM] = eeglab;
pop_editoptions( 'option_storedisk', 1, 'option_savetwofiles', 1, 'option_saveica', 0, 'option_single', 1, 'option_memmapdata', 0, 'option_computeica', 1, 'option_scaleicarms', 1, 'option_rememberfolder', 1);
ALLEEG = eeg_checkset(ALLEEG); EEG = eeg_retrieve(ALLEEG, CURRENTSET);

①EEGLABの起動
②EEGLABのオプションから一度に一つだけデータをメモリ上に読み込む設定にする
これをしないと，大量にロードしたときその分だけメモリが必要になります.
30人，8条件，電極64極とかだと大変なことになります.

cd(currdir);

カレントディレクトリを設定しています.
なんでここに書いたんだろう.
別にもっと前でもいいはず.

allfiles=dir('*.set');
選択したディレクトリに入っているファイルをすべて参照する.
allfilesという構造体に，
name
date
bytes
isdir
datenum
というフィールドが作成され，
情報が格納されます.
これを使っていろいろできるわけです.

for n=1:length(allfiles) %すべての条件にあうファイルをロードする

プログラム中のコメントアウトにあるように，
すべてのファイルに対する処理であることを示しています.
for文開始で，次のendまで.
nというインデックスを，ファイルの数（length(allfiles)）だけ
一つずつ増加させてすべてのファイルに処理をするようになってます.

loadname=allfiles(n).name; %ファイルの名前を取得してloadnameに格納

loadnameを使っていろいろ情報を引き出します.
なので名前がきちんとついていないとダメ.

dataname=loadname(1:42); %ファイル名の前42文字とする
EEG = pop_loadset( 'filename', loadname, 'filepath', pwd);

loadnameから拡張子を外してdatanameにしています.
拡張しぬきの文字数がここでは42ということです.
EEG = pop_loadset(...
は非常に便利でEEGLABを立ち上げなくても機能します.
これでデータをロードします.
pwdでは現在選択中のディレクトリということです.

EEG.filepath=pwd;
ExpName=loadname(1:9); %ファイル名の1-9文字で実験名を命名する
EEG.subject=loadname(11:15); %ファイルの前11-15文字で参加者名を命名する
EEG.condition=loadname(37:42);%ファイルの後ろから37-42文字で条件名を命名する

ここでloadnameから情報を引き出しています.
だからファイル名をきちんと書いておく必要があるんですね.
バラバラだとここでつまづきます.
というかまず一括管理は無理ですね.

[ALLEEG, EEG, CURRENTSET] = eeg_store( ALLEEG, EEG, 0 );
end

というわけでここまででfor文が終了.
データを読み上げて，対応する構造体に格納し，
endまできたら戻って次のデータ，
ということを全部のデータにやると，この次に進みます.

ALLEEG = eeg_checkset(ALLEEG); EEG = eeg_retrieve(ALLEEG, CURRENTSET);
[STUDY, ALLEEG] = std_checkset(STUDY, ALLEEG);
eeglab redraw

ここはほとんど大人の事情.
EEGLAB内でデータの調整をしています.
eeglab redraw
はこの後エラーが出るので，強制的にredrawで止めてます.
ここでredrawされるとポップアップが出るので，
ここで「Remove」を選択すると続きを処理します.

[STUDY ALLEEG] = std_editset( STUDY, ALLEEG, 'name', ExpName, 'filename', ExpName, 'updatedat', 'on', 'savedat', 'on' );
CURRENTSTUDY = 1; EEG = ALLEEG; CURRENTSET = [1:length(EEG)];

eeglab redraw

最後にstudyを保存して終了.
実験名がStudyファイルの名前になるようにしています.
自分でやりたい場合は，ExpNameを直接変更するとできます.

これでフォルダに実験名.studyというファイルができます.

これを作った良さは，また今度.

STUDYファイル②

Fri, 26 Feb 2010 20:25:43 +0900

①からの続きです．
①を踏まえての内容です．

versionはv7.1.3.11.b

さてこの状態でフォルダにすべてのデータファイルが入っています．
たとえば20人，3条件だったら，60個のデータがあるはずです．
EEGLABでは

データの情報を.setに
データの内容を.fdtに

分けて保存するので，実際にフォルダ内のファイルの個数は倍の120個です．
エラーが出たりおかしかったら，ここから確認です．

さてここまでできたらMATLABで以下を実行です．

currdir='E:\FaceEmoCh\FECStudy';

[ALLEEG EEG CURRENTSET ALLCOM] = eeglab;
pop_editoptions( 'option_storedisk', 1, 'option_savetwofiles', 1, 'option_saveica', 0, 'option_single', 1, 'option_memmapdata', 0, 'option_computeica', 1, 'option_scaleicarms', 1, 'option_rememberfolder', 1);
ALLEEG = eeg_checkset(ALLEEG); EEG = eeg_retrieve(ALLEEG, CURRENTSET);
cd(currdir);
allfiles=dir('*.set');
for n=1:length(allfiles) %すべての条件にあうファイルをロードする
loadname=allfiles(n).name; %ファイルの名前を取得してloadnameに格納
dataname=loadname(1:42); %ファイル名の前42文字とする
ExpName=loadname(1:9); %ファイル名の1-9文字で実験名を命名する
EEG = pop_loadset( 'filename', loadname, 'filepath', pwd);
EEG.filepath=pwd;
EEG.subject=loadname(11:15); %ファイルの前11-15文字で参加者名を命名する
EEG.condition=loadname(37:42);%ファイルの後ろから37-42文字で条件名を命名する
[ALLEEG, EEG, CURRENTSET] = eeg_store( ALLEEG, EEG, 0 );
end
ALLEEG = eeg_checkset(ALLEEG); EEG = eeg_retrieve(ALLEEG, CURRENTSET);
[STUDY, ALLEEG] = std_checkset(STUDY, ALLEEG);
eeglab redraw

[STUDY ALLEEG] = std_editset( STUDY, ALLEEG, 'name', ExpName, 'filename', ExpName, 'updatedat', 'on', 'savedat', 'on' );
CURRENTSTUDY = 1; EEG = ALLEEG; CURRENTSET = [1:length(EEG)];

eeglab redraw

この中で
青文字は自分の環境に合わせて変更します．

①自分がデータを入れたフォルダ
②拡張子を除いたファイル名の文字数
③実験名が何文字目から何文字目までか
④参加者名が何文字目から何文字目までか
⑤条件名が何文字目から何文字目までか

途中（一回目のeeglab redraw）でポップアップが出ますので，
（出ないこともあるかも）

Removeを選ぶとその後は自動で続行，STUDYファイルを作成してくれます．

詳しい解説は次で

STUDYファイル①

Fri, 26 Feb 2010 20:09:51 +0900

EEGLABで処理を行う際にSTUDYファイルを作成すると
データの管理に非常に便利です．

市販のソフトなどだとこういうものは勝手に作ってくれるのですが，
EEGLABでは自分できちんと管理しておかなくてはいけません．
ただ逆に言えば自分で管理できる分自由度は上がります．

EEGLAB内ではもちろんいろいろな制約がありますが，
やりたいことはMATLAB上でやってしまってもよいわけです．

なのでこの管理ファイルだけとりあえず作ってしまうと，
データがどこにあるか，どんなデータかなどを
STUDYという一つの構造体から参照できます．

まずは作る前に心構えから．
ファイル名をきれいに管理することが第一歩です．
これはEEGLABやSTUDYファイルに限らず，
大量のデータファイルを扱う場合にはいつでも言えることですね．

個人的には以下の構成をお勧めしています．

「[実験名] [参加者名] [表示しておきたい処理] ... [条件].set」

「」内の半角スペースはそのままファイル名に使います．
区切りをしっかりつけるということですね．
デリミタは半角スペースでなくてもよいのですが，
EEGLAB（MATLABの仕様？）で扱う場合，
半角アンダーバーは「直後の文字の”下付文字”化」
をあらわすようでちょっといやな感じになる時があります．

参加者名は多分「Sub01」とか「s01」でよいと思います．
01にするのは参加者が二桁になる場合，
三桁になる予定なら001として文字数がずれないようにするのがよいです．

表示しておきたい処理にはたとえばいくつのフィルタをかけているとか，
どんなICAをかけたかとか，結果をいろいろ比較したい場合に便利です．
なくても構いません．

条件は実験条件ですね．これも条件間で文字数を一致させるのがよいです．

たとえば私の最近のデータは

「FaceEmoCh Sub01 fil01 ICArej posich.set」

のような感じです．

このような
参加者ごと，条件ごと
のデータができたら，
それを一つのフォルダに入れてしまいます．

これでSTUDYファイルを作る準備ができました．

もちろんGUIから作ることもできます.

ただ一つ一つ情報を打ち込んでいくのはめんどいですので，
コマンドラインかららくらくやる方法を紹介しています.

続く