R/qbpareto2.R

In CaDENCE: Conditional Density Estimation Network Construction and Evaluation

qbpareto2 <-
function(p, prob, scale, shape) 
{
    if (any(p <= 0) || any(p >= 1)) stop("bad input for argument 'p'")
    if (length(prob) == 1) 
        prob <- rep(prob, length(p))
    if (length(scale) == 1) 
        scale <- rep(scale, length(p))
    if (length(shape) == 1) 
        shape <- rep(shape, length(p))
    q <- rep(0, length(p))
    cases <- p > (1 - prob)
    q[cases] <- qpareto2((prob[cases] + p[cases] - 1)/prob[cases], 
        scale = scale[cases], shape = shape[cases])
    q
}