Cara memplot distribusi t di r

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

• dt(x, df) untuk membuat fungsi kepadatan probabilitas
• curve(function, from = NULL, to = NULL) untuk memplot fungsi kepadatan probabilitas

Untuk memplot fungsi kepadatan probabilitas, kita perlu menentukan df (derajat kebebasan) dalam fungsi dt() serta nilai dari dan ke dalam fungsi curve() .

Misalnya, kode berikut mengilustrasikan cara memplot fungsi kepadatan probabilitas untuk distribusi dengan 10 derajat kebebasan dengan sumbu x plot antara -4 dan 4:

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

Mirip dengan distribusi normal , distribusi t simetris di sekitar mean 0.

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

 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 

Kita juga dapat menambahkan beberapa kurva pada grafik untuk membandingkan t distribusi dengan derajat kebebasan yang berbeda. Misalnya, kode berikut membuat t plot distribusi dengan df = 6, df = 10, dan df = 30:

 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)

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