Question

我正在考虑使用Rcpp包将一些C ++用于我的R包的简单部分。我是C ++新手（但很想学习！）。我使用优秀的Rcpp实现了一些简单的cpp程序 - 实际上这个程序包激励我学习C ++ ......

无论如何，我遇到了一个简单的问题，如果我能解决这个问题就会有所帮助。我有一个NumericVector我想要子集然后排序。下面的代码对整个向量进行排序（并且还将处理NAs，这是我需要的）。

我的问题是，我想要提取此向量的一部分，排序并将其用于其他处理 - 我该怎么做？例如，对于长度为10的向量，如何提取和排序元素5:10？

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
RcppExport SEXP rollP(SEXP x) {
  NumericVector A(x); // the data  
  A = sort_unique(A);  
  return A;
}

我从R打电话：

sourceCpp( "rollP.cpp")
rollP(10:1)
# [1]  1  2  3  4  5  6  7  8  9 10

Answer 1

以下是3种变体：

include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector rollP(NumericVector A, int start, int end) {
  NumericVector B(end-start+1) ;
  std::copy( A.begin() + start-1, A.begin() + end, B.begin() ) ;
  return B.sort() ;
}

// [[Rcpp::export]]
NumericVector rollP2(NumericVector A, int start, int end) {
  NumericVector B( A.begin() + start-1, A.begin() + end ) ;
  return B.sort() ;
}

// [[Rcpp::export]]
NumericVector rollP3(NumericVector A, int start, int end) {
  NumericVector B = A[seq(start-1, end-1)] ;
  return B.sort() ;
}

start和end表示基于1的索引，就像您从A[start:end]传递R一样。

Answer 2

你需要研究C ++索引，迭代器和整个位。至少，您需要更改您的界面（vector，fromInd，toInd）并找出您想要返回的内容。

对您的问题的一种解释是将子集从[fromInd, toInd)复制到新的向量中，对其进行排序并将其返回。所有这些都是标准的C ++票价，而优秀的（和免费!!）C++ Annotations这样的好文本将会有所帮助。它也有一个非常强大的STL部分。

Answer 3

您可以在std::slice上使用std::valarray。但是，如果您想特别使用std::vector，则可以使用std::copy提取向量的一部分，然后使用std::sort对提取的向量切片进行排序。

Answer 4

通过使用接收两个迭代器的std::sort实现，您可以非常轻松地完成此任务：

#include <vector>
#include <cinttypes>
#include <algorithm>

template <typename SeqContainer>
SeqContainer slicesort(SeqContainer const& sq, size_t begin, size_t end) {
  auto const b = std::begin(sq)+begin;
  auto const e = std::begin(sq)+end;
  if (b <= std::end(sq) && e <= std::end(sq)) {
    SeqContainer copy(b,e);
    std::sort(copy.begin(),copy.end());
    return copy;
  }
  return SeqContainer();
}

可以像

一样调用

  std::vector<int> v = {3,1,7,3,6,-2,-8,-7,-1,-4,2,3,9};
  std::vector<int> v2 = slicesort(v,5,10);

子集一个矢量并对其进行排序

4 个答案: