【问题标题】:Exporting R regression results to tables for publication将 R 回归结果导出到表中以供发布
【发布时间】:2015-02-18 07:03:36
【问题描述】:

我想知道 R 是否具有存储多个回归结果的功能,决定您想要保留的内容(例如,对置信区间不感兴趣),然后将这些结果导出到 Word 或 .rtf 文档中的方便表格中可以插入学术工作。

我正在考虑Stata 中esttab 命令所允许的功能。

因此,例如在 Stata 中,您可以使用每次执行回归后存储结果

est sto a1 

然后,一旦你存储了各种版本,你就可以使用以下命令

esttab a1 a2 a3 a4 a5 a6 a7 using "C:\...\table_name.rtf", replace compress onecell nogaps stats(N N_g r2) star b(3) se(3) starlevels(* 0.05 ** 0.01 *** 0.001) title(Table 1: This is a Stata Table ) order(order of your variables in the produced table) keep(which variables to keep if there are some you don't want - useful for when you have time series year dummies for instance)

有类似的吗?

响应者提出问题(问题为负 2 ???)使用我的一些代码进行更新。这是非常基本的,Stata 代码准确地捕获了我想要做的事情,所以我对负面评分有点困惑。

假设您有以下回归:

    ctrl <- vv10 ~ ro_rnd + ro_int + ro_col + c_obj + ro_org2 + sc_org + sc_funding + sc_ties + sc_friends + sc_know + r_jobtime + r_wj + sen_match1 + search + lm.res + scenprop

lm.base <- lm(ctrl, data=d)
  temp1 <- summary(lm.base)

lm.main <- update(lm.base, ~. + r_netw_kn + r_brmean6 + r_conf3 + I(r_conf3^2) + r_wtwws, data=d) 
  temp2 <- summary(lm.main) # model 2

lm.intp1 <- update(lm.main, ~. + r_brmean6:r_conf3 + r_brmean6:I(r_conf3^2))
  temp3 <- summary(lm.intp1)

三个回归的摘要存储在三个名为temp1, temp2, and temp3的临时列表中

我想要从这三个列表中导出特定部分(即摘要中的列),以便我可以保留 estimatestandard errorthe stars,模型适合(在通过temp1[9] 进行 lm 回归的情况,理想情况下重新排序变量,将标准误差放在括号内,并减少逗号后的数字数量,这样我就可以直接在表格中使用它们来写论文。 在 Stata 中,这可以通过前面提到的 esttab 命令实现。使用knitr 并不能真正让我灵活地删除我不需要的东西并创建一个包含多个逐步回归的表,其中每一步都添加新变量(这是表示管理结果的常用方法) .

希望这个澄清会有所帮助。

【问题讨论】:

  • 试试这个:library(knitr)
  • 我想您的原始帖子存在两个主要问题:缺少 R 代码和假设 Stata 代码对于您正在处理的 R 人员来说是不言而喻的或熟悉的(在绝大多数情况下不太可能)例)。所以人们经常对问题进行轻微的抨击,或者更经常地忽略问题,而没有持续尝试自己编写代码。 “做这件事的代码是什么?”形式的问题。并不真正符合论坛的目标,尽管有时会被放纵。
  • 在这种情况下,显示来自 esttab 的结果可能比显示语法提供更多信息,但显示您尝试过 R 代码更重要。
  • 谢谢尼克。不知道如何上传 Stata 结果表 - 也不确定这是否有效 - 但我希望我能更好地解释它。我认为通过标记 r 和 stata 我本来可以避免仇恨者:)

标签: r export export-to-excel regression stata


【解决方案1】:

你能用你的一些代码更新这个问题吗?从那里可以确定如何提取。

作为一个没有看到你的模型和R代码的例子,看看提问者如何从线性模型sec.reg这里提取结果:XTS dates from different sources. Using R to calculate beta

提问者将模型的特定部分(alpha)提取到一个表中,循环中包含以下代码部分:
summary[n,1] &lt;- coef(sec.reg)[1]

尝试使用包xtable 获取出版质量表。如果您不熟悉,您可能还想查看包裹 knitr

【讨论】:

  • 嗨,这看起来是朝着这个方向迈出了很好的一步,但我看不出它是如何提取星数和标准错误的。
  • 你能用你的一些代码更新这个问题吗?从那里可以确定如何提取。在名为 sec.reg 的线性模型中,提问者将模型的特定部分 (alpha) 提取到代码部分包含 [n,1] 的表中。我在手机上,因此草率回答。
  • 另见观星包
猜你喜欢
  • 1970-01-01
  • 2021-07-20
  • 1970-01-01
  • 1970-01-01
  • 1970-01-01
  • 2014-01-17
  • 2021-12-18
  • 1970-01-01
  • 1970-01-01
相关资源
最近更新 更多