diff --git a/matlab/prior_posterior_statistics.m b/matlab/prior_posterior_statistics.m
index d63c96c6a65ee7f3d33fb3504738b658a3756dd3..e5d65ef9c762852403928f983402528d41b8562c 100644
--- a/matlab/prior_posterior_statistics.m
+++ b/matlab/prior_posterior_statistics.m
@@ -123,7 +123,7 @@ end
 
 varlist = options_.varlist;
 if isempty(varlist)
-    varlist = M_.endo_names(1:M_.orig_endo_nbr, :);
+    varlist = char(sort(cellstr(M_.endo_names(1:M_.orig_endo_nbr,:))));
 end
 nvar = size(varlist,1);
 SelecVariables = [];
@@ -300,19 +300,19 @@ end
 
 if options_.smoother
     pm3(endo_nbr,gend,ifil(1),B,'Smoothed variables',...
-        '',M_.endo_names(1:M_.orig_endo_nbr, :),M_.endo_names_tex,M_.endo_names,...
+        '',varlist,M_.endo_names_tex,M_.endo_names,...
         varlist,'SmoothedVariables',DirectoryName,'_smooth');
     pm3(exo_nbr,gend,ifil(2),B,'Smoothed shocks',...
         '',M_.exo_names,M_.exo_names_tex,M_.exo_names,...
         M_.exo_names,'SmoothedShocks',DirectoryName,'_inno');
     pm3(endo_nbr,1,ifil(9),B,'Trend_coefficients',...
-    '',M_.endo_names(1:M_.orig_endo_nbr,:),M_.endo_names_tex,M_.endo_names,...
+    '',varlist,M_.endo_names_tex,M_.endo_names,...
         varlist,'TrendCoeff',DirectoryName,'_trend_coeff');
     pm3(endo_nbr,gend,ifil(10),B,'Smoothed constant',...
-        '',M_.endo_names(1:M_.orig_endo_nbr, :),M_.endo_names_tex,M_.endo_names,...
+        '',varlist,M_.endo_names_tex,M_.endo_names,...
         varlist,'Constant',DirectoryName,'_smoothed_constant');
     pm3(endo_nbr,gend,ifil(11),B,'Smoothed trend',...
-        '',M_.endo_names(1:M_.orig_endo_nbr, :),M_.endo_names_tex,M_.endo_names,...
+        '',varlist,M_.endo_names_tex,M_.endo_names,...
         varlist,'Trend',DirectoryName,'_smoothed_trend');
     pm3(endo_nbr,gend,ifil(1),B,'Updated Variables',...
         '',varlist,M_.endo_names_tex,M_.endo_names,...