stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

initialize mode field in slice sampler options

time-shift
Marco Ratto 2016-05-14 14:21:00 +02:00 committed by Johannes Pfeifer
parent e5ab7ed992
commit 1d6961cc4c
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
matlab/check_posterior_sampler_options.m
View File

@ -304,10 +304,10 @@ if init,
% options_.mh_posterior_mode_estimation = 0; % options_.mh_posterior_mode_estimation = 0;
% this is done below, but perhaps preprocessing should do this? % this is done below, but perhaps preprocessing should do this?
if ~isfield(posterior_sampler_options,'mode'), if ~isempty(posterior_sampler_options.mode)
posterior_sampler_options.mode = []; % multimodal case
else % multimodal case
posterior_sampler_options.rotated = 1; posterior_sampler_options.rotated = 1;
posterior_sampler_options.WR=[];
end end
% posterior_sampler_options = set_default_option(posterior_sampler_options,'mode_files',[]); % posterior_sampler_options = set_default_option(posterior_sampler_options,'mode_files',[]);

1
matlab/global_initialization.m
View File

@ -476,6 +476,7 @@ options_.posterior_sampler_options.slice.slice_initialize_with_mode=0;
options_.posterior_sampler_options.slice.use_mh_covariance_matrix=0; options_.posterior_sampler_options.slice.use_mh_covariance_matrix=0;
options_.posterior_sampler_options.slice.WR=[]; options_.posterior_sampler_options.slice.WR=[];
options_.posterior_sampler_options.slice.mode_files=[]; options_.posterior_sampler_options.slice.mode_files=[];
options_.posterior_sampler_options.slice.mode=[];
options_.posterior_sampler_options.slice.initial_step_size=0.8; options_.posterior_sampler_options.slice.initial_step_size=0.8;
% Independent Metropolis-Hastings % Independent Metropolis-Hastings
options_.posterior_sampler_options.imh.proposal_distribution = 'rand_multivariate_normal'; options_.posterior_sampler_options.imh.proposal_distribution = 'rand_multivariate_normal';