bryan/hartman-server
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Fix stats_bg dataframe grouping

Browse Source
This commit is contained in:
Bryan Roessler
2024-09-01 18:58:01 -04:00
parent 79175e7821
commit 3fbde2de24
1 changed files with 2 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
workflow/apps/r/calculate_interaction_zscores5.R
View File

@@ -629,7 +629,7 @@ main <- function() {
write.csv(stats_bg, write.csv(stats_bg,
file = file.path(out_dir, paste0("SummaryStats_BackgroundStrains_", strain, ".csv")), file = file.path(out_dir, paste0("SummaryStats_BackgroundStrains_", strain, ".csv")),
row.names = FALSE) 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 # Filter L values within and outside 2SD of K
results_2sd <- calculate_l_2sd_of_k(df_bg, stats_by_k_bg) 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) write.csv(zscores_interactions, file = file.path(out_dir, "ZScores_Interaction.csv"), row.names = FALSE)
# Define conditions for enhancers and suppressors # Define conditions for enhancers and suppressors
# TODO Add to study config file?
enhancer_condition_L <- zscores_interactions$Avg_Zscore_L >= 2 enhancer_condition_L <- zscores_interactions$Avg_Zscore_L >= 2
suppressor_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 enhancer_condition_K <- zscores_interactions$Avg_Zscore_K >= 2