Comment tracer une distribution de Weibull dans R

Pour tracer la fonction de densité de probabilité pour une distribution de Weibull dans R, nous pouvons utiliser les fonctions suivantes :

• dweibull(x, shape, scale = 1) pour créer la fonction de densité de probabilité.
• courbe(fonction, de = NULL, à = NULL) pour tracer la fonction de densité de probabilité.

Pour tracer la fonction de densité de probabilité, nous devons spécifier la valeur du paramètre de forme et d’échelle dans la fonction dweibull ainsi que les valeurs de et vers dans la fonction courbe() .

Par exemple, le code suivant illustre comment tracer une fonction de densité de probabilité pour une distribution de Weibull avec les paramètres shape = 2 et scale = 1 où l’axe des x du tracé va de 0 à 4 :

curve(dweibull(x, shape=2, scale = 1), from=0, to=4)

Nous pouvons ajouter un titre, modifier l’étiquette de l’axe Y, augmenter la largeur de la ligne et même changer la couleur de la ligne pour rendre l’intrigue plus esthétique :

curve(dweibull(x, shape=2, scale = 1), from=0, to=4, 
    main = 'Weibull Distribution (shape = 2, scale = 1)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

Nous pouvons également ajouter plusieurs courbes au graphique pour comparer les distributions de Weibull avec différents paramètres de forme et d’échelle :

curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

Nous pouvons ajouter une légende au tracé en utilisant la fonction legend() , qui adopte la syntaxe suivante :

légende(x, y=NULL, légende, remplissage, col, bg, lty, cex)

où:

• x, y : les coordonnées x et y utilisées pour positionner la légende
• légende : le texte à mettre dans la légende
• fill : couleur de remplissage à l’intérieur de la légende
• col : la liste des couleurs à utiliser pour les lignes à l’intérieur de la légende
• bg : la couleur de fond de la légende
• lty : style de ligne
• cex : taille du texte dans la légende

Dans notre exemple nous utiliserons la syntaxe suivante pour créer une légende :

#create density plots
curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

#add legend
legend(2, .7, legend=c("shape=2, scale=1", "shape=1.5, scale=1"),
       col=c("red", "blue"), lty=1, cex=1.2)