r - 如何在 Rcpp 上使用 qnorm？

Question

require(inline)
func <- cxxfunction(, 'return Rcpp::wrap( qnorm(0.95,0.0,1.0) );' ,plugin="Rcpp")

错误:没有匹配的函数来调用“qnorm5(double, int, int)”

require(inline)
func <- cxxfunction(, 'return Rcpp::wrap( qnorm(0.95, 0.0, 1.0, 1, 0) );' 
                   ,plugin="Rcpp")

错误:没有匹配的函数来调用“qnorm5(double, double, double, int, int)”

require(inline)
code <-'
double a = qnorm(0.95, 0.0, 1.0);
return  Rcpp::wrap( a );
' 
func <- 
cxxfunction(, code ,plugin="Rcpp")
func()

错误:没有匹配的函数来调用“qnorm5(double, double, double)”

如何在 Rcpp 上使用 qnorm？

Answer 1

通过将 mean 和 sd 参数设为 double，如错误消息所示 - 所以请尝试这是一个完整的示例

library(inline)
f <- cxxfunction(signature(xs="numeric", plugin="Rcpp", body='
     Rcpp::NumericVector x(xs);
     return Rcpp::wrap(Rcpp::qnorm(x, 1.0, 0.0));
')

并查看示例和单元测试 - 我刚刚在单元测试文件 runit.stats.R 中查找了这一点，其中有很多针对这些统计“Rcpp 糖”的测试用例' 函数。

编辑于 2012 年 11 月 14 日: 随着今天发布的 Rcpp 0.10.0，您可以调用 do 签名 R::qnorm(double, double , double, int, int) 如果您想使用针对 Rmath.h 编写的 C 风格代码。 Rcpp Sugar 仍然为您提供矢量化版本。