r 在使用 rcpp 时中止

Question

我正在写一些代码，当运行时，它中止了。 r 版本是 3.5.1。我认为我的 rcpp 代码有问题，但我找不到它。它只是显示 R 会话中止。

################我认为这部分没有任何问题。

#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())];
  }

# ################这部分可能有问题。

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用于求一个向量的所有组合，输出一个矩阵，Subset和vector[-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;
}

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

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

Error in r_test(a, comp, fct, tempf, combnPrim, Fv, vec, Subset) :
Evaluation error: Index out of bounds: [index=10; extent=8]..

这意味着您正在尝试访问不存在的矢量元素。我还可以告诉您发生这种情况的基本原因：在 logical_vector() 中，当 idx 的元素大于 n 时，这将发生。因此，在 r_test() 中的某个时刻，您使用 idx 调用 Subset()，其中的元素大于 x.size().

_{1 你是对的，当用另一个向量子集时，尝试使用 () 进行向量子集化会出错。但是，这在访问向量的 one 元素时有效。}

r 在使用 rcpp 时中止

r aborted when using rcpp

r

rcpp