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: [EMAIL PROTECTED]

The information transmitted in this e-mail is intended only ...{{dropped}}

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Reply via email to