From 08baf6dca5a56a9a1c98cf4ce35ad3bff1a55a97 Mon Sep 17 00:00:00 2001
From: michel <michel@ac1d8469-bf42-47a9-8791-bf33cf982152>
Date: Wed, 22 Apr 2009 21:34:45 +0000
Subject: [PATCH] v4.1: fixing bugs in table for posterior mean
git-svn-id: https://www.dynare.org/svn/dynare/trunk@2615 ac1d8469-bf42-47a9-8791-bf33cf982152
---
matlab/GetPosteriorParametersStatistics.m | 11 +++++------
matlab/row_header_width.m | 6 +++---
2 files changed, 8 insertions(+), 9 deletions(-)
diff --git a/matlab/GetPosteriorParametersStatistics.m b/matlab/GetPosteriorParametersStatistics.m
index eeca2d1d7b..7ac1777081 100644
--- a/matlab/GetPosteriorParametersStatistics.m
+++ b/matlab/GetPosteriorParametersStatistics.m
@@ -58,17 +58,17 @@ NumberOfDraws = TotalNumberOfMhDraws-floor(options_.mh_drop*TotalNumberOfMhDraws
clear record;
pnames=[' ';'beta ';'gamm ';'norm ';'invg ';'unif ';'invg2'];
-tit2 = sprintf('%10s %7s %10s %14s %4s %6s\n',' ','prior mean','post. mean','conf. interval','prior','pstdev');
-pformat = '%-*s %7.3f %8.4f %7.4f %7.4f %4s %6.4f';
+header_width = row_header_width(M_,estim_params_,bayestopt_);
+tit2 = sprintf('%-*s %10s %10s %16s %6s %10s\n',header_width+2,' ','prior mean','post. mean','conf. interval','prior','pstdev');
+pformat = '%-*s %10.3f %10.4f %10.4f %8.4f %6s %10.4f';
disp(' ');disp(' ');disp('ESTIMATION RESULTS');disp(' ');
try
disp(sprintf('Log data density is %f.',oo_.MarginalDensity.ModifiedHarmonicMean))
catch
- [marginal,oo_] = marginal_density(M_, options_, estim_params_, oo_)
+ [marginal,oo_] = marginal_density(M_, options_, estim_params_, oo_);
disp(sprintf('Log data density is %f.',oo_.MarginalDensity.ModifiedHarmonicMean))
end
-header_width = row_header_width(M_,estim_params,bayestopt_);
if np
type = 'parameters';
if TeX
@@ -97,8 +97,7 @@ if np
oo_ = Filloo(oo_,name,type,post_mean,hpd_interval,post_median,post_var,post_deciles,density);
end
end
- disp(sprintf(pformat,name,bayestopt_.p1(ip),...
- header_width, ...
+ disp(sprintf(pformat,header_width,name,bayestopt_.p1(ip),...
post_mean, ...
hpd_interval, ...
pnames(bayestopt_.pshape(ip)+1,:), ...
diff --git a/matlab/row_header_width.m b/matlab/row_header_width.m
index ef3762f110..b013007e5b 100644
--- a/matlab/row_header_width.m
+++ b/matlab/row_header_width.m
@@ -44,13 +44,13 @@ function w=row_header_width(M_,estim_params_,bayestopt_)
end
if nvx
for i=1:nvx
- k = estim_params_.var_exo(i,1)
+ k = estim_params_.var_exo(i,1);
w = max(w,length(deblank(M_.exo_names(k,:))));
end
end
if nvn
for i=1:nvn
- k = estim_params_.var_endo(i,1)
+ k = estim_params_.var_endo(i,1);
w = max(w,length(deblank(M_.endo_names(k,:))));
end
end
@@ -68,7 +68,7 @@ function w=row_header_width(M_,estim_params_,bayestopt_)
k1 = estim_params_.corrn(i,1);
k2 = estim_params_.corrn(i,2);
w = max(w,length(deblank(M_.endo_names(k1,:)))...
- +length(deblank(M_.endo_names(k2,:))))
+ +length(deblank(M_.endo_names(k2,:))));
end
end
--
GitLab