Hi,

We are working with a third order approximation of a model. We have found that the IRFs are highly unstable. When investigating further, we found that they are unstable even if we simulate the model using only the first order terms. This is puzzling because the model is well behaved when doing a first order approximation. We compared oo_.dr.ghx and oo_.dr.ghu that we get from setting order = 1 to the oo_.ghx and oo_.dr.ghu that we get from setting order = 3, and they are very different, which is very odd. Could this be accounting for the unstable behavior? If yes, what can we do? Is it “kosher” to use the matrices from order=1 mixed up with the higher order ones calculated in order=3?

Thanks a lot