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

Hardcode group vars to simplify

Browse Source
This commit is contained in:
Bryan Roessler
2024-09-17 23:30:32 -04:00
parent fdcd97d7f1
commit 219dc15a45
1 changed files with 7 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
qhtcp-workflow/apps/r/calculate_interaction_zscores.R
View File

@@ -186,7 +186,10 @@ calculate_summary_stats <- function(df, variables, group_vars) {
return(list(summary_stats = summary_stats, df_with_stats = df_with_stats)) return(list(summary_stats = summary_stats, df_with_stats = df_with_stats))
} }
calculate_interaction_scores <- function(df, max_conc, variables, group_vars) { calculate_interaction_scores <- function(df, max_conc) {
variables <- c("L", "K", "r", "AUC")
group_vars <- c("OrfRep", "Gene", "num")
# Calculate total concentration variables # Calculate total concentration variables
total_conc_num <- length(unique(df$conc_num)) total_conc_num <- length(unique(df$conc_num))
@@ -206,6 +209,7 @@ calculate_interaction_scores <- function(df, max_conc, variables, group_vars) {
AUC = df %>% filter(conc_num == 0) %>% pull(sd_AUC) %>% first() AUC = df %>% filter(conc_num == 0) %>% pull(sd_AUC) %>% first()
) )
# Calculate per spot
stats <- calculate_summary_stats(df, stats <- calculate_summary_stats(df,
variables = variables, variables = variables,
group_vars = c("OrfRep", "Gene", "num", "conc_num", "conc_num_factor") group_vars = c("OrfRep", "Gene", "num", "conc_num", "conc_num_factor")
@@ -1282,14 +1286,14 @@ main <- function() {
ungroup() ungroup()
message("Calculating reference strain interaction scores") message("Calculating reference strain interaction scores")
reference_results <- calculate_interaction_scores(reference_strain, max_conc, group_vars = c("OrfRep", "Gene", "num")) reference_results <- calculate_interaction_scores(reference_strain, max_conc)
zscores_calculations_reference <- reference_results$calculations zscores_calculations_reference <- reference_results$calculations
zscores_interactions_reference <- reference_results$interactions zscores_interactions_reference <- reference_results$interactions
zscores_calculations_reference_joined <- reference_results$calculations_joined zscores_calculations_reference_joined <- reference_results$calculations_joined
zscores_interactions_reference_joined <- reference_results$interactions_joined zscores_interactions_reference_joined <- reference_results$interactions_joined
message("Calculating deletion strain(s) interactions scores") message("Calculating deletion strain(s) interactions scores")
deletion_results <- calculate_interaction_scores(deletion_strains, max_conc, group_vars = c("OrfRep", "Gene", "num")) deletion_results <- calculate_interaction_scores(deletion_strains, max_conc)
zscores_calculations <- deletion_results$calculations zscores_calculations <- deletion_results$calculations
zscores_interactions <- deletion_results$interactions zscores_interactions <- deletion_results$interactions
zscores_calculations_joined <- deletion_results$calculations_joined zscores_calculations_joined <- deletion_results$calculations_joined