6

私はRの学習を始めたばかりなので、些細な問題に悩まされています。Rで値がどのように出力されるかを理解しようとしています.C++では単に戻り変数を使用しますが、Rの場合はそうではないようです.4つの入力引数を取り、それらの引数を必要な計算を行う c++ 関数です。この myfun を R にロードして 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で変数を返す背後にある基本的な考え方がわかりません。非常に混乱しているようです。誰かがこれらのことを説明する適切なリソースを教えてくれたら本当にありがたいです.

事前にご協力いただきありがとうございます。

4

2 に答える 2

9

すでに C++ に精通している場合は、R と C++ をシームレスに統合するためのRcppプロジェクトが役立つことがあります。これには多くのドキュメントと例が含まれており、さらに CRAN にはそれを使用する 80 を超えるパッケージがあり、別の大規模な例のコーパスを提供しています。

私たちのJSS ペーパー[rcpp]に対応する「紹介」ビネットから始めて、ここのタグを熟読することもできます。

それ以外の場合は、R API に戻ります。R API は、より必要最小限で、C のみに基づいています。

于 2012-11-11T18:09:35.807 に答える
5

以前の質問に答えているときに言ったように、 Rcppを試してみて、コード内のすべての煩わしさを取り除いてください。

#include <Rcpp.h>
using namespace Rcpp ;

void yourCode( std::vector< std::vector<double> >& someVal, 
               std::vector<double>& someVal2 ){ ... }

extern "C" SEXP myfun() { 
    NumericMatrix rate(10, 2); 
    NumericVector dir(10) ;

    std::vector< std::vector<double> > someVal ;
    std::vector< double > someVal2 ;

    // call whatever function that fills someVal and someVal2
    yourCode( someVal, someVal2 ) ;

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

    // structure the output
    List result = List::create( 
        _["rate"] = rate,
        _["dir"] = dir
    ) ;

    return result ;

}
于 2012-11-12T08:40:53.383 に答える