在 R 中使用 RcppArmadillo 包和 rowvec 时出现编译错误

Question

我正在尝试编译以下代码。请参阅下面我到目前为止所尝试的内容。有什么我想念的吗。任何帮助，将不胜感激。

#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]
using namespace Rcpp;

// [[Rcpp::export]]
List beta(const arma::rowvec beta,
          const int n, 
          const int L1,
          const int p,
          const arma::mat YWeight1,
          const arma::mat z){

    double S0=0;

    for(int i = 0;i < n; ++i){
        arma::rowvec zr = z.rows(i,i);
        S0 +=  exp(arma::as_scalar(beta*zr.t()));
    }

    List res;
    res["1"] = S0;
    return(res);
}

我无法复制错误，但这就是我得到的。

no match for call to '(Rcpp::traits::input_parameter<const arma::Row<double> 

等等……

Answer 1

有一个rowvec 转换器。这里的问题是：

filece5923f317b2.cpp:39:34: 错误: type 'Rcpp::traits::input_parameter::type' (aka 'ConstInputParameter >') 不提供调用操作符 rcpp_result_gen = Rcpp::wrap(beta(beta, n, L1, p, YWeight1, z));

几个想法：1. 已经有一个名为 beta() 的函数，2. 有一个名为 beta 的变量，它可能会对 Rcpp 属性造成严重破坏。

解决办法：

• 删除using namespace Rcpp;
• 将函数从 beta() 重命名为 beta_estimator()。
• 指定Rcpp::List的长度
• 通过数字索引而不是字符串访问。

更正代码：

#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]

// [[Rcpp::export]]
Rcpp::List beta_estimator( // renamed function
          const arma::rowvec beta,
          const int n, 
          const int L1,
          const int p,
          const arma::mat YWeight1,
          const arma::mat z){

    double S0 = 0;

    for(int i = 0;i < n; ++i){
        arma::rowvec zr = z.rows(i,i);
        S0 += exp(arma::as_scalar(beta*zr.t()));
    }

    // Specify list length
    Rcpp::List res(1);
    res[0] = S0;
    return(res);
}