*Fieldtrip: Planarをcombineする方法
<data_import>
dataをインポートするときにセンサーの位置情報(ft_read_sens)を保持する必要がある。dataをインポートした時にVariablesに'grad'があるか確認する。
% data import
cfg = [];
cfg.dataset = 'filename';
% data segment
cfg.trialdef.eventtype = 'STI101';
cfg.trialdef.prestim = time;
cfg.trialdef.poststim = time;
cfg.trialdef.eventvalue = eve;
cfg=ft_definetrial(cfg);
% channel select
cfg.channel = {'MEG'};
% baseline correction
cfg.demean = 'yes';
cfg.baselinewindow = [time time]
% filter
cfg.lpfilter = 'yes';
cfg.lpfreq = 30;
% senseor position
grad = ft_read_sens(filename);
% Preprocessing
data1 = ft_preprocessing(cfg);
<data_average>
個人のデータが複数ファイルに分割されている場合、個人内でグランドアベレージをする。この時にoptionで'within'を選択する。
avdata1 = ft_timelockanalysis(cfg,data1); %※file数分繰返す
%
cfg = [];
cfg.method = 'within'
grandavedata = ft_timelockgrandaverage(cfg,avdata1,avdata2,…);
<combine planar>
センサーの位置情報を保持したAverageFileを用いてプラナー電極をコンバインする
% combine planar
cfg = [];
grandavedata _cmb = ft_combineplanar(cfg, grandavedata);
Planar_combine後はレイアウトは'neuromag306cmb.lay'を使う