การแจกแจงปัวซง
| การแจกแจงปัวซง | |||
|---|---|---|---|
ฟังก์ชันมวลความน่าจะเป็น แกนแนวนอนคือดัชนีkซึ่งเป็นจำนวนครั้งที่เกิดขึ้นλคืออัตราการเกิดเหตุการณ์ที่คาดหวัง แกนแนวตั้งคือความน่าจะเป็นของการเกิดเหตุการณ์k ครั้ง เมื่อกำหนดค่า λแล้ว ฟังก์ชันนี้กำหนดได้เฉพาะที่ค่าk เป็นจำนวนเต็ม เท่านั้น เส้นเชื่อมระหว่างแกนเป็นเพียงเส้นนำสายตาเท่านั้น | |||
ฟังก์ชันการกระจายสะสม แกนแนวนอนคือดัชนีkซึ่งเป็นจำนวนครั้งที่เกิดขึ้น ฟังก์ชันการกระจายสะสม (CDF) จะไม่ต่อเนื่องที่จำนวนเต็มkและจะราบเรียบที่จุดอื่นๆ เนื่องจากตัวแปรที่มีการกระจายแบบปัวซงจะมีค่าเป็นจำนวนเต็มเท่านั้น | |||
| สัญกรณ์ | |||
| พารามิเตอร์ | (ประเมิน) | ||
| สนับสนุน | ( จำนวนธรรมชาติเริ่มต้นจาก 0) | ||
| พีเอ็มเอฟ | |||
| ซีดีเอฟ | !}},} หรือหรือ (สำหรับที่ไหนคือฟังก์ชันแกมมาไม่สมบูรณ์บนคือฟังก์ชันพื้นและคือฟังก์ชันแกมมาแบบปรับค่าแล้ว ) | ||
| หมายถึง | |||
| ค่ามัธยฐาน | |||
| โหมด | |||
| ความแปรปรวน | |||
| ความเบี่ยงเบน | |||
| ความโค้งส่วนเกิน | |||
| เอนโทรปี | หรือสำหรับขนาดใหญ่ | ||
| เอ็มจีเอฟ | |||
| ซีเอฟ | |||
| พีจีเอฟ | |||
| ข้อมูลของฟิชเชอร์ | |||
ในทฤษฎีความน่าจะเป็นและสถิติการแจกแจงปัวซง ( / ˈ p w ɑː s ɒ n / ) เป็นการแจกแจงความน่าจะเป็นแบบไม่ต่อเนื่อง ที่แสดงความน่าจะเป็นของการเกิด เหตุการณ์จำนวนหนึ่งในช่วงเวลาที่กำหนด หากเหตุการณ์เหล่านี้เกิดขึ้นด้วยอัตราเฉลี่ยคงที่ที่ทราบ และเป็นอิสระจากเวลาตั้งแต่เหตุการณ์สุดท้าย[ 1 ]นอกจากนี้ยังสามารถใช้สำหรับจำนวนเหตุการณ์ในช่วงเวลาประเภทอื่นที่ไม่ใช่เวลา และในมิติที่มากกว่า 1 (เช่น จำนวนเหตุการณ์ในพื้นที่หรือปริมาตรที่กำหนด) การแจกแจงปัวซงตั้งชื่อตามนักคณิตศาสตร์ชาวฝรั่งเศสSiméon Denis Poissonมีบทบาทสำคัญในการแจกแจงแบบเสถียรไม่ต่อเนื่อง
ภายใต้การแจกแจงปัวซงที่คาดหวังเหตุการณ์λในช่วงเวลาที่กำหนด ความน่าจะเป็นของ เหตุการณ์ kในช่วงเวลาเดียวกันคือ: [ 2 ] : 60 ตัวอย่างเช่น พิจารณาศูนย์บริการลูกค้าที่รับสายโดยเฉลี่ยλ = 3สายต่อนาที ตลอดทั้งวัน หากจำนวนสายที่ได้รับในช่วงเวลาสองช่วงใดๆ ที่ไม่ซ้ำกันเป็นอิสระต่อกัน จำนวนสายkที่ได้รับในนาทีใดๆ จะมีการแจกแจงความน่าจะเป็นแบบปัวซง การรับสาย k = 1 ถึง 4สาย จะมีความน่าจะเป็นประมาณ 0.77 ในขณะที่การรับสาย 0 หรืออย่างน้อย 5 สาย จะมีความน่าจะเป็นประมาณ 0.23
ตัวอย่างคลาสสิกที่ใช้เพื่อกระตุ้นให้เกิดการแจกแจงแบบปัวซงคือจำนวน เหตุการณ์ การสลายตัวของกัมมันตรังสีในช่วงระยะเวลาการสังเกตที่กำหนด[ 3 ]
ประวัติศาสตร์
การนำเสนอการแจกแจงปัวซงได้รับการยกย่องให้แก่นักคณิตศาสตร์และนักฟิสิกส์ชาวฝรั่งเศสSiméon Denis Poisson (1781–1840) ซึ่งตีพิมพ์พร้อมกับทฤษฎีความน่าจะเป็นของเขาในRecherches sur la probabilité des jugements en matière criminelle et en matière civile (1837) [ 4 ] : 205-207งานนี้ตั้งทฤษฎีเกี่ยวกับจำนวนการตัดสินผิดพลาดในประเทศหนึ่งๆ โดยมุ่งเน้นไปที่ตัวแปรสุ่มN บางตัว ที่นับจำนวนเหตุการณ์ที่เกิดขึ้นในช่วงเวลาที่กำหนด อย่างไรก็ตาม ผลลัพธ์ที่คล้ายกันนี้เคยมีมาก่อนแล้วใน ปี 1711 โดยAbraham de MoivreในDe Mensura Sortis seu; de Probabilitate Eventuum in Ludis a Casu Fortuito Pendentibus [ 5 ] : 219 [ 6 ] : 14-15 [ 7 ] : 193 [ 8 ] : 157สิ่งนี้ทำให้เป็นตัวอย่างของกฎของ Stiglerและกระตุ้นให้ผู้เขียนบางคนโต้แย้งว่าการแจกแจงปัวซงควรใช้ชื่อว่า de Moivre [ 9 ] [ 10 ]
ใน ปี พ.ศ. 2403 Simon Newcombได้ปรับการแจกแจงแบบปัวซงให้เข้ากับจำนวนดาวที่พบในหน่วยของพื้นที่[ 11 ] การประยุกต์ใช้ในทางปฏิบัติเพิ่มเติมเกิดขึ้นโดยLadislaus Bortkiewiczใน ปี พ.ศ. 2441 Bortkiewicz แสดงให้เห็นว่าความถี่ที่ทหารในกองทัพปรัสเซียถูกม้าเตะจนเสียชีวิตโดยอุบัติเหตุ นั้น สามารถจำลองได้อย่างดีด้วยการแจกแจงแบบปัวซง[ 12 ] : 23-25
คำจำกัดความ
ฟังก์ชันมวลความน่าจะเป็น
ตัวแปรสุ่ม แบบ ไม่ต่อเนื่อง Xกล่าวได้ว่ามีการแจกแจงแบบปัวซง (Poisson distribution) ที่มีพารามิเตอร์หากมีฟังก์ชันมวลความน่าจะเป็นที่กำหนดโดย: [ 2 ] : 60 ที่ไหน
- kคือจำนวนครั้งที่เกิดขึ้น ()
- eคือจำนวนของออยเลอร์ ()
- k ! = k ( k– 1) ··· (3)(2)(1)คือแฟกทอเรียล
จำนวนจริง บวกλเท่ากับค่าที่คาดหวังของXและเท่ากับความแปรปรวน ของมัน ด้วย[ 13 ]
การแจกแจงปัวซงสามารถนำไปใช้กับระบบที่มีเหตุการณ์ที่เป็นไปได้จำนวนมาก ซึ่งแต่ละเหตุการณ์เกิดขึ้นได้ยากจำนวนเหตุการณ์ดังกล่าวที่เกิดขึ้นในช่วงเวลาที่กำหนดนั้น ภายใต้เงื่อนไขที่เหมาะสม จะเป็นจำนวนสุ่มที่มีการแจกแจงแบบปัวซง
สมการสามารถปรับเปลี่ยนได้ หากแทนที่จะใช้จำนวนเหตุการณ์เฉลี่ยเราได้รับอัตราเฉลี่ยซึ่งเหตุการณ์ต่างๆ เกิดขึ้น จากนั้นและ: [ 14 ]
ตัวอย่าง
การแจกแจงแบบปัวซงอาจมีประโยชน์ในการจำลองเหตุการณ์ต่างๆ เช่น:
- จำนวนอุกกาบาตที่มีเส้นผ่านศูนย์กลางมากกว่าหนึ่งเมตรที่พุ่งชนโลกในหนึ่งปี
- จำนวนโฟตอนเลเซอร์ที่ตกกระทบตัวตรวจจับในช่วงเวลาที่กำหนด
- จำนวนนักเรียนที่ได้คะแนนต่ำและคะแนนสูงในการสอบ และ
- ตำแหน่งของข้อบกพร่องและการเคลื่อนตัวในวัสดุ
ตัวอย่างของการเกิดจุดสุ่มในอวกาศ ได้แก่ ตำแหน่งที่ดาวเคราะห์น้อยพุ่งชนโลก (2 มิติ) ตำแหน่งของความไม่สมบูรณ์ในวัสดุ (3 มิติ) และตำแหน่งของต้นไม้ในป่า (2 มิติ) [ 15 ]
ข้อสมมติฐานและความถูกต้อง
การแจกแจงปัวซงเป็นแบบจำลองที่เหมาะสมหากข้อสมมติฐานต่อไปนี้เป็นจริง:
- kเป็นจำนวนเต็มที่ไม่เป็นลบ ซึ่งระบุจำนวนครั้งที่เหตุการณ์เกิดขึ้นในช่วงเวลาหนึ่ง
- การเกิดเหตุการณ์หนึ่งไม่ได้ส่งผลต่อโอกาสที่จะเกิดเหตุการณ์ที่สอง
- อัตราเฉลี่ยของการเกิดเหตุการณ์นั้นไม่ขึ้นอยู่กับการเกิดเหตุการณ์ใดๆ
- เหตุการณ์สองเหตุการณ์ไม่สามารถเกิดขึ้นพร้อมกันในเวลาเดียวกันได้อย่างแม่นยำ
ถ้าเงื่อนไขเหล่านี้เป็นจริง แสดงว่าkเป็นตัวแปรสุ่มปัวซง และการแจกแจงของkก็เป็นการแจกแจงปัวซงเช่นกัน
การแจกแจงปัวซงยังเป็นลิมิตของการแจกแจงทวินามด้วย ซึ่งความน่าจะเป็นของความสำเร็จในแต่ละครั้งของการทดลองคือ, ที่ไหนคือความคาดหวังและคือจำนวนครั้งของการทดลอง ในขีดจำกัดที่กับคงที่ [ 16 ] [ 17 ] (ดูการแจกแจงที่เกี่ยวข้อง ):
การแจกแจงปัวซงอาจได้มาจากสมการเชิงอนุพันธ์ด้วยเช่นกัน[ 18 ] [ 19 ] [ 20 ]
โดยมีเงื่อนไขเริ่มต้นและได้รับการประเมินที่
ตัวอย่างความน่าจะเป็นสำหรับการแจกแจงปัวซง
ในแม่น้ำสายหนึ่ง น้ำท่วมล้นตลิ่งเกิดขึ้นโดยเฉลี่ย ทุกๆ 100 ปี จงคำนวณความน่าจะเป็นของการเกิดน้ำท่วมล้นตลิ่งจำนวน k = 0, 1, 2, 3, 4, 5 หรือ6 ครั้ง ในช่วงเวลา 100 ปี โดยสมมติว่าแบบจำลองปัวซงเหมาะสม เนื่องจากอัตราการเกิดเหตุการณ์โดยเฉลี่ยคือน้ำท่วมล้นหนึ่งครั้งต่อ 100 ปี ดังนั้นλ = 1 |
ความน่าจะเป็นของ การเกิดน้ำท่วมล้นตลิ่ง 0 ถึง 6 ครั้ง ใน ช่วงเวลา 100 ปี |
ในตัวอย่างนี้ มีรายงานว่าจำนวนประตูเฉลี่ยในการแข่งขันฟุตบอลโลกอยู่ที่ประมาณ 2.5 และแบบจำลองปัวซงมีความเหมาะสม[ 21 ] เนื่องจากอัตราการเกิดเหตุการณ์เฉลี่ยอยู่ที่ 2.5 ประตูต่อแมตช์ ดังนั้นλ = 2.5 |
โอกาสที่จะเกิดประตูตั้งแต่ 0 ถึง 7 ประตูในหนึ่งแมตช์ |
ตัวอย่างที่ละเมิดข้อสมมติฐานของปัวซง
จำนวนนักเรียนที่มาถึงอาคารสโมสรนักศึกษาต่อนาทีนั้น อาจจะไม่เป็นไปตามการแจกแจงแบบปัวซง เนื่องจากอัตราการมาถึงไม่คงที่ (อัตราต่ำในช่วงเวลาเรียน อัตราสูงระหว่างช่วงเวลาเรียน) และการมาถึงของนักเรียนแต่ละคนก็ไม่เป็นอิสระต่อกัน (นักเรียนมักจะมาเป็นกลุ่ม) อัตราการมาถึงที่ไม่คงที่นี้อาจจำลองได้ด้วยการแจกแจงแบบปัวซงผสมและการมาถึงของกลุ่มนักเรียนแทนที่จะเป็นนักเรียนแต่ละคนอาจ จำลองได้ด้วย การแจกแจงแบบปัวซงเชิงประกอบ
จำนวน แผ่นดินไหวขนาด 5 ริกเตอร์ต่อปีในประเทศหนึ่ง อาจไม่เป็นไปตามการกระจายแบบปัวซง หากแผ่นดินไหวขนาดใหญ่ครั้งหนึ่งเพิ่มความน่าจะเป็นของการเกิดแผ่นดินไหวตามมาที่มีขนาดใกล้เคียงกัน
ตัวอย่างที่รับประกันว่าจะมีเหตุการณ์อย่างน้อยหนึ่งเหตุการณ์เกิดขึ้นนั้น ไม่ได้มีการแจกแจงแบบปัวซง แต่สามารถจำลองได้โดยใช้การแจกแจงแบบปัวซงที่ตัดศูนย์ออก
การแจกแจงจำนวนนับที่จำนวนช่วงเวลาที่มีเหตุการณ์เป็นศูนย์สูงกว่าที่คาดการณ์โดยแบบจำลองปัวซง อาจสร้างแบบจำลองได้โดยใช้แบบจำลองที่มีค่าศูนย์มากเกินไป (zero-inflated model )
คุณสมบัติ
สถิติเชิงพรรณนา
- ค่าที่คาดหวังของตัวแปรสุ่มปัวซงคือλ
- ค่าความแปรปรวนของตัวแปรสุ่มปัวซงก็คือλ เช่น กัน
- ค่าสัมประสิทธิ์ความแปรผันคือในขณะที่ดัชนีการกระจายตัวคือ 1 [ 8 ] : 163
- ค่าเบี่ยงเบนสัมบูรณ์เฉลี่ยรอบค่าเฉลี่ยคือ[ 8 ] : 163 !}}.}
- ค่าฐานนิยม ของตัวแปรสุ่มที่มีการแจกแจงแบบปัวซง โดยที่ λไม่ใช่จำนวนเต็มมีค่าเท่ากับซึ่งเป็นจำนวนเต็มที่มากที่สุดที่น้อยกว่าหรือเท่ากับλเขียนได้อีกแบบว่าfloor ( λ )เมื่อλเป็นจำนวนเต็มบวก ค่าฐานนิยมคือλและλ − 1
- ค่าสะสมทั้งหมดของการกระจายปัวซงมีค่าเท่ากับค่าคาดหวังλ โมเมนต์แฟกทอเรียลลำดับที่n ของการกระจายปัวซงคือλ n
- ค่าที่คาดหวังของกระบวนการปัวซงบางครั้งจะถูกแยกออกเป็นผลคูณของความเข้มและการสัมผัส (หรือโดยทั่วไปแล้วจะแสดงเป็นปริพันธ์ของ "ฟังก์ชันความเข้ม" ในช่วงเวลาหรือพื้นที่ ซึ่งบางครั้งเรียกว่า "การสัมผัส") [ 22 ]
ค่ามัธยฐาน
ขอบเขตสำหรับค่ามัธยฐาน () ของการกระจายตัวเป็นที่ทราบและมีความคมชัด : [ 23 ]
ช่วงเวลาที่สูงกว่า
โมเมนต์ ที่ไม่เป็นศูนย์กลางที่สูงกว่าm ของการแจกแจงปัวซงคือพหุนาม Touchardในλ : โดยที่วงเล็บปีกกา { } หมายถึงจำนวนสเตอร์ลิงชนิดที่สอง [ 24 ] [ 1 ] : 6กล่าวอีกนัยหนึ่ง เมื่อกำหนดค่าที่คาดหวังเป็นλ = 1 สูตรของโดบินสกิจะบ่งชี้ว่า โมเมนต์ลำดับที่ nเท่ากับจำนวน พาร์ ติชันของเซตที่มีขนาดn
ขอบเขตบนแบบง่ายคือ: [ 25 ]
ผลรวมของตัวแปรสุ่มที่มีการแจกแจงแบบปัวซง
ถ้าสำหรับถ้า เป็นอิสระต่อกันก็แสดงว่า...[ 26 ] : 65บทกลับคือทฤษฎีบทของ Raikovซึ่งกล่าวว่า ถ้าผลรวมของตัวแปรสุ่มอิสระสองตัวมีการแจกแจงแบบปัวซง ตัวแปรสุ่มอิสระทั้งสองตัวนั้นก็จะมีการแจกแจงแบบปัวซงเช่นกัน [ 27 ] [ 28 ]
เอนโทรปีสูงสุด
It is a maximum-entropy distribution among the set of generalized binomial distributions with mean and ,[29] where a generalized binomial distribution is defined as a distribution of the sum of N independent but not identically distributed Bernoulli variables.
Other properties
- The Poisson distributions are infinitely divisible probability distributions.[30]:233[8]:164
- The directed Kullback–Leibler divergence of from is given by
- If is an integer, then satisfies and [31]
- Bounds for the tail probabilities of a Poisson random variable can be derived using a Chernoff bound argument.[32]:97-98
- The upper tail probability can be tightened (by a factor of at least two) as follows:[33] where is the Kullback–Leibler divergence of from .
- Inequalities that relate the cumulative distribution function of a Poisson random variable to the cumulative distribution function of the standard normal distribution are as follows:[34]
where is the Kullback–Leibler divergence of from and is the Kullback–Leibler divergence of from .
Poisson races
Let and be independent random variables, with then we have that
The upper bound is proved using a standard Chernoff bound.
The lower bound can be proved by noting that is the probability that where which is bounded below by where is relative entropy (See the entry on bounds on tails of binomial distributions for details). Further noting that and computing a lower bound on the unconditional probability gives the result. More details can be found in the appendix of Kamath et al.[35]
Related distributions
As a Binomial distribution with infinitesimal time-steps
การแจกแจงปัวซงสามารถหาได้จากกรณีจำกัดของการแจกแจงทวินามเมื่อจำนวนการทดลองเข้าสู่ค่าอนันต์ และจำนวนความสำเร็จที่คาดหวัง ยังคงคงที่ — ดู หลักการของเหตุการณ์หายากด้านล่าง ดังนั้นจึงสามารถใช้เป็นการประมาณค่าของการแจกแจงทวินามได้หากnมีขนาดใหญ่พอและpมีขนาดเล็กพอ การแจกแจงปัวซงเป็นการประมาณค่าที่ดีของการแจกแจงทวินามหากnมีค่าอย่างน้อย 20 และpมีค่าน้อยกว่าหรือเท่ากับ 0.05 และเป็นการประมาณค่าที่ดีเยี่ยมหากn ≥ 100และnp ≤ 10 [ 36 ]ให้และโดยที่ เป็น ฟังก์ชันความหนาแน่นสะสมของ1การแจกแจงทวินามและการแจกแจงปัวซง ตามลำดับ จะได้ว่า: การอนุมานอย่างหนึ่งของสิ่งนี้ใช้ฟังก์ชันสร้างความน่าจะเป็น [ 37 ] พิจารณาการทดลองแบบเบอร์นูลลี (การโยนเหรียญ) ซึ่งความน่าจะเป็นของความสำเร็จหนึ่งครั้ง (หรือจำนวนความสำเร็จที่คาดหวัง) คือภายในช่วงเวลาที่กำหนด แบ่งช่วงเวลานั้นออกเป็นnส่วน และทำการทดลองในแต่ละช่วงย่อยด้วยความน่าจะเป็นความน่าจะเป็นของ ความสำเร็จ kครั้งจาก การทดลอง nครั้งตลอดช่วงเวลาทั้งหมด จะกำหนดโดยการแจกแจงทวินาม ซึ่งฟังก์ชันก่อกำเนิดคือ: เมื่อพิจารณาลิมิตเมื่อnเพิ่มขึ้นเป็นอนันต์ (โดยที่xคงที่) และใช้คำนิยามลิมิตผลคูณของฟังก์ชันเลขชี้กำลังจะได้เป็นฟังก์ชันก่อกำเนิดของการแจกแจงปัวซง:
ทั่วไป
- ถ้าและเป็นอิสระต่อกัน ดังนั้นความแตกต่างก็คือเป็นไปตามการแจกแจงแบบสเกลแลม (Skellam distribution )
- ถ้าและหากเป็นอิสระต่อกัน การกระจายตัวของโดยมีเงื่อนไขว่าเป็นการแจกแจงแบบทวินามโดยเฉพาะอย่างยิ่ง ถ้าแล้วโดยทั่วไปแล้ว ถ้าX , X , ..., X เป็นตัวแปรสุ่มปัวซงอิสระที่มีพารามิเตอร์λ , λ , ..., λ แล้วที่ให้ไว้ดังนั้นจึงสรุปได้ว่าในความเป็นจริง,
- ถ้าและการกระจายตัวของโดยมีเงื่อนไขว่าX = kเป็นการแจกแจงแบบทวินามดังนั้นการกระจายตัวของ Y จึงเป็นไปตามการกระจายแบบปัวซงในความเป็นจริง ถ้าหากพิจารณาจากเงื่อนไขบางประการเป็นไปตามการแจกแจงแบบพหุนามจากนั้นแต่ละเป็นไปตามการแจกแจงปัวซงแบบอิสระ
- การแจกแจงปัวซงเป็นกรณีพิเศษของการแจกแจงปัวซงแบบผสมแบบไม่ต่อเนื่อง (หรือการแจกแจงปัวซงแบบกระตุก) ที่มีพารามิเตอร์เพียงตัวเดียว[ 38 ] [ 39 ]การแจกแจงปัวซงแบบผสมแบบไม่ต่อเนื่องสามารถอนุมานได้จากการแจกแจงจำกัดของการแจกแจงพหุนามแบบเอกตัวแปร นอกจากนี้ยังเป็นกรณีพิเศษของการแจกแจงปัวซงแบบผสมอีก ด้วย
- สำหรับค่าλ ที่มีขนาดใหญ่เพียงพอ (เช่นλ > 1000 ) การแจกแจงแบบปกติจะมีค่าเฉลี่ยλและความแปรปรวนλ (ส่วนเบี่ยงเบนมาตรฐาน)) เป็นการประมาณที่ดีเยี่ยมของค่าการแจกแจงปัวซงหากλมากกว่าประมาณ10 การแจกแจงปกติจะเป็นการประมาณที่ดีหาก มี การแก้ไขความต่อเนื่อง ที่เหมาะสม กล่าวคือ หากP( X ≤ x )โดยที่xเป็นจำนวนเต็มที่ไม่เป็นลบ ถูกแทนที่ด้วยP( X ≤ x + 0.5 )
- การแปลงค่าความแปรปรวนให้คงที่ : ถ้าจากนั้น[ 8 ] : 168และ[ 40 ] : 196ภายใต้การแปลงนี้ การลู่เข้าสู่ภาวะปกติ (ดังที่)การเพิ่มขึ้น) นั้นเร็วกว่าตัวแปรที่ไม่ได้รับการแปลงมากการแปลงที่ทำให้ความแปรปรวนคงที่ที่ซับซ้อนกว่าเล็กน้อยก็มีอยู่เช่นกัน[ 8 ] : 168หนึ่งในนั้นคือ การ แปลงAnscombe [ 41 ]ดูการแปลงข้อมูล (สถิติ)สำหรับการใช้งานการแปลงทั่วไปเพิ่มเติม
- ถ้าสำหรับทุกt > 0จำนวนการมาถึงในช่วงเวลา[ 0, t ]เป็นไปตามการแจกแจงปัวซงที่มีค่าเฉลี่ยλtแล้วลำดับของเวลาระหว่างการมาถึงจะเป็นตัวแปรสุ่มเอกซ์โพเนนเชียลที่เป็นอิสระและมีการแจกแจงเหมือนกัน โดยมีค่าเฉลี่ย1/ λ [ 42 ] : 317–319
- ฟังก์ชันการกระจายสะสมของการกระจายแบบปัวซงและไคกำลังสองมีความสัมพันธ์กันในลักษณะต่อไปนี้: [ 8 ] : 167และ[ 8 ] : 158
การประมาณค่าปัวซง
สมมติที่ไหนจากนั้น[ 43 ]มีการแจกแจงแบบพหุนามโดยมีเงื่อนไขว่า
หมายความว่า[ 32 ] : 101-102ในบรรดาสิ่งอื่นๆ สำหรับฟังก์ชันที่ไม่เป็นลบใดๆ ถ้าถ้ามีการแจกแจงแบบพหุนามแล้ว ที่ไหน
ปัจจัยของสามารถแทนที่ด้วย 2 ได้หากนอกจากนี้ยังสันนิษฐานว่ามีแนวโน้มเพิ่มขึ้นหรือลดลงอย่างต่อเนื่อง
การแจกแจงปัวซงแบบสองตัวแปร
การแจกแจงนี้ได้รับการขยายไปยังกรณีสองตัวแปร[ 44 ]ฟังก์ชันก่อกำเนิดสำหรับการแจกแจงนี้คือ กับ
การแจกแจงแบบมาร์จินัลคือการแจกแจงปัวซง ( θ )และการแจกแจงปัวซง ( θ )และค่าสัมประสิทธิ์สหสัมพันธ์ถูกจำกัดอยู่ในช่วง
วิธีง่ายๆ ในการสร้างการแจกแจงปัวซงแบบสองตัวแปรคือการนำการแจกแจงปัวซงอิสระสามแบบมาใช้ด้วยวิธีการแล้วจึงตั้งค่าฟังก์ชันความน่าจะเป็นของการแจกแจงปัวซงแบบสองตัวแปรคือ
การแจกแจงปัวซงแบบอิสระ
การแจกแจงปัวซงแบบอิสระ[ 45 ]ที่มีขนาดกระโดดและให้คะแนนเกิดขึ้นใน ทฤษฎี ความน่าจะเป็นอิสระในฐานะลิมิตของการสังเคราะห์อิสระ ซ้ำๆ เมื่อN → ∞
กล่าวอีกนัยหนึ่งคือ ให้เป็นตัวแปรสุ่มเพื่อให้มีค่าด้วยความน่าจะเป็นและค่า 0 ด้วยความน่าจะเป็นที่เหลืออยู่ สมมติด้วยว่าครอบครัวเป็นอิสระอย่างอิสระจากนั้นขีดจำกัดคือของกฎหมายของกำหนดโดยกฎปัวซงอิสระที่มีพารามิเตอร์
คำจำกัดความนี้คล้ายคลึงกับวิธีการหนึ่งในการได้มาซึ่งการแจกแจงปัวซงแบบคลาสสิกจากกระบวนการปัวซง (แบบคลาสสิก)
มาตรการที่เกี่ยวข้องกับกฎปัวซงอิสระจะได้รับจาก[ 46 ] ที่ไหน และได้รับการสนับสนุน
กฎนี้ยังปรากฏในทฤษฎีเมทริกซ์สุ่ม ในรูปของ กฎมาร์เชนโก-ปาสตูร์ค่าคุมูลันต์อิสระของกฎนี้เท่ากับ
การแปลงบางส่วนของกฎนี้
เราให้ค่าของการแปลงที่สำคัญบางอย่างของกฎปัวซงอิสระ การคำนวณสามารถพบได้ในหนังสือLectures on the Combinatorics of Free Probabilityโดย A. Nica และ R. Speicher [ 47 ]
การแปลง R ของกฎปัวซงอิสระกำหนดโดย
การแปลงโคชี (ซึ่งเป็นการแปลงกลับของการแปลงสติลต์เจส ) กำหนดโดย
การแปลง S กำหนดโดย ในกรณีที่
การอนุมานทางสถิติ
การประมาณค่าพารามิเตอร์
กำหนดให้มีค่าที่วัดได้จำนวนnค่า เป็นตัวอย่างสำหรับi = 1, ..., nเราต้องการประมาณค่าของพารามิเตอร์λของประชากรปัวซงที่สุ่มตัวอย่างมา การประมาณค่า ความน่าจะเป็นสูงสุดคือ[ 48 ]
เนื่องจากแต่ละการสังเกตมีค่าเฉลี่ยλดังนั้นค่าเฉลี่ยของตัวอย่างจึงมีค่าเฉลี่ย λ เช่นกัน ด้วยเหตุนี้ การประมาณค่าความน่าจะเป็นสูงสุดจึงเป็นตัวประมาณค่าที่ไม่เอนเอียงของλนอกจากนี้ยังเป็นตัวประมาณค่าที่มีประสิทธิภาพ เนื่องจากความแปรปรวนของมันบรรลุขอบเขตล่างของ Cramér–Rao (CRLB) [ 49 ]ดังนั้นจึงเป็นตัวประมาณค่าที่ไม่เอนเอียงที่มีความแปรปรวนต่ำสุด นอกจากนี้ยังสามารถพิสูจน์ได้ว่าผลรวม (และด้วยเหตุนี้ค่าเฉลี่ย ของตัวอย่าง เนื่องจากเป็นฟังก์ชันหนึ่งต่อหนึ่งของผลรวม) เป็นสถิติที่สมบูรณ์และเพียงพอสำหรับλ
เพื่อพิสูจน์ความเพียงพอ เราอาจใช้ทฤษฎีบทการแยกตัวประกอบพิจารณาการแบ่งฟังก์ชันความน่าจะเป็นมวลของการแจกแจงปัวซงร่วมสำหรับตัวอย่างออกเป็นสองส่วน: ส่วนหนึ่งที่ขึ้นอยู่กับตัวอย่างเพียงอย่างเดียว, เรียกว่าและเป็นสิ่งที่ขึ้นอยู่กับพารามิเตอร์และตัวอย่างเฉพาะผ่านฟังก์ชันเท่านั้นแล้วเป็นสถิติที่เพียงพอสำหรับ
เทอมแรกขึ้นอยู่กับเท่านั้นเทอมที่สองขึ้นอยู่กับตัวอย่างเท่านั้นผ่านทางดังนั้น,เพียงพอแล้ว
ในการหาค่าพารามิเตอร์λที่ทำให้ฟังก์ชันความน่าจะเป็นสำหรับประชากรแบบปัวซงมีค่าสูงสุด เราสามารถใช้ลอการิทึมของฟังก์ชันความน่าจะเป็นได้:
เราหาอนุพันธ์ของโดยสัมพันธ์กับλและเปรียบเทียบกับศูนย์:
การแก้สมการหาค่าλจะได้จุดนิ่ง
ดังนั้นλคือค่าเฉลี่ยของ ค่า k การหาเครื่องหมายของอนุพันธ์อันดับสองของLที่จุดนิ่งจะช่วยกำหนดว่าλเป็น ค่าสุดขั้วประเภทใด
การประเมินอนุพันธ์อันดับสองณ จุดวิเคราะห์จะได้:
ซึ่งเป็นค่าลบของnคูณด้วยส่วนกลับของค่าเฉลี่ยของk นิพจน์นี้จะมีค่าเป็นลบเมื่อค่าเฉลี่ยเป็นบวก ถ้าเงื่อนไขนี้เป็นจริง จุดนิ่งจะทำให้ฟังก์ชันความน่าจะเป็นมีค่าสูงสุด
เพื่อให้ครบถ้วนสมบูรณ์ตระกูลของการแจกแจงจะเรียกว่าสมบูรณ์ก็ต่อเมื่อหมายความว่าสำหรับทุกคนหากบุคคลนั้นเป็น iidแล้วเมื่อทราบการกระจายตัวที่เราต้องการตรวจสอบแล้ว ก็จะเห็นได้ง่ายว่าสถิตินั้นสมบูรณ์
เพื่อให้ความเสมอภาคนี้คงอยู่ต้องเป็น 0 ซึ่งเป็นผลมาจากข้อเท็จจริงที่ว่าไม่มีเงื่อนไขอื่นใดที่จะเป็น 0 สำหรับทุกกรณีในผลรวมและสำหรับค่าที่เป็นไปได้ทั้งหมดของเพราะฉะนั้น,สำหรับทุกคนหมายความว่าและสถิติดังกล่าวได้รับการพิสูจน์แล้วว่าสมบูรณ์
ช่วงความเชื่อมั่น
ช่วงความเชื่อมั่นสำหรับค่าเฉลี่ยของการแจกแจงปัวซงสามารถแสดงได้โดยใช้ความสัมพันธ์ระหว่างฟังก์ชันการแจกแจงสะสมของการแจกแจงปัวซงและการแจกแจงไคกำลังสองการแจกแจงไคกำลังสองมีความสัมพันธ์อย่างใกล้ชิดกับการแจกแจงแกมมาและนี่นำไปสู่การแสดงออกอีกแบบหนึ่ง เมื่อกำหนดค่าสังเกตkจากการแจกแจงปัวซงที่มีค่าเฉลี่ยμช่วงความเชื่อมั่นสำหรับμที่ระดับความเชื่อมั่น1 – αคือ
หรือเทียบเท่า
ที่ไหนคือฟังก์ชันควอนไทล์ (ซึ่งสอดคล้องกับพื้นที่หางด้านล่างp ) ของการแจกแจงไคกำลังสองที่มีnองศาอิสระ และคือฟังก์ชันควอนไทล์ของการแจกแจงแกมมาที่มีพารามิเตอร์รูปร่าง n และพารามิเตอร์มาตราส่วน 1 [ 8 ] : 176-178 [ 50 ]ช่วงนี้เป็น ' ช่วงที่แน่นอน ' ในแง่ที่ว่าความน่าจะเป็นของการครอบคลุมจะไม่น้อยกว่าค่าที่กำหนด1 – α
เมื่อไม่มีควอนไทล์ของการแจกแจงแกมมา ได้มีการเสนอการประมาณค่าที่แม่นยำสำหรับช่วงเวลาที่แน่นอนนี้ (โดยอิงจากการแปลงวิลสัน-ฮิลเฟอร์ตี้ ) ดังนี้: [ 51 ] ที่ไหนหมายถึงค่าเบี่ยงเบนมาตรฐานปกติที่มีพื้นที่หางด้านบนα / 2
สำหรับการประยุกต์ใช้สูตรเหล่านี้ในบริบทเดียวกันกับข้างต้น (โดยกำหนดให้มี ค่าที่วัดได้ nค่าk แต่ละค่าสุ่มมาจากการแจกแจงแบบปัวซงที่มีค่าเฉลี่ยλ ) จะต้องตั้งค่าดังนี้
คำนวณ ช่วงสำหรับμ = nλจากนั้นหาช่วงสำหรับλ
การอนุมานแบบเบย์เซียน
ในการอนุมานแบบเบย์เซียนไพรเออร์คู่ควบสำหรับพารามิเตอร์อัตราλของการแจกแจงปัวซงคือการแจกแจงแกมมา[ 52 ]ให้
แสดงว่าλมีการกระจายตามความหนาแน่น แกมมา gซึ่งกำหนดโดยพารามิเตอร์รูปร่างαและพารามิเตอร์มาตราส่วน ผกผัน β :
จากนั้น เมื่อกำหนดตัวอย่าง ค่าที่วัดได้ nค่าk เช่นเดียว กับก่อนหน้านี้และค่าความน่าจะเป็นก่อนหน้าเป็นGamma( α , β )ค่าความน่าจะเป็นภายหลังจะเป็นดังนี้
โปรดทราบว่าค่าเฉลี่ยภายหลังเป็นแบบเชิงเส้นและกำหนดโดย สามารถแสดงได้ว่าการแจกแจงแกมมาเป็นไพรเออร์เดียวที่ทำให้ค่าเฉลี่ยแบบมีเงื่อนไขเป็นเชิงเส้น ยิ่งไปกว่านั้น ยังมีผลลัพธ์ตรงกันข้ามที่ระบุว่า ถ้าค่าเฉลี่ยแบบมีเงื่อนไขใกล้เคียงกับฟังก์ชันเชิงเส้นในระยะทางที่มากกว่าการแจกแจงก่อนหน้าของλจะต้องใกล้เคียงกับการแจกแจงแกมมาใน ระยะ ทางของ Levy [ 53 ]
ค่าเฉลี่ยภายหลังE[ λ ]เข้าใกล้ค่าประมาณความน่าจะเป็นสูงสุดในขีดจำกัดเมื่อซึ่งเป็น ผลสืบเนื่องโดยตรงจากนิพจน์ทั่วไปของค่าเฉลี่ยของการแจกแจงแกมมา
การแจกแจงการทำนายภายหลัง สำหรับการ สังเกตเพิ่มเติมเพียงครั้งเดียวคือการแจกแจงทวินามเชิงลบ[ 54 ] : 53บางครั้งเรียกว่าการแจกแจงแกมมา-ปัวซง
การประมาณค่าเฉลี่ยปัวซงหลายค่าพร้อมกัน
สมมติคือชุดของตัวแปรสุ่มอิสระจากชุดของการแจกแจงปัวซง โดยแต่ละแบบมีพารามิเตอร์และเราต้องการประมาณค่าพารามิเตอร์เหล่านี้ จากนั้น Clevenson และ Zidek แสดงให้เห็นว่าภายใต้การสูญเสียข้อผิดพลาดกำลังสองแบบนอร์มาไลซ์เมื่อไรจากนั้น เช่นเดียวกับในตัวอย่างของสไตน์สำหรับค่าเฉลี่ยปกติ ตัวประมาณค่า MLEไม่สามารถยอมรับได้[ 55 ]
ในกรณีนี้ จะมี การกำหนดตระกูลของตัวประมาณค่าแบบมินิแม็กซ์สำหรับทุก ๆและเช่น[ 56 ]
การเกิดขึ้นและการประยุกต์ใช้
การประยุกต์ใช้การแจกแจงปัวซงกับข้อมูลการนับ (จำนวนเหตุการณ์) บางประการ: [ 57 ]
- โทรคมนาคม : สายโทรศัพท์ที่เข้ามาในระบบ
- ดาราศาสตร์ : โฟตอนที่เดินทางมาถึงกล้องโทรทรรศน์
- เคมี : การกระจายมวลโมลาร์ของ การพอลิเม อไรเซชันที่มีชีวิต [ 58 ]
- ชีววิทยา : จำนวนการกลายพันธุ์บนสายดีเอ็นเอต่อหน่วยความยาว
- การจัดการ : ลูกค้าที่มาถึงเคาน์เตอร์หรือศูนย์บริการ[ 59 ]
- ด้านการเงินและการประกันภัย : จำนวนความเสียหายหรือการเรียกร้องค่าสินไหมทดแทนที่เกิดขึ้นในช่วงเวลาที่กำหนด
- แผ่นดินไหววิทยา : แบบจำลองปัวซงเชิงอะซิมโทติกของความเสี่ยงสำหรับแผ่นดินไหวขนาดใหญ่[ 60 ]
- กัมมันตภาพรังสี : สลายตัวในช่วงเวลาที่กำหนดในตัวอย่างกัมมันตภาพรังสี[ 61 ]
- ทัศนศาสตร์ : จำนวนโฟตอนที่ปล่อยออกมาในพัลส์เลเซอร์เดียว (จุดอ่อนสำคัญของ โปรโตคอล การกระจายกุญแจควอนตัมซึ่งรู้จักกันในชื่อการแบ่งจำนวนโฟตอน)
ตัวอย่างเพิ่มเติมของการนับเหตุการณ์ที่สามารถจำลองได้ด้วยกระบวนการปัวซง ได้แก่:
- ทหารที่เสียชีวิตจากการถูกม้าเตะในแต่ละปีในแต่ละกอง ทหารม้า ของปรัสเซียตัวอย่างนี้ใช้ในหนังสือของLadislaus Bortkiewicz (1868–1931) [ 12 ] : 23-25
- เซลล์ยีสต์ที่ใช้ในการผลิต เบียร์ กินเนสส์ตัวอย่างนี้ใช้โดยWilliam Sealy Gosset (1876–1937) [ 62 ] [ 63 ]
- สายโทรศัพท์ที่มาถึงศูนย์บริการภายในหนึ่งนาที ตัวอย่างนี้ได้รับการอธิบายโดยAK Erlang (1878–1929) [ 64 ]
- เป้าหมายในกีฬาที่เกี่ยวข้องกับทีมแข่งขันสองทีม[ 65 ]
- จำนวนผู้เสียชีวิตต่อปีในกลุ่มอายุที่กำหนด[ 66 ]
- การเปลี่ยนแปลงอย่างฉับพลันของราคาหุ้นในช่วงเวลาที่กำหนด
- จำนวนครั้ง ที่มีการเข้าถึง เว็บเซิร์ฟเวอร์ต่อนาที (ภายใต้สมมติฐานว่ามีความสม่ำเสมอ )
- การกลายพันธุ์ในดีเอ็นเอ ส่วนหนึ่ง หลังจากได้รับรังสีในปริมาณที่กำหนด
- เซลล์ที่ติดเชื้อที่ระดับความเข้มข้นของการติดเชื้อที่ กำหนด
- แบคทีเรียในของเหลวปริมาณหนึ่ง[ 67 ]
- โฟตอนที่มาถึงวงจรพิกเซลภายใต้ความสว่างที่กำหนดในช่วงเวลาที่กำหนด[ 68 ]
- การลงจอดของระเบิดบิน V-1ในลอนดอนระหว่างสงครามโลกครั้งที่ 2 ได้รับการตรวจสอบโดย RD Clarke ใน ปี พ.ศ. 2489 [ 69 ]
ในทฤษฎีจำนวนเชิงความน่าจะเป็นแกลลาเกอร์แสดงให้เห็นใน ปี 1976 ว่า หากสมมติฐานr-tuple ของจำนวนเฉพาะ ที่ยังไม่ได้ รับ การพิสูจน์บางเวอร์ชันเป็นจริง [ 70 ] แล้วจำนวนเฉพาะในช่วงเวลาสั้นๆ จะเป็นไปตามการแจกแจงปัวซง[ 71 ]
กฎของเหตุการณ์หายาก

