Nothing
beta.mom <- function(x, lower = 0.01, upper = 100) {
x.bar <- mean (x)
n <- length(x)
v <- var(x) * (n-1) / n
R <- 1 / x.bar - 1
f <- function(a){ R * x.bar^2 / (a / x.bar + 1) - v }
u <- uniroot(f, c(lower, upper))
return(c(shape1 = u$root, shape2 = u$root * R))
}
x <- rbeta(50, 2, 5); beta.mom(x)
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.