r - Stargazer 中的 p 值不精确

Question

我希望观星者的回归输出中的星星与“正常输出”中的星星相同。

我生产数据

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)

并获得一些模型估计值，其中 x 的系数的 p 值为 0.1023

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

我想在 LaTeX 中得到这些结果。基于相同的函数，我计算异方差一致的标准估计，并让 stargazer 使用它们。

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

当 alpha=10% 时，观星仪输出在系数处有一个星号，表示显着性。我希望观星者给出与 coeftest 相同的结果。 (由于与 Stata 的可比性，其中 reg L_y x, vce(hc3) 准确地给出了 coeftest 结果。)

我尝试了观星器选项 p.auto、t.auto，但没有帮助。当我执行“stargazer”时，我无法查看底层代码，因为在其他情况下是可能的。该怎么办？

<小时/>

理查兹的回答对我有帮助。我指出了我用来给出多个回归的步骤(比如说 ols_a 和 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)

Answer 1

您需要提供与您的coeftest相关的p值。来自手册页。

p a list of numeric vectors that will replace the default p-values for each model. Matched by element names. These will form the basis of decisions about significance stars

以下内容应该有效。

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

这提供了以下内容。

===============================================
                        Dependent variable:    
                    ---------------------------
                              log(y)           
-----------------------------------------------
x                            -0.00005          


Constant                     6.956***          
                              (0.003)          

-----------------------------------------------
Observations                    100            
R2                             0.026           
Adjusted R2                    0.016           
Residual Std. Error       0.027 (df = 98)      
F Statistic             2.620 (df = 1; 98)     
===============================================
Note:               *p<0.1; **p<0.05; ***p<0.01