Hoe een weibull-verdeling in r te plotten

Om de kansdichtheidsfunctie voor een Weibull-verdeling in R te plotten, kunnen we de volgende functies gebruiken:

• dweibull(x, vorm, schaal = 1) om de kansdichtheidsfunctie te creëren.
• curve(functie, van = NULL, tot = NULL) om de kansdichtheidsfunctie uit te zetten.

Om de waarschijnlijkheidsdichtheidsfunctie te plotten, moeten we de waarde van de vorm- en schaalparameter in de dweibull- functie specificeren, evenals de from- en to- waarden in de curve()- functie.

De volgende code illustreert bijvoorbeeld hoe u een waarschijnlijkheidsdichtheidsfunctie plot voor een Weibull-verdeling met parameters vorm = 2 en schaal = 1, waarbij de x-as van de grafiek van 0 naar 4 gaat:

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

We kunnen een titel toevoegen, het Y-aslabel wijzigen, de rijbreedte vergroten en zelfs de rijkleur wijzigen om de plot esthetisch aantrekkelijker te maken:

 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 

We kunnen ook verschillende curven aan de grafiek toevoegen om Weibull-verdelingen met verschillende vorm- en schaalparameters te vergelijken:

 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)

We kunnen een legenda aan de plot toevoegen met behulp van de legend()- functie, die de volgende syntaxis heeft:

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

Goud:

• x, y: de x- en y-coördinaten die worden gebruikt om de legenda te positioneren
• legenda: de tekst die in de legenda moet worden geplaatst
• vulling: vulkleur binnen de legenda
• col: de lijst met kleuren die moeten worden gebruikt voor de lijnen in de legenda
• bg: de achtergrondkleur van de legenda
• lty: lijnstijl
• cex: grootte van de tekst in de legenda

In ons voorbeeld gebruiken we de volgende syntaxis om een legenda te maken:

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