面板数据的滚动窗口回归

Question

我想在36个月的时间内对面板数据进行滚动窗口回归，并获得每月的截距作为输出。 我的数据有1397个基金（ID），每个月都有252个月收益。

    ï..ID Period Return         RMRF         SMB5          HML          RMW
1-1     1      1   <NA>  0.027131614 -0.000206798 -0.021403548  0.017474395
1-2     1      2   <NA>  0.009564262  0.025552733 -0.011379760  0.022345710
1-3     1      3   <NA>  0.014315746 -0.002847203  0.014032642 -0.007219692
1-4     1      4   <NA>  0.048363609  0.008673535  0.019081777  0.003571456
1-5     1      5   <NA> -0.025923904  0.044332726  0.001115927  0.001927510
1-6     1      6   <NA>  0.026734862  0.001516872  0.005258489  0.003943867
             CMA
1-1 -0.036913366
1-2  0.018724877
1-3  0.016066358
1-4 -0.003265331
1-5  0.018159172
1-6 -0.014966470

         ï..ID Period       Return         RMRF         SMB5          HML
1397-247  1397    247  0.005945991 -0.002882396 -0.013501132  0.000409125
1397-248  1397    248 -0.015762253 -0.012668465 -0.006875024 -0.019310729
1397-249  1397    249   0.03603303  0.032361027 -0.006767910  0.031112181
1397-250  1397    250  0.001932212  0.009555836  0.008110662 -0.004983660
1397-251  1397    251  0.019284853  0.030316686  0.019022510 -0.023171037
1397-252  1397    252  0.013243989  0.022730605  0.026619563  0.003929087
                  RMW          CMA
1397-247 -0.011966912  0.000511407
1397-248  0.008897090 -0.000202207
1397-249 -0.004747638  0.012626697
1397-250  0.012410290 -0.015537292
1397-251  0.013356231 -0.016695288
1397-252 -0.003929087  0.001375180

```

到目前为止，我已经尝试了以下代码：

#Rolling Panel Regression
Roll_Alphas <- roll_regres(formula = plm(Returns ~ RMRF,data=pdata, model = "within"), width= 36)

但收到错误消息：

error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :  0 (non-NA) cases

我认为错误是由于Return变量中的NA值引起的，但是我尝试用随机数替换它们（只是为了查看它是否可以工作），并且仍然遇到相同的错误。

我又如何使用健壮的标准误差创建此回归？