# VIZ <- p_dat
# 
# ggplot(VIZ, aes(x = year, y = m)) +
#   geom_shadowpoint() +
#   scale_y_continuous(expand = expansion(mult = .1))

# VIZ_noref <- bind_rows(rl) %>% 
#   filter(type == "alt")

VIZ <- bind_rows(rl) %>% 
  filter(type == "alt")

ggplot(VIZ, aes(x = year, y = n, col = interaction(sex, dom), shape = type)) +
  geom_vline(xintercept = 2024.5) +
  geom_shadowpoint()

VIZ <- sh %>% 
  filter(sex == "f", type == "alt", age == 64, nat == "ch")
