tanu/LSHTM_analysis
1
0
Fork 0
Code Issues Pull requests Releases Wiki Activity

saving before commiting

Browse source
This commit is contained in:
Tanushree Tunstall 2022-01-31 17:36:00 +00:00
parent ea931b59f3
commit a287b801f7
2 changed files with 12 additions and 12 deletions
Download patch file Download diff file Expand all files Collapse all files

9
scripts/functions/corr_plot_df.R
View file

@ -7,6 +7,7 @@
# merged_df3 or merged_df2!? # merged_df3 or merged_df2!?
# gene: [sanity check] # gene: [sanity check]
# drug: relates to a column name that will need to extracted # drug: relates to a column name that will need to extracted
# ligand_dist_colname = LigDist_colname (variable from plotting_globals()
#colnames_to_extract = c("mutationinformation" #colnames_to_extract = c("mutationinformation"
# , "duet_affinity_change") # , "duet_affinity_change")
@ -18,11 +19,11 @@
# TO DO: SHINY # TO DO: SHINY
#1) Corr type? #1) Corr type?
#2) #2)
################################################################## ##################################################################
corr_data_extract <- function(corr_plot_df corr_data_extract <- function(corr_plot_df
#, gene_name = gene #, gene_name = gene
, drug_name = drug , drug_name = drug
, ligand_dist_colname = LigDist_colname
, colnames_to_extract , colnames_to_extract
, colnames_display_key , colnames_display_key
, extract_scaled_cols = F){ , extract_scaled_cols = F){
@ -44,7 +45,7 @@ corr_data_extract <- function(corr_plot_df
, "duet_stability_change" , "duet_stability_change"
, "ligand_affinity_change" , "ligand_affinity_change"
#, "ligand_distance" #, "ligand_distance"
, LigDist_colname , ligand_dist_colname
, "ddg_foldx" , "ddg_foldx"
, "deepddg" , "deepddg"
, "asa" , "asa"
@ -78,7 +79,7 @@ corr_data_extract <- function(corr_plot_df
colnames_display_key = c(duet_stability_change = "DUET" colnames_display_key = c(duet_stability_change = "DUET"
, ligand_affinity_change = "mCSM-lig" , ligand_affinity_change = "mCSM-lig"
#, ligand_distance = "ligand_distance" #, ligand_distance = "ligand_distance"
#, LigDist_colname = "ligand_distance" #, ligand_dist_colname = "ligand_distance"
, ddg_foldx = "FoldX" , ddg_foldx = "FoldX"
, deepddg = "DeepDDG" , deepddg = "DeepDDG"
, asa = "ASA" , asa = "ASA"
@ -115,7 +116,7 @@ corr_data_extract <- function(corr_plot_df
,"\nRenaming successful") ,"\nRenaming successful")
cat("\nSneak peak...") cat("\nSneak peak...")
#print(head(corr_df)) print(head(corr_df))
# Move drug column to the end # Move drug column to the end
last_col = colnames(corr_df[ncol(corr_df)]) last_col = colnames(corr_df[ncol(corr_df)])

15
scripts/plotting/get_plotting_dfs.R
View file

@ -127,19 +127,19 @@ cat(s1)
#source("other_plots_data.R") #source("other_plots_data.R")
# FIXME # FIXME
source(paste0(plot_script_path, "dm_om_data.R")) #source(paste0(plot_script_path, "dm_om_data.R"))
s2 = c("\nSuccessfully sourced other_plots_data.R") #s2 = c("\nSuccessfully sourced other_plots_data.R")
cat(s2) #cat(s2)
#################################################################### ####################################################################
# Data for Lineage barplots: WF and LF dfs # Data for Lineage barplots: WF and LF dfs
#################################################################### ####################################################################
source(paste0(plot_script_path, "lineage_data.R")) #source(paste0(plot_script_path, "lineage_data.R"))
s3 = c("\nSuccessfully sourced lineage_data.R") #s3 = c("\nSuccessfully sourced lineage_data.R")
cat(s3) #cat(s3)
#################################################################### ####################################################################
# Data for corr plots: # Data for corr plots:
@ -184,8 +184,7 @@ vars1 = ls(envir = .GlobalEnv)[grepl("^cols_to*", ls(envir = .GlobalEnv))]
vars2 = ls(envir = .GlobalEnv)[grepl("pivot_cols_*", ls(envir = .GlobalEnv))] vars2 = ls(envir = .GlobalEnv)[grepl("pivot_cols_*", ls(envir = .GlobalEnv))]
vars3 = ls(envir = .GlobalEnv)[grepl("expected_*", ls(envir = .GlobalEnv))] vars3 = ls(envir = .GlobalEnv)[grepl("expected_*", ls(envir = .GlobalEnv))]
rm(c1 rm( fact_cols
, fact_cols
, infile_metadata , infile_metadata
, infile_params , infile_params
, vars0 , vars0