如何使用 PLM 运行州级固定效果

Question

我正在尝试运行仅具有州级固定效应而不是时间固定效应的回归。

我在尝试：

lm1 <- plm(lnwage ~ age + age^2 + education, data = cps, index = "state", model = "within")

但我没有任何运气，我只是在网上找到同时使用状态和时间固定效应的人的信息。

Answer 1

为了计算任何固定效应，我们可以为相关变量添加一个虚拟变量，好吗？ 让我们考虑plm()文档中的示例，并首先使用状态固定效果的状态虚拟对象执行普通lm() 。

data("Produc", package = "plm")

fe.lm <- lm(log(gsp) ~ 0 + log(pcap) + log(pc) + log(emp) + unemp + 
              factor(state), data=Produc)

> summary(fe.lm)$coef
                                Estimate   Std. Error    t value
log(pcap)                   -0.026149654 0.0290015755 -0.9016632
log(pc)                      0.292006925 0.0251196728 11.6246309
log(emp)                     0.768159473 0.0300917394 25.5272539
unemp                       -0.005297741 0.0009887257 -5.3581508
factor(state)ALABAMA         2.201617056 0.1760038727 12.5089126
factor(state)ARIZONA         2.368088138 0.1751884949 13.5173725
factor(state)ARKANSAS        2.263015801 0.1671716685 13.5370773
...

现在我们使用plm() ，我们还必须添加虚拟对象，但是它没有显示在输出中。

library(plm)
fe.plm <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp + factor(state),
              data=Produc, index=c("state","year"), model="within")

> summary(fe.plm)$coef
              Estimate   Std. Error    t-value      Pr(>|t|)
log(pcap) -0.026149654 0.0290015755 -0.9016632  3.675200e-01
log(pc)    0.292006925 0.0251196728 11.6246309  7.075069e-29
log(emp)   0.768159473 0.0300917394 25.5272539 2.021455e-104
unemp     -0.005297741 0.0009887257 -5.3581508  1.113946e-07

另一种选择是felm() 。

library(lfe)
fe.felm <- felm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp | state | 0, 
                data=Produc) 

> summary(fe.felm)$coef
              Estimate   Std. Error    t value      Pr(>|t|)
log(pcap) -0.026149654 0.0290015755 -0.9016632  3.675200e-01
log(pc)    0.292006925 0.0251196728 11.6246309  7.075069e-29
log(emp)   0.768159473 0.0300917394 25.5272539 2.021455e-104
unemp     -0.005297741 0.0009887257 -5.3581508  1.113946e-07

正如人们所见，一切都会产生相同的值。