C++ 和 R 接口(interface)，获取输出

Question

我刚开始学习 R，因此遇到了看似微不足道的问题。我试图弄清楚 R 中的值是如何输出的。在 C++ 中，我们只使用返回变量，但 R 似乎不是这种情况。假设我有一个函数接受 4 个输入参数，将这些参数传递给c++ 函数执行所需的计算，现在如果我想在 R 中加载这个 myfun 并从 c++ 函数中获取输出，我需要做什么？以下是我尝试使用的模板。

extern "C" {
   SEXP myfun(SEXP S, SEXP A, SEXP B, SEXP C) { //will call this function from R.
     SEXP rate, dir, list, list_names; //declare variables

     PROTECT( rate = allocMatrix(REALSXP, 10, 2) ); //allocate 10x2 matrix of double type?
     PROTECT( dir = allocVector(INTSXP, 10) ); //allocated vector(10) of int type?

     double* p_rate = REAL(rate); //why do I need pointers? 
     int* p_dir = INTEGER(dir);

//这里我调用了一个计算vector<vector<double> > someVal的C++函数和 vector<int> someVal2 .

现在我想将这些值传递给 rate 和 dir。

 for(int i =0; i < 10; i++){
      rate[i][0] = someVal1[i][0];
      rate[i][1] = someVal1[i][1];
      dir[i] = someVal2[i]; 
  }

但是编译器不喜欢它。我不确定在这种情况下如何使用 p_rate 和 p_dir 指针。

另外，我不太明白为什么下面的代码是这样的。它似乎将字符串、速率和目录传递给 list_names，将值传递给速率和目录。为什么我们不能简单地进行 cout 并打印所需的值？

   char *tag[2] = {"rate","dir"};
   PROTECT(list_names = allocVector(STRSXP ,2));
   SET_STRING_ELT(list_names, 0, mkChar(tag[0]));
   SET_STRING_ELT(list_names, 1, mkChar(tag[1]));
   PROTECT(list = allocVector(VECSXP ,2));
   SET_VECTOR_ELT(list, 0, rate);
   SET_VECTOR_ELT(list, 1, dir);
   setAttrib(list, R_NamesSymbol , list_names);

   UNPROTECT(8);
   return (list);
}

我不明白在 R 中返回变量背后的基本思想，这看起来很困惑。如果有人能给我指出解释这些内容的适当资源，我将不胜感激。

提前感谢您的帮助!

Answer 1

如果您已经熟悉 C++，您可能会发现我们的 Rcpp无缝 R 和 C++ 集成的项目很有用。它包含大量文档和示例，此外 CRAN 现在有超过 80 个包使用它，这提供了另一个大型示例语料库。

您可以从与我们的 JSS paper 相对应的“简介”小插图开始，并在此处仔细阅读 [rcpp] 标记。

否则，您将回到更简单且仅基于 C 的 R API。