回归时间序列事件r

Question

我的时间序列如下：

  year         a        b    event         x
1 1950  0.033222        2        0  1055.461
2 1950  0.033222      345        0  1055.461
3 1950  0.033222      290        0  1055.461
4 1950  0.033222      355        0  1055.461
5 1950  0.033222      710        1  1055.461
6 1950  0.033222       36        0  1055.461

我想知道当虚拟变量“event”变为1时x会发生什么。此外，我想知道30年后的效果。所以我的问题是，可能会有一些其他1将在12或21年后出现，例如。如何控制一个单一事件或考虑30年内可能出现的所有事件？

我希望有类似的东西：

lm(x~event, data=data)

其中x是在t年，并且发生了前30年发生的所有事件

我想使用VAR包，但我不知道如何在他们的功能中添加假人。