R/bprior.R

In TKoscik/ez.combat: Easy ComBat Harmonization

Defines functions bprior

Documented in bprior

bprior <- function(gamma.hat){
  m <- mean(gamma_hat)
  v <- var(gamma_hat)
  output <- (m * v + m^3) / v
  return(output)
}