R/bprior.R

bprior <- function(gamma.hat) {
    m = mean(gamma.hat)
    s2 = var(gamma.hat)
    (m * s2 + m^3)/s2
}
yufree/xMSanalyzer documentation built on May 4, 2019, 6:35 p.m.