// [[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只好不停地来发帖了。估计我说的这些东西很多人都会。。。。

鼓励楼主的创新，