diff --git a/matlab/McMCDiagnostics.m b/matlab/McMCDiagnostics.m
index 5763d8ed1353cc85a21aaafe45a381404b528e5a..997422bc765a1083876c5736f515a08cd01e915f 100644
--- a/matlab/McMCDiagnostics.m
+++ b/matlab/McMCDiagnostics.m
@@ -91,6 +91,7 @@ localVars.nblck = nblck;
 localVars.NumberOfMcFilesPerBlock = NumberOfMcFilesPerBlock;
 localVars.Origin = Origin;
 localVars.StepSize = StepSize;
+localVars.mh_drop = options_.mh_drop;
 localVars.NumberOfDraws = NumberOfDraws;
 localVars.NumberOfLines = NumberOfLines;
 localVars.time = time;
diff --git a/matlab/McMCDiagnostics_core.m b/matlab/McMCDiagnostics_core.m
index 9729f2f1b15cef2de978f67f15861c7cbf1dfc86..c566646cf5f486d97ffd761c85e8335fb725fdb6 100644
--- a/matlab/McMCDiagnostics_core.m
+++ b/matlab/McMCDiagnostics_core.m
@@ -65,7 +65,7 @@ for j=fpar:npar,
     ligne   = 0;
     for iter  = Origin:StepSize:NumberOfDraws
         ligne = ligne+1;
-        linea = ceil(0.5*iter);
+        linea = ceil(mh_drop*iter);
         n     = iter-linea+1;
         cinf  = round(n*ALPHA/2);
         csup  = round(n*(1-ALPHA/2));