กลับไปหน้าบทความ

อ่าน 49 นาที

การแจกแจงปัวซง

ตัวเลขสุ่มไม่เท่ากัน

Γ(⌊เค+1⌋,λ)⌊เค⌋!,{\displaystyle {\frac {\Gamma (\lfloor k+1\rfloor ,\lambda )}{\lfloor k\rfloor !}},} หรืออี−λ∑เจ=0⌊เค⌋λเจเจ!

การแจกแจงปัวซง

การแจกแจงปัวซง
ฟังก์ชันมวลความน่าจะเป็น
แกนแนวนอนคือดัชนีkซึ่งเป็นจำนวนครั้งที่เกิดขึ้นλคืออัตราการเกิดเหตุการณ์ที่คาดหวัง แกนแนวตั้งคือความน่าจะเป็นของการเกิดเหตุการณ์k ครั้ง เมื่อกำหนดค่า λแล้ว ฟังก์ชันนี้กำหนดได้เฉพาะที่ค่าk เป็นจำนวนเต็ม เท่านั้น เส้นเชื่อมระหว่างแกนเป็นเพียงเส้นนำสายตาเท่านั้น
ฟังก์ชันการกระจายสะสม
แกนแนวนอนคือดัชนีkซึ่งเป็นจำนวนครั้งที่เกิดขึ้น ฟังก์ชันการกระจายสะสม (CDF) จะไม่ต่อเนื่องที่จำนวนเต็มkและจะราบเรียบที่จุดอื่นๆ เนื่องจากตัวแปรที่มีการกระจายแบบปัวซงจะมีค่าเป็นจำนวนเต็มเท่านั้น
สัญกรณ์ปัวส์(λ){\displaystyle \operatorname {Pois} (\lambda )}
พารามิเตอร์λ(0,){\displaystyle \lambda \in (0,\infty )} (ประเมิน)
สนับสนุนเคเอ็น0{\displaystyle k\in \mathbb {N} _{0}}( จำนวนธรรมชาติเริ่มต้นจาก 0)
พีเอ็มเอฟλเคอีλเค!{\displaystyle {\frac {\แลมบ์ดา ^{k}e^{-\แลมบ์ดา }}{k!}}}
ซีดีเอฟ

Γ(เค+1,λ)เค!,{\displaystyle {\frac {\Gamma (\lfloor k+1\rfloor ,\lambda )}{\lfloor k\rfloor !}},} หรืออีλเจ=0เคλเจเจ!,{\displaystyle e^{-\lambda }\sum _{j=0}^{\lfloor k\rfloor }{\frac {\lambda ^{j}}{j!}},}หรือคิว(เค+1,λ){\displaystyle Q(\lfloor k+1\rfloor ,\lambda )}

(สำหรับเค0,{\displaystyle k\geq 0,}ที่ไหนΓ(x,y){\displaystyle \Gamma (x,y)}คือฟังก์ชันแกมมาไม่สมบูรณ์บนเค{\displaystyle \lfloor k\rfloor }คือฟังก์ชันพื้นและคิว{\displaystyle Q}คือฟังก์ชันแกมมาแบบปรับค่าแล้ว )
หมายถึงλ{\displaystyle \lambda }
ค่ามัธยฐานλ+13150λ{\displaystyle \approx \left\lfloor \lambda +{\frac {1}{3}}-{\frac {1}{50\lambda }}\right\rfloor }
โหมดλ1,λ{\displaystyle \left\lceil \lambda \right\rceil -1,\left\lfloor \lambda \right\rfloor }
ความแปรปรวนλ{\displaystyle \lambda }
ความเบี่ยงเบน1λ{\displaystyle {\frac {1}{\sqrt {\lambda }}}}
ความโค้งส่วนเกิน1λ{\displaystyle {\frac {1}{\lambda }}}
เอนโทรปี

λ[1ln(λ)]+อีλเค=0λเคln(เค!)เค!{\displaystyle \lambda {\Bigl [}1-\ln(\lambda ){\Bigr ]}+e^{-\lambda }\sum _{k=0}^{\infty }{\frac {\lambda ^{k}\ln(k!)}{k!}}}หรือสำหรับขนาดใหญ่λ{\displaystyle \lambda }

12ln(2πอีλ)112λ124λ219360λ3+โอ(1λ4){\displaystyle {\begin{aligned}\approx {\frac {1}{2}}\ln \left(2\pi e\lambda \right)-{\frac {1}{12\lambda }}-{\frac {1}{24\lambda ^{2}}}\\-{\frac {19}{360\lambda ^{3}}}+{\mathcal {O}}\left({\frac {1}{\lambda ^{4}}}\right)\end{aligned}}}
เอ็มจีเอฟเอ็กซ์[λ(อีที1)]{\displaystyle \exp \left[\lambda \left(e^{t}-1\right)\right]}
ซีเอฟเอ็กซ์[λ(อีฉันที1)]{\displaystyle \exp \left[\lambda \left(e^{it}-1\right)\right]}
พีจีเอฟเอ็กซ์[λ(z1)]{\displaystyle \exp \left[\lambda \left(z-1\right)\right]}
ข้อมูลของฟิชเชอร์1λ{\displaystyle {\frac {1}{\lambda }}}

ในทฤษฎีความน่าจะเป็นและสถิติการแจกแจงปัวซง ( / ˈ p w ɑː s ɒ n / ) เป็นการแจกแจงความน่าจะเป็นแบบไม่ต่อเนื่อง ที่แสดงความน่าจะเป็นของการเกิด เหตุการณ์จำนวนหนึ่งในช่วงเวลาที่กำหนด หากเหตุการณ์เหล่านี้เกิดขึ้นด้วยอัตราเฉลี่ยคงที่ที่ทราบ และเป็นอิสระจากเวลาตั้งแต่เหตุการณ์สุดท้าย[ 1 ]นอกจากนี้ยังสามารถใช้สำหรับจำนวนเหตุการณ์ในช่วงเวลาประเภทอื่นที่ไม่ใช่เวลา และในมิติที่มากกว่า 1 (เช่น จำนวนเหตุการณ์ในพื้นที่หรือปริมาตรที่กำหนด) การแจกแจงปัวซงตั้งชื่อตามนักคณิตศาสตร์ชาวฝรั่งเศสSiméon Denis Poissonมีบทบาทสำคัญในการแจกแจงแบบเสถียรไม่ต่อเนื่อง

ภายใต้การแจกแจงปัวซงที่คาดหวังเหตุการณ์λในช่วงเวลาที่กำหนด ความน่าจะเป็นของ เหตุการณ์ kในช่วงเวลาเดียวกันคือ: [ 2 ] : 60λเคอีλเค!.{\displaystyle {\frac {\lambda ^{k}e^{-\lambda }}{k!}}.} ตัวอย่างเช่น พิจารณาศูนย์บริการลูกค้าที่รับสายโดยเฉลี่ยλ = 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) ที่มีพารามิเตอร์λ>0{\displaystyle \lambda >0}หากมีฟังก์ชันมวลความน่าจะเป็นที่กำหนดโดย: [ 2 ] : 60เอฟ(เค;λ)=ปร.(X=เค)=λเคอีλเค!,{\displaystyle f(k;\lambda )=\Pr(X{=}k)={\frac {\lambda ^{k}e^{-\lambda }}{k!}},} ที่ไหน

จำนวนจริง บวกλเท่ากับค่าที่คาดหวังของXและเท่ากับความแปรปรวน ของมัน ด้วย[ 13 ]λ=อี(X)=วาร์(X).{\displaystyle \lambda =\operatorname {E} (X)=\operatorname {Var} (X).}

การแจกแจงปัวซงสามารถนำไปใช้กับระบบที่มีเหตุการณ์ที่เป็นไปได้จำนวนมาก ซึ่งแต่ละเหตุการณ์เกิดขึ้นได้ยากจำนวนเหตุการณ์ดังกล่าวที่เกิดขึ้นในช่วงเวลาที่กำหนดนั้น ภายใต้เงื่อนไขที่เหมาะสม จะเป็นจำนวนสุ่มที่มีการแจกแจงแบบปัวซง

สมการสามารถปรับเปลี่ยนได้ หากแทนที่จะใช้จำนวนเหตุการณ์เฉลี่ยλ,{\displaystyle \lambda ,}เราได้รับอัตราเฉลี่ย{\displaystyle r}ซึ่งเหตุการณ์ต่างๆ เกิดขึ้น จากนั้นλ=ที,{\displaystyle \lambda =rt,}และ: [ 14 ]พี(เค เหตุการณ์ในช่วงเวลา ที)=(ที)เคอีทีเค!.{\displaystyle P(k{\text{ events in interval }}t)={\frac {(rt)^{k}e^{-rt}}{k!}}.}

ตัวอย่าง

หมากฝรั่งถูกทิ้งไว้บนทางเท้าในเมืองเรคยาวิก
หมากฝรั่งบนทางเท้า จำนวนชิ้นหมากฝรั่งบนกระเบื้องแต่ละแผ่นมีลักษณะการกระจายแบบปัวซงโดยประมาณ

การแจกแจงแบบปัวซงอาจมีประโยชน์ในการจำลองเหตุการณ์ต่างๆ เช่น:

  • จำนวนอุกกาบาตที่มีเส้นผ่านศูนย์กลางมากกว่าหนึ่งเมตรที่พุ่งชนโลกในหนึ่งปี
  • จำนวนโฟตอนเลเซอร์ที่ตกกระทบตัวตรวจจับในช่วงเวลาที่กำหนด
  • จำนวนนักเรียนที่ได้คะแนนต่ำและคะแนนสูงในการสอบ และ
  • ตำแหน่งของข้อบกพร่องและการเคลื่อนตัวในวัสดุ

ตัวอย่างของการเกิดจุดสุ่มในอวกาศ ได้แก่ ตำแหน่งที่ดาวเคราะห์น้อยพุ่งชนโลก (2 มิติ) ตำแหน่งของความไม่สมบูรณ์ในวัสดุ (3 มิติ) และตำแหน่งของต้นไม้ในป่า (2 มิติ) [ 15 ]

ข้อสมมติฐานและความถูกต้อง

การแจกแจงปัวซงเป็นแบบจำลองที่เหมาะสมหากข้อสมมติฐานต่อไปนี้เป็นจริง:

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

