r - R中混合分布的模拟

Question

我想模拟数据集，该数据集从 3 个 weibull 分布中具有递减然后增加的风险，但我想让这个风险函数更接近于零，我怎样才能得到大约 0.1 或更少。我怎样才能修复我的代码来拥有这个？

Answer 1

不必使用pexpand rexp，因为当形状参数 = 1 时，指数是 Weibull 分布的退化情况。小于 1 的形状在原点附近将更加“独家”。你可以用两个 Weibull 得到一个漂亮的浴缸，但是因为你有一个三个，所以我只是做了一些小的改动：

 hazmix = function(w1, w2, x) {w1*dweibull(x,.8, 1)/(1-pweibull(x, .8, 1))+ 
    + w2*dweibull(x,1,.5) + 
    ((1-(w1+w2))* dweibull(x, 6, 10)/(1-pweibull(x, 6, 10)))}

 png();plot(hazmix(.2,.4, seq(0, 10, by=.1)), ylim=c(0,1) ); dev.off()

我不保证它是标准化的，如果你需要一个标准化的分布，那么只使用两个带有一个混合参数的 Weibull 会容易得多。