diff --git a/matlab/PosteriorIRF.m b/matlab/PosteriorIRF.m
index c78733e383101727c1ef345a16edc1ca3bdfcaa3..e74c19a49c2fb07c62d55a48e80b02571ce8e596 100644
--- a/matlab/PosteriorIRF.m
+++ b/matlab/PosteriorIRF.m
@@ -311,7 +311,7 @@ for i = 1:M_.exo_nbr
         eval(['oo_.PosteriorIRF.dsge.Mean.' name ' = MeanIRF(:,j,i);']);
         eval(['oo_.PosteriorIRF.dsge.Median.' name ' = MedianIRF(:,j,i);']);
         eval(['oo_.PosteriorIRF.dsge.Var.' name ' = VarIRF(:,j,i);']);
-        eval(['oo_.PosteriorIRF.dsge.Distribution.' name ' = DistribIRF(:,:,j,i);']);
+        eval(['oo_.PosteriorIRF.dsge.deciles.' name ' = DistribIRF(:,:,j,i);']);
         eval(['oo_.PosteriorIRF.dsge.HPDinf.' name ' = HPDIRF(:,1,j,i);']);
         eval(['oo_.PosteriorIRF.dsge.HPDsup.' name ' = HPDIRF(:,2,j,i);']);
     end
@@ -348,7 +348,7 @@ if MAX_nirfs_dsgevar
             eval(['oo_.PosteriorIRF.bvardsge.Mean.' name ' = MeanIRFdsgevar(:,j,i);']);
             eval(['oo_.PosteriorIRF.bvardsge.Median.' name ' = MedianIRFdsgevar(:,j,i);']);
             eval(['oo_.PosteriorIRF.bvardsge.Var.' name ' = VarIRFdsgevar(:,j,i);']);
-            eval(['oo_.PosteriorIRF.bvardsge.Distribution.' name ' = DistribIRFdsgevar(:,:,j,i);']);
+            eval(['oo_.PosteriorIRF.bvardsge.deciles.' name ' = DistribIRFdsgevar(:,:,j,i);']);
             eval(['oo_.PosteriorIRF.bvardsge.HPDinf.' name ' = HPDIRFdsgevar(:,1,j,i);']);
             eval(['oo_.PosteriorIRF.bvardsge.HPDsup.' name ' = HPDIRFdsgevar(:,2,j,i);']);
         end
diff --git a/matlab/dyn_saveas.m b/matlab/dyn_saveas.m
index b691d77727a8dd92189604140eba3a455f986b20..3b00f0aeef634ba2d314bbe5cffd87dda863d9dd 100644
--- a/matlab/dyn_saveas.m
+++ b/matlab/dyn_saveas.m
@@ -32,16 +32,16 @@ function dyn_saveas(h,fname,DynareOptions)
 
 if any(strcmp('eps',cellstr(DynareOptions.graph_format)))
     if exist('OCTAVE_VERSION')
-        eval(['print -depsc2 ' fname '.eps']);
+        eval(['print -depsc2 ' fname '.eps']); % still need to be fixed to take handle into account
     else
-        eval(['print -depsc2 ' fname]);
+        print(h,'-depsc2',[fname,'.eps']) 
     end
 end
 if any(strcmp('pdf',cellstr(DynareOptions.graph_format)))
     if exist('OCTAVE_VERSION')
         warning('Octave cannot create pdf files!')
     else
-        eval(['print -dpdf ' fname]);
+        print(h,'-dpdf',[fname,'.pdf']) 
     end
 end
 if any(strcmp('fig',cellstr(DynareOptions.graph_format)))
diff --git a/matlab/plot_identification.m b/matlab/plot_identification.m
index d439eb6a6877d61c0bb00967d6b0bbd30f6b75b8..4cc2d3a88510ace72037398c7ed66b067265119c 100644
--- a/matlab/plot_identification.m
+++ b/matlab/plot_identification.m
@@ -108,6 +108,7 @@ if SampleSize == 1,
         for ip=1:nparam,
             text(ip,dy(1),name{is(ip)},'rotation',90,'HorizontalAlignment','right','interpreter','none')
         end
+        legend('Moments','Model','LRE model','Location','Best')
         title('Sensitivity bars using derivatives (log-scale)')
         dyn_saveas(hh,[IdentifDirectoryName '/' M_.fname '_sensitivity_' tittxt1 ],options_);
         
diff --git a/matlab/pm3.m b/matlab/pm3.m
index dbc81e7c7ec47158f688d4a41ecd2b1cbe24faba..91ec283ed64ce53af2a8076a0e7ef4ca4591280f 100644
--- a/matlab/pm3.m
+++ b/matlab/pm3.m
@@ -84,7 +84,7 @@ for i = 1:nvar
     eval(['oo_.' name3 '.Mean.' name ' = Mean(:,i);']);
     eval(['oo_.' name3 '.Median.' name ' = Median(:,i);']);
     eval(['oo_.' name3 '.Var.' name ' = Var(:,i);']);
-    eval(['oo_.' name3 '.Distribution.' name ' = Distrib(:,:,i);']);
+    eval(['oo_.' name3 '.deciles.' name ' = Distrib(:,:,i);']);
     eval(['oo_.' name3 '.HPDinf.' name ' = HPD(1,:,i);']);
     eval(['oo_.' name3 '.HPDsup.' name ' = HPD(2,:,i);']);
 end