ถ้าเงื่อนไขเหล่านี้เป็นจริง แสดงว่าkเป็นตัวแปรสุ่มปัวซง และการแจกแจงของkก็เป็นการแจกแจงปัวซงเช่นกัน

การแจกแจงปัวซงยังเป็นลิมิตของการแจกแจงทวินามด้วย ซึ่งความน่าจะเป็นของความสำเร็จในแต่ละครั้งของการทดลองคือพี=λn{\displaystyle p={\frac {\lambda }{n}}}, ที่ไหนλ{\displaystyle \lambda }คือความคาดหวังและn{\displaystyle n}คือจำนวนครั้งของการทดลอง ในขีดจำกัดที่n{\displaystyle n\to \infty }กับλ{\displaystyle \lambda }คงที่ [ 16 ] [ 17 ] (ดูการแจกแจงที่เกี่ยวข้อง ):

ลิมn(nเค)(λn)เค(1λn)nเค=λเคเค!อีλ{\displaystyle \lim _{n\to \infty }{\dbinom {n}{k}}\left({\frac {\lambda }{n}}\right)^{k}\,\left(1-{\frac {\lambda }{n}}\right)^{n-k}={\frac {\lambda ^{k}}{k!}}\,e^{-\lambda }}

การแจกแจงปัวซงอาจได้มาจากสมการเชิงอนุพันธ์ด้วยเช่นกัน[ 18 ] [ 19 ] [ 20 ]

พีเค(ที)ที=λ(พีเค1(ที)พีเค(ที)){\displaystyle {\frac {d\,P_{k}(t)}{dt}}=\lambda \,{\Big (}P_{k-1}(t)-P_{k}(t){\Big )}}

โดยมีเงื่อนไขเริ่มต้นพีเค(0)=δเค0{\displaystyle P_{k}(0)=\delta _{k0}}และได้รับการประเมินที่ที=1{\displaystyle t=1}

ตัวอย่างความน่าจะเป็นสำหรับการแจกแจงปัวซง

ตัวอย่างที่ละเมิดข้อสมมติฐานของปัวซง

จำนวนนักเรียนที่มาถึงอาคารสโมสรนักศึกษาต่อนาทีนั้น อาจจะไม่เป็นไปตามการแจกแจงแบบปัวซง เนื่องจากอัตราการมาถึงไม่คงที่ (อัตราต่ำในช่วงเวลาเรียน อัตราสูงระหว่างช่วงเวลาเรียน) และการมาถึงของนักเรียนแต่ละคนก็ไม่เป็นอิสระต่อกัน (นักเรียนมักจะมาเป็นกลุ่ม) อัตราการมาถึงที่ไม่คงที่นี้อาจจำลองได้ด้วยการแจกแจงแบบปัวซงผสมและการมาถึงของกลุ่มนักเรียนแทนที่จะเป็นนักเรียนแต่ละคนอาจ จำลองได้ด้วย การแจกแจงแบบปัวซงเชิงประกอบ

จำนวน แผ่นดินไหวขนาด 5 ริกเตอร์ต่อปีในประเทศหนึ่ง อาจไม่เป็นไปตามการกระจายแบบปัวซง หากแผ่นดินไหวขนาดใหญ่ครั้งหนึ่งเพิ่มความน่าจะเป็นของการเกิดแผ่นดินไหวตามมาที่มีขนาดใกล้เคียงกัน

ตัวอย่างที่รับประกันว่าจะมีเหตุการณ์อย่างน้อยหนึ่งเหตุการณ์เกิดขึ้นนั้น ไม่ได้มีการแจกแจงแบบปัวซง แต่สามารถจำลองได้โดยใช้การแจกแจงแบบปัวซงที่ตัดศูนย์ออก

การแจกแจงจำนวนนับที่จำนวนช่วงเวลาที่มีเหตุการณ์เป็นศูนย์สูงกว่าที่คาดการณ์โดยแบบจำลองปัวซง อาจสร้างแบบจำลองได้โดยใช้แบบจำลองที่มีค่าศูนย์มากเกินไป (zero-inflated model )

คุณสมบัติ

สถิติเชิงพรรณนา

ค่ามัธยฐาน

ขอบเขตสำหรับค่ามัธยฐาน (ν{\displaystyle \nu }) ของการกระจายตัวเป็นที่ทราบและมีความคมชัด : [ 23 ]λln2ν<λ+13.{\displaystyle \lambda -\ln 2\leq \nu <\lambda +{\frac {1}{3}}.}

ช่วงเวลาที่สูงกว่า

โมเมนต์ ที่ไม่เป็นศูนย์กลางที่สูงกว่าm ของการแจกแจงปัวซงคือพหุนาม Touchardในλ : เค=ฉัน=0เคλฉัน{เคฉัน},{\displaystyle m_{k}=\sum _{i=0}^{k}\lambda ^{i}{\begin{Bmatrix}k\\i\end{Bmatrix}},}โดยที่วงเล็บปีกกา { } หมายถึงจำนวนสเตอร์ลิงชนิดที่สอง [ 24 ] [ 1 ] : 6กล่าวอีกนัยหนึ่ง อี[X]=λ,อี[X(X1)]=λ2,อี[X(X1)(X2)]=λ3,{\displaystyle E[X]=\lambda ,\quad E[X(X-1)]=\lambda ^{2},\quad E[X(X-1)(X-2)]=\lambda ^{3},\cdots } เมื่อกำหนดค่าที่คาดหวังเป็นλ = 1 สูตรของโดบินสกิจะบ่งชี้ว่า โมเมนต์ลำดับที่ nเท่ากับจำนวน พาร์ ติชันของเซตที่มีขนาดn

ขอบเขตบนแบบง่ายคือ: [ 25 ]เค=อี[Xเค](เคบันทึก(เค/λ+1))เคλเคเอ็กซ์(เค22λ).{\displaystyle m_{k}=E[X^{k}]\leq \left({\frac {k}{\log(k/\lambda +1)}}\right)^{k}\leq \lambda ^{k}\exp \left({\frac {k^{2}}{2\lambda }}\right).}

ผลรวมของตัวแปรสุ่มที่มีการแจกแจงแบบปัวซง

ถ้าXฉัน~ปัวส์(λฉัน){\displaystyle X_{i}\sim \operatorname {Pois} (\lambda _{i})}สำหรับฉัน=1,,n{\displaystyle i=1,\dotsc ,n}ถ้า เป็นอิสระต่อกันก็แสดงว่า...ฉัน=1nXฉัน~ปัวส์(ฉัน=1nλฉัน).{\textstyle \sum _{i=1}^{n}X_{i}\sim \operatorname {Pois} \left(\sum _{i=1}^{n}\lambda _{i}\right).}[ 26 ] : 65บทกลับคือทฤษฎีบทของ Raikovซึ่งกล่าวว่า ถ้าผลรวมของตัวแปรสุ่มอิสระสองตัวมีการแจกแจงแบบปัวซง ตัวแปรสุ่มอิสระทั้งสองตัวนั้นก็จะมีการแจกแจงแบบปัวซงเช่นกัน [ 27 ] [ 28 ]

เอนโทรปีสูงสุด

It is a maximum-entropy distribution among the set of generalized binomial distributions Bn(λ){\displaystyle B_{n}(\lambda )} with mean λ{\displaystyle \lambda } and n{\displaystyle n\to \infty },[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 P=Pois(λ){\displaystyle P=\operatorname {Pois} (\lambda )} from P0=Pois(λ0){\displaystyle P_{0}=\operatorname {Pois} (\lambda _{0})} is given byDKL(PP0)=λ0λ+λlogλλ0.{\displaystyle \operatorname {D} _{\text{KL}}(P\parallel P_{0})=\lambda _{0}-\lambda +\lambda \log {\frac {\lambda }{\lambda _{0}}}.}
  • If λ1{\displaystyle \lambda \geq 1} is an integer, then YPois(λ){\displaystyle Y\sim \operatorname {Pois} (\lambda )} satisfies Pr(YE[Y])12{\displaystyle \Pr(Y\geq E[Y])\geq {\frac {1}{2}}} and Pr(YE[Y])12.{\displaystyle \Pr(Y\leq E[Y])\geq {\frac {1}{2}}.}[31]
  • Bounds for the tail probabilities of a Poisson random variable XPois(λ){\displaystyle X\sim \operatorname {Pois} (\lambda )} can be derived using a Chernoff bound argument.[32]:97-98P(Xx)(eλ)xeλxx, for x>λ,P(Xx)(eλ)xeλxx, for x<λ.{\displaystyle {\begin{aligned}P(X\geq x)&\leq {\frac {\left(e\lambda \right)^{x}e^{-\lambda }}{x^{x}}},&{\text{ for }}x>\lambda ,\\[1ex]P(X\leq x)&\leq {\frac {\left(e\lambda \right)^{x}e^{-\lambda }}{x^{x}}},&{\text{ for }}x<\lambda .\end{aligned}}}
  • The upper tail probability can be tightened (by a factor of at least two) as follows:[33]P(Xx)eDKL(QP)max(2,4πDKL(QP)), for x>λ,{\displaystyle P(X\geq x)\leq {\frac {e^{-\operatorname {D} _{\text{KL}}(Q\parallel P)}}{\max {(2,{\sqrt {4\pi \operatorname {D} _{\text{KL}}(Q\parallel P)}}})}},{\text{ for }}x>\lambda ,} where DKL(QP){\displaystyle \operatorname {D} _{\text{KL}}(Q\parallel P)} is the Kullback–Leibler divergence of Q=Pois(x){\displaystyle Q=\operatorname {Pois} (x)} from P=Pois(λ){\displaystyle P=\operatorname {Pois} (\lambda )}.
  • Inequalities that relate the cumulative distribution function of a Poisson random variable XPois(λ){\displaystyle X\sim \operatorname {Pois} (\lambda )} to the cumulative distribution function Φ{\displaystyle \Phi } of the standard normal distribution are as follows:[34]

Φ(sign(kλ)2DKL(QP))<P(Xk)<Φ(sign(k+1λ)2DKL(Q+P)), for k>0,{\displaystyle \Phi {\left(\operatorname {sign} (k-\lambda ){\sqrt {2\operatorname {D} _{\text{KL}}(Q_{-}\parallel P)}}\right)}<P(X\leq k)<\Phi {\left(\operatorname {sign} (k+1-\lambda ){\sqrt {2\operatorname {D} _{\text{KL}}(Q_{+}\parallel P)}}\right)},{\text{ for }}k>0,} where DKL(QP){\displaystyle \operatorname {D} _{\text{KL}}(Q_{-}\parallel P)} is the Kullback–Leibler divergence of Q=Pois(k){\displaystyle Q_{-}=\operatorname {Pois} (k)} from P=Pois(λ){\displaystyle P=\operatorname {Pois} (\lambda )} and DKL(Q+P){\displaystyle \operatorname {D} _{\text{KL}}(Q_{+}\parallel P)} is the Kullback–Leibler divergence of Q+=Pois(k+1){\displaystyle Q_{+}=\operatorname {Pois} (k+1)} from P{\displaystyle P}.

Poisson races

Let XPois(λ){\displaystyle X\sim \operatorname {Pois} (\lambda )} and YPois(μ){\displaystyle Y\sim \operatorname {Pois} (\mu )} be independent random variables, with λ<μ,{\displaystyle \lambda <\mu ,} then we have that e(μλ)2(λ+μ)2e(λ+μ)2λμe(λ+μ)4λμP(XY0)e(μλ)2{\displaystyle {\frac {e^{-({\sqrt {\mu }}-{\sqrt {\lambda }})^{2}}}{(\lambda +\mu )^{2}}}-{\frac {e^{-(\lambda +\mu )}}{2{\sqrt {\lambda \mu }}}}-{\frac {e^{-(\lambda +\mu )}}{4\lambda \mu }}\leq P(X-Y\geq 0)\leq e^{-({\sqrt {\mu }}-{\sqrt {\lambda }})^{2}}}

The upper bound is proved using a standard Chernoff bound.

The lower bound can be proved by noting that P(XY0X+Y=i){\displaystyle P(X-Y\geq 0\mid X+Y=i)} is the probability that Zi2,{\textstyle Z\geq {\frac {i}{2}},} where ZBin(i,λλ+μ),{\textstyle Z\sim \operatorname {Bin} \left(i,{\frac {\lambda }{\lambda +\mu }}\right),} which is bounded below by 1(i+1)2eiD(0.5λλ+μ),{\textstyle {\frac {1}{(i+1)^{2}}}e^{-iD\left(0.5\|{\frac {\lambda }{\lambda +\mu }}\right)},} where D{\displaystyle D} is relative entropy (See the entry on bounds on tails of binomial distributions for details). Further noting that X+YPois(λ+μ),{\displaystyle X+Y\sim \operatorname {Pois} (\lambda +\mu ),} 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]

