r - 从拉普拉斯分布生成随机数

Question

我一直在尝试从双指数(拉普拉斯)分布生成随机数。我现在可以再写代码了。任何帮助，将不胜感激。下面的代码是我写的。

rlaplace = function(u,a,b){
    u = c(runif(ns))
    for(i in 1:ns){
        if(u[i] <= 0.5){
            X = a+b*log(2*u)
        } else{
            X = a-b*log(2*(1-u))
        }
    }
    X
}
z1 = rlaplace(u,a,b)

Answer 1

来自Probability distributions CRAN Task View , 有几个包已经实现了拉普拉斯分布，特别是 distr和 Runuran .

因此，例如，您应该能够安装 distr，然后执行以下操作:

library(distr)
D <- DExp(rate = 1) 
r(D)(1)

代码取自 DExp-class 帮助页面的示例。