用模型总结法计算扫帚后的边际效应

Question

我试图1.在broom中运行嵌套回归后计算边际效果，然后将输出转换为模型摘要可读的格式，以便使用readable将输出编成Word。

library(tidyverse)
library(broom)
library(mfx)
library(modelsummary)


year <- rep(2014:2015, length.out=10000)
group <- sample(c(0,1,2,3,4,5,6), replace=TRUE, size=10000)
female <- sample(c(0,1), replace=TRUE, size=10000)
smoker <- sample(c(0,1), replace=TRUE, size=10000)
dta <- tibble(year = year, group = group, female = female, smoker = smoker)

mods <- dta %>%
  nest(data = c(-year)) %>%
  mutate(model = map(data, ~ glm(smoker ~ female*group, data = .,
                                 family = binomial(link = "probit"))),  # this is the model itself
         mfx=probitmfx(model, data=data))%>%                             # trying to convert to marginal effects 
  tidy()                                                                 # convert to readable format

modelsummary(mods)

我每走一步都会收到各种各样的错误信息。

Answer 1

我想出来了：

mods <- dta %>%
  nest(data = c(-year)) %>%
  mutate(model = map(data, ~ glm(smoker ~ female*group, data = .,
                                 family = binomial(link = "probit"))),
         mfx=map(data, ~ probitmfx(smoker ~ female*group, data = .,
                             )))
mfx <-mods[[4]]

modelsummary(mfx)