条件概率描述的是:已经知道一件事发生之后,另一件事的可能性怎样变化。统计、金融、医学里到处都是这种“带着已知信息再估概率”的问题。

本章先给出条件概率的定义和乘法公式,再讲独立性、贝叶斯定理,以及把条件概率本身当成一种新的概率测度来用。

条件概率

条件概率就是给原有概率再加约束。直观上,这相当于把样本空间压缩到已知事件里,再在这块更小的空间上重新度量。

基本条件概率

定义条件概率

设 AA、BB 是同一概率空间中的事件,且 P(B)>0P(B)>0。事件 AA 在条件 BB 下的概率定义为

P(A∣B)=P(A∩B)P(B).P(A \mid B) = \frac{P(A \cap B)}{P(B)}.
注定义的前提

P(B)=0P(B)=0 时这个比值没有定义。直观上,P(A∣B)P(A\mid B) 就是在得知 BB 已经发生之后,对 AA 的可能性所作的更新。

证明

需要说明 P(E∣F)=P(EF)P(F)P(E\mid F)=\frac{P(EF)}{P(F)} 与概率公理相容。

设样本空间为 Ω\Omega,F⊆ΩF\subseteq\Omega 且 P(F)>0P(F)>0。EFEF 表示 EE 与 FF 同时发生。

在频率解释下,重复大量试验后,P(F)P(F) 是 FF 出现的比例,P(EF)P(EF) 是两者同时出现的比例。一旦已知 FF 发生,新的样本空间就限制在 FF 上。此时 EE 的概率,就是 FF 里同时属于 EE 的那一部分所占的比例,也就是 EFEF 在 FF 中的份额:

P(E∣F)=P(EF)P(F).P(E \mid F) = \frac{P(EF)}{P(F)}.

因此这个定义是把原概率测度限制到条件空间 FF 上之后的自然延拓。

下面用掷骰子把这件事写清楚。

例两次掷骰,已知第一枚是 4

公平六面骰连掷两次。求:在第一枚已经是 4 的条件下,两枚点数之和大于 8 的概率。

两次掷骰的样本空间可用格子图表示:每个点 (x,y)(x,y) 对应第一枚为 xx、第二枚为 yy,一共 3636 个等可能结果。红点是点数之和大于 88 的结果。

两次掷骰的格子图,标出点数之和大于 8 的结果。

记 AA 为“点数之和大于 88”,BB 为“第一枚是 44”。条件 BB 把注意力限制在格子图的第四列。A∩BA\cap B 里只剩 (4,5)(4,5) 和 (4,6)(4,6)。

第一枚为 4 的那一列。

P(A∩B)=236=118,P(B)=636=16.P(A\cap B)=\frac{2}{36}=\frac{1}{18}, \qquad P(B)=\frac{6}{36}=\frac{1}{6}.

因此

P(A∣B)=P(A∩B)P(B)=1/181/6=13.P(A\mid B) = \frac{P(A\cap B)}{P(B)} = \frac{1/18}{1/6} = \frac{1}{3}.

在子空间 BB 里直接数,也是 2/6=1/32/6=1/3。两种算法等价:一种用原来的测度,一种直接在压缩后的样本空间上数。

注等可能只是特例

公平骰子每个面概率相同。现实里多数事件并不等可能。要处理有偏的情形,需要先把定义改写成乘法公式。

推论乘法公式

两边同乘 P(F)P(F),得到

P(EF)=P(E∣F)P(F).P(EF)=P(E\mid F)P(F).
例有偏四面骰掷两次

一枚有偏的四面骰(点数 11 到 44)连掷两次,各面概率为

  • P(1)=0.1P(1)=0.1
  • P(2)=0.2P(2)=0.2
  • P(3)=0.3P(3)=0.3
  • P(4)=0.4P(4)=0.4

求:

  1. 已知第一枚是 22,点数之和为 44 的条件概率。
  2. 已知点数之和为 44,第一枚是 22 的条件概率。
解有偏四面骰

记 AA 为“两点数之和为 44”,BB 为“第一枚是 22”。用树图标出每条分支的概率会更清楚。

有偏四面骰掷两次的树图。

1. P(A∣B)P(A\mid B)。

条件 BB 把第一枚钉在 22。要让和为 44,第二枚必须也是 22。因此

P(A∣B)=P(第二枚是 2)=0.2.P(A\mid B)=P(\text{第二枚是 }2)=0.2.

用定义验算也一样:

P(A∣B)=0.2×0.20.2=0.2.P(A\mid B) = \frac{0.2\times 0.2}{0.2} = 0.2.

这里条件并不改变第二枚仍是 22 的机会,因为两次投掷彼此独立。

2. P(B∣A)P(B\mid A)。

和为 44 的组合有 (1,3)(1,3)、(3,1)(3,1)、(2,2)(2,2)。于是

P(A)=P(1,3)+P(3,1)+P(2,2)=0.1⋅0.3+0.3⋅0.1+0.2⋅0.2=0.10.P(A) = P(1,3)+P(3,1)+P(2,2) = 0.1\cdot 0.3+0.3\cdot 0.1+0.2\cdot 0.2 = 0.10.

树图每条边上的权重都是条件概率。各分支互不相交,所以用可加性把它们加起来。于是

P(B∣A)=P(B∩A)P(A)=P(2,2)P(A)=0.040.10=0.4.P(B\mid A) = \frac{P(B\cap A)}{P(A)} = \frac{P(2,2)}{P(A)} = \frac{0.04}{0.10} = 0.4.

概率不再等可能,做法却没有变:先把样本空间限制到已知事件上,再在这块子空间里定位目标。

格子图适合结果等可能、阶段清楚的实验,例如掷硬币、掷骰子。树图更适合各分支概率不同、过程本身有层次的情形,贝叶斯更新尤其常用。问题规模一大,图就会画不下去,最终还是要靠事件之间的关系来算。

例抛硬币选课

