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

Added the possibility to save the plots generated by the shock_decomposition command.

time-shift
Stéphane Adjemian (Charybdis) 2012-08-24 15:03:35 +02:00
parent 150256268f
commit fe98a0875d
2 changed files with 7 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/graph_decomp.m
View File

@ -1,4 +1,4 @@
function []=graph_decomp(z,shock_names,endo_names,i_var,initial_date)
function []=graph_decomp(z,shock_names,endo_names,i_var,initial_date,DynareModel,DynareOptions)
%function []=graph_decomp(z,varlist,initial_period,freq)
% Copyright (C) 2010-2011 Dynare Team
@ -39,7 +39,7 @@ for j=1:nvar
if ymax-ymin < 1e-6
continue
end
figure('Name',endo_names(i_var(j),:));
fhandle = dyn_figure(DynareOptions,'Name',endo_names(i_var(j),:));
ax=axes('Position',[0.1 0.1 0.6 0.8]);
axis(ax,[xmin xmax ymin ymax]);
plot(ax,x(2:end),z1(end,:),'k-','LineWidth',2)
@ -48,7 +48,7 @@ for j=1:nvar
i_1 = i-1;
yp = 0;
ym = 0;
for k = 1:comp_nbr
for k = 1:comp_nbr
zz = z1(k,i);
if zz > 0
fill([x(i) x(i) x(i+1) x(i+1)],[yp yp+zz yp+zz yp],k);
@ -70,7 +70,7 @@ for j=1:nvar
y1 = 0;
height = 1/comp_nbr;
labels = char(shock_names,'Initial values');
for i=1:comp_nbr
fill([0 0 0.2 0.2],[y1 y1+0.7*height y1+0.7*height y1],i);
hold on
@ -78,5 +78,7 @@ for j=1:nvar
hold on
y1 = y1 + height;
end
dyn_saveas(fhandle,[DynareModel.fname '_shock_decomposition_' endo_names(i_var(j),:)],DynareOptions);
hold off
end

2
matlab/shock_decomposition.m
View File

@ -107,4 +107,4 @@ end
oo_.shock_decomposition = z;
graph_decomp(z,M_.exo_names,M_.endo_names,i_var,options_.initial_date)
graph_decomp(z,M_.exo_names,M_.endo_names,i_var,options_.initial_date,M_,options_)