Rcpp Bindings for Ziggurat RNGs
The RcppZiggurat package regroups a number of implementations of the ‘Ziggurat’ random number generator by Marsaglia and Tsang (JSS, 2000), both in the original version and via the update by Leong, Zhang, Lee, Luk and Villasenor (JSS, 2005).
The purpose of the package is to both study and time these implementations in order to provide a fast N(0,1) generator for use by R. The generators included in R are excellent and very diligently implemented, but not generally the fastest. This package aims to close the performance gap.
You should know what you are doing when deploying this package. If in doubt, stick with the default generators.
The package contains a detailed vignette summarizing the methods, the implementation as well as the statistical properties (examined via three different tests) and the generator speed.
The package currently contains six different implementations. Two are standard from the literature: the original Marsaglia and Tsang approach (which should not be used as the resulting generator has be improved as suggested by Leong et al.), the improved Leong et al variant, a version from the GNU GSL, a version from the GNU Gretl econometric program as well as a version from QuantLib.
The package contains several demo scripts. Try
R> demo(package="RcppZiggurat")More testing never hurts, …
More generators as e.g. the modified Ziggurat implementation suggested by Doornik (though his code is not under an open source licence)
Dirk Eddelbuettel
GPL (>= 2)