diff --git a/matlab/dynare_estimation_1.m b/matlab/dynare_estimation_1.m
index 3fe420a2b..c39d698e8 100644
--- a/matlab/dynare_estimation_1.m
+++ b/matlab/dynare_estimation_1.m
@@ -610,7 +610,7 @@ if (~((any(bayestopt_.pshape > 0) && options_.mh_replic) || (any(bayestopt_.psha
         number_of_plots_to_draw = 0;
         index = [];
         for obs_iter=1:n_varobs
-            if max(abs(measurement_error(obs_iter,:))) > 0.0000001
+            if max(abs(measurement_error(obs_iter,:))) > options_.ME_plot_tol;
                 number_of_plots_to_draw = number_of_plots_to_draw + 1;
                 index = cat(1,index,obs_iter);
             end
diff --git a/matlab/global_initialization.m b/matlab/global_initialization.m
index d26b488c8..762a23281 100644
--- a/matlab/global_initialization.m
+++ b/matlab/global_initialization.m
@@ -467,6 +467,7 @@ options_.prior_trunc = 1e-10;
 options_.smoother = 0;
 options_.posterior_max_subsample_draws = 1200;
 options_.sub_draws = [];
+options_.ME_plot_tol=1e-6;
 % options_.use_mh_covariance_matrix = 0;
 options_.gradient_method = 2; %used by csminwel and newrat
 options_.gradient_epsilon = 1e-6; %used by csminwel and newrat