使用 coeftest 在 stargazer 中添加信息

Question

我正在使用调查设计和稳健误差来估计模型，并尝试使用观星仪显示结果。但是，我无法添加基本信息，例如观察数、R2、调整后的 R2、残差标准。误差和 F 统计量，在表末尾。

我的代码是这样的（我缩写了模型）：

# Svy design.
svy <- svydesign(id = ~conglomerado, 
                 strata = ~estrato, 
                 check.strata = TRUE,
                 weights = ~fact_cal_esi, 
                 data = subset(completos, ocup_ref==1))

# Estimation.
short_model <- svyglm(ln_w ~ sex + age + factor(level),
                      design=svy)
short_model_r <- coeftest(short_model , vcov = vcovHC(short_model , type = "HC1"))

# Table.
stargazer(models,type="text", 
          column.labels = c("Short model"),
          title = 'Table 1: Short model')

但我不知道如何添加这些统计数据。

Answer 1

问题是

svyglm

不会产生像R平方这样的统计数据，然后可以将其传递给

stargazer

。看来您必须自己计算这些额外感兴趣的统计数据，然后使用

stargazer

的

add.lines

参数将它们传递给

stargazer

。您可以在此处的

add.lines

文档中阅读有关

stargazer

的信息：https://cran.r-project.org/web/packages/stargazer/stargazer.pdf.

这些其他帖子提出了有关从

svyglm()

计算这些统计数据的类似问题，并且可能对您这样做有所帮助：

Survey-package：如何从 svyglm-object 获取 R 平方？， https://stats.stackexchange.com/questions/618248/how-to-calculate-f-statistic-from-a-svyglm-model-in-r。