如何通过C ++调用R函数并传递参数

时间:2018-06-26 11:00:41

标签: c++ c r

我正在尝试从C ++程序调用R函数。

rtest = function(input ,output) {
  a <- input
  b <- output 
  outpath <- a+b
  print(a+b)
  return(outpath)
}

这是我的R函数。我需要找到一种通过传递参数从C调用此函数的方法。 calling a R function from python code with passing arguments。在这里,我做了类似的从python调用R的方法。因此,我需要指定R脚本的路径和函数名称,还需要通过python传递参数。我正在C中寻找类似的方法。但是没有得到结果。这可能很简单。任何帮助表示赞赏。

1 个答案:

答案 0 :(得分:3)

这个问题有多种解释,这就是为什么我以前没有尝试回答。这里有几种可能的解释的解决方案:

使用Rcpp定义的C ++函数,该函数从R调用并使用用户定义的R函数http://gallery.rcpp.org/articles/r-function-from-c++/之后:

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector callFunction(NumericVector x, NumericVector y, Function f) {
  NumericVector res = f(x, y);
  return res;
}

/*** R
set.seed(42)
x <- rnorm(1e5)
y <- rnorm(1e5)

rtest <- function(x, y) {
  x + y
}

head(callFunction(x, y, rtest))
head(x + y)
*/

R函数rtest在R中定义,并与它的两个参数一起传递给C ++函数callFunctionRcpp::sourceCpp()的部分结果:

> head(callFunction(x, y, rtest))
[1]  0.95642325 -0.57197358 -1.45084989 -0.18220091  0.07592864  0.56367202

> head(rtest(x, y))
[1]  0.95642325 -0.57197358 -1.45084989 -0.18220091  0.07592864  0.56367202

在R中以及通过C ++调用该函数会得到相同的结果。

使用RInside的C ++程序对存在于C ++中的数据调用用户定义的R函数。在这里,我们有两种可能性:将数据传输到R并在其中调用该函数或将函数移至C ++然后像上面一样在C ++中调用R函数:

#include <RInside.h>

int main(int argc, char *argv[]) {
    // define two vectors in C++
    std::vector<double> x({1.23, 2.34, 3.45});
    std::vector<double> y({2.34, 3.45, 1.23});
    // start R
    RInside R(argc, argv);
    // define a function in R
    R.parseEvalQ("rtest <- function(x, y) {x + y}");
    // transfer the vectors to R
    R["x"] = x;
    R["y"] = y;
    // call the function in R and return the result
    std::vector<double> z = R.parseEval("rtest(x, y)");
    std::cout << z[0] << std::endl;

    // move R function to C++
    Rcpp::Function rtest((SEXP) R.parseEval("rtest"));
    // call the R function from C++
    z = Rcpp::as<std::vector<double> >(rtest(x, y));
    std::cout << z[0] << std::endl;
    exit(0);
}

为了对此进行编译,我使用RInside中的示例随附的GNUmakefile。结果:

$ make -k run
ccache g++ -I/usr/share/R/include -I/usr/local/lib/R/site-library/Rcpp/include -I/usr/local/lib/R/site-library/RInside/include -g -O2 -fdebug-prefix-map=/home/jranke/git/r-backports/stretch/r-base-3.5.0=. -fstack-protector-strong -Wformat -Werror=format-security -Wdate-time -D_FORTIFY_SOURCE=2 -g -Wno-ignored-attributes -Wall    call_function.cpp  -Wl,--export-dynamic -fopenmp -Wl,-z,relro -L/usr/lib/R/lib -lR -lpcre -llzma -lbz2 -lz -lrt -ldl -lm -licuuc -licui18n  -lblas -llapack  -L/usr/local/lib/R/site-library/RInside/lib -lRInside -Wl,-rpath,/usr/local/lib/R/site-library/RInside/lib -o call_function

Running call_function:
3.57
3.57