Cara memplot distribusi weibull di r

Untuk memplot fungsi kepadatan probabilitas untuk distribusi Weibull di R, kita dapat menggunakan fungsi berikut:

• dweibull(x, bentuk, skala = 1) untuk membuat fungsi kepadatan probabilitas.
• curve(function, from = NULL, to = NULL) untuk memplot fungsi kepadatan probabilitas.

Untuk memplot fungsi kepadatan probabilitas, kita perlu menentukan nilai parameter bentuk dan skala dalam fungsi dweibull serta nilai dari dan ke dalam fungsi curve() .

Misalnya, kode berikut mengilustrasikan cara memplot fungsi kepadatan probabilitas untuk distribusi Weibull dengan parameter bentuk = 2 dan skala = 1 dengan sumbu x plot berubah dari 0 hingga 4:

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

Kita dapat menambahkan judul, mengubah label sumbu Y, menambah lebar baris, dan bahkan mengubah warna baris untuk membuat plot lebih estetis:

 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 

Kita juga dapat menambahkan beberapa kurva pada grafik untuk membandingkan distribusi Weibull dengan parameter bentuk dan skala yang berbeda:

 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)

Kita dapat menambahkan legenda ke plot menggunakan fungsi legenda() , yang menggunakan sintaks berikut:

legenda(x, y=NULL, legenda, isi, kolom, bg, lty, cex)

Emas:

• x, y: koordinat x dan y yang digunakan untuk memposisikan legenda
• legenda: teks yang akan dimasukkan ke dalam legenda
• fill: mengisi warna di dalam legenda
• col: daftar warna yang digunakan untuk garis di dalam legenda
• bg : warna latar legenda
• lty: gaya garis
• cex: ukuran teks dalam legenda

Dalam contoh kita, kita akan menggunakan sintaks berikut untuk membuat 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)