R/aprior.R

aprior <- function(gamma.hat) {
    m = mean(gamma.hat)
    s2 = var(gamma.hat)
    (2 * s2 + m^2)/s2
}
BowenNCSU/xMSanalyzer documentation built on May 24, 2019, 7:36 a.m.