วิธีการพล็อตการแจกแจงแบบทวินามใน r

ในการพล็อตฟังก์ชันมวลของความน่าจะเป็นสำหรับ การแจกแจงแบบทวินาม ใน R เราสามารถใช้ฟังก์ชันต่อไปนี้:

• dbinom(x, size, prob) เพื่อสร้างฟังก์ชันมวลความน่าจะเป็น
• plot(x, y, type = ‘h’) เพื่อพล็อตฟังก์ชันมวลความน่าจะเป็น โดยระบุว่าพล็อตนั้นเป็นฮิสโตแกรม (type=’h’)

หากต้องการพล็อตฟังก์ชันมวลของความน่าจะเป็น เพียงระบุ ขนาด   (เช่น จำนวนการทดลอง) และ ปัญหา (เช่น ความน่าจะเป็นที่จะประสบความสำเร็จในการทดลองที่กำหนด) ในฟังก์ชัน dbinom()

ตัวอย่างเช่น โค้ดต่อไปนี้แสดงวิธีการพล็อตฟังก์ชันมวลความน่าจะเป็นสำหรับการแจกแจงแบบทวินามที่มีขนาด = 20 และ prob = 0.3:

 success <- 0:20

plot(success, dbinom(success, size=20, prob=.3), type='h')

แกน x แสดงจำนวนความสำเร็จ และแกน y แสดงความน่าจะเป็นที่จะประสบความสำเร็จตามจำนวนนั้นในการทดลอง 20 ครั้ง

เราสามารถเพิ่มชื่อเรื่อง เปลี่ยนป้ายกำกับแกน และเพิ่มความกว้างของเส้นเพื่อทำให้โครงเรื่องดูสวยงามยิ่งขึ้น:

 success <- 0:20

plot(success,dbinom(success,size=20,prob=.3),
     type='h',
     main='Binomial Distribution (n=20, p=0.3)',
     ylab='Probability',
     xlab ='# Successes',
     lwd=3)

คุณสามารถใช้โค้ดต่อไปนี้เพื่อดูความน่าจะเป็นจริงสำหรับความสำเร็จแต่ละจำนวนที่แสดงในแผนภูมิ:

 #prevent R from displaying numbers in scientific notation
options(scipen=999) 

#define range of successes
success <- 0:20

#display probability of success for each number of trials
dbinom(success, size=20, prob=.3)

[1] 0.00079792266297612 0.00683933711122388 0.02784587252426865
[4] 0.07160367220526231 0.13042097437387065 0.17886305056987975
[7] 0.19163898275344257 0.16426198521723651 0.11439673970486122
[10] 0.06536956554563482 0.03081708090008504 0.01200665489613703
[13] 0.00385928193090119 0.00101783259716075 0.00021810698510587
[16] 0.00003738976887529 0.00000500755833151 0.00000050496386536
[19] 0.00000003606884753 0.00000000162716605 0.00000000003486784

แหล่งข้อมูลเพิ่มเติม

ความรู้เบื้องต้นเกี่ยวกับการแจกแจงแบบทวินาม
ทำความเข้าใจรูปร่างของการแจกแจงแบบทวินาม