Nothing
ddirichlet <- function (x, alpha)
{
if (length(x) != length(alpha))
stop("Mismatch between dimensions of x and alpha in ddirichlet().\n")
logD <- sum(lgamma(alpha)) - lgamma(sum(alpha))
s <- sum((alpha - 1) * log(x))
pd <- exp(sum(s) - logD)
pd[any(x < 0 | x > 1)] <- 0
if(sum(x) != 1) pd <- 0
return(pd)
}
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.