diff --git a/matlab/dynare_identification.m b/matlab/dynare_identification.m
index 89f462ff92c4969fb5322c954429d82b7f20266f..8873fd307902e248a37e3b7434c6de85950ca695 100644
--- a/matlab/dynare_identification.m
+++ b/matlab/dynare_identification.m
@@ -156,7 +156,7 @@ if iload <=0,
     if SampleSize==1,
         BurninSampleSize=0;
     else
-        BurninSampleSize=50;
+        BurninSampleSize=min(SampleSize,50);
     end
     loop_indx = 0;
     file_index = 0;
@@ -350,6 +350,7 @@ if iload <=0,
                     else
                         normaliz1 = ones(1,nparam);
                     end
+                    replic=max(replic,3*length(indJJ));
                     cmm = simulated_moment_uncertainty(indJJ, periods, replic);
                     %                 Jinv=(siJ(:,indok)'*siJ(:,indok))\siJ(:,indok)';
                     %                 MIM=inv(Jinv*cmm*Jinv');