通常,标准误差应在括号中报告。但是,如果你真的想用 t 值替换它们,你可以这样做。
来自lm 文档的示例回归分析:
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
lm.D90 <- lm(weight ~ group - 1)
解决方案 1:
screenreg(list(lm.D9, lm.D90), override.se =
list(c(0.5, 0.5), c(0.5, 0.5)))
这会覆盖标准错误并用您提供的自定义值替换它们。结果表:
=================================
Model 1 Model 2
---------------------------------
(Intercept) 5.03 ***
(0.50)
groupTrt -0.37 4.66 ***
(0.50) (0.50)
groupCtl 5.03 ***
(0.50)
---------------------------------
R^2 0.07 0.98
Adj. R^2 0.02 0.98
Num. obs. 20 20
=================================
*** p < 0.001, ** p < 0.01, * p < 0.05
解决方案 2:
您可以提取系数等并将其保存到对象中,操作该对象,然后将操作后的对象交给 texreg,而不是将模型交给 texreg。
tr1 <- extract(lm.D9)
tr1@se <- c(0.5, 0.5) # enter new values, e.g., t values
tr2 <- extract(lm.D90)
tr2@se <- c(0.5, 0.5) # enter new values, e.g., t values
screenreg(list(tr1, tr2))
目前无法用方括号代替圆括号。