R/ebeta.R

ebeta <-
function (asy, a = 1, b = 1) 
{
    asy[asy > 1 | asy < 0] = NA
    zz = 0 * asy
    lower = rep(0, length(asy))
    upper = rep(1, length(asy))
    diff = 1
    index = 1
    while (diff > 1e-10 && index < 1000) {
        root = pebeta(zz, a, b) - asy
        root[is.na(root)] = 0
        lower[root < 0] = zz[root < 0]
        upper[root > 0] = zz[root > 0]
        zz = (upper + lower)/2
        diff = max(abs(root), na.rm = T)
        index = index + 1
    }
    zz[is.na(asy)] = NA
    return(zz)
}

Try the expectreg package in your browser

Any scripts or data that you put into this service are public.

expectreg documentation built on March 18, 2022, 5:57 p.m.