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.Ta0       = 10;
0026 bayes_parm.Ta0_act = 10; 
0027 
0028 bayes_parm.v0     = 0.01; 
0029 bayes_parm.Tv0    = 0; 
0030 
0031 % --- Input file(relative path from proj_root)
0032 bayes_parm.brainfile = [];
0033 bayes_parm.areafile  = [];
0034 bayes_parm.actfile   = [];
0035 bayes_parm.megfile   = [];
0036 bayes_parm.basisfile = [];
0037 
0038 % --- Output file(relative path from proj_root)
0039 bayes_parm.bayesfile = [];
0040 
0041 % --- Specify time window by sample numbar
0042 bayes_parm.twin_meg  = []; % Time window for analysis
0043 bayes_parm.Tperiod   = []; % Time period for current estimation
0044 bayes_parm.Tnext     = []; % Time step for next period
0045 
0046 % --- Current area and fMRI activity
0047 bayes_parm.area_key  = 'Cortex';
0048 bayes_parm.act_key   = {'Uniform'};
0049 
0050 %% --- Spatial smoothing radius & max radius
0051 bayes_parm.Rfilt  = [6e-3]; % [m]
0052 
0053 %% --- Global window
0054 bayes_parm.reduce_global = 1;
0055 bayes_parm.Rfilt_global  = bayes_parm.Rfilt;
0056 bayes_parm.area_key_global = 'Cortex';
0057 
0058 %% --- Soft normal constraint
0059 bayes_parm.soft_mode = 0;    % = 1: Soft normal constraint
0060 % tangential variance coefficient
0061 bayes_parm.tan_var = 1 ;
0062 % tangential variance max
0063 bayes_parm.var_max = 1/sqrt(2) ;
0064 % ÀÜÀþÊý¸þʬ»¶ = sin( tan_var * Ê¿¶ÑË¡ÀþÊý¸þ¤«¤é¤Î³ÑÅÙÊ¿¶Ñ) <= var_max
0065 
0066 % --- Estimation model
0067 bayes_parm.forward_model    = 'focal';
0068 bayes_parm.noise_estimation_model = 1;
0069 bayes_parm.noise_model = 3;
0070 bayes_parm.update_sx = ON;
0071 bayes_parm.update_v  = OFF;
0072 bayes_parm.flag_baseline  = 0;
0073 bayes_parm.flag_noise_set = 0;
0074 %bayes_parm.megnoisefile    = []; % relative path from proj_root
0075 bayes_parm.megfile_baseline = []; % relative path from proj_root
0076 bayes_parm.twin_noise      = [];
0077 %bayes_parm.twin_global     = [];
0078 bayes_parm.twin_baseline    = [];
0079 bayes_parm.basisfile_global = [];% relative path from proj_root
0080 
0081 %% --- Flag
0082 bayes_parm.patch_norm      = OFF;
0083 bayes_parm.cont_pr         = OFF; 
0084 bayes_parm.temporal_filter = OFF;
0085 bayes_parm.trial_average   = OFF; 
0086 bayes_parm.expand_spatial_filter = ON;     
0087 bayes_parm.remove_crossed_area   = ON;
0088 bayes_parm.variance_orientation  = OFF;  
0089 
0090 %% --- Optinal parameters
0091 bayes_parm.reduce = 1;
0092 bayes_parm.cosval = cos(pi);
0093 bayes_parm.Fdmin = 1e-50;
0094 bayes_parm.a_min = 1e-6;
0095 bayes_parm.a_max = 1e6;
0096 bayes_parm.noise_reg = 0.1;
0097