Celine 在法语课和化学课之间犹豫。她估计法语课拿 A 的概率是 1/21/2,化学课拿 A 的概率是 2/32/3。她决定用一枚公平硬币来选课。问:她选化学课并且拿 A 的概率是多少?

解抛硬币选课

记 CC 为选化学课,AA 为不论哪门课都拿到 A。硬币公平,所以 P(C)=1/2P(C)=1/2。选了化学课后拿 A 的概率是 P(A∣C)=2/3P(A\mid C)=2/3。于是

P(CA)=P(C)P(A∣C)=12⋅23=13.P(CA)=P(C)P(A\mid C)=\frac12\cdot\frac23=\frac13.
注交事件不是条件概率

有人会把“化学课拿 A”听成条件概率。硬币结果并不会改变化学课拿 A 的机会,这里求的是交事件,不是“给定硬币朝向之后的条件概率”。符号语言比口语精确,口语却更容易带进这种混淆。独立性一节会把这件事说透。

例罐中取球

罐中有 88 个红球、44 个白球,不放回地连取两个。

(a) 每次抽取时罐中各球等可能,两个都是红球的概率是多少?

(b) 现在设红球每个重量为 rr,白球每个重量为 ww,下一颗被抽到的概率等于它的重量除以当时罐中所有球的总重量。两个都是红球的概率又是多少?

解罐中取球

(a) 记 R1R_1、R2R_2 分别为第一次、第二次抽到红球。第一次抽到红球后,罐里剩 77 红 44 白,所以 P(R2∣R1)=7/11P(R_2\mid R_1)=7/11。又 P(R1)=8/12=2/3P(R_1)=8/12=2/3,于是

P(R1R2)=P(R1)P(R2∣R1)=23⋅711=1433.P(R_1 R_2)=P(R_1)P(R_2\mid R_1)=\frac{2}{3}\cdot\frac{7}{11}=\frac{14}{33}.

(b) 仍用乘法公式。给红球编号,记 BiB_i 为第一次抽到第 ii 个红球。则

P(R1)=P(⋃i=18Bi)=∑i=18P(Bi)=8r8r+4w.P(R_1) = P\Bigl(\bigcup_{i=1}^8 B_i\Bigr) = \sum_{i=1}^8 P(B_i) = \frac{8r}{8r+4w}.

第一次已经抽走一个红球后,罐里剩 77 红 44 白,同理

P(R2∣R1)=7r7r+4w.P(R_2\mid R_1)=\frac{7r}{7r+4w}.

因此

P(R1R2)=8r8r+4w⋅7r7r+4w=8r⋅7r(8r+4w)(7r+4w).P(R_1 R_2) = \frac{8r}{8r+4w}\cdot\frac{7r}{7r+4w} = \frac{8r\cdot 7r}{(8r+4w)(7r+4w)}.

把乘法公式连用若干次,就得到一般的链式法则。

定理概率的乘法法则

设 E1,E2,…,EnE_1,E_2,\ldots,E_n 为事件。它们依次全部发生的概率为

P(E1E2⋯En)=P(E1) P(E2∣E1) P(E3∣E1E2)⋯P(En∣E1E2⋯En−1).P(E_1 E_2 \cdots E_n) = P(E_1)\, P(E_2\mid E_1)\, P(E_3\mid E_1 E_2) \cdots P(E_n\mid E_1 E_2 \cdots E_{n-1}).
证明

对事件个数 nn 作归纳。

当 n=2n=2 时,乘法法则就是条件概率的定义:

P(E1E2)=P(E1)P(E2∣E1).P(E_1 E_2)=P(E_1)P(E_2\mid E_1).

假设对 n−1n-1 个事件成立:

P(E1⋯En−1)=P(E1)P(E2∣E1)⋯P(En−1∣E1⋯En−2).P(E_1\cdots E_{n-1}) = P(E_1)P(E_2\mid E_1)\cdots P(E_{n-1}\mid E_1\cdots E_{n-2}).

再由条件概率的定义,

P(E1⋯En)=P(E1⋯En−1) P(En∣E1⋯En−1).P(E_1\cdots E_n) = P(E_1\cdots E_{n-1})\, P(E_n\mid E_1\cdots E_{n-1}).

代入归纳假设即得 nn 个事件的公式。

注链式约分

也可以像微分里的链式法则那样直接约分:

P(E1)⋅P(E1E2)P(E1)⋅P(E1E2E3)P(E1E2)⋯P(E1⋯En)P(E1⋯En−1)=P(E1⋯En).P(E_1) \cdot \frac{P(E_1 E_2)}{P(E_1)} \cdot \frac{P(E_1 E_2 E_3)}{P(E_1 E_2)} \cdots \frac{P(E_1\cdots E_n)}{P(E_1\cdots E_{n-1})} = P(E_1\cdots E_n).
例四堆牌各一张 A

一副 5252 张扑克随机分成四堆,每堆 1313 张。求每堆恰好一张 A 的概率。

解四堆牌各一张 A

定义事件:

  • E1E_1:黑桃 A 落在某一堆(这其实就是整个样本空间)
  • E2E_2:黑桃 A 与红心 A 不在同一堆
  • E3E_3:黑桃、红心、方块三张 A 分属三堆
  • E4E_4:四张 A 分属四堆

所求即 P(E1E2E3E4)P(E_1 E_2 E_3 E_4)。由乘法法则,

P(E1E2E3E4)=P(E1) P(E2∣E1) P(E3∣E1E2) P(E4∣E1E2E3).P(E_1 E_2 E_3 E_4) = P(E_1)\, P(E_2\mid E_1)\, P(E_3\mid E_1 E_2)\, P(E_4\mid E_1 E_2 E_3).

P(E1)=1P(E_1)=1。含黑桃 A 的那一堆还要再收 1212 张,这 1212 张等可能地来自剩下 5151 张,红心 A 落在其中的概率是 12/5112/51,所以

P(E2∣E1)=1−1251=3951.P(E_2\mid E_1)=1-\frac{12}{51}=\frac{39}{51}.