As a Binomial distribution with infinitesimal time-steps

การแจกแจงปัวซงสามารถหาได้จากกรณีจำกัดของการแจกแจงทวินามเมื่อจำนวนการทดลองเข้าสู่ค่าอนันต์ และจำนวนความสำเร็จที่คาดหวัง ยังคงคงที่ — ดู หลักการของเหตุการณ์หายากด้านล่าง ดังนั้นจึงสามารถใช้เป็นการประมาณค่าของการแจกแจงทวินามได้หากnมีขนาดใหญ่พอและpมีขนาดเล็กพอ การแจกแจงปัวซงเป็นการประมาณค่าที่ดีของการแจกแจงทวินามหากnมีค่าอย่างน้อย 20 และpมีค่าน้อยกว่าหรือเท่ากับ 0.05 และเป็นการประมาณค่าที่ดีเยี่ยมหากn ≥ 100และnp 10 [ 36 ]ให้เอฟบี{\displaystyle F_{\mathrm {B} }}และเอฟพี{\displaystyle F_{\mathrm {P} }}โดยที่ เป็น ฟังก์ชันความหนาแน่นสะสมของ1การแจกแจงทวินามและการแจกแจงปัวซง ตามลำดับ จะได้ว่า:เอฟบี(เค;n,พี)  เอฟพี(เค;λ=nพี).{\displaystyle F_{\mathrm {B} }(k;n,p)\ \approx \ F_{\mathrm {P} }(k;\lambda =np).} การอนุมานอย่างหนึ่งของสิ่งนี้ใช้ฟังก์ชันสร้างความน่าจะเป็น [ 37 ] พิจารณาการทดลองแบบเบอร์นูลลี (การโยนเหรียญ) ซึ่งความน่าจะเป็นของความสำเร็จหนึ่งครั้ง (หรือจำนวนความสำเร็จที่คาดหวัง) คือλ1{\displaystyle \lambda \leq 1}ภายในช่วงเวลาที่กำหนด แบ่งช่วงเวลานั้นออกเป็นnส่วน และทำการทดลองในแต่ละช่วงย่อยด้วยความน่าจะเป็นλn{\displaystyle {\tfrac {\lambda }{n}}}ความน่าจะเป็นของ ความสำเร็จ kครั้งจาก การทดลอง nครั้งตลอดช่วงเวลาทั้งหมด จะกำหนดโดยการแจกแจงทวินาม พีเค(n)=(nเค)(λn)เค(1λn)nเค,{\displaystyle p_{k}^{(n)}={\binom {n}{k}}\left({\frac {\lambda }{n}}\right)^{\!k}\left(1{-}{\frac {\lambda }{n}}\right)^{\!n-k},} ซึ่งฟังก์ชันก่อกำเนิดคือ: พี(n)(x)=เค=0nพีเค(n)xเค=(1λn+λnx)n.{\displaystyle P^{(n)}(x)=\sum _{k=0}^{n}p_{k}^{(n)}x^{k}=\left(1-{\frac {\lambda }{n}}+{\frac {\lambda }{n}}x\right)^{n}.} เมื่อพิจารณาลิมิตเมื่อnเพิ่มขึ้นเป็นอนันต์ (โดยที่xคงที่) และใช้คำนิยามลิมิตผลคูณของฟังก์ชันเลขชี้กำลังจะได้เป็นฟังก์ชันก่อกำเนิดของการแจกแจงปัวซง: ลิมnพี(n)(x)=ลิมn(1+λ(x1)n)n=อีλ(x1)=เค=0อีλλเคเค!xเค.{\displaystyle \lim _{n\to \infty }P^{(n)}(x)=\lim _{n\to \infty }\left(1{+}{\tfrac {\lambda (x-1)}{n}}\right)^{n}=e^{\lambda (x-1)}=\sum _{k=0}^{\infty }e^{-\lambda }{\frac {\lambda ^{k}}{k!}}x^{k}.}

