Weird stoch_simul

[HelloDynare]

Hi Alisa,

l_pos=strmatch('l',var_list_ ,'exact');

would like to find the position of labor in the var_list_, after we run

stoch_simul(order=1,irf=0,noprint) l y x pi de s;

such that in line 202, dynare will return the variance of l :

variance.l(ii)=oo_.var(l_pos,l_pos);

In some cases, I manually find the position of l in oo_.var in the workspace/variable window, for example, here the 1st postion:

variance.l(ii)=oo_.var(1,1);

I think this command is also okay. BTW, the latest version of MATAB suggests using “string compare” find(strcmp()) to replace strmatch “string match”.

Then

variance.l=NaN(1,length(r_pi_vec))

is about looping over parameters, one hot link is

Loop over parameters

Hi, I know that this problem was broached on this forum many times since I found some topics connected with it but I didn’t find any solution for my problem. I want to look at some statistics of variables in the model when some parameters change. The steady state is calculated by hand and included in ‘initval’. For creating the loop over parameters, I used simple syntax like this: rhos = 0.8:0.05:1.05; for i=1:length(rhos); rho = rhos(i); stoch_simul(order=1); if info; disp('Computation fails…

Note line 182-183:

r_pi_vec=[1.5   1.5 5 1.5];
r_x_vec =[0.125 0   0 1  ];

They seem to be four alternative regimes. Each time, we assign new values to

    set_param_value('r_pi',r_pi_vec(ii));
    set_param_value('r_x',r_x_vec(ii));

r_pi is one-to-one to r_x. The whole loop will be executed 4 times in all.

    [info, oo_, options_] = stoch_simul(M_, options_, oo_, var_list_); %loop over stoch_simul
    if ~info(1)
        %read out current parameter values

If info number is not 1 (error), after running stoch_simul, oo_ and options_ will be updated.

My explanation may contain (many) errors. Please read it with caution.