diff --git a/matlab/PosteriorIRF.m b/matlab/PosteriorIRF.m
index fd37ff48d..8c094a24f 100644
--- a/matlab/PosteriorIRF.m
+++ b/matlab/PosteriorIRF.m
@@ -398,7 +398,7 @@ if options_.TeX && any(strcmp('eps',cellstr(options_.graph_format)))
         figunumber = 0;
 
         for jj=1:nvar
-            if max(abs(MeanIRF(:,jj,ii))) > options_.impulse_responses.plot_threshold
+            if max(abs(MeanIRF(:,jj,ii))) >= options_.impulse_responses.plot_threshold
                 subplotnum = subplotnum+1;
                 
                 if subplotnum == 1 
diff --git a/matlab/PosteriorIRF_core2.m b/matlab/PosteriorIRF_core2.m
index 1fef8cdff..51f4bc08a 100644
--- a/matlab/PosteriorIRF_core2.m
+++ b/matlab/PosteriorIRF_core2.m
@@ -100,7 +100,7 @@ for i=fpar:npar,
     figunumber = 0;
 
     for j=1:nvar
-        if max(abs(MeanIRF(:,j,i))) > options_.impulse_responses.plot_threshold
+        if max(abs(MeanIRF(:,j,i))) >= options_.impulse_responses.plot_threshold
             subplotnum = subplotnum+1;
             if subplotnum == 1 && options_.relative_irf
                 hh = dyn_figure(options_,'Name',['Relative response to orthogonalized shock to ' tit(i,:)]);
diff --git a/matlab/stoch_simul.m b/matlab/stoch_simul.m
index 03f3acc47..d6b6ceadf 100644
--- a/matlab/stoch_simul.m
+++ b/matlab/stoch_simul.m
@@ -215,7 +215,7 @@ if options_.irf
                          y(i_var(j),:)');
                 eval(['oo_.irfs.' deblank(M_.endo_names(i_var(j),:)) '_' ...
                       deblank(M_.exo_names(i,:)) ' = y(i_var(j),:);']);
-                if max(abs(y(i_var(j),:))) > options_.impulse_responses.plot_threshold
+                if max(abs(y(i_var(j),:))) >= options_.impulse_responses.plot_threshold
                     irfs  = cat(1,irfs,y(i_var(j),:));
                     if isempty(mylist)
                         mylist = deblank(var_list(j,:));