Skip to content
Snippets Groups Projects
Commit ab5b922e authored by Marco Ratto's avatar Marco Ratto
Browse files

-) fixed bug introduced in commit 66527df;

-) estimation options for data set must be given also when ppost=1 (fix in testsuite);
parent afc24549
Branches
Tags
No related merge requests found
...@@ -31,6 +31,7 @@ lgy_ = M_.endo_names; ...@@ -31,6 +31,7 @@ lgy_ = M_.endo_names;
x0=[]; x0=[];
options_gsa = set_default_option(options_gsa,'datafile',[]); options_gsa = set_default_option(options_gsa,'datafile',[]);
options_gsa = set_default_option(options_gsa,'rmse',0);
if isfield(options_gsa,'nograph'), if isfield(options_gsa,'nograph'),
options_.nograph=options_gsa.nograph; options_.nograph=options_gsa.nograph;
end end
...@@ -46,7 +47,14 @@ end ...@@ -46,7 +47,14 @@ end
options_.order = 1; options_.order = 1;
if ~isempty(options_gsa.datafile) || isempty(bayestopt_), if ~isempty(options_gsa.datafile) || isempty(bayestopt_) || options_gsa.rmse,
if isempty(options_gsa.datafile) && options_gsa.rmse,
disp('The data file and all relevant estimation options ')
disp('[first_obs, nobs, presample, prefilter, loglinear, lik_init, kalman_algo, ...]')
disp('must be specified for RMSE analysis!');
error('Sensitivity anaysis error!');
end
options_.datafile = options_gsa.datafile; options_.datafile = options_gsa.datafile;
if isfield(options_gsa,'first_obs'), if isfield(options_gsa,'first_obs'),
options_.first_obs=options_gsa.first_obs; options_.first_obs=options_gsa.first_obs;
...@@ -142,7 +150,6 @@ options_gsa = set_default_option(options_gsa,'namendo',[]); ...@@ -142,7 +150,6 @@ options_gsa = set_default_option(options_gsa,'namendo',[]);
options_gsa = set_default_option(options_gsa,'namlagendo',[]); options_gsa = set_default_option(options_gsa,'namlagendo',[]);
options_gsa = set_default_option(options_gsa,'namexo',[]); options_gsa = set_default_option(options_gsa,'namexo',[]);
% RMSE mapping % RMSE mapping
options_gsa = set_default_option(options_gsa,'rmse',0);
options_gsa = set_default_option(options_gsa,'lik_only',0); options_gsa = set_default_option(options_gsa,'lik_only',0);
options_gsa = set_default_option(options_gsa,'var_rmse',options_.varobs); options_gsa = set_default_option(options_gsa,'var_rmse',options_.varobs);
options_gsa = set_default_option(options_gsa,'pfilt_rmse',0.1); options_gsa = set_default_option(options_gsa,'pfilt_rmse',0.1);
......
...@@ -214,5 +214,6 @@ disp('Press ENTER to continue'); pause(5); ...@@ -214,5 +214,6 @@ disp('Press ENTER to continue'); pause(5);
// prefilter=1,mh_jscale=0.5,mh_replic=0, mode_compute=0, mh_drop=0.6); // prefilter=1,mh_jscale=0.5,mh_replic=0, mode_compute=0, mh_drop=0.6);
dynare_sensitivity(nodisplay, graph_format=fig,stab=0, // no need for stability analysis since the posterior sample is surely OK dynare_sensitivity(nodisplay, graph_format=fig,stab=0, // no need for stability analysis since the posterior sample is surely OK
datafile=data_ca1,first_obs=8,nobs=79,prefilter=1,
rmse=1,ppost=1); rmse=1,ppost=1);
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment