Thank you for your responses. I attempted the suggested approach, but it does not seem to work. Specifically, I used the command:

matrix b = irf(0,0,bundle1)

where the estimated SVEC model is stored as bundle1. It appears I may be missing a step or misunderstanding the procedure. I would appreciate further clarification.

The native irf() function (with underlying internal C code) only applies to the fully built-in gretl commands "var" or "vecm". It is understandable that the dual offerings of fully-native var/vecm commands and official but hansl-scripting-based add-on "SVAR" is still a bit confusing.

(BTW, the irf() function takes an optional bundle in the fourth position, not in the third, as duly documented in the function reference.)

For a bundle coming from an estimation with the SVAR addon, you need to access a bundle _member_, not invoke a function on it, because the irfs are already computed. Something like:

matrix b = bundle1.IRFs

where "IRFs" is the name that must be used exactly. Section 2.2 of the SVAR help describes how that matrix is structured, if you want to access certain shocks or targets or horizons. Notice that if you "just" want to produce a plot of an irf, you can use the SVAR functions "IRFplot" or "IRFsave", without having to meddle directly with that matrix.

cheers

sven