3

对于它的功能，我还是个新手Rcpp，更不用说 C++ 本身了，所以对于你们当中的专家来说，这可能看起来微不足道。但是，没有愚蠢的问题，所以无论如何：

我想知道是否有一种方法可以使用索引一次在 C++ 中处理 NumericVector 的多个元素。为了让整个事情更清楚，这里是我正在尝试做的 R 等价物：

# Initial vector
x <- 1:10

# Extract the 2nd, 5th and 8th element of the vector
x[c(2, 5, 8)]
[1] 2 5 8

这是迄今为止我在 R 中使用sourceCpp. 它有效，但对我来说似乎很不方便。有没有更简单的方法来实现我的目标？

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector subsetNumVec(NumericVector x, IntegerVector index) {  
  // Length of the index vector
  int n = index.size();
  // Initialize output vector
  NumericVector out(n);

  // Subtract 1 from index as C++ starts to count at 0
  index = index - 1; 
  // Loop through index vector and extract values of x at the given positions
  for (int i = 0; i < n; i++) {
    out[i] = x[index[i]];
  }

  // Return output
  return out;
}

/*** R
  subsetNumVec(1:10, c(2, 5, 8))
*/
>   subsetNumVec(1:10, c(2, 5, 8))
[1] 2 5 8
4

2 回答 2

1

如果您使用 Armadillo 向量而不是 Rcpp 向量，则可以执行此操作。

Rcpp Gallery有一个包含完整示例的帖子：特别是第二个示例。您的索引条目必须位于（无符号）uvec或umat.

于 2013-07-19T11:39:24.670 回答
0

我认为没有更短的方法！

但是你NumericVector subsetNumVec(NumericVector x, IntegerVector index)很容易出错：

在这条线内

out[i] = x[index[i]];

您无需范围检查即可访问向量。因此，在普通情况下x为空或索引超出范围，您会得到一些未定义的行为。

此外，您的方法可以通过引用调用

NumericVector subsetNumVec(const NumericVector& x, const IntegerVector& index)

没有理由复制两个向量。您只需将减法index = index -1;移至out[i] = x.at(index[i] - 1);

在这里，x.at(index[i] - 1)抛出错误索引。但是你需要一些错误处理（返回空向量或在外部进行处理）。

于 2013-07-19T11:21:16.987 回答