# Simulations with unchanged values

Hi all,

I’m trying to generate some simulated dataset of New Keynesian model for further estimation, but I got the same simulated values even when I changed the parameter values. I have no clue why it happens. Here attached the .mod file. Thank you in advance.

``````var pie R x z xi;
varexo eps_z eps_r;
parameters alpha beta kappa phi_x phi_pie rho_r rho_z sigma sigma_z sigma_r theta >omega; // DECLARATION OF THE DEEP PARAMETERS.
// Parameter Values
alpha   = 0.6;
beta    = 0.99;
phi_x   = 0.125;
phi_pie = 1.5;
rho_r   = 0.75;
rho_z   = 0.9;
sigma   = 1;
sigma_z = 0.3;
sigma_r = 0.2;
theta   = 6;
omega   = 1;
kappa   = ((1-alpha)*(1-alpha*beta)*(omega+theta))/(alpha*sigma*(omega+theta));
// DECLARATION OF THE (LINEAR) DSGE MODEL:
model(linear);
pie         = kappa*x + beta*pie(+1);
x(+1) - x   = sigma*(R - pie(+1) -z);
R           = rho_r*R(-1) + (1-rho_r)*phi_pie*pie + (1-rho_r)*phi_x*x + xi;
z           = rho_z*z(-1) + sigma_z*eps_z;
xi          = sigma_r*eps_r;
end;
shocks;
var eps_z;
stderr 1;
var eps_r;
stderr 1;
end;
stoch_simul(order = 1, periods = 232, simul_replic = 1, nograph);
[sim_array]=get_simul_replications(M_,options_);
y_sim=squeeze(sim_array(strmatch('y',M_.endo_names,'exact'),:,:));``````

Hi,

I do not know the routine `get_simul_replications`, but I suppose it reads data produced by `stoch_simul` and saved on the disk. You have to be aware that Dynare resets the seeds of the random generator between each run, meaning that it will use the same realisation of the sequence of shocks (`eps_z` and `eps_r`) in each run. Are you sure you are changing the values of the parameters? How do you do that?

Best,
Stéphane.

I’m using the following codes to loop over the alphas

``````% loop over values of alpha, get irf using dynare
alphas = 0.7:0.01:0.8;
irfs = zeros(length(alphas), 4*(h+1));
first_time = 1;
for i=1:length(alphas)
global p h;
if first_time;
set_param_value('alpha',alphas(i));
dynare dgploopvar noclearall; % degenerate 50 samples with alpha
first_time = 0;
else
set_param_value('alpha',alphas(i));
info = stoch_simul(var_list_);
if info;
disp(['Computation fails for alpha = ' num2str(alpha)]);
end;
end
% rest codes to extract information I want to know
end
``````

But even when I manually change the parameter values, I get the same simulations. It seems I didn’t solve the model correctly.

That is expected, because you are not handling parameter dependence correctly. `kappa` should be a model-local variable.

Thank you Prof. Pfeifer. It seems that though I run the loop and call Dynare to run the .mod file every time, the inside parameter values are not updated. I’m still trying to figure out why.

I debugged it. It’s because I didn’t update the results in `sim_arry`. It works with following matlab codes.

``````% loop over values of alpha, get irf using dynare
alphas = 0.7:0.01:0.8;
irfs = zeros(length(alphas), 4*(h+1));
first_time = 1;
for i=1:length(alphas)
global p h;
if first_time;
set_param_value('alpha',alphas(i));
dynare dgploopvar noclearall; % degenerate 50 samples with alpha
first_time = 0;
else
set_param_value('alpha',alphas(i));
info = stoch_simul(var_list_);
[sim_array]=get_simul_replications(M_,options_);
y_sim=squeeze(sim_array(strmatch('y',M_.endo_names,'exact'),:,:));
if info;
disp(['Computation fails for alpha = ' num2str(alpha)]);
end;
end
% rest codes to extract information I want to know
end
``````

Thank you both for the advice. I’ll mark this one as solved.