dist irf: 1) 5% shock, 2) plot both ss dist and dist irfs under the t…

…rue parameters
mikkelpm · Jan 7, 2021 · 12a309a · 12a309a
1 parent 9ab6417
commit 12a309a
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/program/hh_model/plot/comput_polfct_distirf.m b/program/hh_model/plot/comput_polfct_distirf.m
@@ -155,7 +155,8 @@
                     'linewidth',1,'edgealpha',alpha_hairline);
             end
 
-            plot(assets_fine_truth,distirf_truth(:,iEpsilon,1)+the_y_wedge,'k-','linewidth',.75); % Truth
+            plot(assets_fine_truth,distirf_truth(:,iEpsilon,1)+the_y_wedge,'k--','linewidth',.5); % Steady state distribution under the true parameters
+            plot(assets_fine_truth,distirf_truth(:,iEpsilon,the_horz+1)+the_y_wedge,'k-','linewidth',.5); % IRF of cross-sec distribution under the true parameters
 
             text(xlim_assets(1)+xloc_text*diff(xlim_assets),the_y_wedge+wedge_text,...
                 ['h = ' num2str(the_horz)],'FontSize',plot_fontsize-2,'FontWeight','bold'); % Horizon label

diff --git a/program/plot_mcmc.m b/program/plot_mcmc.m
@@ -67,7 +67,7 @@
 
 % Specific layout of Asset distribution IRF
 horzs = [0 2 4 8];              % Impulse response horizons to plot
-shock = -10;                    % Shock (unit = 1sd)
+shock = 0.05/0.014;             % Shock (unit = 1sd): here 1sd = 0.014 and we consider a 5% shock
 graph_size_distirf = [6 3.5];   % Graph size
 ylim_distirf = [0 .9];          % y-axis limits
 wedge_param = [.8 -.2];         % Shift distributions vertically by wedge_param(1)+wedge_param(2)*horizon, in order to display on single plot