Compute odds.ratio and p.value for logit model report in stargazer

gistfile1.txt

sg_or <- function(mod){
  OR <- lapply(mod, function(x) as.data.frame(questionr::odds.ratio(x)))
  re <- list(
    or = lapply(OR, `[[`, "OR"),
    ci = lapply(OR, `[`, 2:3),
    p  = lapply(OR, `[[`, "p")
  )
  return(re)
}

# mod <- glm(..., family = "binomial")
# custom_coef <- sg_or(mod)
# stargazer(mod, ci = TRUE, 
#           coef = custom_coef$or, 
#           ci.custom = custom_coef$ci, 
#           p = custom_coef$p)