ทั่วไป

  • ถ้าX1~พีโอฉัน(λ1){\displaystyle X_{1}\sim \mathrm {Pois} (\lambda _{1})\,}และX2~พีโอฉัน(λ2){\displaystyle X_{2}\sim \mathrm {Pois} (\lambda _{2})\,}เป็นอิสระต่อกัน ดังนั้นความแตกต่างก็คือวาย=X1X2{\displaystyle Y=X_{1}-X_{2}}เป็นไปตามการแจกแจงแบบสเกลแลม (Skellam distribution )
  • ถ้าX1~พีโอฉัน(λ1){\displaystyle X_{1}\sim \mathrm {Pois} (\lambda _{1})\,}และX2~พีโอฉัน(λ2){\displaystyle X_{2}\sim \mathrm {Pois} (\lambda _{2})\,}หากเป็นอิสระต่อกัน การกระจายตัวของX1{\displaystyle X_{1}}โดยมีเงื่อนไขว่าX1+X2{\displaystyle X_{1}+X_{2}}เป็นการแจกแจงแบบทวินาม
    โดยเฉพาะอย่างยิ่ง ถ้าX1+X2=เค,{\displaystyle X_{1}+X_{2}=k,}แล้วX1|X1+X2=เค~บีฉันnโอ(เค,λ1/(λ1+λ2)).{\displaystyle X_{1}|X_{1}+X_{2}=k\sim \mathrm {Binom} (k,\lambda _{1}/(\lambda _{1}+\lambda _{2})).}
    โดยทั่วไปแล้ว ถ้าX , X , ..., X เป็นตัวแปรสุ่มปัวซงอิสระที่มีพารามิเตอร์λ , λ , ..., λ แล้ว
    ที่ให้ไว้เจ=1nXเจ=เค,{\displaystyle \sum _{j=1}^{n}X_{j}=k,}ดังนั้นจึงสรุปได้ว่าXฉัน|เจ=1nXเจ=เค~บีฉันnโอ(เค,λฉันเจ=1nλเจ).{\displaystyle X_{i}{\Big |}\sum _{j=1}^{n}X_{j}=k\sim \mathrm {Binom} \left(k,{\frac {\lambda _{i}}{\sum _{j=1}^{n}\lambda _{j}}}\right).}ในความเป็นจริง,{Xฉัน}~เอ็มคุณทีฉันnโอ(เค,{λฉันเจ=1nλเจ}).{\displaystyle \{X_{i}\}\sim \mathrm {Multinom} \left(k,\left\{{\frac {\lambda _{i}}{\sum _{j=1}^{n}\lambda _{j}}}\right\}\right).}
  • ถ้าX~พีโอฉัน(λ){\displaystyle X\sim \mathrm {Pois} (\lambda )\,}และการกระจายตัวของวาย{\displaystyle Y}โดยมีเงื่อนไขว่าX = kเป็นการแจกแจงแบบทวินามวาย(X=เค)~บีฉันnโอ(เค,พี),{\displaystyle Y\mid (X=k)\sim \mathrm {Binom} (k,p),}ดังนั้นการกระจายตัวของ Y จึงเป็นไปตามการกระจายแบบปัวซงวาย~พีโอฉัน(λพี).{\displaystyle Y\sim \mathrm {Pois} (\lambda \cdot p).}ในความเป็นจริง ถ้าหากพิจารณาจากเงื่อนไขบางประการ{X=เค},{\displaystyle \{X=k\},}{วายฉัน}{\displaystyle \{Y_{i}\}}เป็นไปตามการแจกแจงแบบพหุนาม{วายฉัน}(X=เค)~เอ็มคุณทีฉันnโอ(เค,พีฉัน),{\displaystyle \{Y_{i}\}\mid (X=k)\sim \mathrm {Multinom} \left(k,p_{i}\right),}จากนั้นแต่ละวายฉัน{\displaystyle Y_{i}}เป็นไปตามการแจกแจงปัวซงแบบอิสระวายฉัน~พีโอฉัน(λพีฉัน),ρ(วายฉัน,วายเจ)=0.{\displaystyle Y_{i}\sim \mathrm {Pois} (\lambda \cdot p_{i}),\rho (Y_{i},Y_{j})=0.}
  • การแจกแจงปัวซงเป็นกรณีพิเศษของการแจกแจงปัวซงแบบผสมแบบไม่ต่อเนื่อง (หรือการแจกแจงปัวซงแบบกระตุก) ที่มีพารามิเตอร์เพียงตัวเดียว[ 38 ] [ 39 ]การแจกแจงปัวซงแบบผสมแบบไม่ต่อเนื่องสามารถอนุมานได้จากการแจกแจงจำกัดของการแจกแจงพหุนามแบบเอกตัวแปร นอกจากนี้ยังเป็นกรณีพิเศษของการแจกแจงปัวซงแบบผสมอีก ด้วย
  • สำหรับค่าλ ที่มีขนาดใหญ่เพียงพอ (เช่นλ > 1000 ) การแจกแจงแบบปกติจะมีค่าเฉลี่ยλและความแปรปรวนλ (ส่วนเบี่ยงเบนมาตรฐาน)λ{\displaystyle {\sqrt {\lambda }}}) เป็นการประมาณที่ดีเยี่ยมของค่าการแจกแจงปัวซงหากλมากกว่าประมาณ10 การแจกแจงปกติจะเป็นการประมาณที่ดีหาก มี การแก้ไขความต่อเนื่อง ที่เหมาะสม กล่าวคือ หากP( Xx )โดยที่xเป็นจำนวนเต็มที่ไม่เป็นลบ ถูกแทนที่ด้วยP( Xx + 0.5 ) เอฟพีโอฉันโอn(x;λ)เอฟnโอเอ(x;μ=λ,σ2=λ){\displaystyle F_{\mathrm {Poisson} }(x;\lambda )\approx F_{\mathrm {normal} }(x;\mu =\lambda ,\sigma ^{2}=\lambda )}
  • การแปลงค่าความแปรปรวนให้คงที่ : ถ้าX~พีโอฉัน(λ),{\displaystyle X\sim \mathrm {Pois} (\lambda ),}จากนั้น[ 8 ] : 168วาย=2Xเอ็น(2λ;1),{\displaystyle Y=2{\sqrt {X}}\approx {\mathcal {N}}(2{\sqrt {\lambda }};1),}และ[ 40 ] : 196วาย=Xเอ็น(λ;1/4).{\displaystyle Y={\sqrt {X}}\approx {\mathcal {N}}({\sqrt {\lambda }};1/4).}ภายใต้การแปลงนี้ การลู่เข้าสู่ภาวะปกติ (ดังที่)λ{\displaystyle \lambda }การเพิ่มขึ้น) นั้นเร็วกว่าตัวแปรที่ไม่ได้รับการแปลงมากการแปลงที่ทำให้ความแปรปรวนคงที่ที่ซับซ้อนกว่าเล็กน้อยก็มีอยู่เช่นกัน[ 8 ] : 168หนึ่งในนั้นคือ การ แปลงAnscombe [ 41 ]ดูการแปลงข้อมูล (สถิติ)สำหรับการใช้งานการแปลงทั่วไปเพิ่มเติม
  • ถ้าสำหรับทุกt > 0จำนวนการมาถึงในช่วงเวลา[ 0, t ]เป็นไปตามการแจกแจงปัวซงที่มีค่าเฉลี่ยλtแล้วลำดับของเวลาระหว่างการมาถึงจะเป็นตัวแปรสุ่มเอกซ์โพเนนเชียลที่เป็นอิสระและมีการแจกแจงเหมือนกัน โดยมีค่าเฉลี่ย1/ λ [ 42 ] : 317–319 
  • ฟังก์ชันการกระจายสะสมของการกระจายแบบปัวซงและไคกำลังสองมีความสัมพันธ์กันในลักษณะต่อไปนี้: [ 8 ] : 167เอฟปัวซง(เค;λ)=1เอฟχ2(2λ;2(เค+1)) จำนวนเต็ม เค,{\displaystyle F_{\text{Poisson}}(k;\lambda )=1-F_{\chi ^{2}}(2\lambda ;2(k+1))\quad \quad {\text{ integer }}k,}และ[ 8 ] : 158พี(X=เค)=เอฟχ2(2λ;2(เค+1))เอฟχ2(2λ;2เค).{\displaystyle P(X=k)=F_{\chi ^{2}}(2\lambda ;2(k+1))-F_{\chi ^{2}}(2\lambda ;2k).}

การประมาณค่าปัวซง

สมมติX1~ปัวส์(λ1),X2~ปัวส์(λ2),,Xn~ปัวส์(λn){\displaystyle X_{1}\sim \operatorname {Pois} (\lambda _{1}),X_{2}\sim \operatorname {Pois} (\lambda _{2}),\dots ,X_{n}\sim \operatorname {Pois} (\lambda _{n})}ที่ไหนλ1+λ2++λn=1,{\displaystyle \lambda _{1}+\lambda _{2}+\dots +\lambda _{n}=1,}จากนั้น[ 43 ](X1,X2,,Xn){\displaystyle (X_{1},X_{2},\dots ,X_{n})}มีการแจกแจงแบบพหุนาม(X1,X2,,Xn)~มัลท์(เอ็น,λ1,λ2,,λn){\displaystyle (X_{1},X_{2},\dots ,X_{n})\sim \operatorname {Mult} (N,\lambda _{1},\lambda _{2},\dots ,\lambda _{n})}โดยมีเงื่อนไขว่าเอ็น=X1+X2+Xn.{\displaystyle N=X_{1}+X_{2}+\dots X_{n}.}

หมายความว่า[ 32 ] : 101-102ในบรรดาสิ่งอื่นๆ สำหรับฟังก์ชันที่ไม่เป็นลบใดๆเอฟ(x1,x2,,xn),{\displaystyle f(x_{1},x_{2},\dots ,x_{n}),} ถ้า(วาย1,วาย2,,วายn)~มัลท์(,พี){\displaystyle (Y_{1},Y_{2},\dots ,Y_{n})\sim \operatorname {Mult} (m,\mathbf {p} )}ถ้ามีการแจกแจงแบบพหุนามแล้ว อี[เอฟ(วาย1,วาย2,,วายn)]อีอี[เอฟ(X1,X2,,Xn)]{\displaystyle \operatorname {E} [f(Y_{1},Y_{2},\dots ,Y_{n})]\leq e{\sqrt {m}}\operatorname {E} [f(X_{1},X_{2},\dots ,X_{n})]} ที่ไหน(X1,X2,,Xn)~ปัวส์(พี).{\displaystyle (X_{1},X_{2},\dots ,X_{n})\sim \operatorname {Pois} (\mathbf {p} ).}

ปัจจัยของอี{\displaystyle e{\sqrt {m}}}สามารถแทนที่ด้วย 2 ได้หากเอฟ{\displaystyle f}นอกจากนี้ยังสันนิษฐานว่ามีแนวโน้มเพิ่มขึ้นหรือลดลงอย่างต่อเนื่อง

การแจกแจงปัวซงแบบสองตัวแปร

การแจกแจงนี้ได้รับการขยายไปยังกรณีสองตัวแปร[ 44 ]ฟังก์ชันก่อกำเนิดสำหรับการแจกแจงนี้คือ จี(คุณ,วี)=เอ็กซ์[(θ1θ12)(คุณ1)+(θ2θ12)(วี1)+θ12(คุณวี1)]{\displaystyle g(u,v)=\exp[(\theta _{1}-\theta _{12})(u-1)+(\theta _{2}-\theta _{12})(v-1)+\theta _{12}(uv-1)]} กับ θ1,θ2>θ12>0{\displaystyle \theta _{1},\theta _{2}>\theta _{12}>0}

การแจกแจงแบบมาร์จินัลคือการแจกแจงปัวซง ( θ )และการแจกแจงปัวซง ( θ )และค่าสัมประสิทธิ์สหสัมพันธ์ถูกจำกัดอยู่ในช่วง 0ρนาที{θ1θ2,θ2θ1}{\displaystyle 0\leq \rho \leq \min \left\{{\sqrt {\frac {\theta _{1}}{\theta _{2}}}},{\sqrt {\frac {\theta _{2}}{\theta _{1}}}}\right\}}

วิธีง่ายๆ ในการสร้างการแจกแจงปัวซงแบบสองตัวแปรX1,X2{\displaystyle X_{1},X_{2}}คือการนำการแจกแจงปัวซงอิสระสามแบบมาใช้วาย1,วาย2,วาย3{\displaystyle Y_{1},Y_{2},Y_{3}}ด้วยวิธีการλ1,λ2,λ3{\displaystyle \lambda _{1},\lambda _{2},\lambda _{3}}แล้วจึงตั้งค่าX1=วาย1+วาย3,X2=วาย2+วาย3.{\displaystyle X_{1}=Y_{1}+Y_{3},X_{2}=Y_{2}+Y_{3}.}ฟังก์ชันความน่าจะเป็นของการแจกแจงปัวซงแบบสองตัวแปรคือ ปร.(X1=เค1,X2=เค2)=เอ็กซ์(λ1λ2λ3)λ1เค1เค1!λ2เค2เค2!เค=0นาที(เค1,เค2)(เค1เค)(เค2เค)เค!(λ3λ1λ2)เค{\displaystyle \Pr(X_{1}=k_{1},X_{2}=k_{2})=\exp \left(-\lambda _{1}-\lambda _{2}-\lambda _{3}\right){\frac {\lambda _{1}^{k_{1}}}{k_{1}!}}{\frac {\lambda _{2}^{k_{2}}}{k_{2}!}}\sum _{k=0}^{\min(k_{1},k_{2})}{\binom {k_{1}}{k}}{\binom {k_{2}}{k}}k!\left({\frac {\lambda _{3}}{\lambda _{1}\lambda _{2}}}\right)^{k}}

การแจกแจงปัวซงแบบอิสระ

