#' @export
rvngamma <- function (n=1, shape, rate = 1, scale = 1/rate) {
if (any(shape < 0, na.rm=TRUE)) {
stop("shape parameter must be nonnegative")
}
shape <- (1/3 + shape)
rvvapply(rgamma, n.=n, shape=shape, scale=scale)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.