R/dCBellEE.R

In BGFD: Bell-G and Complementary Bell-G Family of Distributions

dCBellEE<-function (x, alpha,beta,lambda, log = FALSE)
{
  G=(1-exp(-alpha*x))^beta
g=alpha*exp(-alpha*x)*beta*(1-exp(-alpha*x))^(beta-1)

  pdf <- x
  pdf[log == FALSE] <- lambda*g*exp(lambda*G)*exp(exp(lambda*G)-1)/(exp(exp(lambda)-1)-1)
  pdf[log == TRUE] <- log(lambda)+log(g)+(lambda*G)+(exp(lambda*G)-1)-log((exp(exp(lambda)-1)-1))
  return(pdf)
}