R中的多面板回归表

Question

我正在尝试创建一个回归表，其中包含两个面板，这两个面板有部分重叠但不完全相同的因变量。对于每个面板，我都想要一个标题，用于指示每个面板的名称以及每个回归的因变量。

# rm(list = ls())

data("mtcars")

# Regressions for panel 1
reg1 <- lm(mpg ~ cyl + disp, data = mtcars)
reg2 <- lm(cyl ~ cyl + disp, data = mtcars)

# Regressions for panel 2
reg3 <- lm(mpg ~ hp + drat, data = mtcars)
reg4 <- lm(cyl ~ hp + drat, data = mtcars)
reg5 <- lm(disp ~ hp + qsec, data = mtcars)

# Panel 1 output
panel1 <- stargazer::stargazer(
  reg1, reg2,
  float = TRUE,
  header = FALSE,
  model.numbers = FALSE,
  omit.table.layout = "n",
  multicolumn = FALSE,
  dep.var.caption = "",
  type = "latex",
  align = TRUE, 
  digits = 2,
  df = FALSE, 
  digits.extra = 2,
  nobs = TRUE,
  omit.stat = c("rsq", "adj.rsq", "ser") 
)

# Panel 2 output
panel2 <- 
  stargazer::stargazer(
  reg3, reg4, reg5,
  float = TRUE,
  header = FALSE,
  model.numbers = FALSE,
  omit.table.layout = "n",
  multicolumn = FALSE,
  dep.var.caption = "",
  type = "latex",
  align = TRUE, 
  digits = 2,
  df = FALSE, 
  digits.extra = 2,
  nobs = TRUE,
  omit.stat = c("rsq", "adj.rsq", "ser")
)

# Plot panels together
table <- 
  starpolishr::star_panel(
  panel1, 
  panel2,
  panel.names = c(
    "Panel 1", 
    "Panel 2"
  ), 
  same.summary.stats = FALSE, 
  same.lhs.vars = FALSE )

# Save as a .tex file
starpolishr::star_tex_write(
  starlist = table, 
  file = paste0(here::here(), "/panel_mwes.tex"), 
  headers = FALSE
)

此代码输出如下: 我需要改变的是：1)将“Panel A: Panel 1”移到因变量名称之上；2)在“Panel B: Panel 2”下添加因变量名称。是否有好的方法可以实现这一点（在stargazer / starpolishr或其他方式中），而不会太过hacky？

- Thomas J. Brailey

---

请勿在您的代码中使用 rm(list = ls())。如果有人运行您的代码，可能会导致大量数据丢失。请注意避免此类操作。 - Martin Gal

@Martin Gal，谢谢你的提醒，我忘记注释掉它了。 - Thomas J. Brailey

这两个表格需要垂直连接成一个单独的表格吗？还是将它们作为两个单独的表格或水平连接也可以？如果是后者，texreg 中的 custom.header 参数是否能够实现你所要求的功能？请参见此处。 - Philip Leifeld

谢谢@PhilipLeifeld--不幸的是，这些面板需要垂直连接。 - Thomas J. Brailey

Answer 1

所以这是我做事的笨拙方式（我不会接受这个作为答案，因为肯定有更好的方法）。

# 1. Flip the panel name and column names around
col_names  <- table[8]
panel_name <- table[10]

table[8]  <- panel_name
table[10] <- col_names

# 2. Insert custom column names below panel 2 header
table <- 
  starpolishr::star_insert_row(
    table, 
    c("\\hline \\\\[-1.8ex] & \\multicolumn{1}{c}{mpg} & \\multicolumn{1}{c}{cyl} & \\multicolumn{1}{c}{disp} \\\\ \\hline \\\\[-1.8ex]"),
    insert.after = c(24)
  )

# Save output
starpolishr::star_tex_write(
  starlist = table, 
  file = paste0(here::here(), "/panel_mwe2.tex"), 
  headers = FALSE
)

以下是我想要实现的大致输出：

- Thomas J. Brailey