R - การแจกแจงแบบทวินาม

แบบจำลองการแจกแจงแบบทวินามเกี่ยวข้องกับการค้นหาความน่าจะเป็นของความสำเร็จของเหตุการณ์ซึ่งมีผลลัพธ์ที่เป็นไปได้เพียงสองรายการในชุดการทดลอง ตัวอย่างเช่นการโยนเหรียญจะให้หัวหรือก้อยเสมอ ความน่าจะเป็นที่จะพบว่า 3 หัวในการโยนเหรียญซ้ำ ๆ เป็นเวลา 10 ครั้งนั้นถูกประมาณในระหว่างการแจกแจงแบบทวินาม

R มีฟังก์ชันในตัวสี่ฟังก์ชันเพื่อสร้างการแจกแจงแบบทวินาม มีอธิบายไว้ด้านล่าง

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

ต่อไปนี้เป็นคำอธิบายของพารามิเตอร์ที่ใช้ -

• x เป็นเวกเตอร์ของตัวเลข

• p คือเวกเตอร์ของความน่าจะเป็น

• n คือจำนวนการสังเกต

• size คือจำนวนการทดลอง

• prob คือความน่าจะเป็นของความสำเร็จของการทดลองแต่ละครั้ง

dbinom ()

ฟังก์ชันนี้ให้การแจกแจงความหนาแน่นของความน่าจะเป็นในแต่ละจุด

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

เมื่อเรารันโค้ดด้านบนจะให้ผลลัพธ์ดังนี้ -

pbinom ()

ฟังก์ชันนี้ให้ความน่าจะเป็นสะสมของเหตุการณ์ เป็นค่าเดียวที่แสดงถึงความน่าจะเป็น

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

print(x)

เมื่อเรารันโค้ดด้านบนจะให้ผลลัพธ์ดังนี้ -

[1] 0.610116

qbinom ()

ฟังก์ชันนี้รับค่าความน่าจะเป็นและให้ตัวเลขที่มีค่าสะสมตรงกับค่าความน่าจะเป็น

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

เมื่อเรารันโค้ดด้านบนจะให้ผลลัพธ์ดังนี้ -

[1] 23

rbinom ()

ฟังก์ชันนี้สร้างจำนวนค่าสุ่มที่ต้องการของความน่าจะเป็นที่กำหนดจากตัวอย่างที่กำหนด

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

print(x)

เมื่อเรารันโค้ดด้านบนจะให้ผลลัพธ์ดังนี้ -

[1] 58 61 59 66 55 60 61 67