Nothing
rST <-
function(n, df=10, skew=1)
{
zstar <- rt(n=n,df=df)
weight <- skew/(skew + 1/skew)
z <- runif(n, -weight, 1 - weight)
signz <- sign(z)
epsilon <- skew^signz
result <- -abs(zstar)/epsilon * signz
return(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.