Update color scheme for interaction plots (#650)

* Use project "divergent" color scheme

* Update final figure name

Co-authored-by: Jaclyn Taroni <jaclyn.n.taroni@gmail.com>

analyses/interaction-plots/scripts/03-plot_interactions.R

@@ -87,6 +87,10 @@ if (!is.na(opts$disease_table)){
 cooccur_file <- opts$infile
 plot_file <- opts$outfile
 
+# get root directory
+root_dir <- rprojroot::find_root(rprojroot::has_dir(".git"))
+
+
 cooccur_df <-
   readr::read_tsv(cooccur_file, col_types = readr::cols()) %>%
   dplyr::mutate(
@@ -114,6 +118,23 @@ cooccur_df <- cooccur_df %>%
   )
 
 
+# Get color palettes
+
+palette_dir <- file.path(root_dir, "figures", "palettes")
+divergent_palette <- readr::read_tsv(file.path(palette_dir, "divergent_color_palette.tsv"), 
+                                     col_types = readr::cols())
+divergent_colors <- divergent_palette %>%
+  dplyr::filter(color_names != "na_color") %>%
+  dplyr::pull(hex_codes)
+na_color <- divergent_palette %>%
+  dplyr::filter(color_names == "na_color") %>%
+  dplyr::pull(hex_codes)
+# not currently using "histologies_color_palette.tsv" as the current 
+# implementation uses `integrated_diagnosis``, which is not covered in that color scheme
+# Since we are only showing the most common subset of diseases, we don't need all colors,
+# but this should probably be updated in the future
+
+
 # create scales for consistent sizing
 # The scales will need to have opts$plotsize elements, 
 # so after getting the unique list, we concatenate on extra elements.
@@ -145,9 +166,9 @@ cooccur_plot <- ggplot(
     limits = yscale,
     breaks = unique(cooccur_df$label2)
   ) +
-  scale_fill_distiller(
-    type = "div",
-    palette = 5,
+  scale_fill_gradientn(
+    colors = divergent_colors,
+    na.value = na_color,
     limits = c(-10, 10),
     oob = scales::squish,
   ) +
@@ -213,6 +234,7 @@ disease_plot <- ggplot(
     y = "Samples with mutations",
     fill = "Diagnosis"
   ) + 
+  # TODO: update to project color scheme (currently requires translation of disease to short_histology)
   colorblindr::scale_fill_OkabeIto() + 
   scale_x_discrete(
     limits = xscale2,