- Are you doing stochastic simulations? Then you can use the
simult_
-function to run joint shocks. - Full geometric decay always takes forever. But people often use the half-life. See e.g. Persistency of IRFs: Definition - #15 by jpfeifer
- That seems like an occasionally binding constraint you could implement with Occbin.