R/summaryOR.R

In IRBLleida/UdBRpackage: Miscellaneous R Functions

#' Funció que donat un model logístic retorna una taula amb:
#' les estimacions dels paràmetres,
#' l'error estàndar,
#' l'estadístic de contrast del paràmetre,
#' el valor p associat,
#' la OR, i
#' l'interval de confiança del OR
#' 
#' @param model a logistic model (an object of class "glm", usually, a result of a call to glm).
#' @param level confidence level of the interval.
#' @param exact logical. Should exact interval be calculated? If FALSE confidence intervals are obtained based on asymptotic normality
#' @return A data.table with estimates from model.
#' @export
summaryOR <- function(model, level = 0.95, exact = T){
  if(exact==T) x <- round(cbind(coef(summary(model)), OR = exp(coef(model)), exp(confint(model, level = level))), 4)
  if(exact==F) x <- round(cbind(coef(summary(model)), OR = exp(coef(model)), exp(confint.default(model, level = level))), 4)
  colnames(x) <- c("Estimate", "SE", "z-value", "p-value", "OR", paste0("LB(", round(100*level, 0), "% CI)"), paste0("UB(", round(100*level, 0), "% CI)"))
  x[which(rownames(x) == "(Intercept)"), 5:7] <- NA
  return(x)
}