我试图从自定义分布中生成随机数,我已经发现了这个问题: 从(任意)连续概率分布模拟, 但不幸的是它没有帮助我,因为方法建议需要分配函数的公式.我的分布是多个均匀分布的组合,基本上分布函数看起来像直方图.一个例子是:
f(x) = {
0 for x < 1
0.5 for 1 <= x < 2
0.25 for 2 <= x < 4
0 for 4 <= x
}
Run Code Online (Sandbox Code Playgroud) 怎样才能从未知密度的分位数随机抽样数据f(x)对x之间0和4R中?
f = function(x) ((x-1)^2) * exp(-(x^3/3-2*x^2/2+x))
Run Code Online (Sandbox Code Playgroud)