> The main problem for me is the "sample" function of R, which I could not > implement in Rcpp. A "sample" function was discussed here recently, e.g. http://lists.r-forge.r-project.org/pipermail/rcpp-devel/2012-December/005025.html (But it has a dependency on RcppArmadillo, which might be undesirable.) Darren