廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

Rcpp の rexp の引数はレートじゃなくて平均だった。

こうやってやると R の rexp と結果が一致する.

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
NumericVector rexp2(int N, double lambda) {
  NumericVector out(N);
  for(int i=0; i<N ;i++){
    out[i] = R::rexp(1/lambda);
  }
  return out;
}