stargazer3 <- function(model, odds.ratios = F, origin_model = NULL, ...) {
  if(!("list" %in% class(model))) model <- list(model)
  if(!("list" %in% class(origin_model))) origin_model <- list(origin_model) 
  if (odds.ratios) {
    coefOR2 <- lapply(model, function(x) exp(x[,1]))
    #seOR2 <- lapply(model, function(x) exp(x[,1]) * exp(x[,2]))
    p2 <- lapply(model, function(x) x[,4])
    obs <- lapply(origin_model, function(x) nobs(x))
    LL <- lapply(origin_model, function(x) round(logLik(x),3))
    AIC <- lapply(origin_model, function(x) round(AIC(x),3))
    stargazer(model, coef = coefOR2,  p = p2,
              star.char = c("+", "*", "**", "***"),
              star.cutoffs = c(.1, .05, .01, .001),
              notes = "+ p < 0.1; * p < 0.05; ** p < 0.01; *** p < 0.001",
              notes.append = FALSE,
              add.lines=list(c("Observations", unlist(obs)),
                             c("Log Likelihood", unlist(LL)),
                             c("Akaike Inf. Crit.", unlist(AIC))), ...)
  } else {
    stargazer(model, ...)
  }
}