r - 在 R 中的模型汇总输出表中计算和显示因变量的均值

Question

我正在使用 R 中 fixest 包中的 feols 命令来估计固定效应回归，在大多数情况下使用权重。为了输出我的结果，我主要使用 modelsummary 包中的 modelsummary 命令。我已经弄清楚如何使用 add_rows 功能来添加有关模型的更多信息(在我的例子中，我使用它来添加指示器以指示给定模型中是否存在控制变量，在类似于 modelsummary 可以本地指示给定模型中固定效应存在的方式)。

但是，我还想向每个模型添加模型因变量的均值，该模型是根据与给定模型估计的相同观测值进行估计的(即考虑到某些观测值因缺失值而下降，并且这可能因模型而异，具体取决于缺失值的结构)，并且最好仅针对我的治疗变量设置为 == 0 的观察结果。我怎样才能做到这一点？

我的基本回归结构如下(一些规范更简单，例如没有 Municipality.ID):

model = feols(depvar ~ i(Year.factor, Treatment.dummy, ref ='2007') + Year.factor + Treatment.dummy + Control.var, data = subset(data.frame, condition < limit), weights = Weight.var, panel.id = c(Year.factor, Municipality.ID), cluster = ~Municipality.ID)

我将这些估计循环到模型列表 mdls。 modelsummary 输出代码大致如下:

mdl.stats <- list(
  list("raw" = "nobs", "clean" = "Observations", "fmt" = 0),
  list("raw" = "r.squared", "clean" = "R Squared", "fmt" = 3),
  list("raw" = "r2", "clean" = "R Squared", "fmt" = 3),
  list("raw" = "r2.within", "clean" = "R2", "fmt" = 3),
  list("raw" = "within.r.squared", "clean" = "R2 (within)", "fmt" = 3),
  list("raw" = "r.squared.within", "clean" = "R2 (within)", "fmt" = 3),
  list("raw" = "FE: TERYT.f", "clean" = "Municipality FE", "fmt" = 0))
cntrls <- rep(c("", "X", ""), times=length(mdls)/3)
cntrls <- as.data.frame(t(c("Controls", cntrls)))
cntrls <- set_names(cntrls, c("Coefficients", names(mdls)))

modelsummary::modelsummary(mdls, gof_map = mdl.stats, stars = TRUE, output = "latex", booktabs = TRUE, add_rows = cntrls)

Answer 1

抱歉，我无法对 Vincent 的回答发表评论。 Fixest 在“my”下提供了“mean of dependent variable”fitstat，这样你的 glance_custom 函数就可以简化为:

glance_custom.fixest <- function(x, ...) {
  out <- data.frame("Mean(DV)" = as.numeric(fitstat(x, type = "my")))
  return(out)
}

而且无需依赖洞察力。