Reputation: 47551
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.
Upvotes: 9
Views: 2867
Reputation: 368201
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.
Upvotes: 8