R/log_factorial.R

Defines functions log_factorial

log_factorial <- function(x) {
  m <- apply(x, c(1, 2), function(t)
    if (is.na(t)) {
      NA
    } else if (t == 0) {
      1
    } else
      sum(log(1:t)))
  return(m)
}
genevievelrobin/mimi documentation built on May 7, 2019, 10:57 a.m.