Tidyverse approach to fitting multiple models using data frame of formulas

multiple_models_tidy_2.R

library(tidyverse)
library(palmerpenguins) # for `penguins` data set

penguin_model_formulas <- tribble(
  ~name,               ~formula,
  "depth only",        "bill_length_mm ~ bill_depth_mm",
  "species only",      "bill_length_mm ~ species",
  "depth and species", "bill_length_mm ~ bill_depth_mm + species",
)

# fit a model for each formula listed above
penguin_models <- penguin_model_formulas %>%
  rowwise(name, formula) %>%
  summarise(model = list(lm(formula, data = penguins))) %>%
  ungroup()

# regression coefficients
library(broom)
penguin_models %>%
  rowwise(name) %>%
  summarise(tidy(model, conf.int = TRUE)) %>%
  ungroup()

# overall tests
library(emmeans)
penguin_models %>%
  rowwise(name) %>%
  summarise(joint_tests(model)) %>%
  ungroup()