R: устойчивые стандартные ошибки в LATEX - PullRequest
Новая
       14

R: устойчивые стандартные ошибки в LATEX

0 голосов
/ 08 сентября 2018

Я использую некоторые регрессии с фиксированным эффектом и устойчивыми стандартными ошибками. И я хотел бы получить их в формате LATEX с пакетом stargazer. Кто-нибудь знает, как это работает простым эффективным способом?

Мой код для регрессий выглядит следующим образом: 

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)

так что это даст мне только вывод регрессий с обычными стандартными ошибками, но я бы хотел использовать надежные.

1 Ответ

0 голосов
/ 08 сентября 2018

Использование некоторых примеров данных из библиотеки 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]))
...