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