假设我有一个简单的数组,具有相应的概率分布。
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)
有没有办法可以使用相同的分布生成另一组数据。由于操作是概率性的,因此它不再需要与初始分布完全匹配,而只是从中生成。
我确实成功地自己找到了一个简单的解决方案。谢谢!
最佳答案
最好的选择是生成经验累积密度函数,近似其倒数,然后转换输入。
复合表达式看起来像
random.points <- approx(
cumsum(pdf_of_data$y)/sum(pdf_of_data$y),
pdf_of_data$x,
runif(10000)
)$y
产量
hist(random.points, 100)
关于R:从概率密度分布生成数据,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/32871602/