R/pCLGeo.R

In ComRiskModel: Fitting of Complementary Risk Models

#' @export
pCLGeo<-function (x, b, q, lambda , log.p = FALSE, lower.tail = TRUE)
	{
	G=(1-((1+(x/b)))^(-q))
	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)
	}