stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed timing bugs in annualized shock decomposition

time-shift
Marco Ratto 2017-02-13 10:12:38 +01:00 committed by Stéphane Adjemian (Charybdis)
parent d7723877cf
commit 1920fb0bfe
2 changed files with 8 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/annualized_shock_decomposition.m
View File

@ -135,10 +135,10 @@ for j=1:nvar
end
for k =1:nterms,
if isstruct(aux),
aux.y = squeeze(yaux(j,k,t0:end));
aux.y = squeeze(yaux(j,k,min((t0-3):-4:1):end));
end
[za(j,k,:), steady_state_a(j,1), gza(j,k,:), steady_state_ga(j,1)] = ...
quarterly2annual(squeeze(z(j,k,t0:end)),steady_state(j),GYTREND0,var_type,islog,aux);
quarterly2annual(squeeze(z(j,k,min((t0-3):-4:1):end)),steady_state(j),GYTREND0,var_type,islog,aux);
end
ztmp=squeeze(za(j,:,:));
if cumfix==0,
@ -178,7 +178,7 @@ end
% realtime
if realtime_ && isstruct(oo_) && isfield(oo_, 'realtime_shock_decomposition'),
init=1;
for i=t0+4:4:t1,
for i=t0:4:t1,
yr=floor(i/4);
za=[];
gza=[];
@ -206,10 +206,10 @@ end
% if k<nterms
% ztmp = squeeze(sum(z(j,[1:k-1,k+1:end-1],t0-4:end)));
% else
ztmp = squeeze(z(j,k,min(t0:-4:1):end));
ztmp = squeeze(z(j,k,min((t0-3):-4:1):end));
% end
if isstruct(aux),
aux.y = squeeze(yaux(j,k,min(t0:-4:1):end));
aux.y = squeeze(yaux(j,k,min((t0-3):-4:1):end));
end
[za(j,k,:), steady_state_a(j,1), gza(j,k,:), steady_state_ga(j,1)] = ...
quarterly2annual(ztmp,steady_state(j),GYTREND0,var_type,islog,aux);

4
matlab/plot_shock_decomposition.m
View File

@ -130,6 +130,8 @@ if isequal(type,'aoa') && isstruct(q2a) && realtime_
end
end
t0=min(options_.shock_decomp.save_realtime);
ini1 = initial_date+t0-1;
t0=t0+(4-ini1.time(2));
if ~isfield(q2a,'var_type'), % private trap for aoa calls
q2a.var_type=1;
end
@ -225,7 +227,7 @@ switch type
if ~isempty(initial_date),
initial_date = initial_date+3;
else
initial_date = dates('0Q4');
initial_date = dates('1Q4');
end
steady_state = 4*steady_state;