bprior <- function(gamma.hat){ m <- mean(gamma_hat) v <- var(gamma_hat) output <- (m * v + m^3) / v return(output) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.