r-如何访问Rcpp::列表中向量的元素



我很困惑。

以下编译并运行良好:

#include <Rcpp.h>
using namespace Rcpp;
// [[Rcpp::export]]
List test(){
    List l;
    IntegerVector v(5, NA_INTEGER);
    l.push_back(v);
    return l;
}

在R:

R) test()
[[1]]
[1] NA NA NA NA NA

但当我尝试在列表中设置IntegerVector时:

// [[Rcpp::export]]
List test(){
    List l;
    IntegerVector v(5, NA_INTEGER);
    l.push_back(v);
    l[0][1] = 1;
    return l;
}

它不编译:

test.cpp:121:8: error: invalid use of incomplete type 'struct SEXPREC'
C:/PROGRA~1/R/R-30~1.0/include/Rinternals.h:393:16: error: forward declaration of 'struct SEXPREC'

这是因为这行:

l[0][1] = 1;

编译器不知道l是一个整数向量列表。本质上,l[0]为您提供了一个SEXP(所有R对象的泛型类型),而SEXP是指向SEXPREC的不透明指针,我们无法访问其中的te定义(因此是不透明的)。所以当你做[1]时,你试图得到第二个SEXPREC,所以不透明度使它不可能,而且这不是你想要的。

你必须明确你正在提取IntegerVector,所以你可以做这样的事情:

as<IntegerVector>(l[0])[1] = 1;

v[1] = 1 ;

IntegerVector x = l[0] ; x[1] = 1 ;

所有这些选项都适用于相同的底层数据结构。

或者,如果您真的想要语法l[0][1],您可以定义自己的数据结构来表示"整数向量列表"。这是一个草图:

template <class T>
class ListOf {
public:
    ListOf( List data_) : data(data_){}
    T operator[](int i){
        return as<T>( data[i] ) ;
    }
    operator List(){ return data ; }
private:
    List data ;
} ;

你可以使用,例如:

// [[Rcpp::export]]
List test2(){
    ListOf<IntegerVector> l = List::create( IntegerVector(5, NA_INTEGER) ) ; 
    l[0][1] = 1 ;
    return l;
}

还要注意,在Rcpp矢量(包括列表)上使用.push_back需要列表数据的完整副本,这可能会导致速度减慢。只有在没有选择的情况下才使用调整大小功能。

最新更新