I have a question, I would like to make a simulation where for 8 consecutive periods a shock is not anticipated, but after 8 periods, the shock must become anticipated, and assume an already anticipated value, from there up to 40 periods.
I understand, that in this post they explain something similar,
If this post is correct, or similar to what I want to do, then should I enter 8 consecutive shocks this way?
//initialize IRF generation
initial_condition_states = repmat(oo_.dr.ys,1,M_.maximum_lag);
shock_matrix = zeros(options_.irf,M_.exo_nbr); %create shock matrix with number of time periods in columns
// set shocks for anticipated shocks
shock_matrix(1,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
shock_matrix(2,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05; shock_matrix(3,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
shock_matrix(4,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
shock_matrix(5,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
shock_matrix(6,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
shock_matrix(7,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
shock_matrix(8,strmatch(‘ERR_P_S’,M_.exo_names,‘exact’)) = -.05;
y2 = simult_(M_,options_,initial_condition_states,oo_.dr,shock_matrix,1);
y_IRF = y2(:,M_.maximum_lag+1:end)-repmat(oo_.dr.ys,1,options_.irf); %deviation from steady state
After that, how do I introduce or indicate that the shock from there on should be anticipated?