0

我想知道如何使用 modelsummary 来组合模型名称和 DV 名称，就像在 Stata 中的 outreg2 中一样？这是代表：

    url <- 'https://vincentarelbundock.github.io/Rdatasets/csv/HistData/Guerry.csv'
dat <- read.csv(url)
models <- list(
  "OLS 1"     = lm(Donations ~ Literacy + Clergy, data = dat),
  "Poisson 1" = glm(Donations ~ Literacy + Commerce, family = poisson, data = dat),
  "OLS 2"     = lm(Crime_pers ~ Literacy + Clergy, data = dat),
  "Poisson 2" = glm(Crime_pers ~ Literacy + Commerce, family = poisson, data = dat),
  "OLS 3"     = lm(Crime_prop ~ Literacy + Clergy, data = dat)
)
modelsummary(models)

#N: DV names
modelsummary(dvnames(models), output = "flextable", estimate="{estimate}{stars}",
             statistic = 'statistic', stars = c('*' = .1, '**' = .05, '***'=0.01))
#N: Model names
modelsummary(models, output = "flextable", estimate="{estimate}{stars}",
             statistic = 'statistic', stars = c('*' = .1, '**' = .05, '***'=0.01))

以下是 DV 和型号名称组合表outreg2在 Stata 中的外观： 在此处输入图像描述

任何信息或建议将不胜感激！

4

1 回答 1

1

您可以使用该add_rows参数并创建自己的自定义函数来自动化该过程：

library(modelsummary)

insert_row <- function(x) {
    out <- c("DV:", names(dvnames(x)))
    out <- data.frame(as.list(out))
    attr(out, "position") <- 0
    return(out)
}

mod <- list(
    lm(mpg ~ hp, mtcars),
    lm(vs ~ hp, mtcars))

modelsummary(mod, add_rows = insert_row(mod))
模型 1 模型 2
DV： mpg 对比
（截距） 30.099 1.217
(1.634) (0.150)
生命值 -0.068 -0.005
(0.010) (0.001)
数量。 32 32
R2 0.602 0.523
R2 调整。 0.589 0.507
AIC 181.2 28.3
BIC 185.6 32.7
日志。 -87.619 -11.134
F 45.460 32.876
于 2021-12-18T02:42:41.307 回答