Eu tento colocar algumas saídas de regressão 2SLS geradas via ivreg do pacote AER em um documento Latex usando o pacote stargazer. Tenho alguns problemas, no entanto, que não consigo me resolver.
1. Não consigo descobrir como inserir o diagnóstico do modelo, conforme fornecido pelo resumo do ivreg. Nomeadamente testes de instrumentos fracos, Wu-Hausmann e Sargan Test. Eu gostaria de tê-los com as estatísticas geralmente relatadas embaixo da tabela, como número de observações, R ao quadrado e Resid. SE. A função stargazer parece não ter um argumento em que você pode fornecer uma lista com diagnósticos adicionais. Eu não coloquei isso no meu exemplo, porque honestamente não tenho idéia por onde começar.
2. Quero trocar os erros padrão normais por erros padrão robustos e a única maneira de fazer isso é produzir objetos com erros padrão robustos e adicioná-los na função stargazer com se = list (). Eu coloquei isso no exemplo de trabalho mínimo abaixo. Existe talvez uma maneira mais elegante de codificar isso ou reestimar o modelo e salvá-lo com erros padrão robustos? Ajuda apreciada.

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))