From b0e7c47a000a679420970f1d05e4f2ed7d60764c Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@ec.europa.eu>
Date: Fri, 11 Jan 2019 23:33:45 +0100
Subject: [PATCH] fixed bug with non-increasing time declarations of moment
restrictions. Modified example to get more sensible graph of ACF.
---
matlab/gsa/map_calibration.m | 4 +++-
tests/gsa/ls2003a.mod | 2 +-
2 files changed, 4 insertions(+), 2 deletions(-)
diff --git a/matlab/gsa/map_calibration.m b/matlab/gsa/map_calibration.m
index 8f6aa859ab..30d394b354 100644
--- a/matlab/gsa/map_calibration.m
+++ b/matlab/gsa/map_calibration.m
@@ -394,6 +394,7 @@ if ~isempty(indx_moment)
% For single legend search which has maximum nbr of restrictions
maxijv=0;
for ij=1:nbr_moment_restrictions
+ endo_prior_restrictions.moment{ij,3} = sort(endo_prior_restrictions.moment{ij,3});
if length(endo_prior_restrictions.moment{ij,3})>maxijv
maxij=ij;maxijv=length(endo_prior_restrictions.moment{ij,3});
end
@@ -477,6 +478,7 @@ if ~isempty(indx_moment)
% end
end
for ij=1:nbr_moment_couples
+ time_matrix{ij} = sort(time_matrix{ij});
if length(time_matrix{ij})>1
if ~DynareOptions.nograph
itmp = (find(plot_indx==ij));
@@ -505,7 +507,7 @@ if ~isempty(indx_moment)
hold off
axis(a)
box on
- set(gca,'xtick',sort(time_matrix{ij}))
+% set(gca,'xtick',sort(time_matrix{ij}))
itmp = min(itmp);
title([endo_prior_restrictions.moment{itmp,1},' vs ',endo_prior_restrictions.moment{itmp,2}],'interpreter','none'),
end
diff --git a/tests/gsa/ls2003a.mod b/tests/gsa/ls2003a.mod
index 0bca629a71..f72476b0cf 100644
--- a/tests/gsa/ls2003a.mod
+++ b/tests/gsa/ls2003a.mod
@@ -88,7 +88,7 @@ end;
moment_calibration;
//y_obs,y_obs, [0.8 1.1]; //[unconditional variance]
-y_obs,y_obs(-(1:4)), +; //[first year acf]
+y_obs,y_obs(1:4), +; //[first year acf]
//y_obs,pie_obs(-4:4), -; //[ccf]
@#for ilag in -2:2
y_obs,R_obs(@{ilag}), -; //[ccf]
--
GitLab