From e8734b1c4be4b9d99184b23567337703a0a90d79 Mon Sep 17 00:00:00 2001
From: Tanushree Tunstall <tanu@tunstall.in>
Date: Thu, 16 Sep 2021 12:43:36 +0100
Subject: [PATCH] sorted merged_df2 and consequently others by position in
 combining_dfs_plotting.R

---
 scripts/functions/combining_dfs_plotting.R   |  8 ++-
 scripts/functions/tests/test_bp_subcolours.R | 70 ++++++++++++++++++--
 2 files changed, 71 insertions(+), 7 deletions(-)

diff --git a/scripts/functions/combining_dfs_plotting.R b/scripts/functions/combining_dfs_plotting.R
index 848face..107c114 100644
--- a/scripts/functions/combining_dfs_plotting.R
+++ b/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
   #-----------------------
diff --git a/scripts/functions/tests/test_bp_subcolours.R b/scripts/functions/tests/test_bp_subcolours.R
index 2156e49..8866ffe 100644
--- a/scripts/functions/tests/test_bp_subcolours.R
+++ b/scripts/functions/tests/test_bp_subcolours.R
@@ -1,20 +1,78 @@
 #!/usr/bin/env Rscript      
-source("~/git/Misc/shiny/myshiny/gid_data.R")
-setwd("~/git/LSHTM_analysis/scripts/functions/")
-source("bp_subcolours.R")
+#source("~/git/Misc/shiny/myshiny/gid_data.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_outcome_colname = "foldx_outcome"
-                  , p_title = "FoldX" )
+                  , stability_colname = "mcsm_na_scaled"
+                  , stability_outcome_colname = "mcsm_na_outcome"
+                  , 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" )