Hi everyone!

I am running some robustness checks on a model that I am solving using Dynare, but I am getting some relatively weird results.

My problem is the following: I have written the model in two different .mod files. In the first, I have a very parsimonious version of the model, where some variables (that I keep constant) are simply treated as parameters. In the second .mod file, I wrote the same model, but declaring some of the parameters as endogenous variables instead, following Ar(1)s with a well-defined autocorrelation and variance for the shock (both declared as exogenous parameters).

As a sanity check, I decided to see what I would get if I set the variance of those shocks to zero in the second model file. In principle, Dynare should give me the same results for both .mod files (since the only difference is that in the first, variables x and y are parameters, while in the second they are constant endogenous variables). All goes well if I use the stoch_simul command for a 2nd order approximation and do not specify “periods”. The output that is printed (policy and transition functions, theoretical moments, etc.) is exactly the same.

However, I get an empty oo_.endo_simul series. To further check for errors, I decided to add the option “periods” to the stoch_simul function so as to generate series of endogenous variables and make sure that everything is consistent. This is where things start going wrong: the policy and transition functions are still the same, but this is no longer true for the rest of the output. In fact, some of the moments of simulated variables look *considerably* different for both models. I am trying to investigate what is the issue here. I have a suspicion, but would like to know if this sounds plausible: the second version of the model declares 10 exogenous shocks, while the first declares only 6 (in different orders). Is it possible that the order in which the shocks are declared matters so much? I would think it should not, especially if the model is stationary and the draw is long enough! This seems, however, to be the likely culprit, since the policy and transition functions (as well as the theoretical moments) are the same for both models - but still I wouldn’t expect the shock order to matter so much!

Thanks in advance!