Reject noisy sensors and trials, and divide trials according to
experimental conditions in division file. If division file
(divfile) was not specified, MEG data will not be divided.
If parm was not specified, recommended values will be used.
Input parameters (recommended value)
- parm.th_eog_v : Threshold for EOG signal in vertical direction
(0.75)
- parm.th_eog_h : Threshold for EOG signal in horizontal
direction (0.75)
- parm.th_meg : Threshold for noise in MEG data (8)
- parm.meg_coef : Threshold for power of MEG data considered as
a candidate of rejected sensor (5.0)
- parm.Ntrial_th : The number of trials over which MEG time
courses exhibit noise on the same sensor (5)
- parm.rm_trial : It can be used to directly specify trials to
be rejected.
- parm.verify : If it is ON (constant), the script shows
trials that should be rejected and the user
must verify each trial by clicking 'OK' button
on a dialog (ON).
- parm.freq : MEG data frequency.
- parm.driftremovalmode
: 'linear' -- remove drift with linear fitting.
: 'guodra' -- remove drift with quodra fitting.
: 'sin' -- remove drift with sin-cos fitting.
- parm.basiocycle: when parm.driftrimovalmode = 'sin', you need
specify this parameter as vector of
basio-cycles, for example parm.basiocycle =
[4,8,16,32] (second). It is recomended that
basiocycle > trial time length.
- parm.baseline_pretrig
: If this parameter is ON, baseline correction
by using median of MEG data at pretrigger
period for each sensor.
- parm.jobmode : Select jobs and order of that in this
preprocessing.
NOTICE that the result of preprocessing would
be changed if you specify the order of jobs.
The following (default) order is recomended.
1 -- Reject dead sensors
2 -- Remove drift and correct baseline
3 -- Remove trials in which EOG signal exceeds
a threshold.
4 -- Reject sensors on which MEG time courses
of some trials are noisy.
5 -- Reject noisy trials
2004-12-28 Taku Yoshioka
2005-08-16 Taku Yoshioka
2006-06-14 Dai Kawawaki
---
Copyright (C) 2011, ATR All Rights Reserved.
License : New BSD License(see VBMEG_LICENSE.txt)0001 function parm = vb_set_preprocess_meg3_parm 0002 % Reject noisy sensors and trials, and divide trials according to 0003 % experimental conditions in division file. If division file 0004 % (divfile) was not specified, MEG data will not be divided. 0005 % If parm was not specified, recommended values will be used. 0006 % 0007 % Input parameters (recommended value) 0008 % - parm.th_eog_v : Threshold for EOG signal in vertical direction 0009 % (0.75) 0010 % - parm.th_eog_h : Threshold for EOG signal in horizontal 0011 % direction (0.75) 0012 % - parm.th_meg : Threshold for noise in MEG data (8) 0013 % - parm.meg_coef : Threshold for power of MEG data considered as 0014 % a candidate of rejected sensor (5.0) 0015 % - parm.Ntrial_th : The number of trials over which MEG time 0016 % courses exhibit noise on the same sensor (5) 0017 % - parm.rm_trial : It can be used to directly specify trials to 0018 % be rejected. 0019 % - parm.verify : If it is ON (constant), the script shows 0020 % trials that should be rejected and the user 0021 % must verify each trial by clicking 'OK' button 0022 % on a dialog (ON). 0023 % - parm.freq : MEG data frequency. 0024 % - parm.driftremovalmode 0025 % : 'linear' -- remove drift with linear fitting. 0026 % : 'guodra' -- remove drift with quodra fitting. 0027 % : 'sin' -- remove drift with sin-cos fitting. 0028 % - parm.basiocycle: when parm.driftrimovalmode = 'sin', you need 0029 % specify this parameter as vector of 0030 % basio-cycles, for example parm.basiocycle = 0031 % [4,8,16,32] (second). It is recomended that 0032 % basiocycle > trial time length. 0033 % - parm.baseline_pretrig 0034 % : If this parameter is ON, baseline correction 0035 % by using median of MEG data at pretrigger 0036 % period for each sensor. 0037 % - parm.jobmode : Select jobs and order of that in this 0038 % preprocessing. 0039 % NOTICE that the result of preprocessing would 0040 % be changed if you specify the order of jobs. 0041 % The following (default) order is recomended. 0042 % 1 -- Reject dead sensors 0043 % 2 -- Remove drift and correct baseline 0044 % 3 -- Remove trials in which EOG signal exceeds 0045 % a threshold. 0046 % 4 -- Reject sensors on which MEG time courses 0047 % of some trials are noisy. 0048 % 5 -- Reject noisy trials 0049 % 2004-12-28 Taku Yoshioka 0050 % 2005-08-16 Taku Yoshioka 0051 % 2006-06-14 Dai Kawawaki 0052 %--- 0053 % 0054 % Copyright (C) 2011, ATR All Rights Reserved. 0055 % License : New BSD License(see VBMEG_LICENSE.txt) 0056 0057 parm.jobmode = [1,2,3,4,5]; 0058 0059 parm.divfile = ['']; 0060 parm.megfile = ['']; 0061 0062 %parm.num_eog_channels = 2; 0063 %parm.eog_channel_id = [211 212]; 0064 %parm.th_eog = [0.75 0.75]; 0065 %parm.th_meg = 8; 0066 %parm.meg_coef = 5; 0067 %parm.Ntrial_th = 5; 0068 %parm.rm_trial = []; 0069 %parm.verify = ON; 0070 0071 parm.freq = 1000; 0072 parm.driftremovalmode = 'sin'; 0073 parm.basiocycle = [4,8,16,32];%if parm.driftremovalmode = 'sin' 0074 parm.baseline_pretrig = ON; 0075 0076 return 0077 0078 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0079 function parm = set_default_preprocess_parm 0080 0081 parm.num_eog_channels = 2; 0082 parm.eog_channel_id = [211 212]; 0083 parm.th_eog = [0.75 0.75]; 0084 parm.th_eog_v = 0.75; 0085 parm.th_eog_h = 0.75; 0086 parm.th_meg = 8; 0087 parm.meg_coef = 5.0; 0088 parm.Ntrial_th = 5; 0089 parm.rm_trial = []; 0090 parm.verify = ON; 0091 0092 return 0093