修复了St用户在R中的效果

Question

我已经找了一个类似的问题，但还没找到我正在寻找的东西。

在Stata中，我可以使用面板数据（学校和时间固定效果）运行下面的模型：

xi: areg total i.year*group_dummy, absorb(school) cluster (school)

如何使用plm包在R中运行类似的回归？

这就是我的尝试：

model1 <- plm(total ~ group_dummy, 
          na.action = na.omit,
          data   = pdlong, 
          method = "within",
          index = c("school", "year"),
          effect = "twoway")

但是我收到一条错误消息（并且不确定应该如何解释）：

x [，ii]：下标超出范围

时出错