R'de weibull dağılımı nasıl çizilir

R’deki bir Weibull dağılımı için olasılık yoğunluk fonksiyonunu çizmek için aşağıdaki fonksiyonları kullanabiliriz:

• olasılık yoğunluk fonksiyonunu oluşturmak için dweibull(x, şekil, ölçek = 1) .
• Olasılık yoğunluk fonksiyonunu çizmek için eğri(işlev, = NULL’dan = NULL’a) kullanın .

Olasılık yoğunluk fonksiyonunun grafiğini çizmek için, curve() fonksiyonunda from ve to değerlerinin yanı sıra dweibull fonksiyonunda şekil ve ölçek parametresinin değerini belirtmemiz gerekir.

Örneğin, aşağıdaki kod, grafiğin x ekseninin 0’dan 4’e gittiği şekil = 2 ve ölçek = 1 parametreleriyle bir Weibull dağılımı için olasılık yoğunluk fonksiyonunun nasıl çizileceğini gösterir:

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

Grafiği estetik açıdan daha hoş hale getirmek için bir başlık ekleyebilir, Y ekseni etiketini değiştirebilir, satır genişliğini artırabilir ve hatta satır rengini değiştirebiliriz:

 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 

Weibull dağılımlarını farklı şekil ve ölçek parametreleriyle karşılaştırmak için grafiğe birkaç eğri de ekleyebiliriz:

 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)

Aşağıdaki sözdizimini alan legend() işlevini kullanarak çizime bir efsane ekleyebiliriz:

açıklama(x, y=NULL, açıklama, dolgu, sütun, bg, lty, cex)

Altın:

• x, y: açıklamayı konumlandırmak için kullanılan x ve y koordinatları
• efsane: efsaneye eklenecek metin
• dolgu: açıklamanın içindeki rengi doldur
• col: açıklamanın içindeki çizgiler için kullanılacak renklerin listesi
• bg: efsanenin arka plan rengi
• lty: çizgi stili
• cex: efsanedeki metnin boyutu

Örneğimizde bir gösterge oluşturmak için aşağıdaki sözdizimini kullanacağız:

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