GLM plot with categorical and continuous variable

[technocrat]

Great reprex. Thanks.

Try starting with

preds.RAIN %>%
  ggplot(aes(temp, exp(fit))) +
  geom_line() +
  facet_wrap(~ rainfall_binary)
  ...