<div dir="ltr"><div><div>Amina and I continued the conversation off list, but for closure purposes, the "multiple" version was tweaked a bit as there was comparison between signed and unsigned ints which did not throw an error for me, but did for Amina. Final multiple version posted below.<br><br></div>Thanks,<br><br></div>Avi<br><div><div><br><br><span style="font-family:monospace,monospace">#include <RcppArmadillo.h><br>using namespace Rcpp;<br>using namespace RcppArmadillo;<br>using namespace arma;<br>//[[Rcpp::depends(RcppArmadillo)]]<br><br>//[[Rcpp::export]]<br>vec Sample_Exp_Multi(size_t n, vec alpha, vec beta) {<br>  size_t m = alpha.size();<br>  vec answers(n);<br>  uvec2 S;<br>  for (size_t i = 0; i < n; ++i){<br>    do {<br>      S = randi<uvec>(2, distr_param(0, m - 1));<br>    } while (S(0) == S(1));<br>    answers(i) = exp(alpha(S(0)) + beta(S(1)));<br>  }<br>  return(answers);<br>}</span><br></div></div></div>