R/pegamma.R

In expectreg: Expectile and Quantile Regression

pegamma <-
function (e, shape, rate = 1, scale = 1/rate) 
{
    if (scale != 1/rate) 
        rate = 1/scale
    u = gamma(1 + shape)/gamma(shape)/rate * pgamma(e, 1 + shape, 
        rate) - e * pgamma(e, shape, rate)
    asy = u/(2 * u + e - shape/rate)
    return(asy)
}