r - 用 t 值替换标准错误而不影响默认格式

Question

我正在学习使用 stargazer 输出回归表。我一直在尝试用 t 值替换标准错误（括号中），并尝试了以下代码：

f1 <- lm(y1 ~ x1 x2, data=mydata)

stargazer(f1, type="text", out="models.txt", report=('vc*t'))

该report=('vc*t')选项确实将标准误差替换为 t 值。但它也删除了标准错误周围的原始括号，并在 t 值前面添加了“t=”。以下是使用 时的系数和 t 值report=('vc*t')：

-0.026***         
t = -3.681                                                          

0.002***          
t = 3.260    

理想的输出应该是这样的：

-0.026***         
(-3.681)                                                          

0.002***          
(3.260) 

我也试过这个se选项：

stargazer(f1, type="text", out="models.txt", se=list(summary(f1)$coef[c(2:4,1),"t value"]))

但它似乎弄乱了表示显着性水平的“*”。

有没有办法解决这个问题？

Answer 1

我想你几乎在那里。由于您向 stargazer 传递了标准误差的新参数，它将使用这些参数来计算显着性。将你的原始 p 值列表传递给 stargazer，你应该已经准备好了。这是一个完整的例子。

data(iris)
ols.result <- lm(Sepal.Length ~ Sepal.Width + Petal.Length, data = iris)
ols.t.stat <- summary(ols.result)$coef[ , "t value"]
ols.p.val  <- summary(ols.result)$coef[ , "Pr(>|t|)"]
stargazer(ols.result, se = list(ols.t.stat), p = list(ols.p.val))