diff --git a/matlab/metropolis_draw.m b/matlab/metropolis_draw.m
index 722feba9569cb1adc72d1da0141790ce5efab463..ad396b30b75becddf9b4204cfa4d7702e8055108 100644
--- a/matlab/metropolis_draw.m
+++ b/matlab/metropolis_draw.m
@@ -70,7 +70,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, round(NumberOfDraws*mh_nblck));
+        options_.sub_draws = min(options_.posterior_max_subsample_draws, ceil(.25*NumberOfDraws));
     else
         if options_.sub_draws>NumberOfDraws*mh_nblck
             skipline()