R/qCB12Geo.r

In ComRiskModel: Fitting of Complementary Risk Models

#' @export
qCB12Geo<-function(p, a, b, k, lambda,log.p = FALSE, lower.tail = TRUE){
	if (log.p == TRUE)
        p <- exp(p)
    if (lower.tail == FALSE)
        p <- 1 - p
	qf <- rep(NaN, length(p))
	t=(p[p >= 0 & p <= 1]/(p[p >= 0 & p <= 1]*lambda+1-lambda))

	qf[p >= 0 & p <= 1] <- a*(((1-t)^(-1/k) -1)^(1/b))
    return(qf)
	}