ln.series |>
filter(lfst_code == 23,
sexs_code == 1,
race_code == 0,
seasonal == "U",
ages_code %in% c(0, 28, 33),
orig_code == 0,
mari_code == 0,
born_code == 0,
vets_code == 0,
education_code == 0,
periodicity_code == "M") |>
left_join(ln.data.1.AllData, by = "series_id") |>
#filter(series_id %in% c("LNS12300001", "LNS12300049", "LNS12300061")) %>%
left_join(ln.ages, by = "ages_code") |>
month_to_date() |>
arrange(desc(date)) |>
mutate(value = as.numeric(value)) |>
ggplot() + theme_minimal() + xlab("") + ylab("Men Employment / Population Ratio") +
geom_line(aes(x = date, y = value/100, color = ages_text)) +
geom_rect(data = nber_recessions |>
filter(Trough >= as.Date("1947-01-01")),
aes(xmin = Peak, xmax = Trough, ymin = -Inf, ymax = +Inf),
fill = 'grey', alpha = 0.5) +
theme(legend.position = c(0.2, 0.2),
legend.title = element_blank()) +
scale_x_date(breaks = seq(1910, 2100, 5) |> paste0("-01-01") |> as.Date(),
labels = date_format("%Y")) +
scale_y_continuous(breaks = 0.01*c(seq(0, 100, 5), seq(100, 500, 50)),
labels = percent_format(accuracy = 1, prefix = ""))