Einfache Matrix zu Matrix-Funktion mit rcpp

Question

Ich versuche, eine einfache Funktion in r mit rcpp zu programmieren. Ich bin herauszufinden, die Syntax durch http://adv-r.had.co.nz/Rcpp.html und die Zucker Funktionen verwendet, konnte aber die folgenden Arbeit nicht bekommen:

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; 
}') 

Grundsätzlich eine Funktion in jedem Element einer Matrix anzuwenden. Wie kann ich es zum Arbeiten bringen?

Answer 1

Leider wird die Matrixmultiplikation für Rcpp NumericMatrix noch nicht unterstützt. Eine Option ist die Verwendung von Armadillo über 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") 

Hinweis % ist das Element weise Multiplikation