sorted merged_df2 and consequently others by position in combining_dfs_plotting.R

scripts/functions/combining_dfs_plotting.R

@@ -153,7 +153,13 @@ combining_dfs_plotting <- function(  my_df_u
     quit()
   }
 
-  # Quick formatting: pretty labels
+  # Quick formatting: ordering df and pretty labels
+  
+  #------------------------------
+  # sorting by column: position
+  #------------------------------
+  merged_df2 = merged_df2[order(merged_df2$position), ]
+  
   #-----------------------
   # mutation_info_labels
   #-----------------------

scripts/functions/tests/test_bp_subcolours.R

@@ -1,20 +1,78 @@
 #!/usr/bin/env Rscript      
-source("~/git/Misc/shiny/myshiny/gid_data.R")
+#source("~/git/Misc/shiny/myshiny/gid_data.R")
-setwd("~/git/LSHTM_analysis/scripts/functions/")
-source("bp_subcolours.R")
 
+source("~/git/LSHTM_analysis/config/gid.R")
+source("~/git/LSHTM_analysis/scripts/plotting/get_plotting_dfs.R")
+source("~/git/LSHTM_analysis/scripts/functions/bp_subcolours.R")
 
+# p1
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "duet_scaled"
+                  , stability_outcome_colname = "duet_outcome"
+                  , p_title = "DUET" )
+
+# p2
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "foldx_scaled"
+                  , stability_outcome_colname = "foldx_outcome"
+                  , p_title = "FoldX" )
+
+# p3
 bp_stability_hmap(plotdf =  merged_df3
                   , stability_colname = "deepddg_scaled"
                   , stability_outcome_colname = "deepddg_outcome"
                   , p_title = "DeepDDG" )
 
+# p4
 bp_stability_hmap(plotdf =  merged_df3
                   , stability_colname = "ddg_dynamut2_scaled"
                   , stability_outcome_colname = "ddg_dynamut2_outcome"
                   , p_title = "Dynamut2" )
 
+# p5
 bp_stability_hmap(plotdf =  merged_df3
-                  , stability_colname = "foldx_scaled"
+                  , stability_colname = "mcsm_na_scaled"
-                  , stability_outcome_colname = "foldx_outcome"
+                  , stability_outcome_colname = "mcsm_na_outcome"
-                  , p_title = "FoldX" )
+                  , p_title = "mCSM-NA" )
+
+# p6
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "ddg_dynamut_scaled"
+                  , stability_outcome_colname = "ddg_dynamut_outcome"
+                  , p_title = "Dynamut" )
+
+# p7
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "ddg_mcsm_scaled"
+                  , stability_outcome_colname = "ddg_mcsm_outcome"
+                  , p_title = "mCSM" )
+
+# p8
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "ddg_duet_scaled"
+                  , stability_outcome_colname = "ddg_duet_outcome"
+                  , p_title = "DUET-d" )
+
+# p9
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "ddg_sdm_scaled"
+                  , stability_outcome_colname = "ddg_sdm_outcome"
+                  , p_title = "SDM" )
+
+# p10
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "ddg_encom_scaled"
+                  , stability_outcome_colname = "ddg_encom_outcome"
+                  , p_title = "ENCoM-Stability" )
+
+# p11
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "dds_encom_scaled"
+                  , stability_outcome_colname = "dds_encom_outcome"
+                  , p_title = "ENCoM-Flexibility" )
+
+# p12
+bp_stability_hmap(plotdf =  merged_df3
+                  , stability_colname = "affinity_scaled"
+                  , stability_outcome_colname = "ligand_outcome"
+                  , p_title = "mCSM-lig" )