Hello,

Can anyone suggest how to simulate a sunspot shock in a linear model as in Figure 4 of the above-mentioned paper, available here: https://www.crei.cat/wp-content/uploads/users/pages/glv04jmcb.pdf

Thank you!

Hello,

Can anyone suggest how to simulate a sunspot shock in a linear model as in Figure 4 of the above-mentioned paper, available here: https://www.crei.cat/wp-content/uploads/users/pages/glv04jmcb.pdf

Thank you!

You should be able to use Bianchi/Nicolo’s method. See e.g.

Thank you very much for pointing me to Bianchi and Nicolo’s recent paper.

I am unable to match the IRFs exactly as in Figure 4 here: https://www.crei.cat/wp-content/uploads/users/pages/glv04jmcb.pdf

Just wanted to check if I am implementing it correctly. To simulate a sunspot shock to output, I have added the following auxiliary equation to the code of Gali et al (2004):

```
omeg = (1/alph)*omeg(-1) + sunspot - (y - expectation(-1)(y));
```

I’ve changed lambda = 0.85 (which is otherwise not feasible in the standard model) and alph=0.5 to ensure this is explosive.

I would really appreciate it if you could please advise if this is the correct way to implement this. Thank you!

Unfortunately, the original paper is not very explicit in how they did this. So I cannot really answer your question.