Population used in OWID/CDC COVID-19 weekly death rates dataset

_chart_usa_population_vaccination_status.r

library(tidyverse)
library(scales)

sf <- 2
width <- 600 * sf
height <- 335 * sf

options(vsc.dev.args = list(width = width, height = height, res = 72 * sf))

df <- read_csv("/Users/ben/Downloads/a.csv")

a <- df |>
  filter(outcome == "death", age_group == "all_ages", vaccination_status == "vaccinated") |>
  select(mmwr_week, unvaccinated_population, vaccinated_population) |>
  pivot_longer(
    cols = c(unvaccinated_population, vaccinated_population),
    names_to = "vaccination_status",
    values_to = "population"
  ) |>
  mutate(vaccination_status = str_remove(vaccination_status, "_population"))

# Stacked Area Plot
a |> ggplot(aes(x = mmwr_week, y = population, fill = vaccination_status)) +
  geom_area(alpha = 0.7) +
  geom_hline(yintercept = 333e6, linetype = "dashed", color = "red", size = 1) +
  scale_y_continuous(
    labels = scales::label_number(scale_cut = scales::cut_si("M"))
  ) +
  labs(
    x = "MMWR Week",
    y = "Population (in millions)",
    fill = "Vaccination Status",
    title = "Population used in OWID/CDC COVID-19 weekly death rates dataset",
    subtitle = "Including reference line at 333M"
  ) +
  theme_minimal() +
  theme(
    legend.position = "bottom",
    plot.title = element_text(size = 16, face = "bold"),
    plot.subtitle = element_text(size = 12)
  )