Я использовал tidyr::nest
для запуска серии моделей логистической регрессии с различными зависимыми переменными. Я хочу вывести результаты в виде единой таблицы html в RMarkdown с каждой моделью в виде столбца и строками в виде экспоненциальных коэффициентов с 99% доверительных интервалов. Я не могу понять рабочий процесс между nest
, tidy
и пакетом визуализации таблиц, таким как stargazer
, чтобы заставить это работать. Если я unnest
свой tidy
вывод и передаю его stargazer
, или если я просто пытаюсь передать вывод nest
ed (переменная во вложенном фрейме данных, называемая моделью ниже) напрямую stargazer
, я не получаю вывода. Я бы предпочел использовать результат tidy
из-за экспоненциальных коэффициентов и 99% доверительных интервалов. Мне действительно нужна эта виньетка, чтобы сделать еще один шаг вперед. и объясните, как использовать выходные данные nest
и tidy
для создания отформатированных таблиц регрессии. Я также просмотрел эту публикацию SO, но мне трудно поверить, что нет более простого способа сделать это, которого я просто упускаю.
Примеры данных, а также мой общий подход к запуску моделей:
id <- 1:2000
gender <- sample(0:1, 2000, replace = T)
age <- sample(17:64, 2000, replace = T)
race <- sample(0:1, 2000, replace = T)
cond_a <- sample(0:1, 2000, replace = T)
cond_b <- sample(0:1, 2000, replace = T)
cond_c <- sample(0:1, 2000, replace = T)
cond_d <- sample(0:1, 2000, replace = T)
df <- data.frame(id, gender, age, race, cond_a, cond_b, cond_c, cond_d)
df %>% gather(c(cond_a, cond_b, cond_c, cond_d), key = "condition", value = "case") %>%
group_by(condition) %>% nest() %>%
mutate(model = map(data, ~glm(case ~ gender + age + race, family = "binomial", data = .)),
tidy = map(model, tidy, exponentiate = T, conf.int = T, conf.level = 0.99))