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

added pretty colnames to corr_data.R

Browse source
This commit is contained in:
Tanushree Tunstall 2021-09-13 10:24:41 +01:00
parent dc5b5e2f11
commit 1249a08fb9
1 changed files with 53 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

53
scripts/plotting/corr_data.R
View file

@ -32,6 +32,32 @@ corr_cols_select <- c("mutationinformation", drug, "mutation_info_labels"
corr_df_m2 = merged_df2[,colnames(merged_df2)%in%corr_cols_select] corr_df_m2 = merged_df2[,colnames(merged_df2)%in%corr_cols_select]
#-----------------------
# formatting: some cols
# Add pretty colnames
#-----------------------
corr_df_m2_f <- corr_df_m2 %>%
rename(
DUET = duet_stability_change
, 'mCSM-lig' = ligand_affinity_change
, FoldX = ddg_foldx
, DeepDDG = deepddg
, ASA = asa
, RSA = rsa
, KD = kd_values
, RD = rd_values
, MAF = af
, 'Log (OR)' = log10_or_mychisq
, '-Log (P)' = neglog_pval_fisher
, Dynamut = ddg_dynamut
, 'ENCoM-DDG'= ddg_encom
, mCSM = ddg_mcsm
, SDM = ddg_sdm
, 'DUET-d' = ddg_duet
, 'ENCoM-DDS'= dds_encom
, Dynamut2 = ddg_dynamut2
, 'mCSM-NA' = mcsm_na_affinity )
#=========================== #===========================
# Corr data for plots: PS # Corr data for plots: PS
# short_df ps: ~merged_df3 # short_df ps: ~merged_df3
@ -53,6 +79,33 @@ if (nrow(corr_df_m3) == nrow(merged_df3) && nrow(merged_df3_comp) == check1) {
, "\nGot: ", check1) , "\nGot: ", check1)
} }
#-----------------------
# formatting: some cols
# Add pretty colnames
#-----------------------
corr_df_m3_f <- corr_df_m3 %>%
rename(
DUET = duet_stability_change
, 'mCSM-lig' = ligand_affinity_change
, FoldX = ddg_foldx
, DeepDDG = deepddg
, ASA = asa
, RSA = rsa
, KD = kd_values
, RD = rd_values
, MAF = af
, 'Log (OR)' = log10_or_mychisq
, '-Log (P)' = neglog_pval_fisher
, Dynamut = ddg_dynamut
, 'ENCoM-DDG'= ddg_encom
, mCSM = ddg_mcsm
, SDM = ddg_sdm
, 'DUET-d' = ddg_duet
, 'ENCoM-DDS'= dds_encom
, Dynamut2 = ddg_dynamut2
, 'mCSM-NA' = mcsm_na_affinity )
########################################################################
cat("\nCorr Data created:" cat("\nCorr Data created:"
, "\n===================================" , "\n==================================="
, "\ncorr_df_m2: created from merged_df2" , "\ncorr_df_m2: created from merged_df2"