mean = mean(!!sym(response_var), na.rm = TRUE), sem = sd(!!sym(response_var), na.rm = TRUE)/sqrt(n()), cat("\n - Seed Treatment:", round(tukey_seed$statistics$MSD, 3 ...