r - R の混合分布からのシミュレーション

Question

3つのワイブル分布からハザードが減少してから増加するデータセットをシミュレートしたいのですが、このハザード関数をよりゼロに近づけたいのですが、どうすれば0.1以下に近づけることができますか。これを持つようにコードを修正するにはどうすればよいですか?

Answer 1

pexpおよびを使用する必要はありませんrexp。なぜなら、形状パラメーターが 1 の場合、指数はワイブル分布の縮退ケースだからです。1 未満の形状は、原点近くでさらに「スクーピー」になります。2 つの Weibull だけで素敵なバスタブを作ることができますが、1 つが 3 つの 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()

正規化されているとは約束しません。正規化された分布が必要な場合は、1 つの混合パラメーターで 2 つの Weibull を使用する方がはるかに簡単です。