andrewheiss

library(tidyverse)
library(patchwork)
library(latex2exp)

logit_df <- tibble(x = seq(0, 100, length.out = 101),
                   logits = seq(-4, 4, length.out = 101)) |> 
  mutate(odds = exp(logits)) |> 
  mutate(probs = plogis(logits))

p1 <- ggplot(logit_df, aes(x = x, y = probs)) +