Question

使用plm运行带有面板数据的回归后，我遇到一个简单的问题，其数据集类似于下面的数据集：

dataset <- data.frame(id = rep(c(1,2,3,4,5), 2),
                      time = rep(c(0,1), each = 5),
                      group = rep(c(0,1,0,0,1), 2),
                      Y = runif(10,0,1))
model <-plm(Y ~ time*group, method = 'fd', effect = 'twoways', data = dataset,
            index = c('id', 'time'))
summary(model)

stargazer(model)

正如您所看到的，模型summary和stargazer显示的表格都会说我的观察数量是10.但是，说{{1}这样更不正确因为我带走了第一个差异后的时间元素？

Answer 1

你对观察的数量是正确的。但是，您的代码并不是您想要它做的（第一个不同的模型）。

如果您想要第一个差异模型，请将参数method切换为model（并删除参数effect，因为它对第一个差异模型没有意义）：

model <-plm(Y ~ time*group, model = 'fd', data = dataset,
            index = c('id', 'time'))
summary(model)

## Oneway (individual) effect First-Difference Model
## 
## Call:
## plm(formula = Y ~ time * group, data = dataset, model = "fd", 
##     index = c("id", "time"))
## 
## Balanced Panel: n = 5, T = 2, N = 10
## Observations used in estimation: 5
## 
## Residuals:
##       Min.    1st Qu.     Median    3rd Qu.       Max. 
## -0.3067240 -0.0012185  0.0012185  0.1367080  0.1700160 
## [...]

在摘要输出中，您可以看到原始数据中的观测数（N = 10）和FD模型中使用的观测数（5）。

使用具有第一个差异的plm时的观察数

1 个答案: