Is there a way to plot the IRFs without running the model again?

Because I am trying to compare the IRFs for different sets of parameters and it’ll save a lot of computation time if I can find a way to get them without executing the command “dynare (model_name)” over and over again.

You have th values of irf in oo_ irf . Save the results file with a different name. Why don’t you put the irfs for all the shocks in a matrix and run them by setting up an m file

I have tried plotting using oo_.irfs and it works. But now I realized I’m back to the original problem: how do I get the impulse responses after a change of parameters without having to run the model again?

I actually agree with you that if I were to change the parameters, I have to re-run the model again.

However, I thought I happened to chance upon this m file that can enable me to do simulate the IRFs without having to run the model again. I guess I probably overlooked something.

I actually realized that because I’m changing parameters that have nothing to do with the steady state, I don’t have to run the model again every time I change that parameter value.

So I can resolve for the new steady state using the following command “[dr_,~,M_,~,oo_] = resol(0,M1_,options1_,oo1_);”. Now the thing is how do I get the impulse response functions from this command?

What I have now is I run a loop such that if there is more than one change in the parameter, at the second change, I’ll skip the “dynare (my model)” command and just go straight to resolving for the new steady-state.

But now, my problem is how do I get the impulse responses from the new steady-state? Because previously, I can simply do oo_.irfs right? But for this, I definitely can’t do just that.

By the way, I think that using oo_.irfs to plot the IRFs might not be so general a way because, correct me if I’m wrong, but let’s say your shock comes from epsilon_z, then you gotta execute “plot(oo_.irfs.Y_epsilon_z)” but what if I have multiple shocks at one time? How do I plot the impulse response functions more generally?

and then plot the results directly from the structures. With the set_param_value you can change values of parameters (for standard deviations, use a new shocks-block with the reset option).