Як побудувати розподіл вейбулла в r

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

• dweibull(x, форма, масштаб = 1) , щоб створити функцію щільності ймовірності.
• curve(function, from = NULL, to = NULL) , щоб побудувати графік функції щільності ймовірності.

Щоб побудувати графік функції щільності ймовірності, нам потрібно вказати значення параметра форми та масштабу у функції Двейбулла , а також значення 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, легенда, заливка, стовпець, 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)