การแจกแจงปัวซงแบบอิสระ[ 45 ]ที่มีขนาดกระโดดα{\displaystyle \alpha }และให้คะแนนλ{\displaystyle \lambda }เกิดขึ้นใน ทฤษฎี ความน่าจะเป็นอิสระในฐานะลิมิตของการสังเคราะห์อิสระ ซ้ำๆ((1λเอ็น)δ0+λเอ็นδα)เอ็น{\displaystyle \left(\left(1-{\frac {\lambda }{N}}\right)\delta _{0}+{\frac {\lambda }{N}}\delta _{\alpha }\right)^{\boxplus N}} เมื่อN

กล่าวอีกนัยหนึ่งคือ ให้Xเอ็น{\displaystyle X_{N}}เป็นตัวแปรสุ่มเพื่อให้Xเอ็น{\displaystyle X_{N}}มีค่าα{\displaystyle \alpha }ด้วยความน่าจะเป็นλเอ็น{\textstyle {\frac {\lambda }{N}}}และค่า 0 ด้วยความน่าจะเป็นที่เหลืออยู่ สมมติด้วยว่าครอบครัวX1,X2,{\displaystyle X_{1},X_{2},\ldots }เป็นอิสระอย่างอิสระจากนั้นขีดจำกัดคือเอ็น{\displaystyle N\to \infty }ของกฎหมายของX1++Xเอ็น{\displaystyle X_{1}+\cdots +X_{N}}กำหนดโดยกฎปัวซงอิสระที่มีพารามิเตอร์λ,α.{\displaystyle \lambda ,\alpha .}

คำจำกัดความนี้คล้ายคลึงกับวิธีการหนึ่งในการได้มาซึ่งการแจกแจงปัวซงแบบคลาสสิกจากกระบวนการปัวซง (แบบคลาสสิก)

มาตรการที่เกี่ยวข้องกับกฎปัวซงอิสระจะได้รับจาก[ 46 ]μ={(1λ)δ0+ν,ถ้า 0λ1ν,ถ้า λ>1,{\displaystyle \mu ={\begin{cases}(1-\lambda )\delta _{0}+\nu ,&{\text{if }}0\leq \lambda \leq 1\\\nu ,&{\text{if }}\lambda >1,\end{cases}}} ที่ไหน ν=12παที4λα2(ทีα(1+λ))2ที{\displaystyle \nu ={\frac {1}{2\pi \alpha t}}{\sqrt {4\lambda \alpha ^{2}-(t-\alpha (1+\lambda ))^{2}}}\,dt} และได้รับการสนับสนุน[α(1λ)2,α(1+λ)2].{\displaystyle [\alpha (1-{\sqrt {\lambda }})^{2},\alpha (1+{\sqrt {\lambda }})^{2}].}

กฎนี้ยังปรากฏในทฤษฎีเมทริกซ์สุ่ม ในรูปของ กฎมาร์เชนโก-ปาสตูร์ค่าคุมูลันต์อิสระของกฎนี้เท่ากับκn=λαn.{\displaystyle \kappa _{n}=\lambda \alpha ^{n}.}

การแปลงบางส่วนของกฎนี้

เราให้ค่าของการแปลงที่สำคัญบางอย่างของกฎปัวซงอิสระ การคำนวณสามารถพบได้ในหนังสือLectures on the Combinatorics of Free Probabilityโดย A. Nica และ R. Speicher [ 47 ]

การแปลง R ของกฎปัวซงอิสระกำหนดโดย อาร์(z)=λα1αz.{\displaystyle R(z)={\frac {\lambda \alpha }{1-\alpha z}}.}

การแปลงโคชี (ซึ่งเป็นการแปลงกลับของการแปลงสติลต์เจส ) กำหนดโดย จี(z)=z+αλα(zα(1+λ))24λα22αz{\displaystyle G(z)={\frac {z+\alpha -\lambda \alpha -{\sqrt {(z-\alpha (1+\lambda ))^{2}-4\lambda \alpha ^{2}}}}{2\alpha z}}}

การแปลง S กำหนดโดย เอส(z)=1z+λ{\displaystyle S(z)={\frac {1}{z+\lambda }}} ในกรณีที่α=1.{\displaystyle \alpha =1.}

การอนุมานทางสถิติ

การประมาณค่าพารามิเตอร์

กำหนดให้มีค่าที่วัดได้จำนวนnค่า เป็นตัวอย่างเคฉัน{0,1,},{\displaystyle k_{i}\in \{0,1,\dots \},}สำหรับi = 1, ..., nเราต้องการประมาณค่าของพารามิเตอร์λของประชากรปัวซงที่สุ่มตัวอย่างมา การประมาณค่า ความน่าจะเป็นสูงสุดคือ[ 48 ]

λ^เอ็มแอลอี=1nฉัน=1nเคฉัน .{\displaystyle {\widehat {\lambda }}_{\mathrm {MLE} }={\frac {1}{n}}\sum _{i=1}^{n}k_{i}\ .}

เนื่องจากแต่ละการสังเกตมีค่าเฉลี่ยλดังนั้นค่าเฉลี่ยของตัวอย่างจึงมีค่าเฉลี่ย λ เช่นกัน ด้วยเหตุนี้ การประมาณค่าความน่าจะเป็นสูงสุดจึงเป็นตัวประมาณค่าที่ไม่เอนเอียงของλนอกจากนี้ยังเป็นตัวประมาณค่าที่มีประสิทธิภาพ เนื่องจากความแปรปรวนของมันบรรลุขอบเขตล่างของ Cramér–Rao (CRLB) [ 49 ]ดังนั้นจึงเป็นตัวประมาณค่าที่ไม่เอนเอียงที่มีความแปรปรวนต่ำสุด นอกจากนี้ยังสามารถพิสูจน์ได้ว่าผลรวม (และด้วยเหตุนี้ค่าเฉลี่ย ของตัวอย่าง เนื่องจากเป็นฟังก์ชันหนึ่งต่อหนึ่งของผลรวม) เป็นสถิติที่สมบูรณ์และเพียงพอสำหรับλ

เพื่อพิสูจน์ความเพียงพอ เราอาจใช้ทฤษฎีบทการแยกตัวประกอบพิจารณาการแบ่งฟังก์ชันความน่าจะเป็นมวลของการแจกแจงปัวซงร่วมสำหรับตัวอย่างออกเป็นสองส่วน: ส่วนหนึ่งที่ขึ้นอยู่กับตัวอย่างเพียงอย่างเดียวx{\displaystyle \mathbf {x} }, เรียกว่าชม.(x){\displaystyle h(\mathbf {x} )}และเป็นสิ่งที่ขึ้นอยู่กับพารามิเตอร์λ{\displaystyle \lambda }และตัวอย่างx{\displaystyle \mathbf {x} }เฉพาะผ่านฟังก์ชันเท่านั้นที(x).{\displaystyle T(\mathbf {x} ).}แล้วที(x){\displaystyle T(\mathbf {x} )}เป็นสถิติที่เพียงพอสำหรับλ.{\displaystyle \lambda .}

พี(x)=ฉัน=1nλxฉันอีλxฉัน!=1ฉัน=1nxฉัน!×λฉัน=1nxฉันอีnλ{\displaystyle P(\mathbf {x} )=\prod _{i=1}^{n}{\frac {\lambda ^{x_{i}}e^{-\lambda }}{x_{i}!}}={\frac {1}{\prod _{i=1}^{n}x_{i}!}}\times \lambda ^{\sum _{i=1}^{n}x_{i}}e^{-n\lambda }}

เทอมแรกชม.(x){\displaystyle h(\mathbf {x} )}ขึ้นอยู่กับเท่านั้นx{\displaystyle \mathbf {x} }เทอมที่สองจี(ที(x)|λ){\displaystyle g(T(\mathbf {x} )|\lambda )}ขึ้นอยู่กับตัวอย่างเท่านั้นผ่านทางที(x)=ฉัน=1nxฉัน.{\textstyle T(\mathbf {x} )=\sum _{i=1}^{n}x_{i}.}ดังนั้น,ที(x){\displaystyle T(\mathbf {x} )}เพียงพอแล้ว

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

(λ)=lnฉัน=1nเอฟ(เคฉันλ)=ฉัน=1nln(อีλλเคฉันเคฉัน!)=nλ+(ฉัน=1nเคฉัน)ln(λ)ฉัน=1nln(เคฉัน!).{\displaystyle {\begin{aligned}\ell (\lambda )&=\ln \prod _{i=1}^{n}f(k_{i}\mid \lambda )\\&=\sum _{i=1}^{n}\ln \!\left({\frac {e^{-\lambda }\lambda ^{k_{i}}}{k_{i}!}}\right)\\&=-n\lambda +\left(\sum _{i=1}^{n}k_{i}\right)\ln(\lambda )-\sum _{i=1}^{n}\ln(k_{i}!).\end{aligned}}}

เราหาอนุพันธ์ของ{\displaystyle \ell }โดยสัมพันธ์กับλและเปรียบเทียบกับศูนย์:

λ(λ)=0n+(ฉัน=1nเคฉัน)1λ=0.{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} \lambda }}\ell (\lambda )=0\iff -n+\left(\sum _{i=1}^{n}k_{i}\right){\frac {1}{\lambda }}=0.\!}

การแก้สมการหาค่าλจะได้จุดนิ่ง

λ=ฉัน=1nเคฉันn{\displaystyle \lambda ={\frac {\sum _{i=1}^{n}k_{i}}{n}}}

ดังนั้นλคือค่าเฉลี่ยของ ค่า k การหาเครื่องหมายของอนุพันธ์อันดับสองของLที่จุดนิ่งจะช่วยกำหนดว่าλเป็น ค่าสุดขั้วประเภทใด

2λ2=λ2ฉัน=1nเคฉัน{\displaystyle {\frac {\partial ^{2}\ell }{\partial \lambda ^{2}}}=-\lambda ^{-2}\sum _{i=1}^{n}k_{i}}

การประเมินอนุพันธ์อันดับสองณ จุดวิเคราะห์จะได้:

2λ2=n2ฉัน=1nเคฉัน{\displaystyle {\frac {\partial ^{2}\ell }{\partial \lambda ^{2}}}=-{\frac {n^{2}}{\sum _{i=1}^{n}k_{i}}}}