อัตราการเกิดเหตุการณ์มีความสัมพันธ์กับความน่าจะเป็นที่เหตุการณ์นั้นจะเกิดขึ้นในช่วงเวลาย่อยเล็กๆ (ไม่ว่าจะเป็นเวลา สถานที่ หรืออื่นๆ) ในกรณีของการแจกแจงปัวซง เราจะสมมติว่ามีช่วงเวลาย่อยที่เล็กพอซึ่งความน่าจะเป็นที่เหตุการณ์จะเกิดขึ้นสองครั้งนั้น "น้อยมาก" ด้วยสมมติฐานนี้ เราสามารถสร้างการแจกแจงปัวซงจากการแจกแจงทวินามได้ โดยอาศัยเพียงข้อมูลเกี่ยวกับจำนวนเหตุการณ์ทั้งหมดที่คาดหวังในช่วงเวลาทั้งหมดเท่านั้น
ให้ แทนจำนวนเหตุการณ์ทั้งหมดในช่วงเวลาทั้งหมดด้วยแบ่งช่วงเวลาทั้งหมดออกเป็นช่วงย่อยมีขนาดเท่ากัน โดยที่(เนื่องจากเราสนใจเฉพาะส่วนเล็ก ๆ ของช่วงเวลาเท่านั้น สมมติฐานนี้จึงมีความหมาย) ซึ่งหมายความว่าจำนวนเหตุการณ์ที่คาดหวังในแต่ละ ช่วงย่อยทั้ง nช่วงนั้นเท่ากับ
ตอนนี้เราสมมติว่าการเกิดเหตุการณ์ในช่วงเวลาทั้งหมดสามารถมองได้ว่าเป็นลำดับของการทดลองแบบเบอร์นูลลีจำนวน n ครั้ง โดยที่การทดลองแบบเบอร์นูลลีครั้งที่ - สอดคล้องกับการตรวจสอบว่าเหตุการณ์เกิดขึ้นในช่วงเวลาย่อยหรือไม่ด้วยความน่าจะเป็นจำนวนเหตุการณ์ทั้งหมดที่คาดการณ์ไว้ในการพิจารณาคดีดังกล่าวจะเป็นจำนวนเหตุการณ์ทั้งหมดที่คาดว่าจะเกิดขึ้นในช่วงเวลาทั้งหมด ดังนั้น สำหรับแต่ละส่วนย่อยของช่วงเวลา เราได้ประมาณการการเกิดเหตุการณ์เป็นกระบวนการเบอร์นูลลีในรูปแบบดังที่เราได้กล่าวไว้ก่อนหน้านี้ เราต้องการพิจารณาเฉพาะช่วงย่อยที่เล็กมากเท่านั้น ดังนั้นเราจึงใช้ลิมิตเป็นมีค่าเข้าสู่ค่าอนันต์
ในกรณีนี้การแจกแจงทวินามจะลู่เข้าสู่สิ่งที่เรียกว่าการแจกแจงปัวซงโดยทฤษฎีบทขีดจำกัดของปัวซง
ในตัวอย่างข้างต้นหลายๆ ตัวอย่าง เช่น จำนวนการกลายพันธุ์ในลำดับดีเอ็นเอที่กำหนด เหตุการณ์ที่ถูกนับนั้นแท้จริงแล้วเป็นผลลัพธ์ของการทดลองแบบไม่ต่อเนื่อง และจะสามารถจำลองได้อย่างแม่นยำยิ่งขึ้นโดยใช้การแจกแจงแบบทวินามนั่นคือ
ในกรณีเช่นนี้nจะมีค่ามากและpจะมีค่าน้อยมาก (ดังนั้นค่าเฉลี่ยnpจึงมีขนาดปานกลาง) จากนั้นจึงสามารถประมาณการกระจายตัวได้โดยใช้การกระจายแบบปัวซงซึ่งไม่ซับซ้อนกว่า
บางครั้งการประมาณค่านี้เรียกว่ากฎของเหตุการณ์หายาก [ 72 ] : 5 เนื่องจาก เหตุการณ์เบอร์นูลลีแต่ละรายการ จำนวน n รายการเกิดขึ้นได้ยาก
ชื่อ "กฎของเหตุการณ์หายาก" อาจทำให้เข้าใจผิดได้ เพราะจำนวนรวมของเหตุการณ์ที่ประสบความสำเร็จในกระบวนการปัวซงไม่จำเป็นต้องหายากเสมอไป หากค่าพารามิเตอร์npไม่น้อย ตัวอย่างเช่น จำนวนสายโทรศัพท์ที่โทรเข้ามายังตู้รับสายที่ไม่ว่างในหนึ่งชั่วโมงนั้น เป็นไปตามการแจกแจงแบบปัวซง โดยเหตุการณ์เหล่านั้นดูเหมือนจะเกิดขึ้นบ่อยสำหรับผู้ให้บริการ แต่ถือว่าหายากจากมุมมองของคนทั่วไปที่แทบจะไม่โทรเข้ามายังตู้รับสายนั้นในชั่วโมงดังกล่าวเลย
ความแปรปรวนของการแจกแจงแบบทวินามเป็น1 − pเท่าของความแปรปรวนของการแจกแจงแบบปัวซง ดังนั้นจึงเกือบเท่ากันเมื่อpมีค่าเล็กมาก
บางครั้ง คำว่ากฎถูกใช้เป็นคำพ้องความหมายของการแจกแจงความน่าจะเป็นและการบรรจบกันในกฎหมายถึงการบรรจบกันในการแจกแจงดังนั้น การแจกแจงปัวซงจึงบางครั้งเรียกว่า "กฎของจำนวนน้อย" เพราะเป็นการแจกแจงความน่าจะเป็นของจำนวนครั้งของการเกิดเหตุการณ์ที่เกิดขึ้นไม่บ่อยนัก แต่มีโอกาสเกิดขึ้นมากมายกฎของจำนวนน้อยเป็นหนังสือของ Ladislaus Bortkiewicz เกี่ยวกับการแจกแจงปัวซง ตีพิมพ์ใน ปี 1898 [ 12 ] [ 73 ]
กระบวนการจุดปัวซง
การแจกแจงปัวซงเกิดขึ้นจากจำนวนจุดของกระบวนการจุดปัวซงที่อยู่ในบริเวณจำกัดบางแห่ง โดยเฉพาะอย่างยิ่ง ถ้าDเป็นปริภูมิบริเวณ เช่น ปริภูมิยุคลิดR dซึ่ง| D |พื้นที่ ปริมาตร หรือโดยทั่วไปแล้ว การวัดแบบเลเบสของบริเวณนั้นมีค่าจำกัด และถ้าN ( D )แทนจำนวนจุดในDแล้ว
การถดถอยปัวซงและการถดถอยทวินามเชิงลบ
การถดถอยแบบปัวซงและ การถดถอยแบบทวิ นามเชิงลบมีประโยชน์สำหรับการวิเคราะห์ที่ตัวแปรตาม (ตัวแปรตอบสนอง) คือจำนวนนับ(0, 1, 2, ... )ของจำนวนเหตุการณ์หรือการเกิดขึ้นในช่วงเวลาหนึ่ง
ชีววิทยา
การทดลอง ของลูเรีย-เดลบรุคได้ทดสอบกับสมมติฐานของการวิวัฒนาการแบบลามาร์ค ซึ่งควรจะส่งผลให้เกิดการแจกแจงแบบปัวซง
Katz และ Miledi วัดศักย์เยื่อหุ้มเซลล์ทั้งที่มีและไม่มีอะเซทิลโคลีน (ACh) [ 74 ]เมื่อมี ACh อยู่ช่องไอออนบนเยื่อหุ้มเซลล์จะเปิดแบบสุ่มในช่วงเวลาสั้นๆ เนื่องจากมีช่องไอออนจำนวนมากซึ่งแต่ละช่องเปิดในช่วงเวลาสั้นๆ จำนวนช่องไอออนทั้งหมดที่เปิดในแต่ละช่วงเวลาจึงมีการกระจายแบบปัวซง เมื่อไม่มี ACh อยู่ ในทางปฏิบัติแล้วจะไม่มีช่องไอออนเปิด ศักย์เยื่อหุ้มเซลล์คือเมื่อหักลบผลกระทบของสัญญาณรบกวนแล้ว คัตซ์และมิเลดีพบว่าค่าเฉลี่ยและความแปรปรวนของศักย์เยื่อหุ้มเซลล์คือและตามลำดับ โดยให้(หน้า 94-95 [ 75 ] )
ในระหว่างการจำลองเซลล์แต่ละครั้ง จำนวนการกลายพันธุ์จะกระจายตัวแบบปัวซงโดยประมาณ[ 76 ]ตัวอย่างเช่น ไวรัส HIV มี เบสคู่ 10,000 คู่ และมีอัตราการกลายพันธุ์ประมาณ 1 ครั้งต่อ เบสคู่ 30,000 คู่ ซึ่งหมายความว่าจำนวนการกลายพันธุ์ต่อการจำลองแต่ละครั้งจะกระจายตัวดังนี้(หน้า 64 [ 75 ] )
การประยุกต์ใช้ในด้านวิทยาศาสตร์อื่นๆ
ในกระบวนการปัวซง จำนวนเหตุการณ์ที่สังเกตได้จะผันผวนรอบค่าเฉลี่ยλโดยมีค่าเบี่ยงเบนมาตรฐานความผันผวนเหล่านี้เรียกว่าสัญญาณรบกวนปัวซงหรือ (โดยเฉพาะในด้านอิเล็กทรอนิกส์) สัญญาณรบกวนช็อต
ความสัมพันธ์ระหว่างค่าเฉลี่ยและค่าเบี่ยงเบนมาตรฐานในการนับเหตุการณ์ที่เกิดขึ้นแยกกันอย่างอิสระนั้นมีประโยชน์ทางวิทยาศาสตร์ โดยการตรวจสอบว่าความผันผวนเปลี่ยนแปลงไปอย่างไรเมื่อเทียบกับสัญญาณเฉลี่ย เราสามารถประมาณการส่วนร่วมของเหตุการณ์แต่ละครั้งได้แม้ว่าส่วนร่วมนั้นจะเล็กเกินกว่าจะตรวจจับได้โดยตรงก็ตาม ตัวอย่างเช่น ประจุeบนอิเล็กตรอนสามารถประมาณได้โดยการหาความสัมพันธ์ระหว่างขนาดของกระแสไฟฟ้ากับสัญญาณรบกวนแบบสุ่ม (shot noise ) ถ้ามี อิเล็กตรอน Nตัวผ่านจุดหนึ่งในช่วงเวลาtโดยเฉลี่ยกระแสไฟฟ้าเฉลี่ย คือเนื่องจากความผันผวนในปัจจุบันน่าจะอยู่ในระดับประมาณนี้(เช่น ค่าเบี่ยงเบนมาตรฐานของกระบวนการปัวซง ) ประจุสามารถประมาณได้จากอัตราส่วน
ตัวอย่างที่พบได้ทั่วไปคือความหยาบที่ปรากฏขึ้นเมื่อขยายภาพถ่าย ความหยาบนั้นเกิดจากความผันผวนแบบปัวซงของจำนวน อนุภาค เงิน ที่ลดลง ไม่ใช่เกิดจากอนุภาคแต่ละอนุภาคเอง การหาความสัมพันธ์ ระหว่าง ความหยาบกับระดับการขยายจะช่วยให้สามารถประเมินสัดส่วนการมีส่วนร่วมของอนุภาคแต่ละอนุภาคได้ (ซึ่งมีขนาดเล็กเกินกว่าจะมองเห็นได้ด้วยตาเปล่า)
ใน ทฤษฎี เซตเชิงสาเหตุองค์ประกอบที่ไม่ต่อเนื่องของปริภูมิเวลาจะเป็นไปตามการแจกแจงแบบปัวซงในปริมาตร
การแจกแจงแบบปัวซงยังปรากฏในกลศาสตร์ควอนตัมโดยเฉพาะอย่างยิ่งในทัศนศาสตร์ควอนตัมกล่าวคือ สำหรับ ระบบ ตัวสั่นฮาร์มอนิกควอนตัมในสถานะโคherentความน่าจะเป็นของการวัดระดับพลังงานเฉพาะนั้นมีการกระจายแบบปัวซง
วิธีการคำนวณ
การแจกแจงแบบปัวซงก่อให้เกิดภารกิจที่แตกต่างกันสองอย่างสำหรับไลบรารีซอฟต์แวร์เฉพาะทาง ได้แก่การประเมินการแจกแจงและสุ่มเลือกหมายเลขตามการกระจายนั้น
การประเมินการแจกแจงปัวซง
การคำนวณสำหรับที่กำหนดและเป็นงานเล็กน้อยที่สามารถทำได้โดยใช้คำจำกัดความมาตรฐานของในแง่ของฟังก์ชันเลขชี้กำลัง ฟังก์ชันกำลัง และฟังก์ชันแฟกทอเรียล อย่างไรก็ตาม นิยามทั่วไปของการแจกแจงปัวซงประกอบด้วยสองพจน์ที่อาจทำให้เกิดค่าเกินขีดจำกัดในคอมพิวเตอร์ได้ง่าย คือλ kและk !เศษส่วนของλ kต่อk !อาจทำให้เกิดข้อผิดพลาดในการปัดเศษซึ่งมีขนาดใหญ่มากเมื่อเทียบกับe − λและดังนั้นจึงให้ผลลัพธ์ที่ผิดพลาด เพื่อความเสถียรทางตัวเลข ฟังก์ชันความน่าจะเป็นมวลของปัวซงจึงควรได้รับการประเมินดังนี้ ซึ่งเทียบเท่าทางคณิตศาสตร์แต่มีเสถียรภาพเชิงตัวเลขสามารถหา ค่าลอการิทึมธรรมชาติของ ฟังก์ชันแกมมาlgamma ได้โดยใช้ ฟังก์ชันใน ไลบรารีมาตรฐานของภาษา C (เวอร์ชัน C99) หรือRฟังก์ชันgammalnในMATLABหรือSciPyหรือlog_gammaฟังก์ชันในFortran 2008 และเวอร์ชันที่ใหม่กว่า
ภาษาโปรแกรมบางภาษามีฟังก์ชันในตัวสำหรับประเมินการแจกแจงปัวซง ได้แก่
การสร้างตัวแปรสุ่ม
งานที่ยากกว่าคือการสุ่มเลือกตัวแปรสุ่ม จำนวนเต็ม จาก1การแจกแจงปัวซงโดยกำหนดค่าให้
โซลูชันนี้จัดทำโดย:
- R : ฟังก์ชัน
rpois(n, lambda); - ไลบรารีวิทยาศาสตร์ GNU (GSL): ฟังก์ชันgsl_ran_poisson
Knuthได้นำเสนออัลกอริธึมง่ายๆ เพื่อสร้างตัวเลขสุ่มที่มีการแจกแจงแบบปัวซง ( การสุ่มตัวอย่างตัวเลขสุ่มเทียม ) : [ 79 ] : 137-138
algorithmpoisson random number (Knuth): init: Let L ← e−λ, k ← 0 and p ← 1. do: k ← k + 1. Generate uniform random number u in [0,1] and let p ← p × u. while p > L. return k − 1.
The complexity is linear in the returned value k, which is λ on average. There are many other algorithms to improve this. Some are given in Ahrens & Dieter, see § References below.
For large values of λ, the value of L = e−λ may be so small that it is hard to represent. This can be solved by a change to the algorithm which uses an additional parameter STEP such that e−STEP does not underflow:
algorithmpoisson random number (Junhao, based on Knuth): init: LetλLeft ← λ, k ← 0 and p ← 1. do: k ← k + 1. Generate uniform random number u in (0,1) and let p ← p × u. while p < 1 and λLeft > 0: ifλLeft > STEP: p ← p × eSTEPλLeft ← λLeft − STEP else: p ← p × eλLeftλLeft ← 0 while p > 1. return k − 1.
The choice of STEP depends on the threshold of overflow. For double precision floating point format the threshold is near e700, so 500 should be a safe STEP.
Other solutions for large values of λ include rejection sampling and using Gaussian approximation.
Inverse transform sampling is simple and efficient for small values of λ, and requires only one uniform random number u per sample. Cumulative probabilities are examined in turn until one exceeds u.
algorithmPoisson generator based upon the inversion by sequential search:[80]:505init: Let x ← 0, p ← e−λ, s ← p. Generate uniform random number u in [0,1]. ในขณะที่ u > s ทำดังนี้: x ← x + 1. p ← p × λ / x. s ← s + p. ส่งคืนค่า x
ดูเพิ่มเติม
- การแจกแจงทวินาม
- การแจกแจงปัวซงแบบผสม
- การแจกแจงคอนเวย์-แม็กซ์เวลล์-ปัวซง
- การแจกแจงเออร์ลัง
- การแจกแจงแบบเอกซ์โปเนนเชียล
- การแจกแจงแกมมา
- การกระจายเฮอร์ไมต์
- ดัชนีการกระจายตัว
- การแจกแจงทวินามเชิงลบ
- การจับกลุ่มแบบปัวซง
- กระบวนการจุดปัวซง
- การถดถอยปัวซง
- การสุ่มตัวอย่างแบบปัวซง
- เวฟเล็ตปัวซง
- ทฤษฎีการเข้าคิว
- ทฤษฎีการต่ออายุ
- เลมมาของร็อบบินส์
- การจัดจำหน่ายสเกลลัม
- การกระจายตัวของทวีดี้
- แบบจำลองที่มีค่าศูนย์มากเกินไป
- การแจกแจงปัวซงแบบตัดศูนย์
Sources
- Ahrens, Joachim H.; Dieter, Ulrich (1974). "วิธีการทางคอมพิวเตอร์สำหรับการสุ่มตัวอย่างจากการกระจายแบบแกมมา เบต้า ปัวซง และทวินาม". Computing . 12 (3): 223– 246. Bibcode : 1974Compu..12..223A . doi : 10.1007/BF02293108 . S2CID 37484126 .
- Ahrens, Joachim H.; Dieter, Ulrich (1982). "การสร้างค่าเบี่ยงเบนปัวซงด้วยคอมพิวเตอร์" . ACM Transactions on Mathematical Software . 8 (2): 163– 179. doi : 10.1145/355993.355997 . S2CID 12410131 .
- Evans, Ronald J.; Boersma, J.; Blachman, NM; Jagers, AA (1988). "เอนโทรปีของการแจกแจงปัวซง: ปัญหา 87-6" . SIAM Review . 30 (2): 314– 317. doi : 10.1137/1030059 .