R:从概率密度分布生成数据

Question

假设我有一个简单的数组，具有相应的概率分布。

library(stats)    
data <- c(0,0.08,0.15,0.28,0.90)
pdf_of_data <- density(data, from= 0, to=1, bw=0.1)

有没有办法可以使用相同的分布生成另一组数据。由于操作是概率性的，因此它不再需要与初始分布完全匹配，而只是从中生成。

我确实成功地自己找到了一个简单的解决方案。谢谢!

Answer 1

最好的选择是生成经验累积密度函数，近似其倒数，然后转换输入。

复合表达式看起来像

random.points <- approx(
  cumsum(pdf_of_data$y)/sum(pdf_of_data$y),
  pdf_of_data$x,
  runif(10000)
)$y

产量

hist(random.points, 100)