比如你要产生2n个服从Laplace分布f(u,b)的随机数,其中u是位置参数,b是尺度参数。
先产生n个随机数x1,x2,...,xn，服从exp(1/b)的随机数；
再产生n个随机数y1,y2,...,yn，服从exp(1/b)的随机数；
第二组数取相反数-y1,-y2,...,-yn，则x1+u,x2+u,...,xn+u,-y1+u,-y2+u,...,-yn+u服从Laplace分布f(u,b)。
我也是第一次产生这样的随机数，我觉得是对的。你试试。

rskewlap(n, mu = 0, alpha = 1, beta = 1,  param = c(mu, alpha, beta))在R包GeneralizedHyperbolic中

方法其实不唯一 楼上方法目测都可行 给个cran task view的链接 各种分布应有尽有
http://cran.r-project.org/web/views/Distributions.html

误于x[xok] : 类别为'closure'的对象不可以取子集
这个错误是什么意思呢？一下子找不出来错误~
大神求教

误于x[xok] : 类别为'closure'的对象不可以取子集
这个错误是什么意思呢？一下子找不出来错误~
大神求教