已知两张 A 已分属两堆后,这两堆还剩 2424 个空位,等可能地来自剩下 5050 张。方块 A 落在这 2424 张里的概率是 24/5024/50,所以

P(E3∣E1E2)=1−2450=2650.P(E_3\mid E_1 E_2)=1-\frac{24}{50}=\frac{26}{50}.

同理

P(E4∣E1E2E3)=1−3649=1349.P(E_4\mid E_1 E_2 E_3)=1-\frac{36}{49}=\frac{13}{49}.

于是

P(E1E2E3E4)=3951⋅2650⋅1349≈0.105.P(E_1 E_2 E_3 E_4) = \frac{39}{51}\cdot\frac{26}{50}\cdot\frac{13}{49} \approx 0.105.

大约有 10.5%10.5\% 的机会每堆恰好一张 A。

习题

锁定源文件在这一节只有标题,没有正文。

贝叶斯定理

贝叶斯定理与贝叶斯式思考

对任意事件 EE、FF,总有

E=EF∪EFc.E = EF \cup EF^c.

EFEF 与 EFcEF^c 互不相交,由可加性

P(E)=P(EF)+P(EFc)=P(E∣F)P(F)+P(E∣Fc)P(Fc)=P(E∣F)P(F)+P(E∣Fc)(1−P(F)).\begin{aligned} P(E) &= P(EF)+P(EF^c) \\ &= P(E\mid F)P(F)+P(E\mid F^c)P(F^c) \\ &= P(E\mid F)P(F)+P(E\mid F^c)\bigl(1-P(F)\bigr). \end{aligned}

这就是:用 FF 发生与不发生这两种情形,给 EE 做一次加权平均。有时直接求 P(E)P(E) 很难,先对另一个事件及其补事件条件化反而容易。

例事故倾向与出险

保险公司把人分成两类:容易出事故的,和不容易出事故的。一年内出事故的概率分别是 0.40.4 和 0.20.2。假设人口中有 30%30\% 属于前一类。新投保人在买保单后一年内出事故的概率是多少?

解事故倾向与出险

记 A1A_1 为一年内出事故,AA 为投保人属于易出事故的一类。先对是否易出事故条件化:

P(A1)=P(A1∣A)P(A)+P(A1∣Ac)P(Ac)=(0.4)(0.3)+(0.2)(0.7)=0.26.P(A_1) = P(A_1\mid A)P(A)+P(A_1\mid A^c)P(A^c) = (0.4)(0.3)+(0.2)(0.7) = 0.26.

上一例用的是 AA 或 AcA^c 之下 A1A_1 的概率。反过来也可以求“已经出事故之后,属于易出事故一类”的概率。

例出险后再判类型

新投保人在买保单后一年内出了事故。他属于易出事故一类的概率是多少?

解出险后再判类型
P(A∣A1)=P(A∩A1)P(A1)=P(A1∣A)P(A)P(A1)=(0.3)(0.4)0.26=613.P(A\mid A_1) = \frac{P(A\cap A_1)}{P(A_1)} = \frac{P(A_1\mid A)P(A)}{P(A_1)} = \frac{(0.3)(0.4)}{0.26} = \frac{6}{13}.

下图用 1,0001{,}000 名投保人的假想模型展示这一更新,其中 EE 就是上文的第一年出险事件 A1A_1。筛选后留下 120+140=260120+140=260 人,易出事故组的占比因此变成 120/260=6/13120/260=6/13。

贝叶斯条件筛选:300 名易出事故者与 700 名其他投保人,出险后分别留下 120 人与 140 人,后验概率为 6/13。

贝叶斯条件筛选:300 名易出事故者与 700 名其他投保人,出险后分别留下 120 人与 140 人,后验概率为 6/13。

已知 P(E∣F)P(E\mid F) 和 P(E∣Fc)P(E\mid F^c),就可以反求 P(F∣E)P(F\mid E)。这就是贝叶斯定理。它不只是代数变形,也刻画了一种更新信念的方式。

定理贝叶斯定理

设 P(B)>0P(B)>0。则

P(A∣B)=P(B∣A)P(A)P(B),P(A\mid B)=\frac{P(B\mid A)P(A)}{P(B)},

其中:

  • P(A∣B)P(A\mid B) 是后验:已经观察到 BB 之后,AA 的概率
  • P(B∣A)P(B\mid A) 是似然:假定 AA 发生时,观察到 BB 的概率
  • P(A)P(A) 是先验:AA 原先的概率
  • P(B)P(B) 是边缘概率。若 {Ai}\{A_i\} 构成完备组,则
P(B)=∑iP(B∣Ai)P(Ai).P(B)=\sum_i P(B\mid A_i)P(A_i).
证明

还需 P(A)>0P(A)>0,使公式里两个条件概率都有定义。由定义,P(A∣B)P(B)=P(A∩B)=P(B∣A)P(A)P(A\mid B)P(B)=P(A\cap B)=P(B\mid A)P(A);再除以 P(B)>0P(B)>0 即得。若分母用全概率公式展开,AiA_i 必须组成有限或可数个互不相交事件的划分,零概率块应在写条件概率前略去。P(A)=0P(A)=0 时,由单调性直接得到后验为零;但在这里的初等事件定义下,P(B∣A)P(B\mid A) 没有定义。

公式很短,却把“先有信念,再拿新证据去修正”这件事写死了。

  • 先验是已有的看法。
  • 似然是:在某个假定之下,眼前这件事有多像会发生。
  • 边缘概率负责把所有可能的解释归一化。
  • 后验是修正之后的新概率。

地平说变成圆球说,就是先验被一串观测一点点压下去的过程:P(地平)P(\text{地平}) 趋于 00,P(圆球)P(\text{圆球}) 趋于 11,给定全部新观测之后,地平说的后验几乎变成 00。

上面例子里 AA 只有“发生或不发生”两种。一般情形要把样本空间剖成任意多个块。先写出全概率公式。

