結合R中的兩個威布爾分佈

Question

我正在研究一個涉及組合兩個威布爾分佈並因此創建雙峯值曲線的項目。然後，我的目標是使用這個預測。我在網上搜索，我似乎無法找到任何東西，或者如果R有一個函數，允許我結合兩個Weibulls。 下面顯示了我用來創建兩個Weibull分佈的代碼，我希望這兩個分佈組合起來構成一個概率密度函數。

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") 

說你要在所關注的

t1 = 30 

搜索你有X之間的時間的概率，並找到最接近T1，然後返回對應的y

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