วิธีการพล็อตการแจกแจงปัวซองใน r

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

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

หากต้องการพล็อตฟังก์ชันมวลของความน่าจะเป็น เพียงระบุ lambda   (เช่น อัตราการเกิดเหตุการณ์) ในฟังก์ชัน dpois()

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

 #define range of "successes"
success <- 0:20

#create plot of probability mass function
plot(success, dois(success, lambda=5), type='h')

แกน x แสดงจำนวน “ความสำเร็จ” เช่น จำนวนเหตุการณ์ที่เกิดขึ้น และแกน y แสดงความน่าจะเป็นที่จะบรรลุจำนวนความสำเร็จนั้นในการทดลอง 20 ครั้ง

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

 success <- 0:20

plot(success, dois(success, lambda=5),
     type='h',
     main='Fish Distribution (lambda = 5)',
     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
dpois(success, lambda=5)

[1] 0.0067379469991 0.0336897349954 0.0842243374886 0.1403738958143
[5] 0.1754673697679 0.1754673697679 0.1462228081399 0.1044448629571
[9] 0.0652780393482 0.0362655774156 0.0181327887078 0.0082421766854
[13] 0.0034342402856 0.0013208616483 0.0004717363030 0.0001572454343
[17] 0.0000491391982 0.0000144527054 0.0000040146404 0.0000010564843
[21] 0.0000002641211