Nothing
#' @importFrom stats dgamma
discrete_gamma <- function(using_terms, max_size=10, shape=10, rate=4) {
n_terms <- length(using_terms)
dist <- dgamma(seq(max_size), shape=shape, rate=rate)
dens <- log(dist/sum(dist))
function(phi) {
x <- length(phi)
if (x > max_size) -Inf
else dens[x] - lchoose(n_terms, x)
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.