使用Rcpp糖函数计算密度

Question

我想像这样计算Gamma分布的密度：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector dg(double sh, double sc){
  return dgamma(seq_len(10)*0.05, sh, sc);
}

然而，当我得到

时，这不能正常工作

> dg(10,0.1)
[1] 0 0 0 0 0 0 0 0 0 0

我想知道代码中有什么问题？

编辑：我想使用矢量版本Rcpp::dgamma而不是R::命名空间中的标量版本。