Actually it is implemented but undocumented. There is a matlab (in dynare) routine called bvar_irf.m which can be called in the mod file (after the estimation command). The routine admits two inputs: the number of lags (any positive integer) and optionally the identification scheme:
- Cholesky, this is the default.
- SquareRoot, using matlab's sqrtm routine.
The advantage of the second approach is that it is independent of the ordering of the variables contrary to the first identification scheme (this advantage comes at the cost in terms of interpretation). The second argument has to be passed as a string. For instance:
- Code: Select all
bvar_irf(4,'SquareRoot');