Nothing
jll_normal <- function(params, X, y) {
p <- length(params)
beta <- params[-p]
sigma <- exp(params[p])
linpred <- X %*% beta
sum(dnorm(y, mean = linpred, sd = sigma, log = TRUE))
}
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.