This might be a very easy question, but I can't really find it. How can I sample a random integer between 1 and N in Rcpp? I know I can sample a uniform number and round it, or include the sample
function to use in the C++ code, but both don't seem like the correct way of doing this.
AFIACR we do not have a "decent" (ie vectorised) sample()
function in Rcpp-land, but Christian Gunning is working on one based on Armadillo via RcppArmadillo.
You can of course sample from U(0,1), scale to N and to the rest. We have R::runif()
for the scalar case and runif()
in sugar which is vectorised.
Edit: As of 2017, both RcppArmadillo and Rcpp now have a sample()
function.
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With