diff --git a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
index 7ff5a570..0e37d1d6 100644
--- a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
+++ b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
@@ -266,10 +266,10 @@ calculate_interaction_scores <- function(df, max_conc, variables = c("L", "K", "
       Zscore_AUC = Delta_AUC / WT_sd_AUC,
 
       # Fit linear models and store in list-columns
-      gene_lm_L = list(lm(Delta_L ~ conc_num_factor, data = cur_data())),
-      gene_lm_K = list(lm(Delta_K ~ conc_num_factor, data = cur_data())),
-      gene_lm_r = list(lm(Delta_r ~ conc_num_factor, data = cur_data())),
-      gene_lm_AUC = list(lm(Delta_AUC ~ conc_num_factor, data = cur_data())),
+      gene_lm_L = list(lm(Delta_L ~ conc_num_factor, data = pick(everything()))),
+      gene_lm_K = list(lm(Delta_K ~ conc_num_factor, data = pick(everything()))),
+      gene_lm_r = list(lm(Delta_r ~ conc_num_factor, data = pick(everything()))),
+      gene_lm_AUC = list(lm(Delta_AUC ~ conc_num_factor, data = pick(everything()))),
 
       # Extract coefficients using purrr::map_dbl
       lm_intercept_L = map_dbl(gene_lm_L, ~ coef(.x)[1]),