Merge pull request #705 from JohannesPfeifer/get_posterior_bug

Fixes bug where number of posterior draws was computed incorrectly and w...
time-shift
Sébastien Villemot 2014-08-14 12:44:58 +02:00
commit 07d321c159
1 changed files with 1 additions and 1 deletions

View File

@ -75,7 +75,7 @@ catch
[marginal,oo_] = marginal_density(M_, options_, estim_params_, oo_);
disp(sprintf('Log data density is %f.',oo_.MarginalDensity.ModifiedHarmonicMean))
end
num_draws=length(NumberOfDraws*options_.mh_nblck);
num_draws=NumberOfDraws*options_.mh_nblck;
hpd_draws = round((1-options_.mh_conf_sig)*num_draws);
if hpd_draws<2
fprintf('posterior_moments: There are not enough draws computes to compute HPD Intervals. Skipping their computation.\n')