Nothing
#' @importFrom stats rcauchy
# wrapped cauchy distribution with rho in (-1,1)
rwrpcauchy<-function (n, location = 0, rho = exp(-1))
{
if(rho == 0)
result <- runif(n, 0, 2 * pi)
else if(abs(rho)==1)
result <- rep(location - pi * (rho<0), n)%%(2*pi)
else {
if(rho < 0){
rho <- abs(rho)
location <- location - pi
}
scale <- -log(rho)
result <- rcauchy(n, location, scale)%%(2 * pi)
}
result
}
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.