From 3fbde2de241de1961e8cc89238354d669a19762f Mon Sep 17 00:00:00 2001 From: Bryan Roessler Date: Sun, 1 Sep 2024 18:58:01 -0400 Subject: [PATCH] Fix stats_bg dataframe grouping --- workflow/apps/r/calculate_interaction_zscores5.R | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/workflow/apps/r/calculate_interaction_zscores5.R b/workflow/apps/r/calculate_interaction_zscores5.R index 08f431a4..b97c6f34 100644 --- a/workflow/apps/r/calculate_interaction_zscores5.R +++ b/workflow/apps/r/calculate_interaction_zscores5.R @@ -629,7 +629,7 @@ main <- function() { write.csv(stats_bg, file = file.path(out_dir, paste0("SummaryStats_BackgroundStrains_", strain, ".csv")), row.names = FALSE) - stats_bg_joined <- left_join(df_bg, stats_bg, by = c("conc_num", "conc_num_factor")) + stats_bg_joined <- left_join(df_bg, stats_bg, by = c("OrfRep", "Gene", "conc_num", "conc_num_factor")) # Filter L values within and outside 2SD of K results_2sd <- calculate_l_2sd_of_k(df_bg, stats_by_k_bg) @@ -705,6 +705,7 @@ main <- function() { write.csv(zscores_interactions, file = file.path(out_dir, "ZScores_Interaction.csv"), row.names = FALSE) # Define conditions for enhancers and suppressors + # TODO Add to study config file? enhancer_condition_L <- zscores_interactions$Avg_Zscore_L >= 2 suppressor_condition_L <- zscores_interactions$Avg_Zscore_L <= -2 enhancer_condition_K <- zscores_interactions$Avg_Zscore_K >= 2