3

我正在自学 Rcpp 并注意到 Rcpp 糖没有示例功能。所以我决定从 C++ 调用基础库中的示例函数。我有两个问题：

1.关于参数prob的类型，我应该使用NumericVector吗？使用矢量类型是否合法？

2. 关于输出的类型，我应该使用 IntegerVector 吗？使用 NumericVector 类型是否合法？

似乎所有这些类型都很好（参见下面的代码），但我想知道哪种类型更好用。

<!-- language-all: lang-html -->
library(inline)
library(Rcpp)

src1 <- '
   RNGScope scope;

  NumericVector thenum(1),myprob(3);

  myprob[0]=0.1;
  myprob[1]=0.5;
  myprob[2]=0.4;

  Environment base("package:base");
  Function sample = base["sample"];

  thenum = sample(3,Named("size",1),Named("prob",myprob));

  return wrap(thenum);
'


src2 <- '
  RNGScope scope;

  IntegerVector theint(1);
  vector<double> myprob(3);
    myprob[0]=0.1;
  myprob[1]=0.5;
  myprob[2]=0.4;
  Environment base("package:base");
  Function sample = base["sample"];

  theint = sample(3,Named("size",1),Named("prob",myprob));

  return wrap(theint);
'


fun1 <- cxxfunction(signature(),body=src1,plugin="Rcpp")
fun2 <- cxxfunction(signature(),body=src2,include='using namespace std;',plugin="Rcpp")

fun1() ## work!
fun2() ## oh this works too!
4

1 回答 1

3

因为您是sample()从 R 调用的，所以整数和数字都像在 R 本身中一样工作：

R> set.seed(42); sample(seq(1L, 5L), 5, replace=TRUE)
[1] 5 5 2 5 4
R> set.seed(42); sample(seq(1.0, 5.0), 5, replace=TRUE)
[1] 5 5 2 5 4
R>
于 2012-09-01T15:17:00.657 回答