Hello, I try to estimate and draw a chart of IRF for 25 combinations based on this model :

yi,t+k=𝛼ik+𝛾tk+𝛽k∆Ri,t+∑𝜗jk∆Ri,t −j+∑𝜃jkyi,t−j+𝜀tk
that
∑𝜗jk∆Ri,t−j
is from 1 to 5

and also
∑𝜃jkyi,t−j
is from 1 to 5


∆Ri,t
is the shock

I do not know how to solve it and draw the chart for 25 combinations in a loop form ( because it could be lag one versus lag one, lag one versus two ......, lag tow versus lag one, lag two versus lag two ....and so on)

I hope I receive your guidance.


Many thanks in advance.

regards,