从R中的混合物分布进行仿真

Question

我想模拟从3个Weibull分布中有递减风险然后递增的数据集，但是我想使该风险函数更接近零，如何才能得到0.1或更少。 我该如何修复我的代码呢？

Answer 1

不必使用pexp和rexp ，因为当形状参数= 1时，指数是威布尔分布的简并情况。 小于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()

我不保证将其标准化，如果您需要标准化的分布，仅使用两个带有一个混合参数的Weibulls会容易得多。