Nothing
loglik.multinom <- function(theta, x) {
probs <- c(theta, 1 - sum(theta))
if (any (probs < 0)) return(NA)
dmultinom(x, size = 100, prob = probs, log = TRUE)
}
maxLik(loglik.multinom, start = rep(0.25, 3), x = c(10, 20, 30, 40)) -> ml;
coef(ml)
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.