R/rST.R

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)
}

Try the betategarch package in your browser

Any scripts or data that you put into this service are public.

betategarch documentation built on May 2, 2019, 10:15 a.m.