Measuring the Effects of a Shock to Monetary Policy - Humboldt ...

112 Bayesian FAVARs with Agnostic Identification
end
%directly called from DO_RESULTS
d = size(squeeze( (results.ira.finalresponses(1).response(1,:,:))),2);
draws = 1:10:d; %show every xth accepted draw response
elseif k(i) == 5
load FOR_MESH_K5_RESULTS;
d = size(squeeze( (results.ira.finalresponses(1).response(1,:,:))),2);
draws = 1:10:d;
elseif k(i) == 7
end
load FOR_MESH_K7_RESULTS;
d = size(squeeze( (results.ira.finalresponses(1).response(1,:,:))),2);
draws = 1:d;
disp (sprintf(’K%d data ready’,k(i)));
n = input.specification.IRA.nsteps;
srs = 1:size(squeeze( (results.ira.finalresponses(1).response(1,:,draws))),1);
smax = size(squeeze( (results.ira.finalresponses(1).response(1,:,draws))),1);
figure(k(i))
axis normal;
grid on;
axis tight;
subplot(2,2,1);
mesh (squeeze(results.ira.finalresponses(1).response(1,srs,draws)));
ylabel(’Horizon’);
xlabel(’Accepted Response Draws’);
zlabel(’Response Scale’);
title (’FFR’);
subplot(2,2,2);
mesh (squeeze(results.ira.finalresponses(1).response(9,srs,draws)));
zlabel(’Response Scale’);
title (’COMMODITY PRICE INDEX’);
subplot(2,2,3);
mesh (squeeze(results.ira.finalresponses(1).response(14,srs,draws)));
zlabel(’Response Scale’);
title ( ’UNEMPLOYMENT’);
subplot(2,2,4);
mesh (squeeze(results.ira.finalresponses(1).response(10,srs,draws)));
zlabel(’Response Scale’);
title ( ’CAPACITY UTIL RATE’);
saveas(gcf, sprintf ( ’FIG_MESH_K%d’ ,k(i)) , ’fig’);
saveas(gcf, sprintf ( ’FIG_MESH_K%d’ ,k(i)) , ’jpg’);