как получить бета-функцию в rcpp

Question

Мне нужно использовать бета-функцию Beta (a, b) в моих вычислениях в Rcpp.Простой пример кода, который я написал, здесь:

cppFunction('double getbeta(double a, double b){
  double res = beta(a, b);
  return res;
}')

Но я получил ошибку в R, говорящую, что

нет соответствующей функции для вызова 'beta'

Это потому, что мы не можем использовать бета-функцию в Rcpp?Если бы мы могли, можно ли использовать соответствующую функцию

lbeta 

?

Answer 1

Используйте пространство имен R:: для функций со скалярными значениями (и Rcpp:: для функций векторизованного сахара):

R> library(Rcpp)
R> cppFunction('double getbeta(double a, double b){
+   double res = R::beta(a, b);
+   return res;
+ }')
R> getbeta(1.0, 2.0)
[1] 0.5
R>