confidence_interval.R

fun <- function(d, est, se, p, expo = FALSE, dp = "%0.2f") {
  mutate(d,
         ci.lo = {{est}} - {{se}} * qnorm((1 + 0.95)/2),
         ci.hi = {{est}} +  {{se}} * qnorm((1 + 0.95)/2),
         p.stars = case_when({{p}} < 0.001 ~ "***",
                             {{p}} < 0.01 ~ "**",
                             {{p}} < 0.05 ~ "*",
                             TRUE ~ ""),
         across(c({{est}}, "ci.lo", "ci.hi"), ~case_when(expo ~ exp(.), TRUE ~ .)),
         across(c({{est}}, "ci.lo", "ci.hi"), ~sprintf(dp, .)),
         txt = paste0({{est}}, " (95% CI: ", ci.lo, " to ", ci.hi, ")", p.stars))
}