clear all;
clc;

load intbank31_results;

names1 = char('\epsilon_{z}','\epsilon_{h}','\epsilon_{mi}','\epsilon_{me}','\epsilon_{a}','\epsilon_{l}','\epsilon_{qk}','\epsilon_{y}','\epsilon_{bh}',...
              '\epsilon_{be}','\epsilon_{d}','\epsilon_{r}','\epsilon_{G}','\epsilon_{\delta^{d}}', '\epsilon_{kb}','\epsilon_{ib}');
names3 = char('\rho_{z}','\rho_{h}','\rho_{mi}','\rho_{me}','\rho_{a}','\rho_{l}','\rho_{qk}','\rho_{y}','\rho_{bh}',...
              '\rho_{be}','\rho_{d}','\rho_{G}','\rho_{\delta^{d}}','\rho_{kb}','\rho_{ib}');
names5 =char('k_{p}','k_{bh}','k_{be}','k_{d}','k_{i}','k_{w}','k_{kb}', 'h^{p}','h^{I}','h^{E}','\Theta','\iota_{w}','\iota_{p}','\phi_{R}','\phi_{\pi}','\phi_{y}','\rho^{fp}');



%%%POSTERIOR DENSITY%%%%%%%%%

e_z  = oo_.posterior_density.shocks_std.e_z;
e_h  = oo_.posterior_density.shocks_std.e_h;
e_mi = oo_.posterior_density.shocks_std.e_mi;
e_me = oo_.posterior_density.shocks_std.e_me;
e_ae = oo_.posterior_density.shocks_std.e_ae;
e_l = oo_.posterior_density.shocks_std.e_l;
e_qk = oo_.posterior_density.shocks_std.e_qk;
e_y = oo_.posterior_density.shocks_std.e_y;
e_bh = oo_.posterior_density.shocks_std.e_bh;
e_be = oo_.posterior_density.shocks_std.e_be;
e_d = oo_.posterior_density.shocks_std.e_d;
e_r_t = oo_.posterior_density.shocks_std.e_r_t;
e_G = oo_.posterior_density.shocks_std.e_G;
e_delta_d = oo_.posterior_density.shocks_std.e_delta_d;
e_kb = oo_.posterior_density.shocks_std.e_kb;
e_ib = oo_.posterior_density.shocks_std.e_ib;

rho_e_z = oo_.posterior_density.parameters.rho_e_z;
rho_e_h = oo_.posterior_density.parameters.rho_e_h;
rho_e_mi = oo_.posterior_density.parameters.rho_mi;
rho_e_me = oo_.posterior_density.parameters.rho_me;
rho_e_ae = oo_.posterior_density.parameters.rho_e_ae;
rho_e_l = oo_.posterior_density.parameters.rho_e_l;
rho_e_qk = oo_.posterior_density.parameters.rho_e_qk;
rho_e_y = oo_.posterior_density.parameters.rho_e_y;
rho_e_bh = oo_.posterior_density.parameters.rho_e_bh;
rho_e_be = oo_.posterior_density.parameters.rho_e_be;
rho_e_d = oo_.posterior_density.parameters.rho_e_d;
rho_e_G = oo_.posterior_density.parameters.rho_G;
rho_e_delta_d = oo_.posterior_density.parameters.rho_delta_d;
rho_kb = oo_.posterior_density.parameters.rho_kb;
rho_ib = oo_.posterior_density.parameters.rho_ib;

