Plot of model

grid <- expand_grid(
  launch_angle = seq(15, 40, length = 50),
  launch_speed = seq(90, 110, length = 50)
)
hats <- fit_23 |>
  predict(newdata = grid, type = "response")
grid <- grid |>
  mutate(prob = hats)

ggplot(grid) +
  geom_contour_fill(
    aes( x = launch_angle, y = launch_speed, z = prob),
    breaks = c(0, .1, .3, .5, .7, .9, 1),
    linewidth = 1.5
  ) +
  scale_fill_viridis_c(option = "H") + 
  theme(text = element_text(size = 18)) +
  labs(x = "Launch Angle", y = "Launch Speed") +
  guides(fill = guide_legend(title = "Prob (HR)"))