# computes log of Beta-Binomial Prior for model size
computeModelPrior <- function(pgamma, p, alpha, beta) {
lbeta(pgamma + alpha, p - pgamma + beta) - lbeta(alpha, beta)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.