R：到 LATEX 中的稳健标准错误答案

【问题标题】：R: robust standard errors into LATEXR：到 LATEX 中的稳健标准错误
【发布时间】：2018-09-08 08:20:21
【问题描述】：

我正在运行一些具有稳健标准误差的固定效应回归。我想用 stargazer 包把它们变成 LATEX 格式。有谁知道这是如何以一种简单有效的方式工作的？

我的回归代码如下所示：

library(plm)
panel_data <- pdata.frame(data, index = c("country","time"))
FE_model <- plm(Y~X+Z, data = panel_data, model = "within",  effect = "twoway")

我的标准错误代码是这样的：

library(lmtest)
coeftest(FE_model, vcov.=vcovHC(FE_model, method = c("arellano")))

我之前在 stargazer 包中使用的代码如下所示：

stargazer(FE_model_1,FE_model_2,FE_model_3, title="FE_model_Results",
      align=TRUE,
      omit.stat=c("LL","ser","f"), no.space = TRUE)

所以这只会给我正常标准误差的回归输出，但我想使用稳健的。

【问题讨论】：

标签： r latex stargazer

【解决方案1】：

使用来自plm 库的一些示例数据：

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

拟合来自plm 库文档的示例模型：

zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
          data = Produc, index = c("state","year"))

按照您的规定估计标准误：

se <- coeftest(zz, vcov.=vcovHC(zz, method = c("arellano")))

生成stargazer 表：

stargazer(zz, title="FE_model_Results",
          align=TRUE,
          omit.stat=c("LL","f"), no.space = TRUE,
          se = list(se[,2]))

【讨论】：