【问题标题】:how to display t value in texreg如何在texreg中显示t值
【发布时间】:2014-11-22 15:16:35
【问题描述】:

我是 R 的初学者,我有一个关于 texreg 的问题。

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

如果有人能给我一些提示,我将不胜感激,谢谢!

【问题讨论】:

  • 什么是texreg?请说明这是来自哪个包裹。还请发布一些示例数据和代码。
  • 由于你没有提供一个测试用例,你应该在包作者准备的文档的第 18-21 页上进行独立研究:jstatsoft.org/v55/i08

标签: r regression texreg


【解决方案1】:

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

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

【讨论】:

  • 社会学家通常应在括号中提供 t 值。
猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2018-10-01
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2017-11-27
  • 1970-01-01
相关资源
最近更新 更多