R/dgp.R

In gp: Maximum Likelihood Estimation of the Generalized Poisson Distribution

Defines functions dgp

Documented in dgp

dgp <- function(y, theta, lambda, logged = TRUE) {

  den <-  log( theta ) + (y - 1) * log( theta + lambda * y ) - theta -
        y * lambda - lgamma(y + 1)
  if ( !logged )  den <- exp(den)
  den

}

Try the gp package in your browser