ซึ่งเป็นค่าลบของnคูณด้วยส่วนกลับของค่าเฉลี่ยของk นิพจน์นี้จะมีค่าเป็นลบเมื่อค่าเฉลี่ยเป็นบวก ถ้าเงื่อนไขนี้เป็นจริง จุดนิ่งจะทำให้ฟังก์ชันความน่าจะเป็นมีค่าสูงสุด

เพื่อให้ครบถ้วนสมบูรณ์ตระกูลของการแจกแจงจะเรียกว่าสมบูรณ์ก็ต่อเมื่ออี(จี(ที))=0{\displaystyle E(g(T))=0}หมายความว่าพีλ(จี(ที)=0)=1{\displaystyle P_{\lambda }(g(T)=0)=1}สำหรับทุกคนλ.{\displaystyle \lambda .}หากบุคคลนั้นXฉัน{\displaystyle X_{i}}เป็น iidพีโอ(λ),{\displaystyle \mathrm {Po} (\lambda ),}แล้วที(x)=ฉัน=1nXฉัน~พีโอ(nλ).{\textstyle T(\mathbf {x} )=\sum _{i=1}^{n}X_{i}\sim \mathrm {Po} (n\lambda ).}เมื่อทราบการกระจายตัวที่เราต้องการตรวจสอบแล้ว ก็จะเห็นได้ง่ายว่าสถิตินั้นสมบูรณ์

อี(จี(ที))=ที=0จี(ที)(nλ)ทีอีnλที!=0{\displaystyle E(g(T))=\sum _{t=0}^{\infty }g(t){\frac {(n\lambda )^{t}e^{-n\lambda }}{t!}}=0}

เพื่อให้ความเสมอภาคนี้คงอยู่จี(ที){\displaystyle g(t)}ต้องเป็น 0 ซึ่งเป็นผลมาจากข้อเท็จจริงที่ว่าไม่มีเงื่อนไขอื่นใดที่จะเป็น 0 สำหรับทุกกรณีที{\displaystyle t}ในผลรวมและสำหรับค่าที่เป็นไปได้ทั้งหมดของλ.{\displaystyle \lambda .}เพราะฉะนั้น,อี(จี(ที))=0{\displaystyle E(g(T))=0}สำหรับทุกคนλ{\displaystyle \lambda }หมายความว่าพีλ(จี(ที)=0)=1,{\displaystyle P_{\lambda }(g(T)=0)=1,}และสถิติดังกล่าวได้รับการพิสูจน์แล้วว่าสมบูรณ์

ช่วงความเชื่อมั่น

ช่วงความเชื่อมั่นสำหรับค่าเฉลี่ยของการแจกแจงปัวซงสามารถแสดงได้โดยใช้ความสัมพันธ์ระหว่างฟังก์ชันการแจกแจงสะสมของการแจกแจงปัวซงและการแจกแจงไคกำลังสองการแจกแจงไคกำลังสองมีความสัมพันธ์อย่างใกล้ชิดกับการแจกแจงแกมมาและนี่นำไปสู่การแสดงออกอีกแบบหนึ่ง เมื่อกำหนดค่าสังเกตkจากการแจกแจงปัวซงที่มีค่าเฉลี่ยμช่วงความเชื่อมั่นสำหรับμที่ระดับความเชื่อมั่น1 – αคือ

12χ2(α/2;2เค)μ12χ2(1α/2;2เค+2),{\displaystyle {\tfrac {1}{2}}\chi ^{2}(\alpha /2;2k)\leq \mu \leq {\tfrac {1}{2}}\chi ^{2}(1-\alpha /2;2k+2),}

หรือเทียบเท่า

เอฟ1(α/2;เค,1)μเอฟ1(1α/2;เค+1,1),{\displaystyle F^{-1}(\alpha /2;k,1)\leq \mu \leq F^{-1}(1-\alpha /2;k+1,1),}

ที่ไหนχ2(พี;n){\displaystyle \chi ^{2}(p;n)}คือฟังก์ชันควอนไทล์ (ซึ่งสอดคล้องกับพื้นที่หางด้านล่างp ) ของการแจกแจงไคกำลังสองที่มีnองศาอิสระ และเอฟ1(พี;n,1){\displaystyle F^{-1}(p;n,1)}คือฟังก์ชันควอนไทล์ของการแจกแจงแกมมาที่มีพารามิเตอร์รูปร่าง n และพารามิเตอร์มาตราส่วน 1 [ 8 ] : 176-178 [ 50 ]ช่วงนี้เป็น ' ช่วงที่แน่นอน ' ในแง่ที่ว่าความน่าจะเป็นของการครอบคลุมจะไม่น้อยกว่าค่าที่กำหนด1 α

เมื่อไม่มีควอนไทล์ของการแจกแจงแกมมา ได้มีการเสนอการประมาณค่าที่แม่นยำสำหรับช่วงเวลาที่แน่นอนนี้ (โดยอิงจากการแปลงวิลสัน-ฮิลเฟอร์ตี้ ) ดังนี้: [ 51 ]เค(119เคzα/23เค)3μ(เค+1)(119(เค+1)+zα/23เค+1)3,{\displaystyle k\left(1-{\frac {1}{9k}}-{\frac {z_{\alpha /2}}{3{\sqrt {k}}}}\right)^{3}\leq \mu \leq (k+1)\left(1-{\frac {1}{9(k+1)}}+{\frac {z_{\alpha /2}}{3{\sqrt {k+1}}}}\right)^{3},} ที่ไหนzα/2{\displaystyle z_{\alpha /2}}หมายถึงค่าเบี่ยงเบนมาตรฐานปกติที่มีพื้นที่หางด้านบนα / 2

สำหรับการประยุกต์ใช้สูตรเหล่านี้ในบริบทเดียวกันกับข้างต้น (โดยกำหนดให้มี ค่าที่วัดได้ nค่าk แต่ละค่าสุ่มมาจากการแจกแจงแบบปัวซงที่มีค่าเฉลี่ยλ ) จะต้องตั้งค่าดังนี้

เค=ฉัน=1nเคฉัน,{\displaystyle k=\sum _{i=1}^{n}k_{i},}คำนวณ ช่วงสำหรับμ = จากนั้นหาช่วงสำหรับλ

การอนุมานแบบเบย์เซียน

ในการอนุมานแบบเบย์เซียนไพรเออร์คู่ควบสำหรับพารามิเตอร์อัตราλของการแจกแจงปัวซงคือการแจกแจงแกมมา[ 52 ]ให้

λ~จีเอเอ(α,เบต้า){\displaystyle \lambda \sim \mathrm {Gamma} (\alpha ,\beta )}

แสดงว่าλมีการกระจายตามความหนาแน่น แกมมา gซึ่งกำหนดโดยพารามิเตอร์รูปร่างαและพารามิเตอร์มาตราส่วน ผกผัน β :