定理全概率公式

设 B1,B2,…,BnB_1,B_2,\ldots,B_n 是样本空间 SS 的有限或可数剖分:当 i≠ji\neq j 时 Bi∩Bj=∅B_i\cap B_j=\emptyset,且 ⋃i=1nBi=S\bigcup_{i=1}^n B_i=S,并且每个 P(Bi)>0P(B_i)>0。则对 SS 中任意事件 AA,

P(A)=∑i=1nP(A∣Bi)P(Bi).P(A)=\sum_{i=1}^n P(A\mid B_i)P(B_i).
证明

B1,…,BnB_1,\ldots,B_n 互不相交且并起来是整个 SS,所以

A=A∩S=A∩(⋃i=1nBi)=⋃i=1n(A∩Bi).A = A\cap S = A\cap\Bigl(\bigcup_{i=1}^n B_i\Bigr) = \bigcup_{i=1}^n (A\cap B_i).

诸 A∩BiA\cap B_i 仍互不相交,于是

P(A)=∑i=1nP(A∩Bi)=∑i=1nP(A∣Bi)P(Bi).P(A) = \sum_{i=1}^n P(A\cap B_i) = \sum_{i=1}^n P(A\mid B_i)P(B_i).
注从二元情形归纳

也可以用本节开头 n=2n=2 的分解作归纳基础。

证明

贝叶斯定理的证明。 设 {Bi}\{B_i\} 是样本空间的剖分,AA 为任意事件。由全概率公式

P(A)=∑i=1nP(A∩Bi)=∑i=1nP(A∣Bi)P(Bi).P(A)=\sum_{i=1}^n P(A\cap B_i)=\sum_{i=1}^n P(A\mid B_i)P(B_i).

对某个固定的 jj,

P(Bj∣A)=P(A∩Bj)P(A)=P(A∣Bj)P(Bj)∑i=1nP(A∣Bi)P(Bi).P(B_j\mid A) = \frac{P(A\cap B_j)}{P(A)} = \frac{P(A\mid B_j)P(B_j)}{\sum_{i=1}^n P(A\mid B_i)P(B_i)}.

若剖分只有 BB 与 BcB^c 两块,则简化为

P(B∣A)=P(A∣B)P(B)P(A∣B)P(B)+P(A∣Bc)P(Bc).P(B\mid A) = \frac{P(A\mid B)P(B)}{P(A\mid B)P(B)+P(A\mid B^c)P(B^c)}.
例嫌疑与特征

调查进行到某一步,探长有 60%60\% 把握认定某嫌疑人有罪。新证据表明罪犯具有某项特征(例如左撇子、秃顶、棕发),人群中有 20%20\% 具有该特征。若发现嫌疑人恰好具有该特征,探长现在应该有多大把握认定他有罪?

解嫌疑与特征

记 GG 为嫌疑人有罪,CC 为他具有该项特征。由贝叶斯定理,

P(G∣C)=P(C∣G)P(G)P(C).P(G\mid C)=\frac{P(C\mid G)P(G)}{P(C)}.

分母用全概率公式:

P(C)=P(C∣G)P(G)+P(C∣Gc)P(Gc).P(C)=P(C\mid G)P(G)+P(C\mid G^c)P(G^c).

代入 P(G)=0.60P(G)=0.60、P(Gc)=0.40P(G^c)=0.40,并取 P(C∣G)=1P(C\mid G)=1(有罪则必有该特征)、P(C∣Gc)=0.20P(C\mid G^c)=0.20,得到

P(C)=(1.0)(0.60)+(0.20)(0.40)=0.68,P(C)=(1.0)(0.60)+(0.20)(0.40)=0.68,

从而

P(G∣C)=(1.0)(0.60)0.68≈0.882.P(G\mid C)=\frac{(1.0)(0.60)}{0.68}\approx 0.882.

观察到该特征之后,把握大约升到 88.2%88.2\%。

例红黑卡片

三张卡片外形相同:第一张两面都红,第二张两面都黑,第三张一面红一面黑。混在帽子里随机抽出一张扣在地上。若朝上的一面是红色,另一面也是红色的概率是多少?

解红黑卡片

记 RRRR、BBBB、RBRB 分别为抽到全红、全黑、红黑卡,记 RR 为朝上一面为红。所求是朝下一面也红的概率。先算

P(R)=P(R∣RR)P(RR)+P(R∣RB)P(RB)+P(R∣BB)P(BB)=1⋅13+12⋅13+0⋅13=12.\begin{aligned} P(R) &= P(R\mid RR)P(RR)+P(R\mid RB)P(RB)+P(R\mid BB)P(BB) \\ &= 1\cdot\frac13+\frac12\cdot\frac13+0\cdot\frac13 = \frac12. \end{aligned}

全红卡朝上为红的概率是 11 而不是 1/21/2,因为它两面都红。朝下也红,当且仅当抽到的是全红卡,且朝上已经是红,于是

P(朝下为红∩R)=P(RR∩R)=P(RR)=13,P(\text{朝下为红}\cap R)=P(RR\cap R)=P(RR)=\frac13,

从而

P(朝下为红∣R)=1/31/2=23.P(\text{朝下为红}\mid R) = \frac{1/3}{1/2} = \frac23.

另一种数法更干净:把每张卡看成两面,六个等可能结果 R1,R2,B1,B2,R3,B3R_1,R_2,B_1,B_2,R_3,B_3。朝上为红的只有 R1,R2,R3R_1,R_2,R_3,其中朝下仍红的是 R1R_1 和 R2R_2,所以概率是 2/32/3。

有人会猜 1/21/2:看见红面之后,似乎只剩“全红卡”和“红黑卡”两种可能。错在把这两种可能当成等可能。全红卡贡献了两张红面,红黑卡只贡献一张。

应用条件概率前,需要明确观测是怎样抽取的。下一题的答案取决于母亲如何选择同行的孩子。

例带着一个孩子散步

一对夫妇有两个孩子。母亲带着其中一个孩子散步,若这个孩子是女孩,两个孩子都是女孩的概率是多少?

解带着一个孩子散步

假设两个孩子的性别独立且男女等可能,母亲也以相同概率选择任一孩子同行,选择与性别无关。以下结论基于这一抽样模型。记

  • G1G_1:年长的孩子是女孩
  • G2G_2:年幼的孩子是女孩
  • GG:被看到的那个孩子是女孩

再对称地定义 B1B_1、B2B_2、BB。所求为 P(G1G2∣G)P(G_1 G_2\mid G)。

P(G1G2∣G)=P(G1G2∩G)P(G).P(G_1 G_2\mid G) = \frac{P(G_1 G_2\cap G)}{P(G)}.

由全概率公式,四种性别组合等可能,各占 1/41/4,且

P(G∣G1G2)=1,P(G∣G1B2)=P(G∣B1G2)=12,P(G∣B1B2)=0,P(G\mid G_1 G_2)=1, \quad P(G\mid G_1 B_2)=P(G\mid B_1 G_2)=\tfrac12, \quad P(G\mid B_1 B_2)=0,

于是

P(G)=1⋅14+12⋅14+12⋅14+0⋅14=12.P(G) = 1\cdot\frac14+\frac12\cdot\frac14+\frac12\cdot\frac14+0\cdot\frac14 = \frac12.

又 P(G1G2∩G)=P(G1G2)=1/4P(G_1 G_2\cap G)=P(G_1 G_2)=1/4,所以

P(G1G2∣G)=1/41/2=12.P(G_1 G_2\mid G)=\frac{1/4}{1/2}=\frac12.

看见的是女孩时,两个都是女孩的概率是 1/21/2。

习题

贝叶斯定理的典型计算已经包含在保险、嫌疑特征、红黑卡片和两个孩子这几例中。独立性一节会再练“信息到底有没有改变概率”。

事件的独立性

独立性的定义

前面把 EE 在 FF 下的条件概率写成

P(E∣F)=P(EF)P(F).P(E\mid F)=\frac{P(EF)}{P(F)}.

若干例子里会出现 P(E∣F)=P(E)P(E\mid F)=P(E),这等价于 P(EF)=P(E)P(F)P(EF)=P(E)P(F)。这时称 EE 与 FF 独立。

定义两个事件独立

事件 EE 与 FF 独立,当且仅当

P(EF)=P(E)P(F).P(EF)=P(E)P(F).

当 P(F)>0P(F)>0 时,这等价于 P(E∣F)=P(E)P(E\mid F)=P(E)。

例点数之和为 6 与第一枚是 4

掷两枚公平骰子。记

  • E1E_1:点数之和为 66
  • FF:第一枚是 44

二者独立吗?

解点数之和为 6 与第一枚是 4

E1E_1 的有利结果是 (1,5),(2,4),(3,3),(4,2),(5,1)(1,5),(2,4),(3,3),(4,2),(5,1),共 55 种,所以 P(E1)=5/36P(E_1)=5/36。第一枚是 44 的概率 P(F)=1/6P(F)=1/6。同时发生只剩 (4,2)(4,2),所以 P(E1∩F)=1/36P(E_1\cap F)=1/36。但

P(E1)P(F)=536⋅16=5216≠136.P(E_1)P(F)=\frac{5}{36}\cdot\frac{1}{6}=\frac{5}{216}\neq\frac{1}{36}.

因此 E1E_1 与 FF 不独立。

独立还会传到补事件上。

命题独立蕴含与补事件独立

若 EE 与 FF 独立,则 EE 与 FcF^c 也独立。

证明

E=EF∪EFcE=EF\cup EF^c,且这两块不相交,所以

P(E)=P(EF)+P(EFc).P(E)=P(EF)+P(EF^c).

由独立性 P(EF)=P(E)P(F)P(EF)=P(E)P(F),于是

P(EFc)=P(E)−P(E)P(F)=P(E)(1−P(F))=P(E)P(Fc).P(EF^c)=P(E)-P(E)P(F)=P(E)(1-P(F))=P(E)P(F^c).

独立和互斥不是一回事。若 AA 与 BB 互斥且都有正概率,则 P(A∩B)=0P(A\cap B)=0,而 P(A)P(B)>0P(A)P(B)>0,二者不可能独立。

多个事件的独立性

两个事件独立,三个事件未必整体独立。

例和为 7 对单枚独立、对同时发生不独立

掷两枚公平骰子。记 EE 为点数之和为 77,FF 为第一枚是 44,GG 为第二枚是 33。EE 对 FF 独立,对 GG 也独立,但对交事件 FGFG 并不独立,因为 P(E∣FG)=1P(E\mid FG)=1。

解和为 7 对单枚独立、对同时发生不独立

P(E)=6/36=1/6P(E)=6/36=1/6,P(F)=P(G)=1/6P(F)=P(G)=1/6,P(E∩F)=P(4,3)=1/36P(E\cap F)=P(4,3)=1/36,于是

P(E∩F)=P(E)P(F)=136.P(E\cap F)=P(E)P(F)=\frac{1}{36}.

EE 与 GG 同理。但 FGFG 就是结果 (4,3)(4,3),此时和一定是 77,所以 P(E∩FG)=P(FG)=1/36P(E\cap FG)=P(FG)=1/36,而

P(E)P(FG)=16⋅136=1216≠136.P(E)P(FG)=\frac16\cdot\frac{1}{36}=\frac{1}{216}\neq\frac{1}{36}.

对 FF、GG 分别独立,并不保证对它们拼成的事件也独立。

定义三个事件相互独立

三个事件 EE、FF、GG 相互独立,是指下列各式同时成立:

P(EFG)=P(E)P(F)P(G),P(EF)=P(E)P(F),P(EG)=P(E)P(G),P(FG)=P(F)P(G).\begin{aligned} P(EFG)&=P(E)P(F)P(G), \\ P(EF)&=P(E)P(F), \\ P(EG)&=P(E)P(G), \\ P(FG)&=P(F)P(G). \end{aligned}

若三者相互独立,则 EE 与由 FF、GG 做成的任何事件都独立。例如对并:

P(E(F∪G))=P(EF∪EG)=P(EF)+P(EG)−P(EFG)=P(E)P(F)+P(E)P(G)−P(E)P(FG)=P(E) P(F∪G).\begin{aligned} P\bigl(E(F\cup G)\bigr) &= P(EF\cup EG) \\ &= P(EF)+P(EG)-P(EFG) \\ &= P(E)P(F)+P(E)P(G)-P(E)P(FG) \\ &= P(E)\,P(F\cup G). \end{aligned}

这个定义可以直接推广到任意多个、甚至无穷多个事件。

例独立试验中的成功次数

一串独立试验,每次成功概率为 pp,失败概率为 1−p1-p。求:

  1. 前 nn 次里至少成功一次的概率
  2. 前 nn 次里恰好成功 kk 次的概率
  3. 所有试验都成功的概率
解独立试验中的成功次数

(a) 先算补事件:前 nn 次全失败。记 EiE_i 为第 ii 次失败,由独立性

P(E1⋯En)=(1−p)n,P(E_1\cdots E_n)=(1-p)^n,

所以至少成功一次的概率是 1−(1−p)n1-(1-p)^n。

(b) 任意一个“恰好 kk 次成功、n−kn-k 次失败”的具体序列,概率都是 pk(1−p)n−kp^k(1-p)^{n-k},这样的序列有 (nk)\binom{n}{k} 个,于是

P(恰好 k 次成功)=(nk)pk(1−p)n−k.P(\text{恰好 }k\text{ 次成功})=\binom{n}{k}p^k(1-p)^{n-k}.

(c) 本小问改记 EiE_i 为第 ii 次成功。前 nn 次全成功的概率是 pnp^n。无穷多次全成功是这些事件的交。由概率的连续性,

P(⋂i=1∞Ei)=lim⁡n→∞pn,P\Bigl(\bigcap_{i=1}^\infty E_i\Bigr) = \lim_{n\to\infty}p^n,

当 p<1p<1 时极限为 00,当 p=1p=1 时为 11。

例5 先于 7 出现

反复独立地掷一对公平骰子。点数之和为 55 先于之和为 77 出现的概率是多少?

解5 先于 7 出现

记 EnE_n 为:前 n−1n-1 次既没有 55 也没有 77,第 nn 次出现 55。所求为

P(⋃n=1∞En)=∑n=1∞P(En).P\Bigl(\bigcup_{n=1}^\infty E_n\Bigr)=\sum_{n=1}^\infty P(E_n).

一次试验得到 55 的概率是 4/364/36,得到 77 的概率是 6/366/36。由独立性,

P(En)=(1−1036)n−1⋅436,P(E_n)=\Bigl(1-\frac{10}{36}\Bigr)^{n-1}\cdot\frac{4}{36},

于是

∑n=1∞P(En)=19∑n=1∞(1318)n−1=19⋅11−13/18=25.\sum_{n=1}^\infty P(E_n) = \frac19\sum_{n=1}^\infty\Bigl(\frac{13}{18}\Bigr)^{n-1} = \frac19\cdot\frac{1}{1-13/18} = \frac25.

也可以对第一次的结果条件化。记 FF 为第一次得 55,GG 为第一次得 77,HH 为第一次两者都不是。则 P(E∣F)=1P(E\mid F)=1,P(E∣G)=0P(E\mid G)=0,P(E∣H)=P(E)P(E\mid H)=P(E),从而

P(E)=19+1318P(E),P(E)=\frac19+\frac{13}{18}P(E),

解得 P(E)=2/5P(E)=2/5。

注互斥事件谁先出现

答案符合直觉:55 的单次概率是 4/364/36,77 是 6/366/36,所以 55 先出现的概率应是 4/104/10。一般地,若单次试验中 EE、FF 互斥,独立重复该试验时,EE 先于 FF 出现的概率为

P(E)P(E)+P(F).\frac{P(E)}{P(E)+P(F)}.

习题

练习和为 7 与第一枚点数

公平骰子连掷两次。点数之和为 77 是否与第一枚的点数独立?证明或举反例,并解释原因。

证明

记 AA 为点数之和为 77,BiB_i 为第一枚是 ii。要独立,需对每个 ii 都有 P(A∩Bi)=P(A)P(Bi)P(A\cap B_i)=P(A)P(B_i)。

AA 的有利结果是 (1,6),(2,5),(3,4),(4,3),(5,2),(6,1)(1,6),(2,5),(3,4),(4,3),(5,2),(6,1),共 66 种,所以 P(A)=6/36=1/6P(A)=6/36=1/6。对每个 ii,P(Bi)=1/6P(B_i)=1/6。A∩BiA\cap B_i 只对应第一枚为 ii、第二枚为 7−i7-i 这一种结果,所以 P(A∩Bi)=1/36P(A\cap B_i)=1/36。于是

P(A)P(Bi)=16⋅16=136=P(A∩Bi).P(A)P(B_i)=\frac16\cdot\frac16=\frac{1}{36}=P(A\cap B_i).

因此 AA 与每个 BiB_i 都独立:知道第一枚是几,并不能改变“和为 77”的概率,因为第二枚总有恰好一种补法。

练习相互独立的一般定义

证明相互独立的一般定义,并说明它对布尔组合仍然成立。

定理相互独立

事件族 {E1,…,En}\{E_1,\ldots,E_n\} 相互独立,是指对任意子集 {Ei1,…,Eik}\{E_{i_1},\ldots,E_{i_k}\},

P(⋂j=1kEij)=∏j=1kP(Eij).P\Bigl(\bigcap_{j=1}^k E_{i_j}\Bigr) = \prod_{j=1}^k P(E_{i_j}).

这对每个 1≤k≤n1\le k\le n 以及每一组下标都要成立。

证明

设该族相互独立。任取 ii,以及由其余事件经布尔运算做成的事件 AA。AA 可以写成若干形如 ⋂(Ej 或 Ejc)\bigcap(E_j\text{ 或 }E_j^c) 的交之并。由容斥与相互独立,每个这样的交的概率都裂成各因子概率(或 1−P(Ej)1-P(E_j))的乘积。对 Ei∩AE_i\cap A 做同样分解,每个含 EiE_i 的交都会多出一个因子 P(Ei)P(E_i),于是 P(Ei∩A)=P(Ei)P(A)P(E_i\cap A)=P(E_i)P(A)。

进一步的条件概率

条件概率也满足概率公理

条件概率本身仍是概率,原先对 PP 成立的结论,对 P( ⋅∣F)P(\,\cdot\mid F) 也应当成立。

命题条件概率的公理

设 EE、FF 为事件,EiE_i(i=1,2,…i=1,2,\ldots)两两不交。则

  • 0≤P(E∣F)≤10\le P(E\mid F)\le 1
  • P(S∣F)=1P(S\mid F)=1,其中 SS 为样本空间
  • 若诸 EiE_i 两两不交,则
P(⋃i=1∞Ei∣F)=∑i=1∞P(Ei∣F).P\Bigl(\bigcup_{i=1}^\infty E_i\Bigm| F\Bigr) = \sum_{i=1}^\infty P(E_i\mid F).
证明

E∩F⊆FE\cap F\subseteq F,所以 P(E∩F)≤P(F)P(E\cap F)\le P(F),从而 0≤P(E∣F)≤10\le P(E\mid F)\le 1。又 S∩F=FS\cap F=F,所以 P(S∣F)=1P(S\mid F)=1。第三款来自可数可加性:

P(⋃i=1∞Ei∣F)=P(⋃i=1∞(Ei∩F))P(F)=∑i=1∞P(Ei∣F).P\Bigl(\bigcup_{i=1}^\infty E_i\Bigm| F\Bigr) = \frac{P\bigl(\bigcup_{i=1}^\infty(E_i\cap F)\bigr)}{P(F)} = \sum_{i=1}^\infty P(E_i\mid F).
命题条件化后仍是概率测度

定义 Q(E)=P(E∣F)Q(E)=P(E\mid F),则 QQ 是 SS 上的概率函数,先前对概率证明过的命题全部适用。例如

Q(E1∪E2)=Q(E1)+Q(E2)−Q(E1∩E2).Q(E_1\cup E_2)=Q(E_1)+Q(E_2)-Q(E_1\cap E_2).
证明

固定 P(F)>0P(F)>0。前面已证明非负性、归一化与可数可加性,这正是 QQ 的全部概率公理。对所写的容斥式,先对 E1∩FE_1\cap F 与 E2∩FE_2\cap F 使用普通容斥,再除以 P(F)P(F);两者的交为 E1∩E2∩FE_1\cap E_2\cap F,于是得到条件化后的公式。

注先做本节习题 1

这等于说:凡是已经对普通概率证明过的结论,都可以原样搬到条件概率上。下面的例子之前,建议先看本节习题 1。

例第一年出险后再估第二年

保险公司仍把投保人分成易出事故与不易出事故两类,一年内出事故的概率仍是 0.40.4 与 0.20.2,新投保人中易出事故的比例是 3/103/10。已知某人第一年出了事故,第二年再出事故的概率是多少?

解第一年出险后再估第二年

记 AA 为易出事故,A1A_1、A2A_2 分别为第一年、第二年出事故。假设投保人的类型不变,且给定类型后各年是否出事故条件独立;仅凭各年的边缘概率,还不足以确定答案。所求为 P(A2∣A1)P(A_2\mid A_1)。对是否易出事故再条件化:

P(A2∣A1)=P(A2∣A∩A1)P(A∣A1)+P(A2∣Ac∩A1)P(Ac∣A1).P(A_2\mid A_1) = P(A_2\mid A\cap A_1)P(A\mid A_1) + P(A_2\mid A^c\cap A_1)P(A^c\mid A_1).

先用贝叶斯定理算 P(A∣A1)P(A\mid A_1)。已知 P(A1∣A)=0.4P(A_1\mid A)=0.4,P(A)=0.3P(A)=0.3,且先前已得 P(A1)=0.26P(A_1)=0.26,于是

P(A∣A1)=(0.4)(0.3)0.26=613,P(Ac∣A1)=713.P(A\mid A_1)=\frac{(0.4)(0.3)}{0.26}=\frac{6}{13}, \qquad P(A^c\mid A_1)=\frac{7}{13}.

假定各年出事故在给定类型后条件独立,则 P(A2∣A∩A1)=0.4P(A_2\mid A\cap A_1)=0.4,P(A2∣Ac∩A1)=0.2P(A_2\mid A^c\cap A_1)=0.2,从而

P(A2∣A1)=25613+15713=1965≈0.2923.P(A_2\mid A_1) =\frac{2}{5}\frac{6}{13}+\frac{1}{5}\frac{7}{13} =\frac{19}{65} \approx 0.2923.
例相继观测下的后验更新

设 H1,…,HnH_1,\ldots,H_n 是互斥且穷尽的假设。以下假定已观测事件的概率为正;后验为零的假设可从条件表达式中略去。

先观测到 E1E_1,由贝叶斯定理得到

P(Hi∣E1)=P(E1∣Hi)P(Hi)P(E1).P(H_i\mid E_1)=\frac{P(E_1\mid H_i)P(H_i)}{P(E_1)}.

再观测到 E2E_2 时,一般的更新公式必须在新似然中保留 E1E_1:

P(Hi∣E1,E2)=P(E2∣Hi,E1)P(Hi∣E1)P(E2∣E1),P(H_i\mid E_1,E_2) =\frac{P(E_2\mid H_i,E_1)P(H_i\mid E_1)}{P(E_2\mid E_1)},

