|
@@ -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]),
|