vb_set_bayes_default_parameters

 Set default value for bayes prior and estimation parameters

 --- Syntax
 bayes_parm = vb_set_bayes_default_parameters;
 bayes_parm = vb_set_bayes_default_parameters(bayes_parm);

 --- History
 2008-11-05 Taku Yoshioka

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function bayes_parm = vb_set_bayes_default_parameters(bayes_parm)
0002 % Set default value for bayes prior and estimation parameters
0003 %
0004 % --- Syntax
0005 % bayes_parm = vb_set_bayes_default_parameters;
0006 % bayes_parm = vb_set_bayes_default_parameters(bayes_parm);
0007 %
0008 % --- History
0009 % 2008-11-05 Taku Yoshioka
0010 %
0011 % Copyright (C) 2011, ATR All Rights Reserved.
0012 % License : New BSD License(see VBMEG_LICENSE.txt)
0013 
0014 %% --- Iteration number
0015 bayes_parm.Npre_train = 1000;
0016 bayes_parm.Ntrain     = 1000;
0017 bayes_parm.skip       = 50;
0018 
0019 %% --- Prior setting
0020 bayes_parm.a0     = 1;
0021 bayes_parm.a0_act = 100; 
0022 % 暫定対応 a0_act=[] としてa0_act を自動推定する設定にすると
0023 % act パターンが実質的にUniform パターンになる。調査中。
0024 
0025 bayes_parm.prior_weight = 10e-5; % relative influence of prior information on data(from 0 to 1 real)
0026 
0027 bayes_parm.v0     = 0.01; 
0028 bayes_parm.Tv0    = 0; 
0029 
0030 % --- Input file(relative path from proj_root)
0031 bayes_parm.brainfile = [];
0032 bayes_parm.areafile  = [];
0033 bayes_parm.actfile   = [];
0034 bayes_parm.megfile   = [];
0035 bayes_parm.basisfile = [];
0036 
0037 % --- Output file(relative path from proj_root)
0038 bayes_parm.bayesfile = [];
0039 
0040 % --- Specify time window by sample numbar
0041 bayes_parm.twin_meg  = []; % Time window for analysis
0042 bayes_parm.Tperiod   = []; % Time period for current estimation
0043 bayes_parm.Tnext     = []; % Time step for next period
0044 
0045 % --- Current area and fMRI activity
0046 bayes_parm.area_key  = 'Cortex';
0047 bayes_parm.act_key   = {'Uniform'};
0048 
0049 %% --- Spatial smoothing radius & max radius
0050 bayes_parm.Rfilt  = [6e-3]; % [m]
0051 
0052 %% --- Global window
0053 bayes_parm.reduce_global = 1;
0054 bayes_parm.Rfilt_global  = bayes_parm.Rfilt;
0055 bayes_parm.area_key_global = 'Cortex';
0056 
0057 %% --- Soft normal constraint
0058 bayes_parm.soft_mode = 0;    % = 1: Soft normal constraint
0059 % tangential variance coefficient
0060 bayes_parm.tan_var = 1 ;
0061 % tangential variance max
0062 bayes_parm.var_max = 1/sqrt(2) ;
0063 % 接線方向分散 = sin( tan_var * 平均法線方向からの角度平均) <= var_max
0064 
0065 % --- Estimation model
0066 bayes_parm.forward_model    = 'focal';
0067 bayes_parm.noise_estimation_model = 1;
0068 bayes_parm.noise_model = 3;
0069 bayes_parm.update_sx = ON;
0070 bayes_parm.update_v  = OFF;
0071 bayes_parm.flag_baseline  = 0;
0072 bayes_parm.flag_noise_set = 0;
0073 %bayes_parm.megnoisefile    = []; % relative path from proj_root
0074 bayes_parm.megfile_baseline = []; % relative path from proj_root
0075 bayes_parm.twin_noise      = [];
0076 %bayes_parm.twin_global     = [];
0077 bayes_parm.twin_baseline    = [];
0078 bayes_parm.basisfile_global = [];% relative path from proj_root
0079 
0080 %% --- Flag
0081 bayes_parm.patch_norm      = OFF;
0082 bayes_parm.cont_pr         = OFF; 
0083 bayes_parm.temporal_filter = OFF;
0084 bayes_parm.trial_average   = OFF; 
0085 bayes_parm.expand_spatial_filter = ON;     
0086 bayes_parm.spatial_smoother ='subj'; % spatial filter calculate from
0087                              % = 'std': Standard brain.
0088                              % = 'subj' Individual brain.
0089 bayes_parm.remove_crossed_area   = ON;
0090 bayes_parm.variance_orientation  = OFF;  
0091 
0092 %% --- Optinal parameters
0093 bayes_parm.reduce = 1;
0094 bayes_parm.cosval = cos(pi);
0095 bayes_parm.Fdmin = 1e-50;
0096 bayes_parm.a_min = 1e-6;
0097 bayes_parm.a_max = 1e6;
0098 bayes_parm.noise_reg = 0.1;
0099