Berikut ini akan diberikan kode-kode untuk membuat kurva dengan R. Sebagai bahan untuk praktek digunakan contoh fungsi kepadatan peluang, yaitu: \[f(x) = \left\{\begin{matrix} \displaystyle \frac{3}{125} x^2 & 0 \leq x \leq 5\\ \\ 0 & x \text{ lainnya.} \ \end{matrix}\right.\]
Langkah pertama yang kita lakukan untuk membuat kurva dari fungsi kepadatan peluang tersebut adalah membangkitkan (generate) datanya terlebih dahulu. Data yang dibangkitkan adalah data variabel \(x\), yaitu data \(x\) dari 0 sampai dengan 5 yang jarak interval datanya adalah 0,1. Kemudian data tersebut dimasukkan ke dalam fungsi kepadatan peluangnya \((f(x)=y).\) Kode-kodenya adalah sebagai berikut.
x <- seq(0, 5, 0.1)
y <- 3/125*x^2
Lakukan plotting data antara \(x\) dan \(y.\)
plot(x, y)
Hasil dari plotting antara \(x\) dan \(y\) tersebut berupa titik-titik. Selanjutnya hubungkan titik-titik tersebut sehingga membentuk garis.
plot(x, y, type = "l")
Luas area di bawah kurva fungsi kepadatan peluang adalah 1. Area di bawah kurva tersebut dapat kita arsir dengan menggunakan kode berikut.
x1 <- c(0, x, 5)
y1 <- c(0, y, 0)
polygon(x1, y1, col="grey", border=NA)
Jika hanya ingin mengarsir sebagian area saja, misalnya nilai \(x\) antara 2 dan 4, maka kode yang digunakan adalah sebagai berikut.
x1 <- seq(2, 4, 0.1)
y1 <- 3/125*x1^2
x1 <- c(2, x1, 4)
y1 <- c(0, y1, 0)
polygon(x1, y1, col="grey", border=NA)
Warna kurva juga bisa diganti.
plot(x, y, type = "l", col="red")
Garis kurva juga bisa dipertebal sesuai keinginan.
plot(x, y, type = "l", col="red", lwd=2)
Tambahkan judulnya.
plot(x, y, type = "l", col="red", lwd=2,
main="Kurva Fungsi Distribusi Peluang")
Ganti label variabel \(x\) dan \(y.\)
plot(x, y, type = "l", col="red", lwd=2,
main="Kurva Fungsi Distribusi Peluang",
xlab="Peubah acak X",
ylab="Fungsi peluang f(x)")
Kurva lain juga bisa ditambahkan.
z <- 3/125*(5-x)^2
lines(x, z, col="blue", lwd=2)
Legenda juga diperlukan untuk mendefinisikan kurva.
legend(3.75,0.3, c(expression(frac(3*x^2,125)), expression(frac(3*(5-x)^2,125))),
lwd=c(3,3), col=c("red", "blue"))
Hilangkan garis kotak legenda agar tampilannya lebih menarik.
legend(3.75,0.3, c(expression(frac(3*x^2,125)), expression(frac(3*(5-x)^2,125))),
lwd=c(3,3), col=c("red", "blue"), box.col=FALSE)
Sekian kode-kode untuk membuat kurva di R.