timtaylor / ageutils Goto Github PK

ageutils::reaggregate_interval_counts(
    lower_bounds = 0:100,
    upper_bounds = c(1:100, Inf),
    counts = runif(101, 100, 1000),
    max_upper = 100,
    breaks = c(0, 10, 50)
)
#> Warning in ageutils::reaggregate_interval_counts(lower_bounds = 0:100,
#> upper_bounds = c(1:100, : `upper_bounds` greater than `max_upper` (100) have
#> been replaced prior to splitting.
#> Error in ageutils::reaggregate_interval_counts(lower_bounds = 0:100, upper_bounds = c(1:100, : `lower_bounds` must be less than `upper_bounds`.

tmp <- c("[0, 1)", "[10, 15)", "[2, 9)")

strcapture(
    pattern = "([[:digit:]]+),[[:space:]]*([[:digit:]]+)",
    x = tmp,
    proto = data.frame(lower=integer(),upper=integer()),
    perl = TRUE
)
#>   lower upper
#> 1     0     1
#> 2    10    15
#> 3     2     9

lu <- 95
 lst <- qs::qread(file = "~/Downloads/list.qs.gz")

  u_f <- function(v) {
    dplyr::lead(v) |>
      dplyr::coalesce(lu)
  }

  lst$data |>
    dplyr::reframe(ageutils::reaggregate_interval_rates(
      lower_bounds = lower_bound,
      upper_bounds = u_f(lower_bound),
      rates = value,
      breaks = lst$breaks,
    ), .by = c(sex))