r - 在 R 中模拟一个奇怪的分布

Question

大家好，我正在尝试模拟看起来像这样的数据:

有人知道我会怎么做吗？

Answer 1

这可以被认为是一个经典问题，将 runif 生成的均匀随机数转换为经验相对频率的逆变换。这是一种使用大约的方法:

freq <- c(0.39, 0.02, 0.15, 0.18, 0.12, 0.09,0.04, 0.01)
sum(freq) # This is a check. The sum must be 1.0.

r_empirical <- function(n, freq) {
  approx(c(0, cumsum(freq)), 0:(length(freq)),
         runif(n), method="constant", f=0)$y
}

x <- r_empirical(1000, freq)

hist(x, breaks=0:length(freq))

下图演示了基本原理。阶梯显示了累积分布，红色箭头表示均匀随机数如何变换: