rcppを使った行列から行列への単純な関数

Question

rcppを使ってrで簡単な関数をコーディングしようとしています。私はhttp://adv-r.had.co.nz/Rcpp.htmlて構文を考え出すと砂糖の機能を使用しますが、作業に以下を得ることができなかったよ。

cppFunction('NumericMatrix fun(NumericMatrix x) { 
    NumericMatrix dd = (exp(-(pow(x, 2))/2)*(3 - 6*pow(x, 2) + pow(x,4)))/sqrt(2*PI); 
    return dd; 
}') 

基本的に行列の各要素に関数を適用します。どうすればそれを動作させることができますか？

Answer 1

残念ながら、行列乗算はRcpp NumericMatrixではまだサポートされていません。 1つの選択肢はArmadilloからRcppArmadilloまでです。

# install.packages("RcppArmadillo") 
Rcpp::cppFunction('arma::mat fun(const arma::mat& x) { 
    arma::mat dd = (exp(-(pow(x, 2))/2)%(3 - 6*pow(x, 2) + pow(x,4)))/sqrt(2*PI); 
    return dd; 
}', depends = "RcppArmadillo") 

注%は要素毎の乗算