使用broom::augment Panel数据模型

Question

我想问一下，

在面板数据中使用augment函数时存在一个问题，特别是在第一差分估计上。

有没有一种方法可以使用函数augment或者任何一种方法来检索回归表，即使是对于第一差分估计？对于其他方法来说，效果非常好。

library(foreign)
library(plm)

crime <- read.dta("http://fmwww.bc.edu/ec-p/data/wooldridge/crime4.dta")

crime.p <- pdata.frame(crime,index=c("county","year"))

panel1 <- plm(log(crmrte) ~ polpc + prbconv + avgsen + density, data = crime.p, model = "within")
panel2 <- plm(log(crmrte) ~ polpc + prbconv + avgsen + density, data = crime.p, model = "random")
panel3 <- plm(log(crmrte) ~ polpc + prbconv + avgsen + density, data = crime.p, model = "fd")


broom::augment(panel1)
broom::augment(panel2)
broom::augment(panel3)

Answer 1

您可以访问broom::augment函数手动生成的组件。下面的代码可能会给你一些指导。我首先将模型作为panel3列表的一部分进行访问，然后按名称加入残差。在某些情况下，残差丢失了，我怀疑这就是为什么broom::augment似乎不起作用的原因。

data.frame("rownames" = row.names(panel3$model), panel3$model) %>%
      left_join(data.frame("rownames" = names(panel3$residuals), "resid" = panel3$residuals))