[R使用RCPP时中止

Question

我写一些代码，并运行时，它中止。将R版本3.5.1是。我觉得有什么毛病我RCPP代码，但我不能找到它。这只是表明R对话中止。

################I don't think there is anything wrong in this part.

#include <Rcpp.h>

Rcpp::LogicalVector logical_index(Rcpp::IntegerVector idx, R_xlen_t n) {
  bool invert = false; 
  Rcpp::LogicalVector result(n, false);

  for (R_xlen_t i = 0; i < idx.size(); i++) {
    if (!invert && idx[i] < 0) invert = true;
    result[std::abs(idx[i])] = true;
  }

  if (!invert) return result;
  return !result;
}


// [[Rcpp::export]]
Rcpp::NumericVector 
  Subset(Rcpp::NumericVector x, Rcpp::IntegerVector idx) {
    return x[logical_index(idx, x.size())];
  }

# ################There maybe something wrong in this part.

library(Rcpp)
library(gRbase)
vec=0:13
a=combnPrim(vec,4)
cppFunction("
NumericVector r_test(NumericMatrix a, Function comp, Function fct,Function 
tempf, Function combnPrim, NumericVector Fv, NumericVector vec, Function 
Subset){
        NumericVector if1(4);
        NumericVector if2(2);
        NumericVector if3(4);
        NumericVector if4(4);
        NumericVector seq(14);
        NumericMatrix b(2,45);
        NumericMatrix c(4,70);
        for (int i1=0; i1<=a.ncol()-1; i1++){
        if1=a (_, i1);
        b=combnPrim(Subset(vec,-(if1)),2);
        for (int i2=0; i2<=b.ncol()-1; i2++){
        if2=b (_, i2);
        c=combnPrim(Subset(Subset(vec,-(if1)),-(if2)),4);
        for (int i3=0; i3<=c.ncol()-1; i3++){
        if3=c (_, i3);
        if4=Subset(Subset(Subset(vec,-(if1)),-(if2)),-(if3));
        seq=tempf(if1,if2,if3,if4);
        }}}return if1;}")

Fv=7.001327
setwd("D:/STAT 602")

sourceCpp("Subset.cpp")

fct=function(seq,data,fact){
  return(anova(lm(data[seq]~as.factor(fact)))$`F value`[1])
}
tempf=function(x1,x2,x3,x4){
  return(c(x1,x2,x3,x4))
}
comp=function(x,y){
  return(x>y)
}
r_test(a,comp, fct, tempf, combnPrim, Fv, vec, Subset)

#

combnPrim用于查找一个向量的所有组合，并且它输出的矩阵，并且子集是为矢量[-i]相同的用途。单独来看，所有功能可在cppFunction的表演非常出色，但把它们像这样的时候，事情出错。

Answer 1

由于这个代码是非常令人费解，我不会去寻找确切，具体的错误。不过，我可以证实，@RalfStubner是正确的，你的问题是，你试图访问不存在的向量的元素。

如果更改logical_index()到

LogicalVector logical_index(IntegerVector idx, R_xlen_t n) {
    bool invert = false;
    LogicalVector result(n, false);

    for (R_xlen_t i = 0; i < idx.size(); i++) {
        if (!invert && idx(i) < 0) invert = true;
        result(std::abs(idx(i))) = true;
    }

    if (!invert) return result;
    return !result;
}

让元素访问是通过()完成，而不是[]确保边界检查，1而不是仅仅中止，我看到以下内容：

r_test(a,comp, fct, tempf, combnPrim, Fv, vec, Subset)

误差在r_test（一，COMP，FCT，tempf，combnPrim，FV，VEC，子集）： 评价错误：索引超出范围：[指数= 10;程度= 8] ..

这意味着你试图访问不存在的向量的元素。我还可以告诉你为什么发生这种情况的根本原因：在logical_vector()，当idx的元素比n更大，这会发生。因此，在r_test()某些时候，您正在使用具有一个元素比Subset()更大idx调用x.size()。

---

1你是正确的，你会得到一个错误试图通过子集化另一个向量的向量时使用()矢量子集。然而，访问向量的一个元素，当这个作品。