如何在texreg中显示t值

时间:2014-11-22 15:16:36

标签: r regression texreg

我是R的初学者,我对texreg有疑问。

我一直在网上搜索信息好几天,但我找不到多少。我想通过使用texreg在我导出的表中包含t值,并且我希望我的t值位于[]中每个系数的正下方。

如果有人能给我提供一些提示,真的很感激,谢谢!

1 个答案:

答案 0 :(得分:0)

通常,标准错误应在括号中报告。但是,如果您真的想用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))

目前无法用方括号替换圆括号。