r - 使用 R 生成密度对象的随机偏差

Question

我创建了一个密度对象 dd，如下所示:

x1 <- rnorm(1000) 
x2 <- rnorm(1000, 3, 2) 
x <- rbind(x1, x2)
dd <- density(x) 
plot(dd)

这会产生这种非常非高斯分布:

alt text http://www.cerebralmastication.com/wp-content/uploads/2009/09/nongaus.png

我最终希望获得与该分布的随机偏差，类似于 rnorm 与正态分布的偏差。

我试图解决这个问题的方法是获取内核的 CDF，然后如果我向它传递累积概率(逆 CDF)，则让它告诉我变量。这样我就可以将均匀随机变量的向量转换为密度的绘图。

看来我想做的事情应该是其他人在我之前做过的基本事情。有没有简单的方法或简单的功能来做到这一点？我讨厌重新发明轮子。

FWIW我发现this R Help article但我无法理解他们在做什么，而且最终的输出似乎并没有产生我想要的结果。但这可能是我不明白的一步。

我考虑过使用Johnson distribution from the suppdists package但约翰逊不会给我提供我的数据所具有的漂亮的双峰驼峰。

Answer 1

替代方法:

sample(x, n, replace = TRUE)