Rcpp错误：不兼容的类型（从％s到％s）

Question

我最近开始使用Rcpp包将我的R代码段写入C++。给定一个数据矩阵，我有以下Rcpp函数，它可以计算一些观察值的协方差的内核重加估计值。

cppFunction(' 
     NumericVector get_cov_1obs(NumericMatrix cdata, int ID, float radius){ 

     int nrow = cdata.nrow(), ncol = cdata.ncol(); 
     float norm_ = 0; 
     float w; 
     NumericMatrix out(ncol, ncol); 

     NumericMatrix outer_prod(ncol, ncol); 

     for (int i=0; i<ncol;i++){ 
     for (int j=0;j<ncol;j++){ 
     out(i,j) = 0; 
     outer_prod(i,j) = 0; 
     } 
     } 

     for (int i=0; i<nrow;i++){ 
     w = exp(-(i-ID)*(i-ID)/(2*radius)); 
     norm_ += w; 
     for (int j=0; j<ncol;j++){ 
     for (int k=0;k<ncol;k++){ 
     outer_prod(j,k) = cdata(i,j) * cdata(i,k); 
     } 
     } 

     for (int j=0; j<ncol;j++){ 
     for (int k=0;k<ncol;k++){ 
     out(j,k) += outer_prod(j,k)*w; 
     } 
     } 
     } 

     for (int i=0; i<ncol;i++){ 
     for (int j=0;j<ncol;j++){ 
     out(i,j) /= norm_; 
     } 
     } 

     return out; 
     }') 

我想快点估计内核rewieghted协方差matricies对数据集中的所有意见，并将它们存储为一个数组。由于RCPP不处理数组我写了下面的一个R函数：

get_kern_cov_C = function(data, radius){ 
    # data is data for which we wish to estimate covariances 
    # radius is the radius of the gaussian kernel 

    # calculate covariances: 
    kern_cov = array(0, c(ncol(data),ncol(data),nrow(data))) 
    for (i in 1:nrow(data)){ 
    kern_cov[,,i] = get_cov_1obs(cdata=data, ID = i-1, radius=radius) 
    } 
    return(kern_cov) 
} 

这似乎很好地工作（还有很多，很多快于R）但问题是，每一个现在，然后（似乎是随机的）我得到以下形式的错误：

Error in kern_cov[, , i] = get_cov_1obs(cdata = data, ID = i - 1, radius = radius) : 
    incompatible types (from X to Y) 

其中X或者是内置或NULL和ÿ加倍。

我大致了解为什么会发生这种情况（我试图将一个内置/ NULL变量放入一个double），但我不确定是在代码中的错误。我怀疑这可能与内存管理有关，因为它只会一次又一次地发生。

Answer 1

您也可以在C（++）级别测试NULL，在这种情况下可能应该这样做。

至于它为什么会发生：恐怕你需要调试这个。