在 Rcpp 中调用 R 函数 [重复]

Question

我试图在 Rcpp 中调用 sd(x)，它是一个 R 函数。我看过一个在 Rcpp 中调用 R 函数 dbate(x) 的例子，它运行良好。

// dens calls the pdf of beta distribution in R
//[[Rcpp::export]]
double dens(double x, double a, double b)
{
    return R::dbeta(x,a,b,false);
}

但是当我厌倦了将这种方法应用于 sd(x) 如下时，它出错了。

// std calls the sd function in R
//[[Rcpp::export]]
double std(NumericVector x)
{
  return R::sd(x);
}

有谁知道为什么这不起作用？

Answer 1

您的代码存在一些问题。

1. std 与 C++ Standard Library namespace 相关
   • 这是触发：

     错误：将“std”重新定义为不同类型的符号

2. R:: 是处理Rmath functions 的命名空间。其他 R 函数将不会在此范围内找到。
3. 要在 C++ 中直接调用 R 函数，您必须使用 Rcpp::Environment 和 Rcpp::Function，如示例 sd_r_cpp_call() 中所示。
   • 这种方法存在许多问题，包括但不限于速度损失。
4. 最好使用Rcpp sugar 表达式或实现您自己的方法。

说到这里，我们来谈谈代码：

#include <Rcpp.h>

//' @title Accessing R's sd function from Rcpp
// [[Rcpp::export]]
double sd_r_cpp_call(const Rcpp::NumericVector& x){

  // Obtain environment containing function
  Rcpp::Environment base("package:stats"); 

  // Make function callable from C++
  Rcpp::Function sd_r = base["sd"];    

  // Call the function and receive its list output
  Rcpp::NumericVector res = sd_r(Rcpp::_["x"] = x,
                                 Rcpp::_["na.rm"]  = true); // example of additional param

  // Return test object in list structure
  return res[0];
}


// std calls the sd function in R
//[[Rcpp::export]]
double sd_sugar(const Rcpp::NumericVector& x){
  return Rcpp::sd(x); // uses Rcpp sugar
}

/***R
x = 1:5
r = sd(x)
v1 = sd_r_cpp_call(x)
v2 = sd_sugar(x)

all.equal(r,v1)
all.equal(r,v2)
*/