R/pCFGeo.R

In ComRiskModel: Fitting of Complementary Risk Models

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