diff --git a/matlab/estimation/metropolis_draw.m b/matlab/estimation/metropolis_draw.m
index 505d1dbeb6b5deaa7964a366b882f5fabf2cd2a7..7c41680989757ab6f9714748eac43ea152963b9c 100644
--- a/matlab/estimation/metropolis_draw.m
+++ b/matlab/estimation/metropolis_draw.m
@@ -68,7 +68,7 @@ if init
     mh_nblck = options_.mh_nblck;
     % set sub_draws option if empty
     if isempty(options_.sub_draws)
-        options_.sub_draws = min(options_.posterior_max_subsample_draws, ceil(.25*NumberOfDraws));
+        options_.sub_draws = min(options_.posterior_max_subsample_draws, ceil(NumberOfDraws*mh_nblck));
     else
         if options_.sub_draws>NumberOfDraws*mh_nblck
             skipline()