其中

P(E2∣E1)=∑jP(E2∣Hj,E1)P(Hj∣E1).P(E_2\mid E_1)=\sum_j P(E_2\mid H_j,E_1)P(H_j\mid E_1).

如果给定每个假设后 E1E_1 与 E2E_2 条件独立,新似然才可简化为 P(E2∣Hi)P(E_2\mid H_i),从而

P(Hi∣E1,E2)=P(E2∣Hi)P(E1∣Hi)P(Hi)P(E2∣E1)P(E1).P(H_i\mid E_1,E_2) =\frac{P(E_2\mid H_i)P(E_1\mid H_i)P(H_i)}{P(E_2\mid E_1)P(E_1)}.

分母是联合证据 P(E1∩E2)P(E_1\cap E_2),由两个概率相乘得到,而不是相除。例如,从若干枚已知偏率的硬币中选定一枚,但不知道选中了哪枚,然后反复抛掷这同一枚硬币。给定硬币身份后,各次抛掷独立,因此每次更新只需保留当前后验与新观测的似然,不必单独存储此前所有结果。

多重条件概率

多个条件可以先交成一个事件,再套用原来的定义。

定义多重条件概率

事件 AA 在条件 B1,…,BnB_1,\ldots,B_n 下的概率定义为

P(A∣B1,…,Bn)=P(A∩B1∩⋯∩Bn)P(B1∩⋯∩Bn),P(A\mid B_1,\ldots,B_n) = \frac{P(A\cap B_1\cap\cdots\cap B_n)}{P(B_1\cap\cdots\cap B_n)},

前提是分母为正。

注多个条件合成一个事件

令 B=B1∩⋯∩BnB=B_1\cap\cdots\cap B_n。只要 P(B)>0P(B)>0,直接应用条件概率的定义即可:

P(A∣B1,…,Bn)=P(A∣B)=P(A∩B)P(B).P(A\mid B_1,\ldots,B_n) =P(A\mid B) =\frac{P(A\cap B)}{P(B)}.

这里无需归纳证明,也不要求各个 BiB_i 互相独立。

例系统故障的多重条件

某计算机系统的故障可能同时受到软件故障 SS、硬件故障 HH、网络故障 NN、电涌 PP、误操作 UU 影响。已知

  • P(故障∣S∩H∩N∩P∩U)=0.90P(\text{故障}\mid S\cap H\cap N\cap P\cap U)=0.90
  • P(S∩H∩N∩P∩U)=0.01P(S\cap H\cap N\cap P\cap U)=0.01
  • P(故障∣S∩H∩N)=0.75P(\text{故障}\mid S\cap H\cap N)=0.75
  • P(S∩H∩N)=0.05P(S\cap H\cap N)=0.05
  • P(P∩U∣S∩H∩N)=0.20P(P\cap U\mid S\cap H\cap N)=0.20

把这些条件合在一起,系统故障的综合条件概率是多少?

解系统故障的多重条件

先核验联合概率:

P(S∩H∩N∩P∩U)=P(P∩U∣S∩H∩N) P(S∩H∩N)=0.20×0.05=0.01.P(S\cap H\cap N\cap P\cap U) = P(P\cap U\mid S\cap H\cap N)\,P(S\cap H\cap N) = 0.20\times 0.05 = 0.01.

再由定义,

P(故障∩S∩H∩N∩P∩U)=0.90×0.01=0.009,P(\text{故障}\cap S\cap H\cap N\cap P\cap U) = 0.90\times 0.01 = 0.009,

从而

P(故障∣S,H,N,P,U)=0.0090.01=0.90.P(\text{故障}\mid S,H,N,P,U) = \frac{0.009}{0.01} = 0.90.

这与直接给出的条件概率一致:多重条件并没有引入新的机制,只是把已经交在一起的事件当成单个条件来用。

习题

练习条件概率的容斥与全概率

证明把 Q(E)=P(E∣F)Q(E)=P(E\mid F) 看成概率测度后,

Q(E1∪E2)=Q(E1)+Q(E2)−Q(E1E2)Q(E_1\cup E_2)=Q(E_1)+Q(E_2)-Q(E_1 E_2)

成立,并且由此得到

P(E1∣F)=P(E1∣E2F)P(E2∣F)+P(E1∣E2cF)P(E2c∣F).P(E_1\mid F) = P(E_1\mid E_2 F)P(E_2\mid F) + P(E_1\mid E_2^c F)P(E_2^c\mid F).
证明

Q(S)=P(S∣F)=1Q(S)=P(S\mid F)=1;Q(E)≥0Q(E)\ge 0;对两两不交的序列,

Q(⋃i=1∞Ei)=∑i=1∞Q(Ei).Q\Bigl(\bigcup_{i=1}^\infty E_i\Bigr) = \sum_{i=1}^\infty Q(E_i).

因此 QQ 满足概率公理,两个事件的容斥公式对 QQ 成立。

再看第二式。先假设 Q(E2)>0Q(E_2)>0 且 Q(E2c)>0Q(E_2^c)>0,把 QQ 的全概率公式写出来:

Q(E1)=Q(E1∣E2)Q(E2)+Q(E1∣E2c)Q(E2c).Q(E_1)=Q(E_1\mid E_2)Q(E_2)+Q(E_1\mid E_2^c)Q(E_2^c).

而

Q(E1∣E2)=Q(E1E2)Q(E2)=P(E1E2∣F)P(E2∣F)=P(E1∣E2F),Q(E_1\mid E_2) = \frac{Q(E_1 E_2)}{Q(E_2)} = \frac{P(E_1 E_2\mid F)}{P(E_2\mid F)} = P(E_1\mid E_2 F),

同理 Q(E1∣E2c)=P(E1∣E2cF)Q(E_1\mid E_2^c)=P(E_1\mid E_2^c F)。代回即得。若某一分支在 QQ 下概率为零,省略该分支即可,无需对零概率事件定义条件概率。