在不影响默认格式的情况下,用 t 值替换标准错误
Replace standard errors with t values without affecting the default format
我正在学习使用 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="。这是使用 report=('vc*t')
时的系数和 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"]))
但它似乎弄乱了表示显着性水平的“*”。
有没有办法解决这个问题?
我想你快到了。由于您向 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))
我正在学习使用 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="。这是使用 report=('vc*t')
时的系数和 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"]))
但它似乎弄乱了表示显着性水平的“*”。
有没有办法解决这个问题?
我想你快到了。由于您向 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))