SHA |>
filter(HC == "HC65",
MEASURE == "PARPIB",
HP == "HPTOT",
HF == "HFTOT",
LOCATION %in% c("FRA", "DEU", "KOR")) |>
left_join(SHA_var$LOCATION, by = "LOCATION") |>
year_to_date() |>
left_join(colors, by = c("Location" = "country")) |>
mutate(obsValue = obsValue/100) |>
ggplot() + geom_line(aes(x = date, y = obsValue, color = color)) +
scale_color_identity() + theme_minimal() + add_3flags +
scale_x_date(breaks = seq(1920, 2100, 5) |> paste0("-01-01") |> as.Date(),
labels = date_format("%Y")) +
scale_y_continuous(breaks = 0.01*seq(-7, 60, 0.01),
labels = scales::percent_format(accuracy = 0.01),
limits = c(0, 0.0012)) +
ylab("Epidemio. surveillance,risk & disease control (% of GDP)") + xlab("")