R'de bir t dağılımı nasıl çizilir

R’deki at 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 dt(x, df)
• olasılık yoğunluk fonksiyonunu çizmek için eğri(işlev, = NULL’dan = NULL’a)

Olasılık yoğunluk fonksiyonunu çizmek için dt() fonksiyonunda df’yi (serbestlik derecesi) ve curve() fonksiyonunda from ve to değerlerini belirtmemiz gerekir.

Örneğin, aşağıdaki kod, grafiğin x ekseninin -4 ile 4 arasında olduğu 10 serbestlik derecesine sahip bir dağılım için olasılık yoğunluk fonksiyonunun nasıl çizileceğini gösterir:

 curve(dt(x, df=10), from=-4, to=4)

Normal dağılıma benzer şekilde t dağılımı da ortalama 0 civarında simetriktir.

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(dt(x, df=10), from=-4, to=4, 
    main = 't Distribution (df = 10)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

Farklı serbestlik derecelerine sahip t dağılımlarını karşılaştırmak için grafiğe birden fazla eğri de ekleyebiliriz. Örneğin, aşağıdaki kod df = 6, df = 10 ve df = 30 olan t dağılım grafiklerini oluşturur:

 curve(dt(x, df=6), from=-4, to=4, col='blue') 
curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE)
curve(dt(x, df=30), from=-4, to=4, col='green', 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(dt(x, df=6), from=-4, to=4, col='blue') 
curve(dt(x, df=10), from=-4, to=4, col='red', add=TRUE)
curve(dt(x, df=30), from=-4, to=4, col='green', add=TRUE)

#add legend
legend(-4, .3, legend=c("df=6", "df=10", "df=30"),
       col=c("blue", "red", "green"), lty=1, cex=1.2)