k_p = oo_.posterior_density.parameters.k_p;
k_bh = oo_.posterior_density.parameters.k_bh;
k_be = oo_.posterior_density.parameters.k_be;
k_d = oo_.posterior_density.parameters.k_d;
k_i = oo_.posterior_density.parameters.k_i;
k_w = oo_.posterior_density.parameters.k_w;
k_kb = oo_.posterior_density.parameters.k_kb;
hab_p = oo_.posterior_density.parameters.hab_p;
hab_i = oo_.posterior_density.parameters.hab_i;
hab_e = oo_.posterior_density.parameters.hab_e;
chi_sb = oo_.posterior_density.parameters.Dis_sb;
iota_w = oo_.posterior_density.parameters.iota_w;
iota_p = oo_.posterior_density.parameters.iota_p;
phi_R = oo_.posterior_density.parameters.phi_R;
phi_pi = oo_.posterior_density.parameters.phi_pi;
phi_y = oo_.posterior_density.parameters.phi_y;
rho_fp = oo_.posterior_density.parameters.rho_fp;





 Post1_stderr =[e_z(:,2) e_h(:,2)  e_mi(:,2) e_me(:,2)  e_ae(:,2) e_l(:,2) e_qk(:,2) e_y(:,2) e_bh(:,2) e_be(:,2) e_d(:,2) e_r_t(:,2) e_G(:,2) e_delta_d(:,2)  e_kb(:,2) e_ib(:,2)];
 Post1_param  =[rho_e_z(:,2) rho_e_h(:,2) rho_e_mi(:,2) rho_e_me(:,2) rho_e_ae(:,2) rho_e_l(:,2) rho_e_qk(:,2) rho_e_y(:,2) rho_e_bh(:,2) rho_e_be(:,2) rho_e_d(:,2) rho_e_G(:,2) rho_e_delta_d(:,2) rho_kb(:,2) rho_ib(:,2)];  
 Post3_param  = [k_p(:,2) k_bh(:,2) k_be(:,2) k_d(:,2)  k_i(:,2) k_w(:,2) k_kb(:,2) hab_p(:,2) hab_i(:,2)  hab_e(:,2) chi_sb(:,2) iota_w(:,2) iota_p(:,2) phi_R(:,2) phi_pi(:,2) phi_y(:,2) rho_fp(:,2)];


  
Post1_stderr_support =[e_z(:,1) e_h(:,1)  e_mi(:,1) e_me(:,1)  e_z(:,1) e_ae(:,1) e_l(:,1) e_qk(:,1) e_y(:,1) e_bh(:,1) e_be(:,1) e_d(:,1) e_r_t(:,1) e_G(:,1) e_delta_d(:,1)];
Post1_param_support  =[rho_e_z(:,1) rho_e_h(:,1) rho_e_mi(:,1) rho_e_me(:,1) rho_e_ae(:,1) rho_e_l(:,1) rho_e_qk(:,1) rho_e_y(:,1) rho_e_bh(:,1) rho_e_be(:,1) rho_e_d(:,1) rho_e_G(:,1) rho_e_delta_d(:,1) rho_kb(:,1) rho_ib(:,1) ];
Post3_param_support  =[ k_p(:,1) k_bh(:,1) k_be(:,1) k_d(:,1)  k_i(:,1) k_w(:,1)  k_kb(:,1) hab_p(:,1) hab_i(:,1)  hab_e(:,1) chi_sb(:,1) iota_w(:,1) iota_p(:,1) phi_R(:,1) phi_pi(:,1) phi_y(:,1) rho_fp(:,1)];



%%%%%%%%PRIOR DENSITY%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

e_z1  = oo_.prior_density.shocks_std.e_z;
e_h1  = oo_.prior_density.shocks_std.e_h;
e_mi1 = oo_.prior_density.shocks_std.e_mi;
e_me1 = oo_.prior_density.shocks_std.e_me;
e_ae1 = oo_.prior_density.shocks_std.e_ae;
e_l1 = oo_.prior_density.shocks_std.e_l;
e_qk1 = oo_.prior_density.shocks_std.e_qk;
e_y1 = oo_.prior_density.shocks_std.e_y;
e_bh1 = oo_.prior_density.shocks_std.e_bh;
e_be1 = oo_.prior_density.shocks_std.e_be;
e_d1 = oo_.prior_density.shocks_std.e_d;
e_r_t1 = oo_.prior_density.shocks_std.e_r_t;
e_G1 = oo_.prior_density.shocks_std.e_G;
e_delta_d1 = oo_.prior_density.shocks_std.e_delta_d;
e_kb1 = oo_.prior_density.shocks_std.e_kb;
e_ib1 = oo_.prior_density.shocks_std.e_ib;

rho_e_z1 = oo_.prior_density.parameters.rho_e_z;
rho_e_h1 = oo_.prior_density.parameters.rho_e_h;
rho_e_mi1 = oo_.prior_density.parameters.rho_mi;
rho_e_me1 = oo_.prior_density.parameters.rho_me;
rho_e_ae1 = oo_.prior_density.parameters.rho_e_ae;
rho_e_l1 = oo_.prior_density.parameters.rho_e_l;
rho_e_qk1 = oo_.prior_density.parameters.rho_e_qk;
rho_e_y1 = oo_.prior_density.parameters.rho_e_y;
rho_e_bh1 = oo_.prior_density.parameters.rho_e_bh;
rho_e_be1 = oo_.prior_density.parameters.rho_e_be;
rho_e_d1 = oo_.prior_density.parameters.rho_e_d;
rho_e_G1 = oo_.prior_density.parameters.rho_G;
rho_e_delta_d1 = oo_.prior_density.parameters.rho_delta_d;
rho_kb1 = oo_.prior_density.parameters.rho_kb;
rho_ib1 = oo_.prior_density.parameters.rho_ib;

