From 8c24b5210cb680f64464784b9762fe902f15f670 Mon Sep 17 00:00:00 2001
From: Houtan Bastani <houtan@dynare.org>
Date: Wed, 12 Jul 2017 18:12:26 +0200
Subject: [PATCH] sur gibbs: add plots, save parameter values

---
 matlab/surgibbs.m | 17 ++++++++++++++++-
 1 file changed, 16 insertions(+), 1 deletion(-)

diff --git a/matlab/surgibbs.m b/matlab/surgibbs.m
index bd3f62abc..3c7721ff9 100644
--- a/matlab/surgibbs.m
+++ b/matlab/surgibbs.m
@@ -131,4 +131,19 @@ for i = 1:ndraws
     Sigma_upper_chol = chol(Omegabar, 'upper');
     beta = rand_multivariate_normal(betabar', Sigma_upper_chol, nvars)';
     oo_.surgibbs.betadraws(i, :) = beta';
-end
\ No newline at end of file
+end
+
+% save parameter values
+oo_.surgibbs.beta = (sum(oo_.surgibbs.betadraws)/ndraws)';
+figure
+nrows = 5;
+ncols = floor(nvars/nrows);
+if mod(nvars, nrows) ~= 0
+    ncols = ncols + 1;
+end
+for j = 1:length(pnamesall)
+    M_.params(strmatch(pnamesall{j}, M_.param_names, 'exact')) = oo_.surgibbs.beta(j);
+    subplot(nrows, ncols, j)
+    histogram(oo_.surgibbs.betadraws(:, j))
+    title(pnamesall{j}, 'Interpreter', 'none')
+end