负二项分布
不同来源对负二项分布的定义略有差异:随机变量的最小可能取值可能是 (仅计失败的次数,或反之),亦可能是 (总次数,不论成败);参数 可能表示每次试验成功的概率,也可能表示失败的概率;试验的终止条件可能是成功 次或失败 次。[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)
帕斯卡分布概率质量函数
对所有可能
值求和,一定等于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.