R/pCExpGeo.r

In ComRiskModel: Fitting of Complementary Risk Models

#' @export
pCExpGeo<-function (x,  alpha,lambda , log.p = FALSE, lower.tail = TRUE)
	{
	G=(1-exp(-alpha*x))
	cdf <- x
    cdf[log.p == FALSE & lower.tail == TRUE] <-(1-lambda)*G/(1-lambda*G)
	cdf[log.p == TRUE & lower.tail == TRUE] <-(1-lambda)*log(G)-log(1-lambda*G)
    cdf[log.p == FALSE & lower.tail == FALSE] <-(1-G)/(1-lambda*G)
    cdf[log.p == TRUE & lower.tail == FALSE] <-log(1-G)-log(1-lambda*G)
    return(cdf)
	}