require(dplyr, quietly = TRUE)
data(Vanderpump)
Vanderpump %>%
select(vstatus, agegrp, smoker) %>%
tbl_strata(
strata = smoker,
.tbl_fun = ~ .x %>%
tbl_summary(by = agegrp)
) %>%
cosm_sum(bold = TRUE, head_label = " ")
mhor(vstatus ~ agegrp / smoker, data = Vanderpump)
Run the code above in your browser using DataLab