R: Надежная диагностика SE и модели в таблице Stargazer
Я пытаюсь поместить некоторые результаты регрессии 2SLS, сгенерированные через ivreg из пакета AER, в документ Latex, используя пакет stargazer. Однако у меня есть пара проблем, которые я не могу решить сам.
1. Я не могу понять, как вставить диагностику модели в соответствии с кратким описанием ivreg. А именно тесты слабых инструментов, Wu-Hausmann и Sargan Test. Я хотел бы иметь их со статистикой, обычно сообщаемой под таблицей, как число наблюдений, R в квадрате и Resid. SE. Кажется, у функции stargazer нет аргумента, в котором вы можете предоставить список с дополнительной диагностикой. Я не включил это в свой пример, потому что, честно говоря, понятия не имею, с чего начать.
2. Я хочу заменить обычные стандартные ошибки надежными стандартными ошибками, и единственный найденный способ сделать это - создать объекты с устойчивыми стандартными ошибками и добавить их в функцию stargazer с помощью se = list (). Я положил это в минимальный рабочий пример ниже. Может быть, есть более элегантный способ кодирования этого или, возможно, переоценка модели и сохранение ее с надежными стандартными ошибками? Помощь приветствуется.
library(AER)
library(stargazer)
y <- rnorm(100, 5, 10)
x <- rnorm(100, 3, 15)
z <- rnorm(100, 3, 7)
a <- rnorm(100, 1, 7)
b <- rnorm(100, 3, 5)
# Fitting IV models
fit1 <- ivreg(y ~ x + a |
a + z,
model = TRUE)
fit2 <- ivreg(y ~ x + a |
a + b + z,
model = TRUE)
# Here are the se's and the diagnostics i want
summary(fit1, vcov = sandwich, diagnostics=T)
summary(fit2, vcov = sandwich, diagnostics=T)
# Getting robust se's, i think HC0 is the standard
# used with "vcov=sandwich" from the above summary
cov1 <- vcovHC(fit1, type = "HC0")
robust1 <- sqrt(diag(cov1))
cov2 <- vcovHC(fit2, type = "HC0")
robust2 <- sqrt(diag(cov1))
# Create latex table
stargazer(fit1, fit2, type = "latex", se=list(robust1, robust2))