[英]R: robust standard errors into LATEX
我正在运行带有固定标准误差的固定效应回归。 我想通过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)
因此,这只会给我正常标准误差的回归输出,但我想使用可靠的回归。
使用来自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]))
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.