Quiero las mismas estrellas para significaciones en la salida de regresión en stargazer que en la "salida normal".

Produzco datos

library("stargazer"); library("lmtest"); library("sandwich")
set.seed(1234)
df <- data.frame(y=1001:1100)
df$x <- c(1:70,-100:-71) + rnorm(100, 0, 74.8)
model <- lm(log(y) ~ x, data=df)

y obtener algunas estimaciones del modelo donde el coeficiente de x tiene un valor p de 0.1023

coeftest(model, vcov = vcovHC(model, type="HC3"))

Quiero tener estos resultados en LaTeX. Basado en la misma función, calculo estimaciones estándar consistentes de heteroscedasticidad y dejo que Stargazer las use.

stderr_HC3_model <- sqrt(diag(vcovHC(model, type = "HC3")))
stargazer(model, se=list(stderr_HC3_model))

La salida del observador de estrellas tiene una estrella en el coeficiente que indica significancia cuando alfa = 10%. Quiero que stargazer dé lo mismo que el coeftest. (Debido a la comparabilidad con Stata donde reg L_y x, vce (hc3) da exactamente los resultados de coeftest).

Jugué con las opciones de stargazer p.auto, t.auto que no me ayudaron. Cuando ejecuto "stargazer" no puedo ver el código subyacente, ya que es posible en otros casos. ¿Qué hacer?

La respuesta de Richards me ayudó. Indico los pasos que usé para dar más de una regresión (digamos ols_a y ols_b).

ses <- list(coeftest(ols_a, vcov = vcovHC(ols_a, type="HC3"))[,2],
        coeftest(ols_b, vcov = vcovHC(ols_b, type="HC3"))[,2])
pvals <- list(coeftest(ols_a, vcov = vcovHC(ols_a, type="HC3"))[,4],
          coeftest(ols_b, vcov = vcovHC(ols_b, type="HC3"))[,4])
stargazer(ols_a, ols_b, type="text", p=pvals, se=ses)