k_p1 = oo_.prior_density.parameters.k_p;
k_bh1 = oo_.prior_density.parameters.k_bh;
k_be1 = oo_.prior_density.parameters.k_be;
k_d1 = oo_.prior_density.parameters.k_d;
k_i1 = oo_.prior_density.parameters.k_i;
k_w1 = oo_.prior_density.parameters.k_w;
k_kb1 = oo_.prior_density.parameters.k_kb;
hab_p1 = oo_.prior_density.parameters.hab_p;
hab_i1 = oo_.prior_density.parameters.hab_i;
hab_e1 = oo_.prior_density.parameters.hab_e;
chi_sb1 = oo_.prior_density.parameters.Dis_sb;
iota_w1 = oo_.prior_density.parameters.iota_w;
iota_p1 = oo_.prior_density.parameters.iota_p;
phi_R1 = oo_.prior_density.parameters.phi_R;
phi_pi1 = oo_.prior_density.parameters.phi_pi;
phi_y1 = oo_.prior_density.parameters.phi_y;
rho_fp1 = oo_.prior_density.parameters.rho_fp;





 Prior1_stderr =[e_z1(:,2) e_h1(:,2)  e_mi1(:,2) e_me1(:,2)  e_ae1(:,2) e_l1(:,2) e_qk1(:,2) e_y1(:,2) e_bh1(:,2) e_be1(:,2) e_d1(:,2) e_r_t1(:,2) e_G1(:,2) e_delta_d1(:,2)  e_kb1(:,2) e_ib1(:,2)];
 Prior1_param  =[rho_e_z1(:,2) rho_e_h1(:,2) rho_e_mi1(:,2) rho_e_me1(:,2) rho_e_ae1(:,2) rho_e_l1(:,2) rho_e_qk1(:,2) rho_e_y1(:,2) rho_e_bh1(:,2) rho_e_be1(:,2) rho_e_d1(:,2) rho_e_G1(:,2) rho_e_delta_d1(:,2) rho_kb1(:,2) rho_ib1(:,2) ];
 Prior3_param  =[k_p1(:,2) k_bh1(:,2) k_be1(:,2) k_d1(:,2)  k_i1(:,2) k_w1(:,2) k_kb1(:,2) hab_p1(:,2) hab_i1(:,2) hab_e1(:,2) chi_sb1(:,2) iota_w1(:,2) iota_p1(:,2) phi_R1(:,2) phi_pi1(:,2) phi_y1(:,2) rho_fp1(:,2)];


  
Prior1_stderr_support =[e_z1(:,1) e_h1(:,1)  e_mi1(:,1) e_me1(:,1)  e_z1(:,1) e_ae1(:,1) e_l1(:,1) e_qk1(:,1) e_y1(:,1) e_bh1(:,1) e_be1(:,1) e_d1(:,1) e_r_t1(:,1) e_G1(:,1) e_delta_d1(:,1)  e_kb1(:,1) e_ib1(:,1)];
Prior1_param_support  =[rho_e_z1(:,1) rho_e_h1(:,1) rho_e_mi1(:,1) rho_e_me1(:,1) rho_e_ae1(:,1) rho_e_l1(:,1) rho_e_qk1(:,1) rho_e_y1(:,1) rho_e_bh1(:,1)  rho_e_be1(:,1) rho_e_d1(:,1) rho_e_G1(:,1) rho_e_delta_d1(:,1) rho_kb1(:,1) rho_ib1(:,1)]; 
Prior3_param_support  =[k_p1(:,1) k_bh1(:,1) k_be1(:,1) k_d1(:,1)  k_i1(:,1) k_w1(:,1)  k_kb1(:,1) hab_p1(:,1) hab_i1(:,1)   hab_e1(:,1) chi_sb1(:,1) iota_w1(:,1) iota_p1(:,1) phi_R1(:,1) phi_pi1(:,1) phi_y1(:,1) rho_fp1(:,1)];






