From da3b23d50266b5bfc281dcae16cc12107074532b Mon Sep 17 00:00:00 2001 From: Tanushree Tunstall Date: Tue, 6 Oct 2020 19:12:32 +0100 Subject: [PATCH] indicated hardcoded active site residues for pnca --- scripts/data_selective_pressure.R | 1 + 1 file changed, 1 insertion(+) diff --git a/scripts/data_selective_pressure.R b/scripts/data_selective_pressure.R index 95730a8..eae722e 100644 --- a/scripts/data_selective_pressure.R +++ b/scripts/data_selective_pressure.R @@ -38,6 +38,7 @@ dnds_df_all = merged_df2[, cols_to_select] str(dnds_df_all) is.numeric(dnds_df_all$position) +# !!!!BEWARE HARDCODED value for !!! active_site_positions = c(8, 13, 49, 51, 57, 68, 71, 96, 103, 133, 134, 137, 138) aa_site_i = which(dnds_df_all$position%in%active_site_positions)