如何运行滚动窗口面板回归？

Question

我有一个公司年面板数据，其中包含10000家公司的30年数据：

 firm       year        y        x1         x2
 A          1997       150      130        202
 A          1998       200      250        230
 A          1999       350       22         50
 A          2000        30      150         30
 B          1995        50      350         30
 B          1996       350      270         59
 B          1997       789      360        382
 B          1998       360      385         75
 B          1999       362      567        123

我想运行宽度为6的滚动窗口以进行面板回归：

plm(y~ lag(x1, 1)+ lag(log(x2), 1), data= dt, 
    model= "within", effect= "twoways",
    index= c("firm", "year"))

我认为，对于10000家公司来说，滚动窗口应该从1年到30年运行10000次。 我想获取每个公司年观测值的滞后系数（x1，1）。