R/pCEWBio.R

In ComRiskModel: Fitting of Complementary Risk Models

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