R/getKcutoff.R

getKcutoff <-
function (m, n, s = 0.45014) 
{
    mu = (sqrt(m - 0.5) + sqrt(n - 0.5))^2
    sigma = (sqrt(m - 0.5) + sqrt(n - 0.5)) * (1/sqrt(m - 0.5) + 
        1/sqrt(n - 0.5))^(1/3)
    return(sqrt(mu/n + s * sigma/n))
}

Try the ruv package in your browser

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

ruv documentation built on Aug. 31, 2019, 1:04 a.m.