UA MATH564 概率论 计算至少有一个发生的概率:Boole不等式
UA MATH564 概率论 计算至少有一个发生的概率:Boole不等式
- Boole不等式
- Gunias不等式
- 钟开莱-艾尔迪希(Chung-Erdos)不等式
概率作为一种特殊的测度,满足有限可加性:
P(⋃i=1nAi)=∑i=1nP(Ai),Ai∩Aj=ϕ,∀i≠jP(\bigcup_{i=1}^{n} A_i) = \sum_{i=1}^{n} P(A_i),A_i \cap A_j = \phi,\forall i \ne jP(i=1⋃nAi)=i=1∑nP(Ai),Ai∩Aj=ϕ,∀i=j
Ai,AjA_i,A_jAi,Aj作为样本空间σ\sigmaσ-代数中的元素,Ai∩Aj=ϕA_i \cap A_j = \phiAi∩Aj=ϕ的概率解释是Ai,AjA_i,A_jAi,Aj这两个事件互斥。上式左边的含义是互斥事件A1,A2,⋯,AnA_1,A_2,\cdots,A_nA1,A2,⋯,An中至少有一个发生的概率,右边的含义是互斥事件A1,A2,⋯,AnA_1,A_2,\cdots,A_nA1,A2,⋯,An的概率之和。但在具体应用中,事件互斥的假设实在是太强了,我们希望在没有互斥这个条件下,也能有一个计算P(⋃i=1nAi)P(\bigcup_{i=1}^{n} A_i)P(⋃i=1nAi)的方法。
Boole不等式
尝试估计某个量的值之前,我们往往会尝试找到一些上界和下界,然后试图找到更小的上界和更大的下界,逐步逼近准确值。要寻找P(⋃i=1nAi)P(\bigcup_{i=1}^{n} A_i)P(⋃i=1nAi)的上界和下界,最容易想到的是外测度的有限次可加性:
P(⋃i=1nAi)≤∑i=1nP(Ai)P(\bigcup_{i=1}^{n} A_i) \le \sum_{i=1}^{n} P(A_i)P(i=1⋃nAi)≤i=1∑nP(Ai)
当且仅当A1,A2,⋯,AnA_1,A_2,\cdots,A_nA1,A2,⋯,An互斥时取等。这个不等式叫做Boole不等式。要证明这个不等式也非常容易,我们构造一个集列{Bi}i=1n\{B_i\}_{i=1}^n{Bi}i=1n,其中B1=A1,B2=A2−A1,⋯,Bn=An−⋃i=1n−1AnB_1 = A_1,B_2 = A_2-A_1,\cdots,B_n = A_n - \bigcup_{i=1}^{n-1}A_nB1=A1,B2=A2−A1,⋯,Bn=An−⋃i=1n−1An,显然Bi∩Bj=ϕ,∀i≠jB_i \cap B_j = \phi, \forall i \ne jBi∩Bj=ϕ,∀i=j,并且
⋃i=1nBi=⋃i=1nAi\bigcup_{i=1}^n B_i = \bigcup_{i=1}^{n} A_ii=1⋃nBi=i=1⋃nAi
根据概率的可加性和单调性(Bi⊂Ai⇒P(Bi)≤P(Ai)B_i \subset A_i \Rightarrow P(B_i) \le P(A_i)Bi⊂Ai⇒P(Bi)≤P(Ai))
P(⋃i=1nAi)=P(⋃i=1nBi)=∑i=1nP(Bi)≤∑i=1nP(Ai)P(\bigcup_{i=1}^{n} A_i) = P(\bigcup_{i=1}^n B_i ) = \sum_{i=1}^n P(B_i) \le \sum_{i=1}^{n} P(A_i)P(i=1⋃nAi)=P(i=1⋃nBi)=i=1∑nP(Bi)≤i=1∑nP(Ai)
基于de Morgan律很容易就能得到Boole不等式的另一种形式。
P(⋂i=1nAi)=1−P(⋂i=1nAiˉ)=1−P(⋃i=1nAˉi)≥1−∑i=1nP(Aˉi)=1−∑i=1n[1−P(Ai)]=∑i=1nP(Ai)−(n−1)P(\bigcap_{i=1}^n A_i) = 1-P(\bar{\bigcap_{i=1}^n A_i}) = 1-P(\bigcup_{i=1}^n \bar{A}_i) \ge 1 - \sum_{i=1}^{n} P(\bar{A}_i) \\ = 1-\sum_{i=1}^{n} [1-P(A_i)] = \sum_{i=1}^{n} P(A_i) - (n-1)P(i=1⋂nAi)=1−P(i=1⋂nAiˉ)=1−P(i=1⋃nAˉi)≥1−i=1∑nP(Aˉi)=1−i=1∑n[1−P(Ai)]=i=1∑nP(Ai)−(n−1)
Boole不等式给出了非常粗糙的上界和下界,后续的任务就是试图推导出更精确的上界和下界。Gunias不等式和钟开莱-艾尔迪希(Chung-Erdos)不等式分别提供了至少有一个事件发生概率的更精确的上界和下界。
Gunias不等式
因为现在没有A1,A2,⋯,AnA_1,A_2,\cdots,A_nA1,A2,⋯,An互斥的假设,所以他们两两之间交集不一定为0,在这nnn个事件中,总有一个事件与其他事件“重合程度”最大,即∃k∈{1,2,⋯,n}\exists k \in \{1,2,\cdots,n\}∃k∈{1,2,⋯,n},
∑i=1,i≠kP(Ai∩Ak)≥∑i=1,i≠jP(Ai∩Aj),∀j∈{1,2,⋯,n}\sum_{i=1,i \ne k}P(A_i \cap A_k) \ge \sum_{i=1,i \ne j}P(A_i \cap A_j),\forall j \in \{1,2,\cdots,n\}i=1,i=k∑P(Ai∩Ak)≥i=1,i=j∑P(Ai∩Aj),∀j∈{1,2,⋯,n}
下面我们构造两个集列{Bi}i=1n\{B_i\}_{i=1}^n{Bi}i=1n与{Ci}i=1n\{C_i\}_{i=1}^n{Ci}i=1n:
B1=Ak,B2=A1−Ak,⋯,Bk=Ak−1−Ak,Bk+1=Ak+1−Ak,⋯,Bn=An−AkC1=Ak,C2=A1−Ak,C3=A2−A1−Ak,⋯,Ck=Ak−1−⋃j=1,j≠k−1kAjCk+1=Ak+1−⋃j=1kAj,⋯,Cn=An−⋃j=1n−1AjB_1 = A_k,B_2=A_1-A_k,\cdots,B_{k} = A_{k-1} - A_k,B_{k+1} = A_{k+1}-A_{k},\cdots,B_n = A_n - A_k \\ C_1 = A_k,C_2 = A_1-A_k,C_3 = A_2 - A_1 - A_k,\cdots,C_{k} = A_{k-1} - \bigcup_{j=1,j \ne k-1}^kA_j \\ C_{k+1} = A_{k+1} - \bigcup_{j=1}^{k} A_j,\cdots,C_n = A_n - \bigcup_{j=1}^{n-1} A_jB1=Ak,B2=A1−Ak,⋯,Bk=Ak−1−Ak,Bk+1=Ak+1−Ak,⋯,Bn=An−AkC1=Ak,C2=A1−Ak,C3=A2−A1−Ak,⋯,Ck=Ak−1−j=1,j=k−1⋃kAjCk+1=Ak+1−j=1⋃kAj,⋯,Cn=An−j=1⋃n−1Aj
显然Ci⊂BiC_i \subset B_iCi⊂Bi,因此
⋃i=1nCi⊂⋃i=1nBi\bigcup_{i=1}^n C_i \subset \bigcup_{i=1}^n B_ii=1⋃nCi⊂i=1⋃nBi
集列{Ci}i=1n\{C_i\}_{i=1}^n{Ci}i=1n与Boole不等式证明中构造的集列类似,
⋃i=1nCi=⋃i=1nAi\bigcup_{i=1}^n C_i = \bigcup_{i=1}^n A_i i=1⋃nCi=i=1⋃nAi
根据概率的单调性,
P(⋃i=1nAi)=P(⋃i=1nCi)≤P(⋃i=1nBi)P(\bigcup_{i=1}^n A_i) = P(\bigcup_{i=1}^n C_i) \le P(\bigcup_{i=1}^n B_i)P(i=1⋃nAi)=P(i=1⋃nCi)≤P(i=1⋃nBi)
根据Boole不等式
P(⋃i=1nBi)≤∑i=1nP(Bi)=P(Ak)+P(A1−Ak)+⋯+P(Ak−1−Ak)+P(Ak+1−Ak)+⋯+P(An−Ak)=P(Ak)+P(A1)−P(A1∩Ak)+⋯+P(Ak−1)−P(Ak−1∩Ak)+P(Ak+1)−P(Ak+1∩Ak)+⋯+P(An)−P(An∩Ak)=∑i=1nP(Ai)−∑j=1,j≠knP(Aj∩Ak)P(\bigcup_{i=1}^n B_i) \le \sum_{i=1}^n P(B_i) = P(A_k) + P(A_1-A_k) + \cdots \\+ P(A_{k-1}-A_k) + P(A_{k+1}-A_k) + \cdots + P(A_n - A_k) \\ = P(A_k) + P(A_1)-P(A_1\cap A_k) + \cdots + P(A_{k-1})-P(A_{k-1} \cap A_k) \\+ P(A_{k+1})-P(A_{k+1}\cap A_k) + \cdots + P(A_n) - P(A_n \cap A_k) \\ = \sum_{i=1}^n P(A_i) - \sum_{j=1,j \ne k}^{n}P(A_j \cap A_k)P(i=1⋃nBi)≤i=1∑nP(Bi)=P(Ak)+P(A1−Ak)+⋯+P(Ak−1−Ak)+P(Ak+1−Ak)+⋯+P(An−Ak)=P(Ak)+P(A1)−P(A1∩Ak)+⋯+P(Ak−1)−P(Ak−1∩Ak)+P(Ak+1)−P(Ak+1∩Ak)+⋯+P(An)−P(An∩Ak)=i=1∑nP(Ai)−j=1,j=k∑nP(Aj∩Ak)
这个就是Gunias不等式,也可以将其写成:
P(⋃i=1nAi)≤∑i=1nP(Ai)−supk∈{1,⋯,n}∑j=1,j≠knP(Aj∩Ak)P(\bigcup_{i=1}^n A_i) \le \sum_{i=1}^n P(A_i) - \sup_{k \in \{1,\cdots,n\}} \sum_{j=1,j \ne k}^{n}P(A_j \cap A_k)P(i=1⋃nAi)≤i=1∑nP(Ai)−k∈{1,⋯,n}supj=1,j=k∑nP(Aj∩Ak)
钟开莱-艾尔迪希(Chung-Erdos)不等式
Boole不等式给出的至少有一个事件发生的概率的下界非常粗糙,甚至都不能保证下界是非负的,但钟开莱-艾尔迪希不等式给出的下界一定是非负的:
P(⋃i=1nAi)≥(∑i=1nP(Ai))2∑i,j=1nP(Ai∩Aj)P(\bigcup_{i=1}^n A_i) \ge \frac{\left( \sum_{i=1}^n P(A_i) \right)^2}{\sum_{i,j=1}^n P(A_i \cap A_j)}P(i=1⋃nAi)≥∑i,j=1nP(Ai∩Aj)(∑i=1nP(Ai))2
证明
当n=2n=2n=2时,钟开莱-艾尔迪希不等式为
P(A1∪A2)≥P(A1)2+P(A2)2+2P(A1)P(A2)P(A1)+P(A2)+2P(A1∩A2)⇐[P(A1)+P(A2)+2P(A1∩A2)]P(A1∪A2)≥P(A1)2+P(A2)2+2P(A1)P(A2)⇐P(A1)[P(A1∪A2)−P(A1)]+P(A2)[P(A1∪A2)−P(A2)]+2[P(A1∩A2)P(A1∪A2)−P(A1)P(A2)]≥0⇐2P(A1)P(A2)−[P(A1)+P(A2)]P(A1∩A2)+2[P(A1∩A2)P(A1∪A2)−P(A1)P(A2)]≥0⇐[2P(A1∪A2)−(P(A1)+P(A2))]P(A1∩A2)≥0P(A_1 \cup A_2) \ge \frac{P(A_1)^2+ P(A_2)^2 + 2P(A_1)P(A_2)}{P(A_1) + P(A_2) + 2P(A_1 \cap A_2)} \\ \Leftarrow [P(A_1) + P(A_2) + 2P(A_1 \cap A_2)]P(A_1 \cup A_2) \ge P(A_1)^2+ P(A_2)^2 + 2P(A_1)P(A_2) \\ \Leftarrow P(A_1)[P(A_1 \cup A_2)-P(A_1)] + P(A_2)[P(A_1 \cup A_2)-P(A_2)] \\+2[P(A_1 \cap A_2)P(A_1 \cup A_2)-P(A_1)P(A_2)] \ge 0 \\ \Leftarrow 2P(A_1)P(A_2) - [P(A_1)+P(A_2)]P(A_1\cap A_2) \\ +2[P(A_1 \cap A_2)P(A_1 \cup A_2)-P(A_1)P(A_2)] \ge 0 \\ \Leftarrow [2P(A_1\cup A_2)-(P(A_1)+P(A_2))]P(A_1 \cap A_2) \ge 0P(A1∪A2)≥P(A1)+P(A2)+2P(A1∩A2)P(A1)2+P(A2)2+2P(A1)P(A2)⇐[P(A1)+P(A2)+2P(A1∩A2)]P(A1∪A2)≥P(A1)2+P(A2)2+2P(A1)P(A2)⇐P(A1)[P(A1∪A2)−P(A1)]+P(A2)[P(A1∪A2)−P(A2)]+2[P(A1∩A2)P(A1∪A2)−P(A1)P(A2)]≥0⇐2P(A1)P(A2)−[P(A1)+P(A2)]P(A1∩A2)+2[P(A1∩A2)P(A1∪A2)−P(A1)P(A2)]≥0⇐[2P(A1∪A2)−(P(A1)+P(A2))]P(A1∩A2)≥0
显然2P(A1∪A2)−(P(A1)+P(A2))≥02P(A_1\cup A_2)-(P(A_1)+P(A_2)) \ge 02P(A1∪A2)−(P(A1)+P(A2))≥0,因此当n=2n=2n=2时,钟开莱-艾尔迪希不等式成立。
假设n=kn=kn=k时,钟开莱-艾尔迪希不等式成立,则n=k+1n=k+1n=k+1时,
P(⋃i=1k+1Ai)=P(Ak+1∪⋃i=1kAi)≥P(Ak+1)2+P(⋃i=1kAi)2+2P(A1)P(⋃i=1kAi)P(Ak+1)+P(⋃i=1kAi)+2P(Ak+1∩⋃i=1kAi)P(\bigcup_{i=1}^{k+1} A_i) = P(A_{k+1} \cup \bigcup_{i=1}^k A_i ) \ge \frac{P(A_{k+1})^2+ P(\bigcup_{i=1}^k A_i )^2 + 2P(A_1)P(\bigcup_{i=1}^k A_i )}{P(A_{k+1}) + P(\bigcup_{i=1}^k A_i ) + 2P(A_{k+1} \cap \bigcup_{i=1}^k A_i )}P(i=1⋃k+1Ai)=P(Ak+1∪i=1⋃kAi)≥P(Ak+1)+P(⋃i=1kAi)+2P(Ak+1∩⋃i=1kAi)P(Ak+1)2+P(⋃i=1kAi)2+2P(A1)P(⋃i=1kAi)
根据Boole不等式,
P(Ak+1)+P(⋃i=1kAi)+2P(Ak+1∩⋃i=1kAi)≤P(Ak+1)+∑i=1kP(Ai)+2∑i=1kP(Ak+1∩Ai)≤∑i,j=1k+1P(Ai∩Aj)P(A_{k+1}) + P(\bigcup_{i=1}^k A_i ) + 2P(A_{k+1} \cap \bigcup_{i=1}^k A_i ) \\ \le P(A_{k+1}) + \sum_{i=1}^k P(A_i) + 2\sum_{i=1}^k P(A_{k+1} \cap A_i)\le \sum_{i,j=1}^{k+1} P(A_i \cap A_j)P(Ak+1)+P(i=1⋃kAi)+2P(Ak+1∩i=1⋃kAi)≤P(Ak+1)+i=1∑kP(Ai)+2i=1∑kP(Ak+1∩Ai)≤i,j=1∑k+1P(Ai∩Aj)
对P(⋃i=1kAi)P(\bigcup_{i=1}^k A_i )P(⋃i=1kAi)使用钟开莱-艾尔迪希不等式,说明分子大于(∑i=1k+1P(Ai))2\left( \sum_{i=1}^{k+1} P(A_i) \right)^2(∑i=1k+1P(Ai))2,这一步留给读者自己计算。综上,钟开莱-艾尔迪希不等式对所有的自然数都成立。
UA MATH564 概率论 计算至少有一个发生的概率:Boole不等式相关推荐
- UA MATH564 概率论 计算至少有一个发生的概率:容斥原理与庞加莱公式
UA MATH564 概率论 计算至少有一个发生的概率:容斥原理与庞加莱公式 事件的并的Poincare公式 事件的交的Poincare公式 上一讲介绍了P(⋃i=1nAi)P(\bigcup_{i= ...
- UA MATH564 概率论 计算至少有一个发生的概率:Waring公式
UA MATH564 概率论 计算至少有一个发生的概率:Waring公式 基于组合学的证明方法 基于概率论的证明方法 在推出Poincare公式后,计算nnn个事件中至少有一个事件发生的概率的问题就彻 ...
- UA MATH564 概率论 QE练习题 信封问题
UA MATH564 概率论 QE练习题 信封问题 2015年1月的第二题 2015年5月的第一题 这一篇介绍QE理论中出现了两个信封问题相关的题目:2015年1月的第二题和2015年5月的第一题.信 ...
- UA MATH564 概率论 高阶矩的计算:有限差算子方法2
UA MATH564 概率论 高阶矩的计算:有限差算子方法2 二项分布的高阶矩 Poisson分布的高阶矩 负二项分布的高阶矩 这一讲介绍有限差算子方法怎么用来计算常见的离散型随机变量的高阶矩.这里提 ...
- UA MATH564 概率论 高阶矩的计算:有限差算子方法1
UA MATH564 概率论 高阶矩的计算:有限差算子方法1 有限差算子及其性质 0的有限差与Stirling数之间的关系 计算离散型随机变量的n阶矩 有限差算子(finite difference) ...
- UA MATH564 概率论 依概率收敛的一个例题
UA MATH564 概率论 依概率收敛的一个例题 Part (a) Let Yn∼U(−1/n,1/n)Y_n \sim U(-1/n,1/n)Yn∼U(−1/n,1/n), Zn∼U(n,n+1 ...
- UA MATH564 概率论I 离散型随机变量
UA MATH564 概率论I 离散型随机变量 古典概型 排列数与组合数 二项式定理,负二项式定理与多项式定理 第二类Stirling数 独立性与条件概率 全概率公式 贝叶斯公式 Bernoulli试 ...
- UA MATH564 概率论 多项分布
UA MATH564 概率论 多项分布 多项式系数与多项式定理 多项分布的定义 多项分布的性质 多项式系数与多项式定理 下面讨论一个组合问题,假设我们要把nnn个物体分成rrr组,每一组有n1,n2, ...
- UA MATH564 概率论 多元随机变量的变换 理论与应用1
UA MATH564 概率论 多元随机变量的变换 理论与应用1 定义法 Jacobi行列式法 假设X1,⋯,XnX_1,\cdots,X_nX1,⋯,Xn是nnn个一元随机变量,存在一个C1(Rn ...
最新文章
- Rust语言开发基础(六)基础语法
- html调用百度地图语音播报,实现百度地图导航演示的语音播放功能
- [知识库分享系列] 三、Web(高性能Web站点建设)
- ASP.NET Web API 基本操作(CRUD)
- 音视频编解码的国际标准
- python中内置的数据结构有几种?
- PID控制器开发笔记之八:带死区的PID控制器的实现
- 整个互联网真的是呈现出一种勃勃的生机
- 调试器工作原理--CPU软件断点/硬件断点/单步执行标识
- 转件工程--实践者的研究方法阅读笔记1
- Keras网络层之“关于Keras的层(Layer)”
- JS getElementsByClassName
- angular解决跨域问题
- 快速下载720云高清全景图片
- 如何做好信息化和数字化建设,看这一篇就够了--童亚斋
- 九宫格拼图,避免不可还原
- 优化iphone存储空间,以及优化iphone内存。
- 模型解释器-LIME
- 微信文件夹的dat文件怎么打开_微信dat文件如何打开?
- CS01/CS02/CS03 BOM修改 CSAP_MAT_BOM_MAINTAIN
热门文章
- (转载)数据库范式及宽表窄表理解
- Leetcode 231. 2的幂 解题思路及C++实现
- Linux I2C子系统分析-I2C设备驱动
- this和self区别
- 软件设计原则——合成复用原则
- 案例:用户信息列表展示||1. 需求 2. 设计 3. 开发4. 测试 5. 部署运维
- 计算机组成原理——总线结构
- java vector arraylist linkedlist用法与区别
- Linux个人用户安装JDK(与root用户不同版本)
- [YTU]_2907(类重载实现矩阵加法)