R：plm 个体和时间固定效应，但没有其他回归量

Question

我想运行只包括时间和单个固定效应（即没有其他右侧变量）的回归。

我尝试使用plm：

plm(y ~ -1,data=data, effect="twoways", model="within")

但是，语法不正确，仅从模型公式中抑制 -1 也不起作用。

错误信息是：Error in uniqval[as.character(effect), , drop = F] : incorrect number of dimensions

对于 y 仅基于时间和单个固定效应的回归，使用 plm 的正确语法是什么？

谢谢！

Answer 1

在我看来，plm 不支持。请改用lm 方法：

library(plm)
data(Grunfeld)

# not possible with plm
mod <- plm(inv ~ -1, data=Grunfeld, model="within", effect = "twoways")

# use lm instead
mod2 <- lm(inv ~ -1 + factor(firm) + factor(year), data=Grunfeld)