R/barplot_bun_comp_stack_med.R

In utah-osa/hcctools:

barplot_bun_comp_stack_med <- function(df) {
    df <- df%>%
        dplyr::rowwise() %>%
        dplyr::mutate(pr_list = pr_list %>%
                   stringr::str_split("_") %>%
                   unlist() %>%
                   unique() %>%
                   paste0(collapse=","))

    providers <- df %>% dplyr::pull(pr_list)

    vt_meds <- df %>% dplyr::pull(vt_med)

    ## add the pc list here to demarcate different provider/pc bundles.
    pc_lists <- df %>% dplyr::pull(pc_list)

    pr_pc_list <- paste0("Provider(s):",providers,"\n CPT:", pc_lists)

    surg_med <- df$surg_bun_med
    medi_med <- df$medi_bun_med
    radi_med <- df$radi_bun_med
    path_med <- df$path_bun_med
    anes_med <- df$anes_bun_med
    fac_med <- df$fac_bun_med

    component_meds <-
        rbind(surg_med, medi_med, path_med, radi_med, anes_med, fac_med) %>%
        tidyr::replace_na(0) %>%
        as.data.frame()

    component_meds$ID <- as.factor(1:nrow(component_meds))

    # fill_cols <- c(
    #     "Anesthesia" = "#003f5c",
    #     "Facility" = "#444e86",
    #     "Medicine Proc" = "#955196",
    #     "Pathology" = "#dd5182",
    #     "Radiology" = "#ff6e54",
    #     "Surgery" = "#ffa600"
    # )

    melted <- reshape2::melt(component_meds, id = "ID") %>%
        dplyr::mutate(
            ID = dplyr::case_when(
                ID == 1 ~ "Surgery",
                ID == 2 ~ "Medicine Proc",
                ID == 3 ~ "Pathology",
                ID == 4 ~ "Radiology",
                ID == 5 ~ "Anesthesia",
                ID == 6 ~ "Facility"
            )
        )

    melted <-  melted %>%
        dplyr::rowwise()%>%
        dplyr::mutate(variable = pr_pc_list %>%
                          purrr::pluck(
                                    as.numeric(
                                        str_extract(
                                            as.character(variable),
                                            "[[:digit:]]+")
                                        )
                                    )
                      )

    ggplot2::ggplot(melted, ggplot2::aes(variable, value, fill = ID, group = ID)) +
        ggplot2::geom_bar(stat = "identity", position = "stack") +
        ggplot2::scale_fill_manual(values = cust_color_pal1) +
        ggplot2::ggtitle("Bundles (with component) by Provider ") +
        ggplot2::xlab("Component Bundle") +
        ggplot2::ylab("Price ($)") +
        ggplot2::scale_y_continuous(labels = scales::dollar_format(prefix = "$")) +
        ggplot2::coord_flip()+
        ggplot2::theme(

            axis.text.x = ggplot2::element_text(angle = -45, hjust = 1, vjust=1)
            #axis.title.x=element_blank(),
            #     axis.text.x = element_blank(),
            #     axis.ticks.x = element_blank()
        )
}