【问题标题】:Is there anyway to export feols model using stargazer in R?无论如何在R中使用stargazer导出feols模型?
【发布时间】:2021-02-19 22:53:56
【问题描述】:

我使用 feols 模型(最固定的包)运行了一堆模型,但我无法使用 stargazer 将模型导出到表中。关于如何做到这一点的任何建议?

看起来我可以使用“etable”功能,但我想使用 stargazer,因为我想在我的表格中添加几行注释并按照我想要的方式格式化表格(例如使用“table.layout” " 观星器中的函数)。

【问题讨论】:

  • 我不相信stargazer 支持这种模式。但是,modelsummary package 开箱即用地支持它。该软件包允许您添加注释,并且表格可以无限定制。当然,这是一个不同的包,这就是为什么我把它写成评论而不是答案。 (免责声明:我是作者)
  • @Vincent 这应该是一个答案
  • @Jakob 我不确定,因为作者明确声明他们想使用 stargazer 包的一个特定参数。 FWIW,网站上有一个fixest 示例:vincentarelbundock.github.io/modelsummary/articles/…
  • @Vincent 我明白了,但是如果你在谷歌上搜索如何制作最合适的乳胶表,你就会到这里。如果你赶时间,你就不会去看 cmets。
  • @Jakob 好点。我添加了一个答案。

标签: r export stargazer


【解决方案1】:

我不相信 stargazer 支持这种模式。但是,it is supported out-of-the-box by the modelsummary package. 这个包允许你添加注释,并且它生成的表格非常可定制,因为modelsummary 支持多个后端包来创建和自定义表格:kableExtragtflextable、@987654332 @。表格还可以导出为多种格式,包括 HTML、Markdown、LaTeX、JPG、data.frame 或 PDF。

(免责声明:我是modelsummary的作者。)

这是一个简单线性回归模型的示例:

library(fixest)
library(modelsummary)

# create a toy dataset
base <- iris
names(base) <- c("y", "x1", "x_endo_1", "x_inst_1", "fe")
base$x_inst_2 <- 0.2 * base$y + 0.2 * base$x_endo_1 + rnorm(150, sd = 0.5)
base$x_endo_2 <- 0.2 * base$y - 0.2 * base$x_inst_1 + rnorm(150, sd = 0.5)

# estimate
mod <- feols(y ~ sw(x1, x_endo_1, x_inst_1) | fe, data = base)

# table
modelsummary(mod)

您可以使用fixest 提供的各种公式函数,例如逐步包含协变量:

mod <- feols(y ~ sw(x1, x_endo_1, x_inst_1) | fe, data = base)
modelsummary(mod)

并且modelsummary 还支持工具变量估计。这将并排显示两个阶段:

mod <- feols(y ~ x1 | fe | x_endo_1 + x_endo_2 ~ x_inst_1 + x_inst_2, data = base)
modelsummary(summary(mod, stage = 1:2))

【讨论】:

    【解决方案2】:

    您也可以使用 fixest 中的etable 函数来导出输出表:

    library(fixest)
    data("mtcars")
    
    # models
    model1 <- feols(mpg ~ cyl + disp, data=mtcars)
    model2 <- feols(mpg ~ cyl +  hp, data=mtcars)
    
    # data.frame output
    df <- etable(list(model1, model2), tex=FALSE)
    
    # Latex output
    etable(list(model1, model2), tex=TRUE)
    

    您还可以使用file 参数将输出保存在本地。

    etable(list(model1, model2), tex=FALSE, file ='tt.txt')

    【讨论】:

      猜你喜欢
      • 2018-08-13
      • 2014-12-06
      • 2021-12-20
      • 1970-01-01
      • 1970-01-01
      • 2021-12-05
      • 2021-03-24
      • 1970-01-01
      • 2017-06-28
      相关资源
      最近更新 更多