R/sCBellEE.R

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

sCBellEE<-function (x,  alpha, beta, lambda , log.p = FALSE, lower.tail = TRUE) 
	{
	G=(1-exp(-alpha*x))^beta
	cdf <- x
    cdf[log.p == FALSE & lower.tail == TRUE] <-  (exp(exp(lambda*G)-1)-1)/(exp(exp(lambda)-1)-1)        
    cdf[log.p == TRUE & lower.tail == TRUE]  <-  log(exp(exp(lambda*G)-1)-1)-log(exp(exp(lambda)-1)-1)  
    cdf[log.p == FALSE & lower.tail == FALSE] <- (exp(exp(lambda)-1)-exp(exp(lambda*G)-1))/((exp(exp(lambda)-1)-1))  
    cdf[log.p == TRUE & lower.tail == FALSE] <-  log(exp(exp(lambda)-1)-exp(exp(lambda*G)-1))-log((exp(exp(lambda)-1)-1))
	sf<-1-cdf
    return(sf)
	}