Your friendly neighbourhood co-author of RcppArmadillo
here: I can assure you that it does not provide rmultinom
, but Rcpp does. In fact, it simply passes through to R itself as a quick grep
would have told you:
inline void rmultinom(int n, double* prob, int k, int* rn)
{ return ::rmultinom(n, prob, k, rn); }
So I would suggest your first write a five-line C program against the R API to make sure you know how to have rmultinom
do what you want it to do, and then use Rcpp and RcppArmadillo to do the same thing on the data in your vector.