Nothing
regbeta <- function(pa, ly, sly1, x, n) {
phi <- exp(pa[1]) ; b <- pa[-1]
m <- exp( tcrossprod(b, x) )
m <- m / ( 1 + m )
a1 <- m * phi ; a2 <- phi - a1
- n * lgamma(phi) + sum( lgamma(a1) ) + sum( lgamma(a2) ) - sum(a1 * ly) - phi * sly1
}
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.