evals_train %>%
map(.f = ~lm(score ~ .x + <VARIABLES SELECTED>, data = evals_train)) %>%
map_df(.f = ~get_regression_summaries(.x)$adj_r_squared) %>%
select(-ID,
-score,
-<VARIABLE 1 SELECTED>,
-<VARIABLE 2 SELECTED>
) %>%
pivot_longer(cols = everything(),
names_to = "variable",
values_to = "adj_r_sq") %>%
slice_max(adj_r_sq)