負二項分布
不同來源對負二項分布的定義略有差異:隨機變量的最小可能取值可能是 (僅計失敗的次數,或反之),亦可能是 (總次數,不論成敗);參數 可能表示每次試驗成功的概率,也可能表示失敗的概率;試驗的終止條件可能是成功 次或失敗 次。[1] |
機率質量函數
|
參數 |
(實)
(實) |
---|
值域 |
![{\displaystyle k\in \{0,1,2,\ldots \}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eeb23cb33394a6e5787b660d262e699a46d7ac53) |
---|
機率質量函數 |
![{\displaystyle {\frac {\Gamma (r+k)}{k!\,\Gamma (r)}}\,p^{r}\,(1-p)^{k}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/29b36036bd48d80042165f08945251a7cb3aa154) |
---|
累積分布函數 |
![{\displaystyle I_{p}(r,k+1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f46c8a380b0adc968b21fb2f409d302dd9f72cb) |
---|
期望值 |
![{\displaystyle r\,{\frac {1-p}{p}}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/692eb5d9ba406e2347f0e851b67fdc6f9da47c04) |
---|
眾數 |
![{\displaystyle \lfloor (r-1)\,(1-p)/p\rfloor {\text{ if }}r>1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8754b51bae6c798a900c32fb75e18d4c3c0ab87a)
![{\displaystyle 0{\text{ if }}r\leq 1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e17a51ee4bc9e1274904c826d0bcc38f3019006) |
---|
變異數 |
![{\displaystyle r\,{\frac {1-p}{p^{2}}}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/15c426a081d840dd354df29be717c63dcf097c43) |
---|
偏度 |
![{\displaystyle {\frac {2-p}{\sqrt {r\,(1-p)}}}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/132d6d339904324fc484d791b8a262ea2fba6d39) |
---|
峰度 |
![{\displaystyle {\frac {6}{r}}+{\frac {p^{2}}{r\,(1-p)}}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5f2b83ca53011e8be585ffefffc609afdd64ef6a) |
---|
動差母函數 |
![{\displaystyle \left({\frac {p}{1-(1-p)e^{t}}}\right)^{r}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1563a81bf6c916d0cd396cf7f08ddaae6afccd5a) |
---|
特徵函數 |
![{\displaystyle \left({\frac {p}{1-(1-p)e^{i\,t}}}\right)^{r}\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aadf94ebd07b67261b89ec2c7c901f246c7f786f) |
---|
負二項分布(Negative binomial distribution)是統計學上一種描述在一系列獨立同分布的伯努利試驗中,成功次數達到指定次數(記為
)時失敗次數的離散概率分布。比如,如果我們定義擲骰子隨機變量
值為
時成功,所有
為失敗,這時我們反覆擲骰子直到1出現3次(成功次數
),此時非1數字出現次數的概率分布即為負二項分布。
帕斯卡分布(Pascal distribution,來自布萊茲·帕斯卡 (Blaise Pascal))和波利亞分布(Polya distribution,又稱罐子模型,來自喬治·波利亞 (George Pólya))均是負二項分布的特例。在工程、氣候等領域中經常用「負二項分布」或「帕斯卡分布」來描述變量
為整數的情況,而使用「波利亞分布」來描述
取到實數值
的情況。
對於「相關的離散事件」("associated discrete events")的發生,例如龍捲風爆發,相比於泊松分布,波利亞分布由於允許其平均值和方差不同,而能夠給出更精確的模型。在流行病學中,它已被用於模擬傳染病的疾病傳播,其中可能的繼發感染數量可能因個體和環境而異[2]。 更一般地說,由於正協方差項,事件具有正相關的事件導致比獨立事件更大的方差可能是合適的。
「負二項分布」與「二項分布」的區別在於:「二項分布」是固定試驗總次數
的獨立試驗中,成功次數k的分布;而「負二項分布」是所有到r次成功時即終止的獨立試驗中,失敗次數k的分布。
術語「負二項式」可能是因為出現在分布的概率質量函數公式中的某個二項式係數可以用負數更簡單地寫出[3]。
定義
若每次伯努利試驗有兩種可能的結果,分別為成功或者失敗。在每次試驗中,成功的概率為
,失敗的概率為
。反覆進行該伯努利試驗,直到觀察到第
次成功發生。此時試驗失敗次數
的分布即為負二項分布(或稱帕斯卡分布),那麼:
- 若隨機變量
服從參數為
和
的負二項分布,則記為
.
在實際生活中,我們可以使用負二項分布描述某種機器在壞掉前,能夠工作的天數的分布。此時,「成功」的事件可以指機器正常工作一天,「失敗」的事件可以指機器故障的一天。如果我們使用負二項分布來描述運動員在獲取r個獎牌前嘗試的次數的分布,此時,「失敗」的事件指運動員的一次嘗試,「成功」的事件指運動員獲取一枚獎牌。如果使用負二項分布來描述擲一枚硬幣出現r次正面前,出現硬幣反面的次數的分布,「成功」的事件指出現硬幣的正面,「失敗」的事件指出現硬幣的反面。
概率質量函數
帕斯卡分布
當
是整數時的負二項分布又稱帕斯卡分布,其概率質量函數為:
其中
是失敗的次數,
是成功的次數,
是事件成功的概率。在負二項分布的概率質量函數中,由於
次伯努利試驗為獨立同分布,每個成功
次、失敗
次的事件的概率為
。由於第
次成功一定是最後一次試驗,所以應該在
次試驗中選擇
次成功,使用排列組合二項係數獲取所有可能的選擇數。
二項係數與負二項名稱來源
括號中為二項式係數表達式:
![{\displaystyle {\binom {k+r-1}{r-1}}={\frac {(k+r-1)!}{k!\,(r-1)!}}={\frac {(k+r-1)(k+r-2)\dotsm (r)}{k!}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4233f4c3dfbc50adb8020c92513fa80ac89ff142)
該表達式可以寫成帶負值參數的二項係數的形式,如下式所示,解釋了「負二項」名稱的來源:
![{\displaystyle {\begin{aligned}&{\frac {(k+r-1)\dotsm (r)}{k!}}\\[6pt]={}&(-1)^{k}{\frac {(-r)(-r-1)(-r-2)\dotsm (-r-k+1)}{k!}}=(-1)^{k}{\binom {-r}{k}}.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6ffcb8bf581835d46e6607be6a55cc58882ee921)
概率質量函數對所有可能k值求和為1
帕斯卡分布概率質量函數
對所有可能
值求和,一定等於1:
證明如下:
其中第三步用到了二項序列展開。
幾何分布
取
,負二項分布等於幾何分布。其概率質量函數為
。
例子
舉例說,若我們擲骰子,擲到一即視為成功。則每次擲骰的成功率是
。要擲出三次一,所需的擲骰次數屬於集合{ 3, 4, 5, 6, ... }。擲到三次一的擲骰次數是負二項分布的隨機變數。要在第三次擲骰時,擲到第三次一,則之前兩次都要擲到一,其機率為
。注意擲骰是伯努利試驗,之前的結果不影響隨後的結果。
若要在第四次擲骰時,擲到第三次一,則之前三次之中要有剛好兩次擲到一,在三次擲骰中擲到2次1的機率為
。第四次擲骰要擲到一,所以要將前面的機率再乘
:
。
相關分布
幾何分布(在 { 0, 1, 2, 3, ... } 上)是負二項分布的一個特例,其中
![{\displaystyle \operatorname {Geom} (p)=\operatorname {NB} (1,\,1-p).\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ed323ec384df0476b4d752add267a989163533bb)
參見
參考文獻
- ^ DeGroot, Morris H. Probability and Statistics Second. Addison-Wesley. 1986: 258–259. ISBN 0-201-11366-X. LCCN 84006269. OCLC 10605205.
- ^ e.g. J.O. Lloyd-Smith, S.J. Schreiber, P.E. Kopp, and W.M. Getz (2005), Superspreading and the effect of individual variation on disease emergence, Nature, 438, 355–359. doi:10.1038/nature04153
The overdispersion parameter is usually denoted by the letter
in epidemiology, rather than
as here.
- ^ Casella, George; Berger, Roger L. Statistical inference
2nd. Thomson Learning. 2002: 95. ISBN 0-534-24312-6.