age_susceptibility_covid.r

library("readr")
library("dplyr")

df <- readr::read_csv("https://raw.githubusercontent.com/epiforecasts/inc2prev/master/outputs/estimates_age_ab.csv") ## nolint

td <- df |>
  dplyr::filter(name == "infections") |>
  dplyr::select(t_index, date) |>
  dplyr::distinct() |>
  dplyr::arrange(t_index)

susc <- df |>
  dplyr::filter(name == "gen_dab") |>
  dplyr::select(t_index, age_group = variable, mean, q5, q25, q75, q95) |>
  dplyr::left_join(td, by = "t_index") |>
  dplyr::mutate(
    age_group = dplyr::recode(age_group, `2-10` = "02-10")
  ) |>
  dplyr::arrange(date, age_group)

readr::write_csv(susc, "age_susceptibility_covid.csv")