问题描述
我正在尝试 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))
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)