如何使用安装包产生拉普拉斯分布

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有哪位大神告知一下怎么产生拉普拉斯随机数呢？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：拉普拉斯 如何使用 普拉斯 安装包 随机数 拉普拉斯 如何 就业统计 论坛 统计

比如你要产生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

求证1加1 发表于 2014-10-15 22:44
方法其实不唯一 楼上方法目测都可行 给个cran task view的链接 各种分布应有尽有
http://cran.r-project.o ...

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

求证1加1 发表于 2014-10-15 22:44
方法其实不唯一 楼上方法目测都可行 给个cran task view的链接 各种分布应有尽有
http://cran.r-project.o ...

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

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

额 那就把代码贴出来大家讨论讨论咯