3

いくつかの R コードを Rcpp コードに変換しており、平均ベクトルと標準偏差ベクトルが与えられた観測ベクトルの尤度を計算する必要があります。平均が 0 で標準偏差が 1 であると仮定すると、この関数を記述できます (これを実行するには、'inline' および 'Rcpp' パッケージをロードする必要があります)。

dtest1 = cxxfunction(signature( x = "numeric"),
                      'Rcpp::NumericVector xx(x);
                       return::wrap(dnorm(xx, 0.0, 1.0));', 
                       plugin='Rcpp')

そして結果は予想通り。

> dtest1(1:3) 
[1] 0.241970725 0.053990967 0.004431848

しかし、関数を作ろうとすると

dtest2 = cxxfunction(signature( x = "numeric", y="numeric", z="numeric" ),
                  'Rcpp::NumericVector xx(x);
                   Rcpp::NumericVector yy(y);
                   Rcpp::NumericVector zz(z);
                   return::wrap(dnorm(xx, yy, zz));',
                   plugin='Rcpp')

これにより、さまざまな手段を渡すことができ、標準偏差は以下に示すエラーになります。私が作成しようとしている関数を作成する方法はありますか、または正規密度を手動でプログラムする必要がありますか?

エラー

Error in compileCode(f, code, language = language, verbose = verbose) :   
    Compilation ERROR, function(s)/method(s) not created! file31c82bff9d7c.cpp: In function ‘SEXPREC* file31c82bff9d7c(SEXP, SEXP, SEXP)’:
file31c82bff9d7c.cpp:33:53: error: no matching function for call to 
     ‘dnorm4(Rcpp::NumericVector&, Rcpp::NumericVector&, Rcpp::NumericVector&)’
file31c82bff9d7c.cpp:33:53: note: candidates are:
/home/chris/R/x86_64-pc-linux-gnu-library/3.0/Rcpp/include/Rcpp/stats/norm.h:106:1: 
     note: template<int RTYPE, bool NA, class T> Rcpp::stats::D0<RTYPE, NA, T> Rcpp::dnorm4(const Rcpp::VectorBase<RTYPE, NA, VECTOR>&, bool)
/home/chris/R/x86_64-pc-linux-gnu-library/3.0/Rcpp/include/Rcpp/stats/norm.h:107:1: 
     note: template<int RTYPE, bool NA, class T> Rcpp::stats::D1<RTYPE, NA, T> Rcpp::dnorm4(const Rcpp::VectorBase<RTYPE, NA, VECTOR>&, double, bool)
/home/chris/R/x86_64-pc-linux-gnu-library/3.0/Rcpp/include/Rcpp/stats/norm.h:108:1: 
     note: template<int RTYPE, bool NA, class T> Rcpp::stats::D2<RTYPE, NA, T> Rcpp::dnorm4(const Rcpp::VectorBase<RTYPE, NA,
In addition: Warning message:
running command '/usr/lib/R/bin/R CMD SHLIB file31c82bff9d7c.cpp 2> file31c82bff9d7c.cpp.err.txt' had status 1
4

1 に答える 1

7

砂糖dnormは、最初の引数に関してのみベクトル化されます。

簡単にするために (もう少し複雑ですが、ここではまだ気にする必要はありません)、

dnorm(xx, 0.0, 1.0)

オーバーロードを使用する

NumericVector dnorm( NumericVector, double, double )

そして、2番目の呼び出しは次のようなものを使用しようとします

NumericVector dnorm( NumericVector, NumericVector, NumericVector )

これは実装されていません。それを実装することはできますが、優先順位リストで十分に高くする必要があります。

それまでの間、次のような小さなラッパーを書くのは簡単です (これは引数の長さなどを処理しません ...):

NumericVector my_dnorm( NumericVector x, NumericVector means, NumericVector sds){
    int n = x.size() ;
    NumericVector res(n) ;
    for( int i=0; i<n; i++) res[i] = R::dnorm( x[i], means[i], sds[i] ) ;
    return res ;
}
于 2013-08-09T19:39:36.290 に答える