stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

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
time-shift
michel 2009-04-22 21:34:45 +00:00
parent d55470ca88
commit 08baf6dca5
2 changed files with 8 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
matlab/GetPosteriorParametersStatistics.m
View File

@ -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,:), ...

6
matlab/row_header_width.m
View File

@ -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