R/dgpois.R

In jlivsey/LatentGaussCounts: Integer-valued Time Series Models based on Gaussian Copula

#' pdf of generalized poisson
#'
#' @param x series
#' @param lambda 1 parameter
#' @param eta 2 parameter
#'
#' @return pdf
#' @export
#'
#'

dgpois <- function(x,lambda,eta){
  tmp = log(lambda) + (x-1)*log(lambda+eta*x) -lgamma(x+1) - lambda - eta*x
  return(exp(tmp))
}
# # GenPoisson density from VGAM package--I need to check if this is right
# # code below uses lambda instead of eta and theta instead of eta