如何为R中含多个因变量的多变量多重回归模型生成发表级结果表格？

stargazer原生不支持cbind()构造的多因变量回归返回对象（属于mvreg类），因此会触发你遇到的报错，以下是3种可替代的实现方案，均可以输出带显著性星号、符合发表规范的结果表格：

方案1：使用modelsummary包（最推荐） #

modelsummary是目前主流的stargazer替代包，原生支持多因变量回归结果，可输出HTML、LaTeX、Word等多种格式，自定义程度极高，可轻松匹配stargazer的经典表格风格。
示例代码如下：

# 安装包（首次使用运行）
# install.packages("modelsummary")
library(modelsummary)

# 直接输出多因变量回归结果，type参数指定输出格式
# output = "html" 输出网页格式，output = "latex"输出tex代码，output = "docx"输出word文件
modelsummary(m1, output = "html", stars = TRUE)

• 可通过gof_map参数调整需要展示的拟合优度指标，title参数添加表格标题，对齐方式、字体等均可自定义。

方案2：使用texreg包 #

texreg是另一款成熟的发表级表格生成包，对多因变量回归有良好支持，对LaTeX格式的适配度更高。
示例代码如下：

# 安装包（首次使用运行）
# install.packages("texreg")
library(texreg)

# 输出HTML格式结果，可直接复制到论文或导出保存
htmlreg(m1, star.symbol = "*", stars = c(0.01, 0.05, 0.1), doctype = FALSE)

• 调整custom.note参数可以自定义表注内容，默认的显著性标注规则和stargazer完全一致。

方案3：拆分多因变量模型喂给stargazer（无需更换工具） #

如果不想更换惯用的stargazer包，可以先把多因变量回归拆分为多个单因变量回归，再传入stargazer生成表格：

# 批量生成3个单因变量回归模型
model_list <- lapply(list(dv1, dv2, dv3), function(dv) lm(dv ~ iv1 + iv2))

# 用stargazer输出结果
stargazer(model_list, type = "html", dep.var.labels = c("dv1", "dv2", "dv3"))

• 该方案输出的表格格式和你之前用stargazer生成的完全一致，仅需要多一步拆分模型的操作。

内容的提问来源于stack exchange，提问作者stm