h=figure('Position', [600, 0, 1000, 900]);
axes ('position', [0, 0, 1, 1]);


figure(1)
for j = 1:16;
    subplot(4,4,j), plot(Post1_stderr_support(:,j),Post1_stderr(:,j),'k', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
                    plot(Prior1_stderr_support(:,j),Prior1_stderr(:,j),'-.k', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','r'); hold on;
                    %xlabel(Post1_stderr_support(:,j));
                    %ylabel('% dev from SS');
                    grid on
                    title(names1(j,:),'FontSize',10)
axis tight;                   
end;
legend('Posterior', 'Prior',  'Location', 'Best', 'Orientation', 'horizontal')

% 
% h=figure('Position', [600, 0, 1000, 900]);
% axes ('position', [0, 0, 1, 1]);
% 
% 
% figure(2)
% for j = 1:6;
%     subplot(3,2,j), plot(Post2_stderr_support(:,j),Post2_stderr(:,j),'b', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
%                     plot(Prior2_stderr_support(:,j),Prior2_stderr(:,j),'--r', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
%                     %xlabel(Post1_stderr_support(:,j));
%                     %ylabel('% dev from SS');
%                     grid on
%                     title(names2(j,:),'FontSize',10)
% axis tight;                   
% end;
% legend('Posterior', 'Prior',  'Location', 'Best', 'Orientation', 'horizontal')
% 
% h=figure('Position', [600, 0, 1000, 900]);
% axes ('position', [0, 0, 1, 1]);
% 
% 
% figure(3)
% for j = 1:10;
%     subplot(5,2,j), plot(Post1_param_support(:,j),Post1_param(:,j),'b', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
%                     plot(Prior1_param_support(:,j),Prior1_param(:,j),'--r', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','r'); hold on;
%                     %xlabel(Post1_stderr_support(:,j));
%                     %ylabel('% dev from SS');
%                     grid on
%                     title(names3(j,:),'FontSize',10)
% axis tight;                   
% end;
% 
% legend('Posterior', 'Prior',  'Location', 'Best', 'Orientation', 'horizontal')
% h=figure('Position', [600, 0, 1000, 900]);
% axes ('position', [0, 0, 1, 1]);
% % 
% 
% figure(4)
% for j = 1:4;
%     subplot(2,2,j), plot(Post2_param_support(:,j),Post2_param(:,j),'b', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
%                     plot(Prior2_param_support(:,j),Prior2_param(:,j),'--r', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','r'); hold on;
%                     %xlabel(Post1_stderr_support(:,j));
%                     %ylabel('% dev from SS');
%                     grid on
%                     title(names4(j,:),'FontSize',10)
% axis tight;                   
% end;
% legend('Posterior', 'Prior',  'Location', 'Best', 'Orientation', 'horizontal')
% 
% h=figure('Position', [600, 0, 1000, 900]);
% axes ('position', [0, 0, 1, 1]);
% 
% 
% figure(5)
% for j = 1:8;
%     subplot(4,2,j), plot(Post3_param_support(:,j),Post3_param(:,j),'b', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
%                     plot(Prior3_param_support(:,j),Prior3_param(:,j),'--r', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','r'); hold on;
%                     %xlabel(Post1_stderr_support(:,j));
%                     %ylabel('% dev from SS');
%                     grid on
%                     title(names5(j,:),'FontSize',10)
% axis tight;                   
% end;
% legend('Posterior', 'Prior',  'Location', 'Best', 'Orientation', 'horizontal')
% 
% h=figure('Position', [600, 0, 1000, 900]);
% axes ('position', [0, 0, 1, 1]);


% figure(1)
% for j = 1:17;
%     subplot(5,4,j), plot(Post3_param_support(:,j),Post3_param(:,j),'k', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','b'); hold on;
%                     plot(Prior3_param_support(:,j),Prior3_param(:,j),'-.k', 'LineWidth',2, 'MarkerSize', 5,'MarkerEdgeColor','k', 'MarkerFaceColor','r'); hold on;
%                  
%                     %xlabel(Post1_stderr_support(:,j));
%                     %ylabel('% dev from SS');
%                     grid on
%                     title(names5(j,:),'FontSize',10)
% axis tight;                   
% end;
% %legend('Posterior', 'Prior',  'Location', 'Best', 'Orientation', 'horizontal')
% 