จี(λα,เบต้า)=เบต้าαΓ(α)λα1อีเบต้าλ สำหรับ λ>0.{\displaystyle g(\lambda \mid \alpha ,\beta )={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\;\lambda ^{\alpha -1}\;e^{-\beta \,\lambda }\qquad {\text{ for }}\lambda >0\,\!.}

จากนั้น เมื่อกำหนดตัวอย่าง ค่าที่วัดได้ nค่าk เช่นเดียว กับก่อนหน้านี้และค่าความน่าจะเป็นก่อนหน้าเป็นGamma( α , β )ค่าความน่าจะเป็นภายหลังจะเป็นดังนี้

λ~จีเอเอ(α+ฉัน=1nเคฉัน,เบต้า+n).{\displaystyle \lambda \sim \mathrm {Gamma} {\left(\alpha +\sum _{i=1}^{n}k_{i},\beta +n\right)}.}

โปรดทราบว่าค่าเฉลี่ยภายหลังเป็นแบบเชิงเส้นและกำหนดโดย อี[λเค1,,เคn]=α+ฉัน=1nเคฉันเบต้า+n.{\displaystyle E[\lambda \mid k_{1},\ldots ,k_{n}]={\frac {\alpha +\sum _{i=1}^{n}k_{i}}{\beta +n}}.} สามารถแสดงได้ว่าการแจกแจงแกมมาเป็นไพรเออร์เดียวที่ทำให้ค่าเฉลี่ยแบบมีเงื่อนไขเป็นเชิงเส้น ยิ่งไปกว่านั้น ยังมีผลลัพธ์ตรงกันข้ามที่ระบุว่า ถ้าค่าเฉลี่ยแบบมีเงื่อนไขใกล้เคียงกับฟังก์ชันเชิงเส้นในแอล2{\displaystyle L_{2}}ระยะทางที่มากกว่าการแจกแจงก่อนหน้าของλจะต้องใกล้เคียงกับการแจกแจงแกมมาใน ระยะ ทางของ Levy [ 53 ]

ค่าเฉลี่ยภายหลังE[ λ ]เข้าใกล้ค่าประมาณความน่าจะเป็นสูงสุดλ^เอ็มแอลอี{\displaystyle {\widehat {\lambda }}_{\mathrm {MLE} }}ในขีดจำกัดเมื่อα0,เบต้า0,{\displaystyle \alpha \to 0,\beta \to 0,}ซึ่งเป็น ผลสืบเนื่องโดยตรงจากนิพจน์ทั่วไปของค่าเฉลี่ยของการแจกแจงแกมมา

การแจกแจงการทำนายภายหลัง สำหรับการ สังเกตเพิ่มเติมเพียงครั้งเดียวคือการแจกแจงทวินามเชิงลบ[ 54 ] : 53บางครั้งเรียกว่าการแจกแจงแกมมา-ปัวซง

การประมาณค่าเฉลี่ยปัวซงหลายค่าพร้อมกัน

สมมติX1,X2,,Xพี{\displaystyle X_{1},X_{2},\dots ,X_{p}}คือชุดของตัวแปรสุ่มอิสระจากชุดของพี{\displaystyle p}การแจกแจงปัวซง โดยแต่ละแบบมีพารามิเตอร์λฉัน,{\displaystyle \lambda _{i},}ฉัน=1,,พี,{\displaystyle i=1,\dots ,p,}และเราต้องการประมาณค่าพารามิเตอร์เหล่านี้ จากนั้น Clevenson และ Zidek แสดงให้เห็นว่าภายใต้การสูญเสียข้อผิดพลาดกำลังสองแบบนอร์มาไลซ์แอล(λ,λ^)=ฉัน=1พีλฉัน1(λ^ฉันλฉัน)2,{\textstyle L(\lambda ,{\hat {\lambda }})=\sum _{i=1}^{p}\lambda _{i}^{-1}({\hat {\lambda }}_{i}-\lambda _{i})^{2},}เมื่อไรพี>1,{\displaystyle p>1,}จากนั้น เช่นเดียวกับในตัวอย่างของสไตน์สำหรับค่าเฉลี่ยปกติ ตัวประมาณค่า MLEλ^ฉัน=Xฉัน{\displaystyle {\hat {\lambda }}_{i}=X_{i}}ไม่สามารถยอมรับได้[ 55 ]

ในกรณีนี้ จะมี การกำหนดตระกูลของตัวประมาณค่าแบบมินิแม็กซ์สำหรับทุก ๆ0<ซี2(พี1){\displaystyle 0<c\leq 2(p-1)}และ(พี2+พี1){\displaystyle b\geq (p-2+p^{-1})}เช่น[ 56 ]λ^ฉัน=(1ซี+ฉัน=1พีXฉัน)Xฉัน,ฉัน=1,,พี.{\displaystyle {\hat {\lambda }}_{i}=\left(1-{\frac {c}{b+\sum _{i=1}^{p}X_{i}}}\right)X_{i},\qquad i=1,\dots ,p.}

การเกิดขึ้นและการประยุกต์ใช้

การประยุกต์ใช้การแจกแจงปัวซงกับข้อมูลการนับ (จำนวนเหตุการณ์) บางประการ: [ 57 ]

ตัวอย่างเพิ่มเติมของการนับเหตุการณ์ที่สามารถจำลองได้ด้วยกระบวนการปัวซง ได้แก่:

ในทฤษฎีจำนวนเชิงความน่าจะเป็นแกลลาเกอร์แสดงให้เห็นใน ปี 1976 ว่า หากสมมติฐานr-tuple ของจำนวนเฉพาะ ที่ยังไม่ได้ รับ การพิสูจน์บางเวอร์ชันเป็นจริง [ 70 ] แล้วจำนวนเฉพาะในช่วงเวลาสั้นๆ จะเป็นไปตามการแจกแจงปัวซง[ 71 ]

กฎของเหตุการณ์หายาก

การเปรียบเทียบการแจกแจงปัวซง (เส้นสีดำ) และการแจกแจงทวินามโดยที่n = 10 (วงกลมสีแดง), n = 20 (วงกลมสีน้ำเงิน), n = 1000 (วงกลมสีเขียว) การแจกแจงทั้งหมดมีค่าเฉลี่ยเท่ากับ 5 แกนแนวนอนแสดงจำนวนเหตุการณ์kยิ่งnมีค่ามากขึ้น การแจกแจงปัวซงก็ยิ่งเป็นการประมาณค่าที่ดีขึ้นสำหรับการแจกแจงทวินามที่มีค่าเฉลี่ยเท่ากัน 

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

ให้ แทนจำนวนเหตุการณ์ทั้งหมดในช่วงเวลาทั้งหมดด้วยλ.{\displaystyle \lambda .}แบ่งช่วงเวลาทั้งหมดออกเป็นn{\displaystyle n}ช่วงย่อยฉัน1,,ฉันn{\displaystyle I_{1},\dots ,I_{n}}มีขนาดเท่ากัน โดยที่n>λ{\displaystyle n>\lambda }(เนื่องจากเราสนใจเฉพาะส่วนเล็ก ๆ ของช่วงเวลาเท่านั้น สมมติฐานนี้จึงมีความหมาย) ซึ่งหมายความว่าจำนวนเหตุการณ์ที่คาดหวังในแต่ละ ช่วงย่อยทั้ง nช่วงนั้นเท่ากับλ/n.{\displaystyle \lambda /n.}

ตอนนี้เราสมมติว่าการเกิดเหตุการณ์ในช่วงเวลาทั้งหมดสามารถมองได้ว่าเป็นลำดับของการทดลองแบบเบอร์นูลลีจำนวน n ครั้ง โดยที่ฉัน{\displaystyle i}การทดลองแบบเบอร์นูลลีครั้งที่ - สอดคล้องกับการตรวจสอบว่าเหตุการณ์เกิดขึ้นในช่วงเวลาย่อยหรือไม่ฉันฉัน{\displaystyle I_{i}}ด้วยความน่าจะเป็นλ/n.{\displaystyle \lambda /n.}จำนวนเหตุการณ์ทั้งหมดที่คาดการณ์ไว้ในn{\displaystyle n}การพิจารณาคดีดังกล่าวจะเป็นλ,{\displaystyle \lambda ,}จำนวนเหตุการณ์ทั้งหมดที่คาดว่าจะเกิดขึ้นในช่วงเวลาทั้งหมด ดังนั้น สำหรับแต่ละส่วนย่อยของช่วงเวลา เราได้ประมาณการการเกิดเหตุการณ์เป็นกระบวนการเบอร์นูลลีในรูปแบบบี(n,λ/n).{\displaystyle {\textrm {B}}(n,\lambda /n).}ดังที่เราได้กล่าวไว้ก่อนหน้านี้ เราต้องการพิจารณาเฉพาะช่วงย่อยที่เล็กมากเท่านั้น ดังนั้นเราจึงใช้ลิมิตเป็นn{\displaystyle n}มีค่าเข้าสู่ค่าอนันต์

ในกรณีนี้การแจกแจงทวินามจะลู่เข้าสู่สิ่งที่เรียกว่าการแจกแจงปัวซงโดยทฤษฎีบทขีดจำกัดของปัวซ

ในตัวอย่างข้างต้นหลายๆ ตัวอย่าง เช่น จำนวนการกลายพันธุ์ในลำดับดีเอ็นเอที่กำหนด เหตุการณ์ที่ถูกนับนั้นแท้จริงแล้วเป็นผลลัพธ์ของการทดลองแบบไม่ต่อเนื่อง และจะสามารถจำลองได้อย่างแม่นยำยิ่งขึ้นโดยใช้การแจกแจงแบบทวินามนั่นคือ X~บี(n,พี).{\displaystyle X\sim {\textrm {B}}(n,p).}

ในกรณีเช่นนี้nจะมีค่ามากและpจะมีค่าน้อยมาก (ดังนั้นค่าเฉลี่ยnpจึงมีขนาดปานกลาง) จากนั้นจึงสามารถประมาณการกระจายตัวได้โดยใช้การกระจายแบบปัวซงซึ่งไม่ซับซ้อนกว่า X~ปัวส์(nพี).{\displaystyle X\sim {\textrm {Pois}}(np).}

บางครั้งการประมาณค่านี้เรียกว่ากฎของเหตุการณ์หายาก [ 72 ] : 5 เนื่องจาก เหตุการณ์เบอร์นูลลีแต่ละรายการ จำนวน n รายการเกิดขึ้นได้ยาก

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

ความแปรปรวนของการแจกแจงแบบทวินามเป็น1 − pเท่าของความแปรปรวนของการแจกแจงแบบปัวซง ดังนั้นจึงเกือบเท่ากันเมื่อpมีค่าเล็กมาก

บางครั้ง คำว่ากฎถูกใช้เป็นคำพ้องความหมายของการแจกแจงความน่าจะเป็นและการบรรจบกันในกฎหมายถึงการบรรจบกันในการแจกแจงดังนั้น การแจกแจงปัวซงจึงบางครั้งเรียกว่า "กฎของจำนวนน้อย" เพราะเป็นการแจกแจงความน่าจะเป็นของจำนวนครั้งของการเกิดเหตุการณ์ที่เกิดขึ้นไม่บ่อยนัก แต่มีโอกาสเกิดขึ้นมากมายกฎของจำนวนน้อยเป็นหนังสือของ Ladislaus Bortkiewicz เกี่ยวกับการแจกแจงปัวซง ตีพิมพ์ใน ปี 1898 [ 12 ] [ 73 ]

กระบวนการจุดปัวซง

การแจกแจงปัวซงเกิดขึ้นจากจำนวนจุดของกระบวนการจุดปัวซงที่อยู่ในบริเวณจำกัดบางแห่ง โดยเฉพาะอย่างยิ่ง ถ้าDเป็นปริภูมิบริเวณ เช่น ปริภูมิยุคลิดR dซึ่ง| D |พื้นที่ ปริมาตร หรือโดยทั่วไปแล้ว การวัดแบบเลเบสของบริเวณนั้นมีค่าจำกัด และถ้าN ( D )แทนจำนวนจุดในDแล้ว

พี(เอ็น(ดี)=เค)=(λ|ดี|)เคอีλ|ดี|เค!.{\displaystyle P(N(D)=k)={\frac {(\lambda |D|)^{k}e^{-\lambda |D|}}{k!}}.}

การถดถอยปัวซงและการถดถอยทวินามเชิงลบ

การถดถอยแบบปัวซงและ การถดถอยแบบทวิ นามเชิงลบมีประโยชน์สำหรับการวิเคราะห์ที่ตัวแปรตาม (ตัวแปรตอบสนอง) คือจำนวนนับ(0, 1, 2, ... )ของจำนวนเหตุการณ์หรือการเกิดขึ้นในช่วงเวลาหนึ่ง

ชีววิทยา

การทดลอง ของลูเรีย-เดลบรุคได้ทดสอบกับสมมติฐานของการวิวัฒนาการแบบลามาร์ค ซึ่งควรจะส่งผลให้เกิดการแจกแจงแบบปัวซง

Katz และ Miledi วัดศักย์เยื่อหุ้มเซลล์ทั้งที่มีและไม่มีอะเซทิลโคลีน (ACh) [ 74 ]เมื่อมี ACh อยู่ช่องไอออนบนเยื่อหุ้มเซลล์จะเปิดแบบสุ่มในช่วงเวลาสั้นๆ เนื่องจากมีช่องไอออนจำนวนมากซึ่งแต่ละช่องเปิดในช่วงเวลาสั้นๆ จำนวนช่องไอออนทั้งหมดที่เปิดในแต่ละช่วงเวลาจึงมีการกระจายแบบปัวซง เมื่อไม่มี ACh อยู่ ในทางปฏิบัติแล้วจะไม่มีช่องไอออนเปิด ศักย์เยื่อหุ้มเซลล์คือวี=เอ็นเปิดวีไอออน+วี0+วีเสียงรบกวน{\displaystyle V=N_{\text{open}}V_{\text{ion}}+V_{0}+V_{\text{noise}}}เมื่อหักลบผลกระทบของสัญญาณรบกวนแล้ว คัตซ์และมิเลดีพบว่าค่าเฉลี่ยและความแปรปรวนของศักย์เยื่อหุ้มเซลล์คือ8.5×103วี{\displaystyle 8.5\times 10^{-3}\;\mathrm {V} }และ(29.2×106วี)2{\displaystyle (29.2\times 10^{-6}\;\mathrm {V} )^{2}}ตามลำดับ โดยให้วีไอออน=107วี{\displaystyle V_{\text{ion}}=10^{-7}\;\mathrm {V} }(หน้า 94-95 [ 75 ] )

ในระหว่างการจำลองเซลล์แต่ละครั้ง จำนวนการกลายพันธุ์จะกระจายตัวแบบปัวซงโดยประมาณ[ 76 ]ตัวอย่างเช่น ไวรัส HIV มี เบสคู่ 10,000 คู่ และมีอัตราการกลายพันธุ์ประมาณ 1 ครั้งต่อ เบสคู่ 30,000 คู่ ซึ่งหมายความว่าจำนวนการกลายพันธุ์ต่อการจำลองแต่ละครั้งจะกระจายตัวดังนี้พีโอฉัน(1/3){\displaystyle \mathrm {Pois} (1/3)}(หน้า 64 [ 75 ] )

การประยุกต์ใช้ในด้านวิทยาศาสตร์อื่นๆ

ในกระบวนการปัวซง จำนวนเหตุการณ์ที่สังเกตได้จะผันผวนรอบค่าเฉลี่ยλโดยมีค่าเบี่ยงเบนมาตรฐานσเค=λ.{\displaystyle \sigma _{k}={\sqrt {\lambda }}.}ความผันผวนเหล่านี้เรียกว่าสัญญาณรบกวนปัวซงหรือ (โดยเฉพาะในด้านอิเล็กทรอนิกส์) สัญญาณรบกวนช็อ

ความสัมพันธ์ระหว่างค่าเฉลี่ยและค่าเบี่ยงเบนมาตรฐานในการนับเหตุการณ์ที่เกิดขึ้นแยกกันอย่างอิสระนั้นมีประโยชน์ทางวิทยาศาสตร์ โดยการตรวจสอบว่าความผันผวนเปลี่ยนแปลงไปอย่างไรเมื่อเทียบกับสัญญาณเฉลี่ย เราสามารถประมาณการส่วนร่วมของเหตุการณ์แต่ละครั้งได้แม้ว่าส่วนร่วมนั้นจะเล็กเกินกว่าจะตรวจจับได้โดยตรงก็ตาม ตัวอย่างเช่น ประจุeบนอิเล็กตรอนสามารถประมาณได้โดยการหาความสัมพันธ์ระหว่างขนาดของกระแสไฟฟ้ากับสัญญาณรบกวนแบบสุ่ม (shot noise ) ถ้ามี อิเล็กตรอน Nตัวผ่านจุดหนึ่งในช่วงเวลาtโดยเฉลี่ยกระแสไฟฟ้าเฉลี่ย คือฉัน=อีเอ็น/ที{\displaystyle I=eN/t}เนื่องจากความผันผวนในปัจจุบันน่าจะอยู่ในระดับประมาณนี้σฉัน=อีเอ็น/ที{\displaystyle \sigma _{I}=e{\sqrt {N}}/t}(เช่น ค่าเบี่ยงเบนมาตรฐานของกระบวนการปัวซง ) ประจุอี{\displaystyle e}สามารถประมาณได้จากอัตราส่วนทีσฉัน2/ฉัน.{\displaystyle t\sigma _{I}^{2}/I.}

ตัวอย่างที่พบได้ทั่วไปคือความหยาบที่ปรากฏขึ้นเมื่อขยายภาพถ่าย ความหยาบนั้นเกิดจากความผันผวนแบบปัวซงของจำนวน อนุภาค เงิน ที่ลดลง ไม่ใช่เกิดจากอนุภาคแต่ละอนุภาคเอง การหาความสัมพันธ์ ระหว่าง ความหยาบกับระดับการขยายจะช่วยให้สามารถประเมินสัดส่วนการมีส่วนร่วมของอนุภาคแต่ละอนุภาคได้ (ซึ่งมีขนาดเล็กเกินกว่าจะมองเห็นได้ด้วยตาเปล่า)

ใน ทฤษฎี เซตเชิงสาเหตุองค์ประกอบที่ไม่ต่อเนื่องของปริภูมิเวลาจะเป็นไปตามการแจกแจงแบบปัวซงในปริมาตร

การแจกแจงแบบปัวซงยังปรากฏในกลศาสตร์ควอนตัมโดยเฉพาะอย่างยิ่งในทัศนศาสตร์ควอนตัมกล่าวคือ สำหรับ ระบบ ตัวสั่นฮาร์มอนิกควอนตัมในสถานะโคherentความน่าจะเป็นของการวัดระดับพลังงานเฉพาะนั้นมีการกระจายแบบปัวซง

วิธีการคำนวณ

การแจกแจงแบบปัวซงก่อให้เกิดภารกิจที่แตกต่างกันสองอย่างสำหรับไลบรารีซอฟต์แวร์เฉพาะทาง ได้แก่การประเมินการแจกแจงพี(เค;λ){\displaystyle P(k;\lambda )}และสุ่มเลือกหมายเลขตามการกระจายนั้น

การประเมินการแจกแจงปัวซง

การคำนวณพี(เค;λ){\displaystyle P(k;\lambda )}สำหรับที่กำหนดเค{\displaystyle k}และλ{\displaystyle \lambda }เป็นงานเล็กน้อยที่สามารถทำได้โดยใช้คำจำกัดความมาตรฐานของพี(เค;λ){\displaystyle P(k;\lambda )}ในแง่ของฟังก์ชันเลขชี้กำลัง ฟังก์ชันกำลัง และฟังก์ชันแฟกทอเรียล อย่างไรก็ตาม นิยามทั่วไปของการแจกแจงปัวซงประกอบด้วยสองพจน์ที่อาจทำให้เกิดค่าเกินขีดจำกัดในคอมพิวเตอร์ได้ง่าย คือλ kและk !เศษส่วนของλ kต่อk !อาจทำให้เกิดข้อผิดพลาดในการปัดเศษซึ่งมีขนาดใหญ่มากเมื่อเทียบกับe λและดังนั้นจึงให้ผลลัพธ์ที่ผิดพลาด เพื่อความเสถียรทางตัวเลข ฟังก์ชันความน่าจะเป็นมวลของปัวซงจึงควรได้รับการประเมินดังนี้ เอฟ(เค;λ)=เอ็กซ์[เคlnλλlnΓ(เค+1)],{\displaystyle \!f(k;\lambda )=\exp \left[k\ln \lambda -\lambda -\ln \Gamma (k+1)\right],} ซึ่งเทียบเท่าทางคณิตศาสตร์แต่มีเสถียรภาพเชิงตัวเลขสามารถหา ค่าลอการิทึมธรรมชาติของ ฟังก์ชันแกมมาlgamma ได้โดยใช้ ฟังก์ชันใน ไลบรารีมาตรฐานของภาษา C (เวอร์ชัน C99) หรือRฟังก์ชันgammalnในMATLABหรือSciPyหรือlog_gammaฟังก์ชันในFortran  2008 และเวอร์ชันที่ใหม่กว่า

ภาษาโปรแกรมบางภาษามีฟังก์ชันในตัวสำหรับประเมินการแจกแจงปัวซง ได้แก่

  • R : ฟังก์ชันdpois(x, lambda);
  • Excel : ฟังก์ชันPOISSON( x, mean, cumulative)ที่มีแฟล็กสำหรับระบุการกระจายแบบสะสม
  • Mathematica : การแจกแจงปัวซงแบบเอกตัวแปรเป็น, [ 77 ]การแจกแจงปัวซงแบบทวิตัวแปรเป็น,. [ 78 ]PoissonDistribution[λ{\displaystyle \lambda }]MultivariatePoissonDistribution[θ12,{\displaystyle \theta _{12},}{ θ1θ12,{\displaystyle \theta _{1}-\theta _{12},}θ2θ12{\displaystyle \theta _{2}-\theta _{12}}}]

การสร้างตัวแปรสุ่ม

งานที่ยากกว่าคือการสุ่มเลือกตัวแปรสุ่ม จำนวนเต็ม จาก1การแจกแจงปัวซงโดยกำหนดค่าให้λ.{\displaystyle \lambda .}

โซลูชันนี้จัดทำโดย:

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 .
ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Poisson_distribution&oldid=1355788475#Generating_Poisson-distributed_random_variables "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ การแจกแจงปัวซง

Γ(⌊เค+1⌋,λ)⌊เค⌋!,{\displaystyle {\frac {\Gamma (\lfloor k+1\rfloor ,\lambda )}{\lfloor k\rfloor !}},} หรืออี−λ∑เจ=0⌊เค⌋λเจเจ!

ประวัติศาสตร์

การนำเสนอการแจกแจงปัวซงได้รับการยกย่องให้แก่นักคณิตศาสตร์และนักฟิสิกส์ชาวฝรั่งเศส Siméon Denis Poisson (1781–1840) ซึ่งตีพิมพ์พร้อมกับทฤษฎีความน่าจะเป็นของเขาใน Recherches sur la probabilité des jugements en matière criminelle et en matière civile (1837) [ 4...

ฟังก์ชันมวลความน่าจะเป็น

ตัวแปรสุ่ม แบบ ไม่ต่อเนื่อง X กล่าวได้ว่ามีการแจกแจงแบบปัวซง (Poisson distribution) ที่มีพารามิเตอร์ 0"}}"> 0}"> λ > 0 {\displaystyle \lambda >0} 0}"> หากมี ฟังก์ชันมวลความน่าจะเป็น ที่กำหนดโดย: [ 2 ] : 60 เอฟ ( เค ; λ ) = ปร. ( X = เค ) = λ เค อี − λ เค !

ตัวอย่าง

การแจกแจงแบบปัวซงอาจมีประโยชน์ในการจำลองเหตุการณ์ต่างๆ เช่น: