old_scripts/plots/population_pyramids_2001_census.R

In EU-ECDC/HerpesZosterModel: ZM

# Load required packages
library(eurostat) # For downloading data
library(dplyr)
library(ggplot2) # For plotting
library(scales)

# Load data
# Eurostat data "Population by sex, age group and country of citizenship"
cens_01nsctz <- get_eurostat(id = "cens_01nsctz")
cens_01nsctz$age <- factor(cens_01nsctz$age,
                           levels = c("Y_LT5", "Y5-9", "Y10-14", "Y15-19",
                                      "Y20-24", "Y25-29", "Y30-34", "Y35-39",
                                      "Y40-44", "Y45-49", "Y50-54", "Y55-59",
                                      "Y60-64", "Y65-69", "Y70-74", "Y75-79",
                                      "Y80-84", "Y_GE85", "TOTAL", "UNK"))

#cens_01nsctz %>% glimpse()
cens_01nsctz <- cens_01nsctz %>%
  filter(sex != "T", age != "TOTAL", age != "UNK")

theme_set(theme_bw())
p <- ggplot(data = cens_01nsctz, 
            mapping = aes(x = age, fill = sex, 
                          y = ifelse(test = sex == "M", 
                                     yes = -values, no = values))) +
  geom_bar(stat = "identity") +
  scale_y_continuous(labels = abs, limits = max(cens_01nsctz$values) * c(-1, 1)) +
  labs(y = "Population", x = "Age") +
  coord_flip() +
  facet_wrap(. ~ geo, scales = "free_x") +
  theme(axis.text.x  = element_text(angle = 90, vjust = 1)) +
  scale_fill_grey(start = 0.6, end = 0.4) +
  scale_x_discrete(labels = function(x){format(x, scientific = FALSE)})

# Save output
tiff(filename = "./figures/pyramids_2001cens.tif",
     width = 1200, height = 800)
p
dev.off()