Nothing
#' @export
rskellam <- function(n, lambda1, lambda2=lambda1){
# Skellam random variables
if (missing(n)|missing(lambda1)) stop("first 2 arguments are required")
if (length(n)>1) n <- length(n)
lambda1 <- rep(lambda1,length.out=n)
lambda2 <- rep(lambda2,length.out=n)
oops <- !(is.finite(lambda1)&(lambda1>=0)&is.finite(lambda2)&(lambda2>=0))
if(any(oops)) warning("NaNs produced")
ret <- rep(NaN,length.out=n)
n <- n-sum(oops)
ret[!oops] <- stats::rpois(n,lambda1[!oops])-stats::rpois(n,lambda2[!oops])
ret
}
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.