rotation_group_diff.R

library(readabs)
library(tidyverse)

gf_cur <- read_lfs_grossflows("current")
gf_prev <- read_lfs_grossflows("previous")

ur_incoming <- gf_cur |> 
  filter(lfs_previous == "Incoming rotation group") |> 
  group_by(lfs_current, date) |> 
  summarise(value = sum(persons)) |> 
  pivot_wider(names_from = lfs_current) |> 
  mutate(ur = Unemployed / (`Employed full-time` + `Employed part-time` + `Unemployed`)) |> 
  select(date, ur) |> 
  mutate(series = "incoming")

ur_outgoing <- gf_prev |> 
  filter(lfs_current == "Outgoing rotation group") |> 
  group_by(lfs_previous, date) |>
  summarise(value = sum(persons)) |>
  pivot_wider(names_from = lfs_previous) |>
  mutate(ur = Unemployed / (`Employed full-time` + `Employed part-time` + `Unemployed`)) |>
  select(date, ur) |>
  mutate(series = "outgoing") 

ur <- ur_incoming |> 
  bind_rows(ur_outgoing) |> 
  pivot_wider(names_from = series,
              values_from = ur)  |> 
  mutate(diff = incoming - outgoing,
         abs_diff = abs(diff))

ur |> 
  mutate(pctile = percent_rank(abs_diff))  |> 
  filter(date == max(date))