在 R 中使用 stargazer 在回歸 output 中包括 t 統計

Question

我嘗試使用 -stargazer- function 將三個回歸包含到一個表中。 我有以下代碼：

library(Jmisc)
library(tidyverse)
library(sandwich)
library(lmtest)
library(multiwayvcov)
library(stargazer)

set.seed(123)
df <- data.frame(
  x1 = rnorm(10, mean=0, sd=1),
  x2 = rnorm(10, mean=0, sd=1),
  y = rnorm(10, mean=0, sd=1)
)

r1 <- lm(y ~ x1 + x2, df)
cov1 <- vcovHC(r1, type="HC1", cluster="clustervar")
robust.se1 <- sqrt(diag(cov1))
t1 <- coef(r1)/robust.se1

r2 <- lm(y ~ x1, df)
cov2 <- vcovHC(r2, type="HC1", cluster="clustervar")
robust.se2 <- sqrt(diag(cov2))
t2 <- coef(r2)/robust.se2

r3 <- lm(y ~ x2, df)
cov3 <- vcovHC(r3, type="HC1", cluster="clustervar")
robust.se3 <- sqrt(diag(cov3))
t3 <- coef(r3)/robust.se2

stargazer(r1, r2, r3, 
se = NULL, 
t = list(t1, t2, t3), 
align=TRUE, 
type="html", 
nobs=TRUE, 
out="StargazerTest.txt")

生成的表報告標准誤差，而不是我創建的 t 統計量。 這很可能是由於底部的-stargazer- function。 我已經查找了它的目錄，但仍然不明白如何讓它做我想做的事。

Answer 1

如此處所述，您可以使用report指定所需的值（自stargazer 5.0 起）。 在您的情況下，刪除se = NULL和t = list(t1, t2, t3)並輸入：

report = ('c*t')

如：

stargazer(r1, r2, r3, 
          report = ('c*t'), 
          align=TRUE, 
          type="html", 
          nobs=TRUE, 
          out="StargazerTest.txt")

編輯：由於您需要使用穩健標准誤差，您應該使用 function coeftest （庫lmtest ）而不是手動計算穩健標准誤差。 下面是一個關於你的回歸的例子：

library(Jmisc)
library(tidyverse)
library(sandwich)
library(lmtest)
library(multiwayvcov)
library(stargazer)

set.seed(123)
df <- data.frame(
  x1 = rnorm(10, mean=0, sd=1),
  x2 = rnorm(10, mean=0, sd=1),
  y = rnorm(10, mean=0, sd=1)
)

r1 <- lm(y ~ x1 + x2, df)
cov1 <- vcovHC(r1, type="HC1", cluster="clustervar")
robust.se1 <- sqrt(diag(cov1))
t1 <- coef(r1)/robust.se1

foo <- coeftest(r1, vcov = vcovHC(r1, type = "HC1"))

stargazer(foo,  
          report = ('c*t'), 
          align=TRUE, 
          type="html", 
          nobs=TRUE,
          out="StargazerTest.txt")

請注意， foo給出了與t1相同的 t 值，但也顯示了系數、se 等，這使得stargazer可以正常工作