Combinaison de deux distributions de Weibull dans R

Question

Je travaille sur un projet qui consiste à combiner deux distributions de Weibull et ainsi créer une courbe à double pic. Je vise ensuite à faire des prédictions en utilisant cela. J'ai cherché en ligne et je n'arrive pas à trouver quoi que ce soit sur ce sujet ou si R a une fonction qui me permet de combiner deux Weibulls. Ci-dessous montre le code que j'ai utilisé pour créer les deux distributions de Weibull que je souhaite combiner pour créer une seule fonction de densité de probabilité.

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

Toute aide serait incroyable.

Merci!

Answer 1

Serait-il logique de combiner les distributions de probabilité et d'utiliser ensuite l'élément «y» de votre liste finale pour faire des prédictions? Si c'est le cas, cela devrait fonctionner. L'AUC finale est toujours ~ 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") 

que vous voulez la probabilité à un moment d'intérêt

t1 = 30 

Recherche parmi les x que vous avez et trouver le plus proche de t1 puis retourner l'y correspondant

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