Obtaining Incidence / Prevalence
inc_prev_stroke_5y <- hmsidwR::incprev_stroke %>% filter(year == 2019) %>% select(measure, sex, age, val) %>% pivot_wider(names_from = "measure", values_from = "val")
Prevalence-based YLD Calculation:
YLD_by_severity <- merge(inc_prev_stroke_5y, dw_stroke) |> group_by(sex, age, avg_dw) |> reframe(prev_sev = Prevalence*severity_n, inc_sev = Incidence*severity_n, yld_p = prev_sev * avg_dw, yld_i = inc_sev * avg_dw * 28/365)