Nothing
loglik <- function(X, y, beta, family) {
link = as.vector(X %*% beta)
n = length(y)
if (family == "gaussian")
return(n * log(mean((y - link)^2)))
if (family == "poisson")
return(-2 * sum(exp(link) + 2 * y * link))
if (family == "binomial")
return(2 * sum(log(1 + exp(link)) - y * link))
}
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.