Residualize function

Residualize_func.txt

cols_to_res <- colnames(cog)[str_detect(colnames(cog), "_T")]

col_res <- function(df, response_var, predictor_string, model_func = lm){
  # predictor_string <- paste(predictor_vars, collapse = " + ")
  # Construct the full formula dynamically
  formula_string <- paste(response_var, "~", predictor_string)
  print(formula_string)
  
  current_formula <- as.formula(formula_string)
  
  m <- model_func(formula = current_formula, data = df, na.action = na.exclude)
  
  return_df <- data.frame(as.numeric(resid(m)))
  names(return_df) <- paste(response_var)
  return(return_df)
}

cog.r <- cols_to_res %>% 
  map_dfc(~col_res(cog, ., "1 + site_id_l", model_func = lm)) 

cog.r <- cbind(cog[, !..res_col], cog.r)