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

shock_decomposition.m has to rely on plot_shock_decomposition.m, to make behaviour robust.

time-shift
Marco Ratto 2017-01-19 16:52:11 +01:00 committed by Stéphane Adjemian (Charybdis)
parent f66191f59f
commit d304fd26c1
1 changed files with 2 additions and 29 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
matlab/shock_decomposition.m
View File

@ -120,34 +120,7 @@ end
oo_.shock_decomposition = z;
if options_.use_shock_groups
shock_groups = M_.shock_groups.(options_.use_shock_groups);
shock_ind = fieldnames(shock_groups);
ngroups = length(shock_ind);
shock_names = shock_ind;
for i=1:ngroups,
shock_names{i} = (shock_groups.(shock_ind{i}).label);
end
zz = zeros(endo_nbr,ngroups+2,gend);
for i=1:ngroups
for j = shock_groups.(shock_ind{i}).shocks
k = find(strcmp(j,cellstr(M_.exo_names)));
zz(:,i,:) = zz(:,i,:) + z(:,k,:);
z(:,k,:) = 0;
end
end
zothers = sum(z(:,1:nshocks,:),2);
zz(:,ngroups+1,:) = sum(z(:,1:nshocks+1,:),2);
if any(any(zothers)),
shock_names = [shock_names; {'Others + Initial Values'}];
end
zz(:,ngroups+2,:) = z(:,nshocks+2,:);
z = zz;
else
shock_names = M_.exo_names;
end
if ~options_.no_graph.shock_decomposition
graph_decomp(z,shock_names,M_.endo_names,i_var,options_.initial_date,M_,options_)
plot_shock_decomposition(M_,oo_,options_,varlist);
end