Stargazer 回归输出，对于稳健的标准误差显示错误的重要性星号

如何解决Stargazer 回归输出，对于稳健的标准误差显示错误的重要性星号

我正在尝试 1. 获取观星者输出以显示正确的标准错误并显示正确的意义（星星）； 2. 添加稳健的 F 统计量（来自线性假设）。

问题 1： 意义星不正确。如果您将 r_mod1 与 stargazer 的输出进行比较，您会注意到 r_mod1 的已婚变量只有 1 颗星，但 stargazer 输出显示两颗星。

问题 2： 我不知道如何获得稳健的 F 统计量。我读到您必须在某处手动执行此操作，因为 stargazer 没有执行此操作的功能。这是真的？对于像 stargazer 这样开发良好的软件包来说，这似乎是一个基本限制。

任何帮助将不胜感激！

我的代码：

# Libraries

library(tidyverse)
library(Ecdat)
library(lmtest) 
library(sandwich) 
library(stargazer)
library(car)

# Get data from Ecdat
data(Wages)
wages <- Wages
rm(Wages)


# A sample model

mod1 <- lm(lwage ~ ed + exp + I(exp^2) + wks + bluecol + ind + 
             south + smsa + married + sex + union + black,data = wages)


# Robust standard errors
r_mod1 <- coeftest(mod1,vcov=(vcovHC(mod1,"HC1")))

# Robust F-stats
lh1 <- linearHypothesis(mod1,names(coef(mod1))[-1],white.adjust = T)

# Extract robust values
se1 <- r_mod1[,2]
p1 <- r_mod1[,4]
t1 <- r_mod1[,3]


# STARGAZER
stargazer(mod1,type = "text",column.labels=c("Model one (All)"),title = "Regression results",style = "qje",single.row = TRUE,t.auto = T,# Necessary?
          p.auto = T,# Necessary?
          se = list(se1),p = list(p1),t = list(t1))