Чтобы построить функцию плотности вероятности для распределения Вейбулла в R, мы можем использовать следующие функции:
- dweibull(x, shape, scale = 1) для создания функции плотности вероятности.
- кривая (функция, от = NULL, до = NULL) для построения графика функции плотности вероятности.
Чтобы построить график функции плотности вероятности, нам нужно указать значение параметра формы и масштаба в функции dweibull вместе со значениями from и to в функции curve() .
Например, в следующем коде показано, как построить график функции плотности вероятности для распределения Вейбулла с параметрами shape = 2 и scale = 1, где ось x графика находится в диапазоне от 0 до 4:
curve(dweibull(x, shape=2, scale = 1), from=0, to=4)
Мы можем добавить заголовок, изменить метку оси Y, увеличить ширину линии и даже изменить цвет линии, чтобы сделать график более эстетичным:
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
Мы также можем добавить на график более одной кривой, чтобы сравнить распределения Вейбулла с различными параметрами формы и масштаба:
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() , которая имеет следующий синтаксис:
легенда (x, y = NULL, легенда, заливка, col, bg, lty, cex)
куда:
- x, y: координаты x и y, используемые для позиционирования легенды.
- легенда: текст для легенды
- fill: цвет заливки внутри легенды
- col: список цветов, которые будут использоваться для линий внутри легенды.
- bg: цвет фона для легенды
- lty: стиль линии
- cex: размер текста в легенде
В нашем примере мы будем использовать следующий синтаксис для создания легенды:
#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)
