Simon Frost
2006-May-20 19:05 UTC
[R] Generating random numbers from skew power exponential distribution?
Dear R Help List, I'm trying to do some MCMC work, fitting a Bayesian version of the skew exponential power (sometimes called the skew power exponential) distribution, for which I need to draw random numbers. The implementation of rSEP in the gamlss library seems a little clunky (and uses integration); is there a simple/efficient way to generate random numbers from this distribution, say using a rejection scheme, as there is for the power exponential and the skew normal distribution? Thanks! Simon Simon D.W. Frost, M.A., D.Phil. Department of Pathology University of California, San Diego Antiviral Research Center 150 W. Washington St., Suite 100 San Diego, CA 92103 USA Tel: +1 619 543 8080 x275 Fax: +1 619 298 0177 Email: sdfrost at ucsd.edu The information transmitted in this e-mail is intended only ...{{dropped}}