R: plm中的数据过滤

Question

我有一个14年的pdata.frame x 89个观测和10个变量+4个模型。

这些虚拟变量仅用于过滤(必要时)我的数据。在使用Stata时，我只需在代码的末尾添加一个"if VAR==1“。如何在R中与plm软件包一起使用？

示例

Stata码

quietly xtreg y x1 x2 if x3==1, fe

R码

plm( y ~ x1 + x2, data = PANEL, model = "within")

我必须创建单独的面板，已经过滤过的数据，还是可以在运行plm时做到这一点？

Answer 1

您可以在plm (subset=)中使用子集选项，也可以在拟合数据之前对数据进行子集。

使用来自区域==6上的包子集的数据集，

library(plm)
data("Produc", package = "plm")
fit1 = plm(gsp ~ hwy + pc, data = Produc, subset = region == 6)
fit2 = plm(gsp ~ hwy + pc, data = subset(Produc, region == 6))
identical(coefficients(fit1), coefficients(fit2))