Jak wykreślić rozkład weibulla w r

Aby wykreślić funkcję gęstości prawdopodobieństwa dla rozkładu Weibulla w R, możemy skorzystać z następujących funkcji:

• dweibull(x, kształt, skala = 1) , aby utworzyć funkcję gęstości prawdopodobieństwa.
• krzywa(funkcja, od = NULL, do = NULL), aby wykreślić funkcję gęstości prawdopodobieństwa.

Aby wykreślić funkcję gęstości prawdopodobieństwa, należy określić wartość parametru kształtu i skali w funkcji dweibulla oraz wartości od i do w funkcji curve() .

Na przykład poniższy kod ilustruje, jak wykreślić funkcję gęstości prawdopodobieństwa dla rozkładu Weibulla z parametrami kształt = 2 i skala = 1, gdzie oś x wykresu przechodzi od 0 do 4:

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

Możemy dodać tytuł, zmienić etykietę osi Y, zwiększyć szerokość wiersza, a nawet zmienić kolor wiersza, aby wykres był bardziej estetyczny:

 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 

Możemy również dodać do wykresu kilka krzywych, aby porównać rozkłady Weibulla z różnymi parametrami kształtu i skali:

 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)

Legendę do wykresu możemy dodać za pomocą funkcji legend() , która przyjmuje następującą składnię:

legenda(x, y=NULL, legenda, wypełnienie, col, bg, lty, cex)

Złoto:

• x, y: współrzędne x i y użyte do umiejscowienia legendy
• legenda: tekst, który należy umieścić w legendzie
• fill: kolor wypełnienia legendy
• col: lista kolorów linii wewnątrz legendy
• bg: kolor tła legendy
• lty: styl linii
• cex: rozmiar tekstu w legendzie

W naszym przykładzie użyjemy następującej składni do utworzenia legendy:

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