R/bprior.R

In 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)
}

Try the ez.combat package in your browser