R/aprior.R

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