diff --git a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
index 85d3e84d..fcb4f48c 100644
--- a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
+++ b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
@@ -293,10 +293,10 @@ calculate_interaction_scores <- function(df, max_conc, variables, group_vars = c
       lm_Score_K = sapply(lm_K, function(model) coef(model)[2] * max_conc + coef(model)[1]),
       lm_Score_r = sapply(lm_r, function(model) coef(model)[2] * max_conc + coef(model)[1]),
       lm_Score_AUC = sapply(lm_AUC, function(model) coef(model)[2] * max_conc + coef(model)[1]),
-      r_squared_L = sapply(lm_L, function(model) summary(model)$r.squared),
-      r_squared_K = sapply(lm_K, function(model) summary(model)$r.squared),
-      r_squared_r = sapply(lm_r, function(model) summary(model)$r.squared),
-      r_squared_AUC = sapply(lm_AUC, function(model) summary(model)$r.squared),
+      R_Squared_L = sapply(lm_L, function(model) summary(model)$r.squared),
+      R_Squared_K = sapply(lm_K, function(model) summary(model)$r.squared),
+      R_Squared_r = sapply(lm_r, function(model) summary(model)$r.squared),
+      R_Squared_AUC = sapply(lm_AUC, function(model) summary(model)$r.squared),
       Sum_Zscore_L = sum(Zscore_L, na.rm = TRUE),
       Sum_Zscore_K = sum(Zscore_K, na.rm = TRUE),
       Sum_Zscore_r = sum(Zscore_r, na.rm = TRUE),
@@ -324,19 +324,22 @@ calculate_interaction_scores <- function(df, max_conc, variables, group_vars = c
       "se_L", "se_K", "se_r", "se_AUC",
       "Raw_Shift_L", "Raw_Shift_K", "Raw_Shift_r", "Raw_Shift_AUC",
       "Z_Shift_L", "Z_Shift_K", "Z_Shift_r", "Z_Shift_AUC",
-      "WT_L", "WT_K", "WT_r", "WT_AUC", "WT_sd_L", "WT_sd_K", "WT_sd_r", "WT_sd_AUC",
-      "Exp_L", "Exp_K", "Exp_r", "Exp_AUC", "Delta_L", "Delta_K", "Delta_r", "Delta_AUC",
+      "WT_L", "WT_K", "WT_r", "WT_AUC",
+      "WT_sd_L", "WT_sd_K", "WT_sd_r", "WT_sd_AUC",
+      "Exp_L", "Exp_K", "Exp_r", "Exp_AUC",
+      "Delta_L", "Delta_K", "Delta_r", "Delta_AUC",
       "Zscore_L", "Zscore_K", "Zscore_r", "Zscore_AUC",
       "NG", "SM", "DB") %>%
     ungroup()
 
   # Also arrange results by Z_lm_L and NG
   interactions <- stats %>%
-    select("OrfRep", "Gene", "num", "Raw_Shift_L", "Raw_Shift_K",  "Raw_Shift_AUC", "Raw_Shift_r",
+    select("OrfRep", "Gene", "num",
+      "Raw_Shift_L", "Raw_Shift_K", "Raw_Shift_AUC", "Raw_Shift_r",
       "Z_Shift_L", "Z_Shift_K", "Z_Shift_r", "Z_Shift_AUC",
       "lm_Score_L", "lm_Score_K", "lm_Score_AUC", "lm_Score_r",
       "R_Squared_L", "R_Squared_K", "R_Squared_r", "R_Squared_AUC",
-      "Sum_Z_Score_L", "Sum_Z_Score_K", "Sum_Z_Score_r", "Sum_Z_Score_AUC",
+      "Sum_Zscore_L", "Sum_Zscore_K", "Sum_Zscore_r", "Sum_Zscore_AUC",
       "Avg_Zscore_L", "Avg_Zscore_K", "Avg_Zscore_r", "Avg_Zscore_AUC",
       "Z_lm_L", "Z_lm_K", "Z_lm_r", "Z_lm_AUC",
       "NG", "SM", "DB") %>%