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

Merge branch 'gsa_plot' into 'master' 

Fixes to GSA plots

See merge request Dynare/dynare!1631
time-shift
Sébastien Villemot 2018-11-15 09:12:29 +00:00
parent 7b0ef3d080 30361ccb49
commit 6e5465c0ff
3 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
matlab/gsa/map_calibration.m
View File

@ -42,7 +42,6 @@ end
pvalue_ks = DynareOptions.opt_gsa.pvalue_ks; pvalue_ks = DynareOptions.opt_gsa.pvalue_ks;
indx_irf = []; indx_irf = [];
indx_moment = []; indx_moment = [];
DynareOptions.nodisplay = 1;
init = ~DynareOptions.opt_gsa.load_stab; init = ~DynareOptions.opt_gsa.load_stab;
options_mcf.pvalue_ks = DynareOptions.opt_gsa.pvalue_ks; options_mcf.pvalue_ks = DynareOptions.opt_gsa.pvalue_ks;

4
matlab/gsa/mcf_analysis.m
View File

@ -34,6 +34,8 @@ if DynareOptions.TeX
else else
param_names_tex = options_mcf.param_names_tex; param_names_tex = options_mcf.param_names_tex;
end end
else
param_names_tex = strrep(options_mcf.param_names,'_','\_');
end end
amcf_name = options_mcf.amcf_name; amcf_name = options_mcf.amcf_name;
amcf_title = options_mcf.amcf_title; amcf_title = options_mcf.amcf_title;
@ -78,7 +80,7 @@ if ~isempty(indmcf) && ~DynareOptions.nograph
skipline() skipline()
xx=[]; xx=[];
if ~ isempty(xparam1), xx=xparam1(indmcf); end if ~ isempty(xparam1), xx=xparam1(indmcf); end
scatter_mcf(lpmat(ibeha,indmcf),lpmat(inobeha,indmcf), param_names(indmcf), ... scatter_mcf(lpmat(ibeha,indmcf),lpmat(inobeha,indmcf), param_names_tex(indmcf), ...
'.', [fname_,'_',amcf_name], OutputDirectoryName, amcf_title,xx, DynareOptions, ... '.', [fname_,'_',amcf_name], OutputDirectoryName, amcf_title,xx, DynareOptions, ...
beha_title, nobeha_title) beha_title, nobeha_title)
end end

4
matlab/gsa/scatter_mcf.m
View File

@ -143,7 +143,7 @@ for i = 1:p
end end
if i==1 if i==1
if nflag == 1 if nflag == 1
ylabel(vnames(j,:),'Rotation',45,'interpreter','none', ... ylabel(vnames(j,:),'Rotation',45, ...
'HorizontalAlignment','right','VerticalAlignment','middle'); 'HorizontalAlignment','right','VerticalAlignment','middle');
else else
ylabel([num2str(j),' '],'Rotation',90) ylabel([num2str(j),' '],'Rotation',90)
@ -151,7 +151,7 @@ for i = 1:p
end end
if j==1 if j==1
if nflag == 1 if nflag == 1
title(vnames(i,:),'interpreter','none','Rotation',45, ... title(vnames(i,:),'Rotation',45, ...
'HorizontalAlignment','left','VerticalAlignment','bottom') 'HorizontalAlignment','left','VerticalAlignment','bottom')
else else
title(num2str(i)) title(num2str(i))