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

Fix dataframe mutation order

Browse Source
This commit is contained in:
Bryan Roessler
2024-09-01 17:54:08 -04:00
parent e1f57ff7d7
commit f4c9985e40
1 changed files with 9 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -561,11 +561,11 @@ main <- function() {
max_conc <- max(df$conc_num_factor) max_conc <- max(df$conc_num_factor)
# QC # QC
# Filter the df above sd tolerance # Filter the df above sd tolerance for raw L vs K plot
df_above_tolerance <- df %>% filter(DB == 1) df_above_tolerance <- df %>% filter(DB == 1)
# Exclude variables above delta bg tolerance # Set vars above the delta background tolerance to NA
df_na <- df_above_tolerance %>% df_na <- df %>%
mutate( mutate(
L = ifelse(DB == 1, NA, L), L = ifelse(DB == 1, NA, L),
r = ifelse(DB == 1, NA, r), r = ifelse(DB == 1, NA, r),
@@ -577,10 +577,12 @@ main <- function() {
# Generate QC PDFs and HTMLs # Generate QC PDFs and HTMLs
message("Generating QC plots") message("Generating QC plots")
variables <- c("L", "K", "r", "AUC", "delta_bg") variables <- c("L", "K", "r", "AUC", "delta_bg")
# generate_and_save_plots(df, out_dir_qc, "Before_QC", variables, include_qc = TRUE) generate_and_save_plots(df, out_dir_qc, "Before_QC", variables, include_qc = TRUE)
# generate_and_save_plots(df_above_tolerance, out_dir_qc, "Raw_L_vs_K_above_delta_bg_threshold", variables, include_qc = TRUE) generate_and_save_plots(df_above_tolerance, out_dir_qc, "Raw_L_vs_K_above_delta_bg_threshold", variables, include_qc = TRUE)
# generate_and_save_plots(df_na, out_dir_qc, "After_QC", variables) generate_and_save_plots(df_na, out_dir_qc, "After_QC", variables)
# generate_and_save_plots(df_no_zeros, out_dir_qc, "No_Zeros", variables) generate_and_save_plots(df_no_zeros, out_dir_qc, "No_Zeros", variables)
rm(df, df_above_tolerance, df_no_zeros)
# Calculate summary statistics # Calculate summary statistics
message("Calculating summary statistics for all strains") message("Calculating summary statistics for all strains")