R/dCBellL.R

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

dCBellL<-function (x, b, q, lambda, log = FALSE)
{
	G=(1-((1+(x/b)))^(-q))
	g=(q/b)*(1+(x/b))^(-q-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)
}