R/igamma.R

igamma <- function (s, u, log = FALSE) {
  robj <- pgamma(u, shape = s, rate = 1, log.p = TRUE) + lgamma(s)
  if (!log) robj <- exp(robj)
  robj
}
igollini/tailloss documentation built on May 18, 2019, 3:40 a.m.