Wie zeichnet man eine Weibull-Verteilung in R

Um die Wahrscheinlichkeitsdichtefunktion für eine Weibull-Verteilung in R darzustellen, können wir die folgenden Funktionen verwenden:

• dweibull(x, shape, scale = 1), um die Wahrscheinlichkeitsdichtefunktion zu erstellen.
• curve(function, from = NULL, to = NULL) zur Darstellung der Wahrscheinlichkeitsdichtefunktion.

Um die Wahrscheinlichkeitsdichtefunktion plotten, müssen wir den Wert für die Form und Skalenparameter in der dweibull-Funktion spezifizieren zusammen mit der von und bis Werte in der curve()-Funktion.

Der folgende Code zeigt beispielsweise, wie eine Wahrscheinlichkeitsdichtefunktion für eine Weibull-Verteilung mit den Parametern Form = 2 und Skala = 1 dargestellt wird, wobei die x-Achse des Diagramms zwischen 0 und 4 liegt:

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

Wir können einen Titel hinzufügen, die Beschriftung der y-Achse ändern, die Linienbreite erhöhen und sogar die Linienfarbe ändern, um das Diagramm ästhetischer zu gestalten:

curve(dweibull(x, shape=2, scale = 1), from=0, to=4, 
    main = 'Weibull Distribution (shape = 2, scale = 1)', #Title hinzufügen
    ylab = 'Density', #y-Achsentitel ändern
    lwd = 2, #Linienstärke auf 2 setzen
    col = 'steelblue') #Linienfarbe auf steelblue ändern

Wir können dem Diagramm auch mehr als eine Kurve hinzufügen, um Weibull-Verteilungen mit verschiedenen Form- und Skalierungsparametern zu vergleichen:

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)

Wir können dem Plot eine Legende hinzufügen, indem wir die Funktion legend() verwenden, die die folgende Syntax annimmt:

legend(x, y=NULL, legend, fill, col, bg, lty, cex)

wobei:

• x, y: Die x- und y-Koordinaten, mit denen die Legende positioniert wird
• legend: Der Text in der Legende
• fill: Füllfarbe innerhalb der Legende
• col: Die Liste der Farben, die für die Linien innerhalb der Legende verwendet werden sollen
• bg: die Hintergrundfarbe für die Legende
• lty: Linienstil
• cex: Textgröße in der Legende

In unserem Beispiel verwenden wir die folgende Syntax, um eine Legende zu erstellen:

#Dichtediagramme erstellen
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)
#Legende hinzufügen
legend(2, .7, legend=c("shape=2, scale=1", "shape=1.5, scale=1"),
       col=c("red", "blue"), lty=1, cex=1.2)