Impulse responses don't display when solving to second order

[davidamdur]

I am using 4.0.4, and I have a simple Dynare model that behaves as follows. When I do stoch_simul(order=1), everything works fine and I see four sets of impulse response graphs. However, when I do stoch_simul(order=2), Dynare appears to solve the model successfully, but no impulse response graphs appear. Has anyone run into this? Is it (possibly) a bug in 4.0.4? I am new to Dynare, so I may well be missing something.

More details: when I run Dynare on Villaverde’s rbc.mod (from the ‘examples’ section of the website), the impulse response graphs appear both for first order and second order. However, for my model, they appear at first order only. I believe Dynare has solved my model successfully to second order, because I see nonzero coefficients on cross-product terms in the policy and transition functions table. Dynare also displays theoretical moments, variance decomposition, matrix of correlations and coefficients of autocorrelation. Just no impulse response graphs at second order. I can provide my mod file if that would help.

Any help greatly appreciated.

Thanks,
Dave

[davidamdur]

Here’s an update. When I look at oo_.irfs (at second order), I see all NaNs. Also, when I turn the “simul” flag on, all the simulated moments are also NaN. I went a little deeper and stepped through stoch_simul (and dependent functions) in the debugger. I noticed that my simulated endogenous variables grew very quickly in magnitude with each iterated period. Sure enough, after about five periods, all the simulated variables were NaN.

So it looks like I have an explosive model, at least to second order. What confuses me a bit is that the Blanchard-Kahn condition was satisfied, and Dynare did find a unique solution for the policy and transition functions. I also get well-defined theoretical moments. I thought that if Blanchard-Kahn was satisfied, the resulting solution was theoretically non-explosive?

I do suspect the “cuplrit” here is a near-unit-root in the model. I have a two-country international model with incomplete markets; without any frictions, such models have a unit root in the state variable for (relative) wealth. I’ve imposed a portfolio adjustment cost (similar to Schmitt-Grohe and Uribe (2003)), which should make the model theoretically stationary. However, perhaps it is “close enough” to being nonstationary that the Monte Carlo simulations fail.

So, does anyone out there have experience simulating international models with near-unit-roots in Dynare? Any ideas/workarounds for getting Monte Carlo simulations to work in this type of environment? As a newbie to Dynare, I appreciate any thoughts or advice.

Thanks,
Dave

[davidamdur]

I think I fixed my problem. My shock variances were way too high. When I reset the shock variances to something reasonable (0.01^2 instead of 1), everything worked fine. With the variances accidentally that high and the reversion back to steady-state fairly slow, the simulation was probably wandering very far from the steady state, where the dynamics are not very accurate.

It might be nice if MATLAB displayed an error for this case. For example, if the Monte Carlo simulation in simult_ yields (absolute) values for endogenous variables that are above a certain threshold, error out with a message like “Error simulating the model. The model appears explosive using the shock variances supplied. Ensure that the model is stationary and/or consider trying smaller shock variances.”

Anyway, I’m glad to be back in business (for the moment).

Dave

[davidamdur]

Just a note that that I found these closely related threads:

Error in output results - what does it mean?

I am trying to estimate a DSGE model with monopolistic competition and price rigidity according to Rotemberg (1983). The model, in levels, takes Peter Ireland’s framework. I have attached the mod.file I’m getting the follow output “error”. What does it mean? There is no impulse response functions neither. Can someone have a look and tell me where the problem might be? Merci/Thanks, pycs Irelandmoney4.mod (1.71 KB)

2nd order decision rule

Hi, I am solving a relatively big DSGE model using Dynare version 4.0.3: MODEL SUMMARY Number of variables: 260 Number of stochastic shocks: 1 Number of state variables: 86 Number of jumpers: 26 Number of static variables: 154 I did not have problems with the 1st order approximation, when all the irf`s were plotted correctly. Only changing to “order=2”, to run the 2nd approximation, the log file did not report errors, but all the irf 's values are NaN. Since Dynare…

Pruning in Dynare++

I’m getting explosive IRFs for a second-order approximation to a model in dynare++, which is chracteristic of 3rd and 4th order terms getting introduced from multiplication of 1st and 2nd order ones. Of course, the problem could be with my model somewhere, but it’d be nice to rule out the possibility of higher order terms creeping in. So does anyone know if dynare++ already has a pruning algorithm in it, or if there’s any other diagnostic I could use to see if this is the source of the problem…

and especially (for Dynare 3.0.65):

Pruning algorithm for second order approximation

Dear all, due to the ongoing difficulties with explosive simulations with second order approximations, I introduced the option for the pruning algorithm in Dynare 3.065 described in Kim, Kim, Schaumburg and Sims (2005): “Calculating and using second order accurate solutions of discrete time dynamic equilibrium models” In order to use it correctly under Windows, please proceed as follows: Download the simult_pruning.m file in the attachment to the Dynare directory, generally under …\dynare_v3…

Sounds like an implementation of Chris Sim’s pruning algorithm for version 4 would be a crowd pleaser.

[davidamdur]

Stephan has now posted a version of Chris Sims’ pruning algorithm for Dynare 4.0.4 here:

Pruning implementation for Dynare 4.0.4?

Does anyone have (or know of) an implementation of Chris Sims’ pruning algorithm that runs in Dynare 4.0.4? I’m interested in a modification of simult_.m along these lines: https://forum.dynare.org/t/pruning-algorithm-for-second-order-approximation/632/1 The above implementation was for 3.065 but doesn’t seem to run in 4.0.4. Any pointers appreciated. Thanks, Dave

This may be helpful if you’re encountering explosive IRFs and simulations at second-order, even though your model is stable to first-order. It was very helpful for me. Read Kim, Kim, Schaumberg and Sims (JEDC 2008) for details.

[ammm]

Thank you for this post, I have the same issue and reducing the shock variances or using the pruning algorithm helps. I just want to add that this algorithm is included in Dynare versions 4.3 and higher. Simply include the option ‘pruning’ into stoch_simul().