Merge pull request #1027 from JohannesPfeifer/irf_plot_threshold

Make irf_plot_threshold consistent with manual

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 

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,:)]);

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,:));