Come tracciare una distribuzione di weibull in r

Per tracciare la funzione di densità di probabilità per una distribuzione di Weibull in R, possiamo utilizzare le seguenti funzioni:

• dweibull(x, forma, scala = 1) per creare la funzione di densità di probabilità.
• curve(funzione, da = NULL, a = NULL) per tracciare la funzione di densità di probabilità.

Per tracciare la funzione di densità di probabilità, dobbiamo specificare il valore del parametro di forma e scala nella funzione dweibull , nonché i valori da e a nella funzione curve() .

Ad esempio, il codice seguente illustra come tracciare una funzione di densità di probabilità per una distribuzione Weibull con parametri forma = 2 e scala = 1 dove l’asse x del grafico va da 0 a 4:

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

Possiamo aggiungere un titolo, modificare l’etichetta dell’asse Y, aumentare la larghezza della riga e persino cambiare il colore della riga per rendere la trama esteticamente più gradevole:

 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 

Possiamo anche aggiungere diverse curve al grafico per confrontare distribuzioni di Weibull con diversi parametri di forma e scala:

 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)

Possiamo aggiungere una legenda alla trama utilizzando la funzione legend() , che accetta la seguente sintassi:

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

Oro:

• x, y: le coordinate xey utilizzate per posizionare la legenda
• legenda: il testo da inserire nella legenda
• fill: colore di riempimento all’interno della legenda
• col: l’elenco dei colori da utilizzare per le righe all’interno della legenda
• bg: il colore di sfondo della legenda
• lty: stile della linea
• cex: dimensione del testo nella legenda

Nel nostro esempio utilizzeremo la seguente sintassi per creare una legenda:

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