概率统计14 | 您所在的位置:网站首页 › 几何分布的公式是什么 › 概率统计14 |
我家小朋友年方1岁半,家里每天上午都要出去遛小孩。现在小朋友有两项爱好,在家翻垃圾桶,出门捡烟头。 翻垃圾桶可以有效地限制,捡烟头可是防不胜防。 也许烟头能散发出特殊的能量波动,小区的绿化带和草坪上的大部分烟头都能被小朋友准确地发现,如果他在不规则的前进路线中突然停下了,那肯定是看到了新的烟头。 在我的严密监视下,小朋友捡烟头的几率已经从原来的“绝不放过”下降到了现在的20%,如果他直到回家还没有烟头,就能得到一个棒棒糖作为奖励。当然,每个烟头对他来说仍然存在着独特的魅力,在最初遇到的3个烟头中没有经受住考验的概率是多少? 设随机变量X是小朋友最终捡起烟头时所发现的烟头个数,那么: 最初遇到的3个烟头中没有抵抗住诱惑的概率可以用P(X ≤ 3)表示: 每个烟头都是一次考验,我的关注点是,小朋友在第几次考验时会捡起烟头? 质量函数小朋友每次遇到烟头的行为都是一个独立的随机试验,试验只有成功和失败两种结果,且对于每个随机试验来说,成功的概率都是相同的。对于小朋友自己来说,捡起烟头是成功,错过才是失败,只要成功一次就没有棒棒糖,试验结束。现在用p表示成功率(捡起烟头的概率),q = 1 – p表示失败率,随机变量X表示第一次成功时所经历的试验次数,那么在进行了r次试验后才遇到第一次成功的概率(或者说在第r次试验取得成功前,需要经历r-1次失败的概率)可以表示为: 这个分布就是几何分布(Geometric distribution),X服从几何分布,记为X~GE(p)。 下面的python代码展示了r取不同值时的P(X=r)。 1 import numpy as np 2 import matplotlib.pyplot as plt 3 from scipy import stats 4 5 p = 0.2 # 成功率 6 print('X ~ GE({0})'.format(p)) 7 8 rs = np.array(range(1, 11, 1)) # 随机变量的取值 9 # 几何分布 X ~ GE(0.2) 10 ps = stats.geom.pmf(rs, p) # 每个随机变量对应的概率 11 for i, r in enumerate(rs): 12 print('P(X={0})={1}'.format(r, ps[i])) 13 14 plt.bar(left=rs, height=ps, width=0.5) 15 plt.xlabel('r, X=r表示第r次试验才成功') 16 plt.ylabel('P(X=r)') 17 plt.rcParams['font.sans-serif'] = ['SimHei'] # 用来正常显示中文标签 18 plt.show()X ~ GE(0.2) P(X=1)=0.2 P(X=2)=0.16000000000000003 P(X=3)=0.12800000000000003 P(X=4)=0.10240000000000003 P(X=5)=0.08192000000000002 P(X=6)=0.06553600000000002 P(X=7)=0.052428800000000025 P(X=8)=0.041943040000000015 P(X=9)=0.033554432000000016 P(X=10)=0.026843545600000015 对于几何分布来说,r是大于等于1的自然数,并且只有成功的概率处于(0, 1)之间才有意义,因此当r = 1时,P(X = r)达到最大值,随着r的增大,P(X = r)也越来越小。 几何分布是一个典型的长尾分布,第一次就成功的概率是最高的,这看似有违直觉。这里需要注意的是,“第r次才取得成功”和“在r次之内取得成功”是两回事,后者才是概率越来越大。 实际中有不少随机变量服从几何分布,例如某产品的不合格率为0.05,则首次查到不合格品的检查次数X~GE(0.05)。 分布的形状现在回到问题的关注点:小朋友在第几次考验时会捡起烟头?这个问题并不能确切地回答,可以回答的是,在r次考验之内捡起烟头的概率。这个概率可以用分布函数表示: 另一种计算方法是: P(X > r)表示在取得第一次成功时,前r次试验都失败的概率。F(r)更为专业的说法是:为得到1次成功而进行r次伯努利试验,r的概率分布。
下面的代码展示了几何分布的形状: 1 fs = stats.geom.cdf(rs, p) # 每个r对应的分布 2 for i, r in enumerate(rs): 3 print('F({0})=P(X |
CopyRight 2018-2019 实验室设备网 版权所有 |