小编Bro*_*oom的帖子

我是 Rcpp 的新手。

我创建了一个 rcpp 函数，它接受一个包含 2 列的数据帧和一个向量作为输入，并返回一个向量。我的数据如下

set.seed(10)
min= sort(rnorm(1000,800,sd=0.1))
max= min+0.02
k=data.frame(min,max)
explist= sort(rnorm(100,800,sd=0.2))

然后我调用cfilter.cpp

k$output <- cfilter(k,explist)

#include <Rcpp.h> 
using namespace Rcpp;
// [[Rcpp::export]]
CharacterVector cfilter(DataFrame k, NumericVector explist) {

  NumericVector col1 = k["min"];
  NumericVector col2 = k["max"];
  NumericVector exp = explist ;

  int n = col1.size();
  int j = 0;
  CharacterVector out(n);

  for (int i=0; i<n ; i++){
    out[i]=NA_STRING;

    while(exp[j]<= col2[i]){

        if( exp[j]>= col1[i] && exp[j]<= col2[i] ){
          out[i]="Y";
          break;
        }
        else if(exp[j]>col2[i]){

          break;
        }
        else …