*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'を使う