diff --git a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
index 588534c0..3039e5e2 100644
--- a/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
+++ b/qhtcp-workflow/apps/r/calculate_interaction_zscores.R
@@ -825,10 +825,9 @@ generate_rank_plot_configs <- function(df_filtered, variables, is_lm = FALSE) {
         rectangles <- NULL
       }
 
-      print(df_filtered %>% select(all_of(c(x_var, y_var))))
-
       # Check if there is sufficient variation
       if (length(unique(df_filtered[[x_var]])) < 2 || length(unique(df_filtered[[y_var]])) < 2) {
+        print(df_filtered %>% select(all_of(c("scan", "Plate", "Row", "Col", "num", "conc_num", x_var, y_var))))
         message("Not enough variation in ", x_var, " or ", y_var, " to fit linear model. Skipping.")
         next
       }