data-raw/gl_reg_deaths.R

In ONSdigital/Bayesian-demographic-accounts: Bayesian Demographic Accounts

library(dplyr)
library(assertr)

pxweb_nonbirths <- readRDS("data-raw/pxweb_nonbirths.rds")

gl_reg_deaths <- pxweb_nonbirths %>%
    filter(event == "Death") %>%
    mutate(age = time - cohort - (triangle == "Upper")) %>%
    mutate(is_valid_age = !is.na(age) & (age >= 0)) %>%
    verify(is_valid_age | (count == 0L)) %>%
    filter(is_valid_age) %>%
    filter(age <= max(age * (count > 0L))) %>%
    count(cohort, age, gender, time, wt = count, name = "count")

save(gl_reg_deaths,
     file = "data/gl_reg_deaths.rda",
     compress = "xz")