R - Distribusi Binomial

Model distribusi binomial berkaitan dengan menemukan probabilitas keberhasilan suatu peristiwa yang hanya memiliki dua kemungkinan hasil dalam serangkaian percobaan. Misalnya, melempar koin selalu menghasilkan kepala atau ekor. Probabilitas menemukan tepat 3 kepala dalam melempar koin berulang kali selama 10 kali diperkirakan selama distribusi binomial.

R memiliki empat fungsi built-in untuk menghasilkan distribusi binomial. Mereka dijelaskan di bawah ini.

dbinom(x, size, prob)
pbinom(x, size, prob)
qbinom(p, size, prob)
rbinom(n, size, prob)

Berikut ini adalah deskripsi parameter yang digunakan -

x adalah vektor angka.
p adalah vektor probabilitas.
n adalah jumlah pengamatan.
size adalah jumlah percobaan.
prob adalah probabilitas keberhasilan setiap percobaan.

dbinom ()

Fungsi ini memberikan distribusi kepadatan probabilitas di setiap titik.

# Create a sample of 50 numbers which are incremented by 1.
x <- seq(0,50,by = 1)

# Create the binomial distribution.
y <- dbinom(x,50,0.5)

# Give the chart file a name.
png(file = "dbinom.png")

# Plot the graph for this sample.
plot(x,y)

# Save the file.
dev.off()

Ketika kita menjalankan kode di atas, hasilnya adalah sebagai berikut -

pbinom ()

Fungsi ini memberikan probabilitas kumulatif dari suatu peristiwa. Ini adalah nilai tunggal yang mewakili probabilitas.

# Probability of getting 26 or less heads from a 51 tosses of a coin.
x <- pbinom(26,51,0.5)

print(x)

Ketika kita menjalankan kode di atas, hasilnya adalah sebagai berikut -

[1] 0.610116

qbinom ()

Fungsi ini mengambil nilai probabilitas dan memberikan angka yang nilai kumulatifnya cocok dengan nilai probabilitas.

# How many heads will have a probability of 0.25 will come out when a coin
# is tossed 51 times.
x <- qbinom(0.25,51,1/2)

print(x)

Ketika kita menjalankan kode di atas, hasilnya adalah sebagai berikut -

[1] 23

rbinom ()

Fungsi ini menghasilkan jumlah nilai acak yang diperlukan dari probabilitas tertentu dari sampel tertentu.

# Find 8 random values from a sample of 150 with probability of 0.4.
x <- rbinom(8,150,.4)

print(x)

Ketika kita menjalankan kode di atas, hasilnya adalah sebagai berikut -

[1] 58 61 59 66 55 60 61 67