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

saving

Browse source
This commit is contained in:
Tanushree Tunstall 2022-09-06 21:35:59 +01:00
parent 23bb087ea3
commit ade1739753
3 changed files with 27 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

24
scripts/plotting/plotting_thesis/embb/embb_appendix_tables.R
View file

@ -1,5 +1,5 @@
#!/usr/bin/env Rscript
#source("~/git/LSHTM_analysis/config/katg.R")
source("~/git/LSHTM_analysis/config/embb.R")
#source("~/git/LSHTM_analysis/scripts/plotting/get_plotting_dfs.R")
#=======
@ -150,8 +150,12 @@ nrow(Out_df_ligS)
Out_df_ligS$drug_site = ifelse(Out_df_ligS$position%in%aa_pos_drug, "drug", "no")
table(Out_df_ligS$drug_site)
Out_df_ligS$heme_site = ifelse(Out_df_ligS$position%in%aa_pos_hem, "heme", "no")
table(Out_df_ligS$heme_site)
Out_df_ligS$dsl_site = ifelse(Out_df_ligS$position%in%aa_pos_dsl, "dsl", "no")
table(Out_df_ligS$dsl_site)
Out_df_ligS$cdl_site = ifelse(Out_df_ligS$position%in%aa_pos_cdl, "cdl", "no")
table(Out_df_ligS$cdl_site)
#--------------------
# write output file: KS test within grpup
@ -219,8 +223,11 @@ nrow(Out_df_ncaS)
Out_df_ncaS$drug_site = ifelse(Out_df_ncaS$position%in%aa_pos_drug, "drug", "no")
table(Out_df_ncaS$drug_site)
Out_df_ncaS$heme_site = ifelse(Out_df_ncaS$position%in%aa_pos_hem, "heme", "no")
table(Out_df_ncaS$heme_site)
Out_df_ncaS$dsl_site = ifelse(Out_df_ncaS$position%in%aa_pos_dsl, "dsl", "no")
table(Out_df_ncaS$dsl_site)
Out_df_ncaS$cdl_site = ifelse(Out_df_ncaS$position%in%aa_pos_cdl, "cdl", "no")
table(Out_df_ncaS$cdl_site)
#--------------------
# write output file: KS test within grpup
@ -446,8 +453,11 @@ nrow(combined_table)
combined_table$drug_site = ifelse(combined_table$position%in%aa_pos_drug, "drug", "no")
table(combined_table$drug_site)
combined_table$heme_site = ifelse(combined_table$position%in%aa_pos_hem, "heme", "no")
table(combined_table$heme_site)
combined_table$dsl_site = ifelse(combined_table$position%in%aa_pos_dsl, "dsl", "no")
table(combined_table$dsl_site)
combined_table$cdl_site = ifelse(combined_table$position%in%aa_pos_cdl, "cdl", "no")
table(combined_table$cdl_site)
#--------------------
# write output file: KS test within grpup