在犰狳矢量类的RcppArmadillo样本

时间:2017-04-11 17:58:45

标签: c++ r rcpp armadillo

我们一直在使用sample中的RcppArmadillo函数随机抽样NumericVector个对象。但是,我们注意到在Armadillo类型(vecuvec)上使用相同的功能是不可能的。我们已经查看了sample.h文件中的函数定义,它看起来像是一个模板函数,应该可以使用这些类型,但我们还没有能够弄清楚如何使它与Armadillo课程没有与NumericVector库中的IntegerVectorRcpp类型进行大量转换。

例如,我们将此函数编写在名为try.cpp的文件中。

// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
#include <RcppArmadilloExtensions/sample.h>

using namespace arma;
using namespace Rcpp;

// [[Rcpp::export]]
arma::uvec sample_index(const int &size){
    arma::uvec sequence = linspace<uvec>(0, size-1, size);
    arma::uvec out = sample(sequence, size, false);
    return out;
}

运行上面的代码会产生以下错误:

src/try.cpp|11 col 22 error| no matching function for call to 'sample' [cpp/gcc]      

~/Library/R/3.3/library/Rcpp/include/Rcpp/sugar/functions/sample.h|401 col 1 error| note: candidate function not viable: no known conversion from 'arma::uvec' (aka 'Col<unsigned int>') to 'int' for 1st argument [cpp/gcc]

~/Library/R/3.3/library/Rcpp/include/Rcpp/sugar/functions/sample.h|437 col 1 error| note: candidate template ignored: could not match 'Vector' against 'Col' [cpp/gcc]

对此的任何帮助将不胜感激:)

1 个答案:

答案 0 :(得分:1)

如果将来遇到此问题,问题似乎与正在使用的命名空间中的sample函数的多个定义有关。特别是键入定义了所需函数的名称空间可以解决问题。特别是,需要从sample调用Rcpp::RcppArmadillo函数。

以下代码可以根据需要运行。

// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
#include <RcppArmadilloExtensions/sample.h>

// [[Rcpp::export]]
arma::uvec sample_index(const int &size){
    arma::uvec sequence = arma::linspace<arma::uvec>(0, size-1, size);
    arma::uvec out = Rcpp::RcppArmadillo::sample(sequence, size, false);
    return out;
}