VIZ <- 
  left_join(filter(pen_n, type == "alt") %>% 
              dplyr::summarize(n = sum(n), .by = c("year", "nat", "age")),
            filter(pop) %>% 
              dplyr::summarize(pop = sum(pop), .by = c("year", "nat", "age")), 
  by = c("year", "nat", "age")) %>% 
  filter(n > 0, year == 2023, age == 71)
