Rの2つのワイブル分布を組み合わせる

Question

私は2つのワイブル分布を組み合わせてダブルピークカーブを作成するプロジェクトに取り組んでいます。私はそれを使って予測をするつもりです。私はオンラインで検索しましたが、これに関して何かを見つけることができないか、Rに2つのワイブルを組み合わせる機能があると思われます。 以下は、単一の確率密度関数を作成するために結合したい2つのワイブル分布を作成するために使用したコードを示しています。

curve(dweibull(x, scale=30.59898985, shape=2.27136646),from=0, to=70, main="Weibull distribution") 
curve(dweibull(x, scale=19.39743639, shape=1.22800332),from=0, to=70, main="Weibull distribution") 

何か助けが驚くでしょう。

ありがとうございます！

Answer 1

確率分布を結合し、最後のリストの要素 "y"を使用して予測を行うのは意味がありますか？もしそうなら、これはうまくいくはずです。最終的なAUCは依然として〜1です。

dwb1 <- curve(dweibull(x, scale=30.59898985, shape=2.27136646),from=0, to=70, main="Weibull distribution") 
dwb2 <- curve(dweibull(x, scale=19.39743639, shape=1.22800332),from=0, to=70, main="Weibull distribution") 

# combine 
final.dwb <- lapply(c("x", "y"), (function(i){ 
    (dwb1[[i]] + dwb2[[i]])/2 
})) 
names(final.dwb) <- c("x", "y") 

# plot 
plot(final.dwb$y ~ final.dwb$x, xlim=c(0,70), main = "combined Weibull distributions", type = "n", las = 2) 
lines(final.dwb$y ~ final.dwb$x, xlim=c(0,70), main = "combined Weibull distributions") 

あなたが興味のある時点で、あなたが持っているXの中

t1 = 30 

検索をする確率をしたいと言うとT1と、その後

対応するYを返すに最も近いを見つけます

id <- which.min(abs(t1 - final.dwb$x)) 
final.dwb$y[id]