我发现一个让R快几十倍的方法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

用rcpp或者rcpp armadillo来compile c++函数。这样写成的函数比R函数呼叫速度快出几倍乃至几千倍。有兴趣的可以研究一下

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏16 回帖

关键词：COMPILE Rcpp ARMA comp Mad

// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
using namespace Rcpp;
using namespace arma;
// [[Rcpp::export]]
vec samplex(){
  vec x(100000);
  for(int i = 0; i<100000; ++i){
  double u1=R::runif(0,1);
  double a=1,b=1.7,mid=0;
  double temp=0;
  double m=exp(1.7*log(1.7))-1;
  for (int l = 0;  l<50; ++l){
  mid=(a+b)/2;
  temp=(exp(mid*log(mid))-1)/m-u1;
  if(temp<0){
  a=mid;
  }
  if(temp>0){
  b=mid;
  }
  }
  x[i]=mid;
  }
  return x;
}

这是用rcpp写的c++函数，关于从一个分布sample的。用i7处理器sample10万次只要0.61秒。用R实现这个不花个1分钟是不太可能的

而且这个方法上手极其简单。你只要装：rstudio，r和rtool即可。package只要rcpp和rcpparmadillo

哎。为了下一个pdf只好不停地来发帖了。估计我说的这些东西很多人都会。。。。

鼓励楼主的创新，

应该鼓励

学习一下

虽然这个结果很早就有了，但还是鼓励楼主继续探索，后续再介绍一些在R中用C的方法

顶一下