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)