Как построить распределение вейбулла в r

Чтобы построить функцию плотности вероятности для распределения Вейбулла в R, мы можем использовать следующие функции:

• dweibull(x, shape, масштаб = 1) для создания функции плотности вероятности.
• кривая(функция, from = NULL, to = NULL) для построения функции плотности вероятности.

Чтобы построить функцию плотности вероятности, нам нужно указать значение параметра формы и масштаба в функции dweibull , а также значения from и to в функции Curve() .

Например, следующий код показывает, как построить график функции плотности вероятности для распределения Вейбулла с параметрами shape = 2 и масштабом = 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, легенда, заливка, столбец, фон, 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)