Home > vbmeg > functions > template > vb_set_preprocess_meg3_parm.m

vb_set_preprocess_meg3_parm

PURPOSE ^

Reject noisy sensors and trials, and divide trials according to

SYNOPSIS ^

function parm = vb_set_preprocess_meg3_parm

DESCRIPTION ^

 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)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SUBFUNCTIONS ^

SOURCE CODE ^

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 

Generated on Mon 22-May-2023 06:53:56 by m2html © 2005