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

Added legends for the figures.

time-shift
Jukka Heikkonen 2014-11-14 14:53:37 +01:00
parent 90487c03e0
commit 83931dca66
1 changed files with 34 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
matlab/gsa/map_calibration.m
View File

@ -98,6 +98,13 @@ if ~isempty(indx_irf),
if nrow*(nrow-1)>nbr_irf_restrictions,
ncol=nrow-1;
end
% For single legend search which has maximum nbr of restrictions
maxijv=0;
for ij=1:nbr_irf_restrictions
if length(DynareOptions.endogenous_prior_restrictions.irf{ij,3})>maxijv
maxij=ij;maxijv=length(DynareOptions.endogenous_prior_restrictions.irf{ij,3});
end
end
for ij=1:nbr_irf_restrictions,
figure(h1),
mat_irf{ij}=mat_irf{ij}(irestrictions,:);
@ -116,13 +123,23 @@ if ~isempty(indx_irf),
leg = num2str(DynareOptions.endogenous_prior_restrictions.irf{ij,3}(1));
if size(mat_irf{ij},2)>1,
leg = [leg,':' ,num2str(DynareOptions.endogenous_prior_restrictions.irf{ij,3}(end))];
end
title([DynareOptions.endogenous_prior_restrictions.irf{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.irf{ij,2}, '(', leg,')'],'interpreter','none'),
end
title([DynareOptions.endogenous_prior_restrictions.irf{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.irf{ij,2}, '(', leg,')'],'interpreter','none'),
%set(legend_h,'Xlim',[0 1]);
if ij==maxij
leg1 = num2str(DynareOptions.endogenous_prior_restrictions.irf{ij,3}(:));
[legend_h,object_h,plot_h,text_strings]=legend(leg1);
Position=get(legend_h,'Position');Position(1:2)=[-0.055 0.95-Position(4)];
set(legend_h,'Position',Position);
end
% hc = get(h,'Children');
%for i=2:2:length(hc)
%end
indx1 = find(indx_irf(:,ij)==0);
indx2 = find(indx_irf(:,ij)~=0);
atitle=[DynareOptions.endogenous_prior_restrictions.irf{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.irf{ij,2}, '(', leg,')'];
fprintf(['%4.1f%% of the prior support matches IRF ',atitle,' inside [%4.1f, %4.1f]\n'],length(indx1)/length(irestrictions)*100,DynareOptions.endogenous_prior_restrictions.irf{ij,4})
% aname=[type '_irf_calib_',int2str(ij)];
aname=[type '_irf_calib_',endo_prior_restrictions.irf{ij,1},'_VS_',endo_prior_restrictions.irf{ij,2}];
atitle=[type ' IRF Calib: Parameter(s) driving ',DynareOptions.endogenous_prior_restrictions.irf{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.irf{ij,2}, '(', leg,')'];
[proba, dproba] = stab_map_1(xmat, indx1, indx2, aname, 0);
@ -143,6 +160,13 @@ if ~isempty(indx_moment)
if nrow*(nrow-1)>nbr_moment_restrictions,
ncol=nrow-1;
end
% For single legend search which has maximum nbr of restrictions
maxijv=0;
for ij=1:nbr_moment_restrictions
if length(DynareOptions.endogenous_prior_restrictions.moment{ij,3})>maxijv
maxij=ij;maxijv=length(DynareOptions.endogenous_prior_restrictions.moment{ij,3});
end
end
for ij=1:nbr_moment_restrictions,
figure(h2),
mat_moment{ij}=mat_moment{ij}(irestrictions,:);
@ -163,11 +187,17 @@ if ~isempty(indx_moment)
leg = [leg,':' ,num2str(DynareOptions.endogenous_prior_restrictions.moment{ij,3}(end))];
end
title([DynareOptions.endogenous_prior_restrictions.moment{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.moment{ij,2},'(',leg,')'],'interpreter','none'),
if ij==maxij
leg1 = num2str(DynareOptions.endogenous_prior_restrictions.moment{ij,3}(:));
[legend_h,object_h,plot_h,text_strings]=legend(leg1);
Position=get(legend_h,'Position');Position(1:2)=[-0.055 0.95-Position(4)];
set(legend_h,'Position',Position);
end
indx1 = find(indx_moment(:,ij)==0);
indx2 = find(indx_moment(:,ij)~=0);
atitle=[DynareOptions.endogenous_prior_restrictions.moment{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.moment{ij,2}, '(', leg,')'];
fprintf(['%4.1f%% of the prior support matches MOMENT ',atitle,' inside [%4.1f, %4.1f]\n'],length(indx1)/length(irestrictions)*100,DynareOptions.endogenous_prior_restrictions.moment{ij,4})
% aname=[type '_moment_calib_',int2str(ij)];
aname=[type '_moment_calib_',DynareOptions.endogenous_prior_restrictions.moment(ij,1),'_VS_',DynareOptions.endogenous_prior_restrictions.moment(ij,2)];
atitle=[type ' MOMENT Calib: Parameter(s) driving ',DynareOptions.endogenous_prior_restrictions.moment{ij,1},' vs ',DynareOptions.endogenous_prior_restrictions.moment{ij,2}, '(', leg,')'];
[proba, dproba] = stab_map_1(xmat, indx1, indx2, aname, 0);