本文是斯坦福大学CS229机器学习课程的基础材料,原始文件下载
原文作者:Arian Maleki , Tom Do
翻译:石振宇
审核和修改制作:黄海广
备注:请关注github的更新。
CS229 机器学习课程复习材料-概率论
文章目录
概率论复习和参考
概率论是对不确定性的研究。通过这门课,我们将依靠概率论中的概念来推导机器学习算法。这篇笔记试图涵盖适用于CS229的概率论基础。概率论的数学理论非常复杂,并且涉及到“分析”的一个分支:测度论。在这篇笔记中,我们提供了概率的一些基本处理方法,但是不会涉及到这些更复杂的细节。
1. 概率的基本要素
为了定义集合上的概率,我们需要一些基本元素,
-
样本空间$\Omega $:随机实验的所有结果的集合。在这里,每个结果 $w \in \Omega $ 可以被认为是实验结束时现实世界状态的完整描述。
-
事件集(事件空间) F \mathcal{F} F:元素 A ∈ F A \in \mathcal{F} A∈F 的集合(称为事件)是 $\Omega $ 的子集(即每个 A ⊆ Ω A \subseteq \Omega A⊆Ω 是一个实验可能结果的集合)。
备注: F \mathcal{F} F需要满足以下三个条件:
(1) ∅ ∈ F \emptyset \in \mathcal{F} ∅∈F
(2) A ∈ F ⟹ Ω \ A ∈ F A \in \mathcal{F} \Longrightarrow \Omega \backslash A \in \mathcal{F} A∈F⟹Ω\A∈F
(3) A 1 , A 2 , ⋯ A i ∈ F ⟹ ∪ i A i ∈ F A_1,A_2,\cdots A_{i} \in \mathcal{F}\Longrightarrow\cup_{i} A_{i} \in \mathcal{F} A1,A2,⋯Ai∈F⟹∪iAi∈F
-
概率度量 P P P:函数 P P P是一个$ \mathcal{F} \rightarrow \mathbb{R}$的映射,满足以下性质:
-
对于每个 A ∈ F A \in \mathcal{F} A∈F, P ( A ) ≥ 0 P(A) \geq 0 P(A)≥0,
-
P ( Ω ) = 1 P(\Omega) = 1 P(Ω)=1
-
如果 A 1 , A 2 , ⋯ A_1 ,A_2 ,\cdots A1,A2,⋯ 是互不相交的事件 (即 当$ i \neq j 时 , 时, 时,A_{i} \cap A_{j}=\emptyset$ ), 那么:
P ( ∪ i A i ) = ∑ i P ( A i ) P\left(\cup_{i} A_{i}\right)=\sum_{i} P\left(A_{i}\right) P(∪iAi)=i∑P(Ai)
以上三条性质被称为概率公理。
举例:
考虑投掷六面骰子的事件。样本空间为 Ω = { 1 , 2 , 3 , 4 , 5 , 6 } \Omega= \{1,2,3,4,5,6\} Ω={1,2,3,4,5,6}。最简单的事件空间是平凡事件空间 F = { ∅ , Ω } \mathcal{F} = \{\emptyset,\Omega\} F={∅,Ω}.另一个事件空间是 Ω \Omega Ω的所有子集的集合。对于第一个事件空间,满足上述要求的唯一概率度量由 P ( ∅ ) = 0 P(\emptyset) = 0 P(∅)=0, p ( Ω ) = 1 p(\Omega)= 1 p(Ω)=1给出。对于第二个事件空间,一个有效的概率度量是将事件空间中每个事件的概率分配为 i / 6 i/6 i/6,这里 i i i 是这个事件集合中元素的数量;例如 P ( { 1 , 2 , 3 , 4 } ) = 4 / 6 P(\{1,2,3,4\}) =4/6 P({1,2,3,4})=4/6, P ( { 1 , 2 , 3 } ) = 3 / 6 P(\{1,2,3\}) =3/6 P({1,2,3})=3/6。
性质:
- 如果 A ⊆ B A \subseteq B A⊆B,则:$ P(A) \leq P(B)$
- P ( A ∩ B ) ≤ m i n ( P ( A ) , P ( B ) ) P(A \cap B) \leq min(P(A),P(B) ) P(A∩B)≤min(P(A),P(B))
- (布尔不等式): P ( A ∪ B ) ≤ P ( A ) + P ( B ) P(A \cup B) \leq P(A)+P(B) P(A∪B)≤P(A)+P(B)
- P ( Ω ∣ A ) = 1 − P ( A ) P(\Omega |A ) =1-P(A) P(Ω∣A)=1−P(A)
- (全概率定律):如果 A 1 , ⋯ , A k A_1,\cdots,A_k A1,⋯,Ak是一些互不相交的事件并且它们的并集是 Ω \Omega Ω,那么它们的概率之和是1
1.1 条件概率和独立性
假设
B
B
B是一个概率非0的事件,我们定义在给定
B
B
B的条件下
A
A
A 的条件概率为:
P
(
A
∣
B
)
≜
P
(
A
∩
B
)
P
(
B
)
P(A | B) \triangleq \frac{P(A \cap B)}{P(B)}
P(A∣B)≜P(B)P(A∩B)
换句话说,
P
(
A
∣
B
P(A|B
P(A∣B)是度量已经观测到
B
B
B事件发生的情况下
A
A
A事件发生的概率,两个事件被称为独立事件当且仅当
P
(
A
∩
B
)
=
P
(
A
)
P
(
B
)
P(A \cap B) = P(A)P(B)
P(A∩B)=P(A)P(B)(或等价地,
P
(
A
∣
B
)
=
P
(
A
)
P(A|B) = P(A)
P(A∣B)=P(A))。因此,独立性相当于是说观察到事件
B
B
B对于事件
A
A
A的概率没有任何影响。
2. 随机变量
考虑一个实验,我们翻转10枚硬币,我们想知道正面硬币的数量。这里,样本空间 Ω \Omega Ω的元素是长度为10的序列。例如,我们可能有 w 0 = { H , H , T , H , T , H , H , T , T , T } ∈ Ω w_0 = \{H,H,T,H,T,H,H,T,T,T\}\in\Omega w0={H,H,T,H,T,H,H,T,T,T}∈Ω。然而,在实践中,我们通常不关心获得任何特定正反序列的概率。相反,我们通常关心结果的实值函数,比如我们10次投掷中出现的正面数,或者最长的背面长度。在某些技术条件下,这些函数被称为随机变量。
更正式地说,随机变量 X X X是一个的 Ω ⟶ R \Omega \longrightarrow \mathbb{R} Ω⟶R函数。通常,我们将使用大写字母 X ( ω ) X(\omega) X(ω)或更简单的 X X X(其中隐含对随机结果 ω \omega ω的依赖)来表示随机变量。我们将使用小写字母 x x x来表示随机变量的值。
举例:
在我们上面的实验中,假设
X
(
ω
)
X(\omega)
X(ω)是在投掷序列
ω
\omega
ω中出现的正面的数量。假设投掷的硬币只有10枚,那么
X
(
ω
)
X(\omega)
X(ω)只能取有限数量的值,因此它被称为离散随机变量。这里,与随机变量
X
X
X相关联的集合取某个特定值
k
k
k的概率为:
P
(
X
=
k
)
:
=
P
(
{
ω
:
X
(
ω
)
=
k
}
)
P(X=k) :=P(\{\omega : X(\omega) =k\})
P(X=k):=P({ω:X(ω)=k})
举例:
假设
X
(
ω
)
X(\omega)
X(ω)是一个随机变量,表示放射性粒子衰变所需的时间。在这种情况下,
X
(
ω
)
X(\omega)
X(ω)具有无限多的可能值,因此它被称为连续随机变量。我们将
X
X
X在两个实常数
a
a
a和
b
b
b之间取值的概率(其中
a
<
b
a < b
a<b)表示为:
P
(
a
≤
X
≤
b
)
:
=
P
(
{
ω
:
a
≤
X
(
ω
)
≤
b
}
)
P(a \leq X \leq b) :=P(\{\omega : a \leq X(\omega) \leq b\})
P(a≤X≤b):=P({ω:a≤X(ω)≤b})
2.1 累积分布函数
为了指定处理随机变量时使用的概率度量,通常可以方便地指定替代函数(CDF、PDF和PMF),在本节和接下来的两节中,我们将依次描述这些类型的函数。
累积分布函数(CDF)是函数
F
X
:
R
→
[
0
,
1
]
F_{X} : \mathbb{R} \rightarrow[0,1]
FX:R→[0,1],它将概率度量指定为:
F
X
(
x
)
≜
P
(
X
≤
x
)
F_{X}(x) \triangleq P(X \leq x)
FX(x)≜P(X≤x)
通过使用这个函数,我们可以计算任意事件发生的概率。图1显示了一个样本CDF函数。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Tk2AXLQe-1603691202074)(images/c89fd93c5d6dccce89762d57fcd66dac.png)]
- 0 ≤ F X ( x ) ≤ 1 0 \leq F_{X}(x)\leq 1 0≤FX(x)≤1
- lim x → − ∞ F X ( x ) = 0 \lim _{x \rightarrow-\infty} F_{X}(x)=0 limx→−∞FX(x)=0
- lim x → ∞ F X ( x ) = 1 \lim _{x \rightarrow\infty} F_{X}(x)=1 limx→∞FX(x)=1
- x ≤ y ⟹ F X ( x ) ≤ F X ( y ) x \leq y \Longrightarrow F_{X}(x)\leq F_{X}(y) x≤y⟹FX(x)≤FX(y)
2.2 概率质量函数
当随机变量
X
X
X取有限种可能值(即,
X
X
X是离散随机变量)时,表示与随机变量相关联的概率度量的更简单的方法是直接指定随机变量可以假设的每个值的概率。特别地,概率质量函数(PMF)是函数
p
X
:
Ω
→
R
p_{X} : \Omega \rightarrow \mathbb{R}
pX:Ω→R,这样:
p
X
(
x
)
≜
P
(
X
=
x
)
p_{X}(x) \triangleq P(X=x)
pX(x)≜P(X=x)
在离散随机变量的情况下,我们使用符号 V a l ( X ) Val(X) Val(X)表示随机变量 X X X可能假设的一组可能值。例如,如果 X ( ω ) X(\omega) X(ω)是一个随机变量,表示十次投掷硬币中的正面数,那么 V a l ( X ) = { 0 , 1 , 2 , . . . , 10 } Val(X) =\{0,1,2,...,10\} Val(X)={0,1,2,...,10}。
性质:
- 0 ≤ p X ( x ) ≤ 1 0 \leq p_{X}(x)\leq 1 0≤pX(x)≤1
- ∑ x ∈ V al ( X ) p X ( x ) = 1 \sum_{x \in V \text { al }(X)} p_{X}(x)=1 ∑x∈V al (X)pX(x)=1
- ∑ x ∈ A p X ( x ) = P ( X ∈ A ) \sum_{x \in A} p_{X}(x)=P(X \in A) ∑x∈ApX(x)=P(X∈A)
2.3 概率密度函数
对于一些连续随机变量,累积分布函数
F
X
(
x
)
F_X (x)
FX(x)处可微。在这些情况下,我们将**概率密度函数(PDF)**定义为累积分布函数的导数,即:
f
X
(
x
)
≜
d
F
X
(
x
)
d
x
f_{X}(x) \triangleq \frac{d F_{X}(x)}{d x}
fX(x)≜dxdFX(x)
请注意,连续随机变量的概率密度函数可能并不总是存在的(即,如果它不是处处可微)。
根据微分的性质,对于很小的
Δ
x
\Delta x
Δx,
P
(
x
≤
X
≤
x
+
Δ
x
)
≈
f
X
(
x
)
Δ
x
P(x \leq X \leq x+\Delta x) \approx f_{X}(x) \Delta x
P(x≤X≤x+Δx)≈fX(x)Δx
CDF和**PDF(当它们存在时!)都可用于计算不同事件的概率。但是应该强调的是,任意给定点的概率密度函数(PDF)**的值不是该事件的概率,即
f
X
(
x
)
≠
P
(
X
=
x
)
f _X (x) \not = P(X = x)
fX(x)=P(X=x)。例如,
f
X
(
x
)
f _X (x)
fX(x)可以取大于1的值(但是
f
X
(
x
)
f _X (x)
fX(x)在
R
\mathbb{R}
R的任何子集上的积分最多为1)。
性质:
- f X ( x ) ≥ 0 f_X(x)\geq 0 fX(x)≥0
- ∫ − ∞ ∞ f X ( x ) = 1 \int_{-\infty}^{\infty} f_{X}(x)=1 ∫−∞∞fX(x)=1
- ∫ x ∈ A f X ( x ) d x = P ( X ∈ A ) \int_{x \in A} f_{X}(x) d x=P(X \in A) ∫x∈AfX(x)dx=P(X∈A)
2.4 期望
假设
X
X
X是一个离散随机变量,其PMF为
p
X
(
x
)
p_X (x)
pX(x),
g
:
R
⟶
R
g : \mathbb{R} \longrightarrow \mathbb{R}
g:R⟶R是一个任意函数。在这种情况下,
g
(
X
)
g(X)
g(X)可以被视为随机变量,我们将
g
(
X
)
g(X)
g(X)的期望值定义为:
E
[
g
(
X
)
]
≜
∑
x
∈
V
a
l
(
X
)
g
(
x
)
p
X
(
x
)
E[g(X)] \triangleq \sum_{x \in V a l(X)} g(x) p_{X}(x)
E[g(X)]≜x∈Val(X)∑g(x)pX(x)
如果
X
X
X是一个连续的随机变量,其PDF 为
f
X
(
x
)
f _X (x)
fX(x),那么
g
(
X
)
g(X)
g(X)的期望值被定义为:
E
[
g
(
X
)
]
≜
∫
−
∞
∞
g
(
x
)
f
X
(
x
)
d
x
E[g(X)] \triangleq \int_{-\infty}^{\infty} g(x) f_{X}(x) d x
E[g(X)]≜∫−∞∞g(x)fX(x)dx
直觉上, g ( X ) g(X) g(X)的期望值可以被认为是 g ( x ) g(x) g(x)对于不同的 x x x值可以取的值的“加权平均值”,其中权重由 p X ( x ) p_X(x) pX(x)或 f X ( x ) f_X(x) fX(x)给出。作为上述情况的特例,请注意,随机变量本身的期望值,是通过令 g ( x ) = x g(x) = x g(x)=x得到的,这也被称为随机变量的平均值。
性质:
- 对于任意常数 a ∈ R a \in \mathbb{R} a∈R, E [ a ] = a E[a]=a E[a]=a
- 对于任意常数 a ∈ R a \in \mathbb{R} a∈R, E [ a f ( X ) ] = a E [ f ( X ) ] E[af(X)]=aE[f(X)] E[af(X)]=aE[f(X)]
- (线性期望): E [ f ( X ) + g ( X ) ] = E [ f ( X ) ] + E [ g ( X ) ] E[f(X)+g(X)]=E[f(X)]+E[g(X)] E[f(X)+g(X)]=E[f(X)]+E[g(X)]
- 对于一个离散随机变量 X X X, E [ 1 { X = k } ] = P ( X = k ) E[1\{X=k\}]=P(X=k) E[1{X=k}]=P(X=k)
2.5 方差
随机变量
X
X
X的方差是随机变量
X
X
X的分布围绕其平均值集中程度的度量。形式上,随机变量
X
X
X的方差定义为:
Var
[
X
]
≜
E
[
(
X
−
E
(
X
)
)
2
]
\operatorname{Var}[X] \triangleq E\left[(X-E(X))^{2}\right]
Var[X]≜E[(X−E(X))2]
使用上一节中的性质,我们可以导出方差的替代表达式:
E
[
(
X
−
E
[
X
]
)
2
]
=
E
[
X
2
−
2
E
[
X
]
X
+
E
[
X
]
2
]
=
E
[
X
2
]
−
2
E
[
X
]
E
[
X
]
+
E
[
X
]
2
=
E
[
X
2
]
−
E
[
X
]
2
\begin{aligned} E\left[(X-E[X])^{2}\right] &=E\left[X^{2}-2 E[X] X+E[X]^{2}\right] \\ &=E\left[X^{2}\right]-2 E[X] E[X]+E[X]^{2} \\ &=E\left[X^{2}\right]-E[X]^{2} \end{aligned}
E[(X−E[X])2]=E[X2−2E[X]X+E[X]2]=E[X2]−2E[X]E[X]+E[X]2=E[X2]−E[X]2
其中第二个等式来自期望的线性,以及 E [ X ] E[X] E[X]相对于外层期望实际上是常数的事实。
性质:
- 对于任意常数 a ∈ R a \in \mathbb{R} a∈R, V a l [ a ] = 0 Val[a]=0 Val[a]=0
- 对于任意常数 a ∈ R a \in \mathbb{R} a∈R, V a r [ a f ( X ) ] = a 2 V a r [ f ( X ) ] Var[af(X)]=a^2Var[f(X)] Var[af(X)]=a2Var[f(X)]
举例:
计算均匀随机变量
X
X
X的平均值和方差,任意
x
∈
[
0
,
1
]
x \in [0,1]
x∈[0,1],其PDF为
p
X
(
x
)
=
1
p_X(x)= 1
pX(x)=1,其他地方为0。
E
[
X
]
=
∫
−
∞
∞
x
f
X
(
x
)
d
x
=
∫
0
1
x
d
x
=
1
2
E[X]=\int_{-\infty}^{\infty} x f_{X}(x) d x=\int_{0}^{1} x d x=\frac{1}{2}
E[X]=∫−∞∞xfX(x)dx=∫01xdx=21
E
[
X
2
]
=
∫
−
∞
∞
x
2
f
X
(
x
)
d
x
=
∫
0
1
x
2
d
x
=
1
3
E\left[X^{2}\right]=\int_{-\infty}^{\infty} x^{2} f_{X}(x) d x=\int_{0}^{1} x^{2} d x=\frac{1}{3}
E[X2]=∫−∞∞x2fX(x)dx=∫01x2dx=31
V
a
r
[
X
]
=
E
[
X
2
]
−
E
[
X
]
2
=
1
3
−
1
4
=
1
12
Var[X]=E[X^2]-E[X]^2=\frac{1}{3}-\frac{1}{4}=\frac{1}{12}
Var[X]=E[X2]−E[X]2=31−41=121
举例:
假设对于一些子集 A ⊆ Ω A \subseteq \Omega A⊆Ω,有 g ( x ) = 1 { x ∈ A } g(x) = 1\{x \in A\} g(x)=1{x∈A},计算 E [ g ( X ) ] E[g(X)] E[g(X)]?
离散情况:
E
[
g
(
X
)
]
=
∑
x
∈
V
a
l
(
X
)
1
{
x
∈
A
}
P
X
(
x
)
d
x
=
∑
x
∈
A
P
X
(
x
)
d
x
=
P
(
x
∈
A
)
E[g(X)]=\sum_{x \in V a l(X)} 1\{x \in A\} P_{X}(x) d x=\sum_{x \in A} P_{X}(x) d x=P(x \in A)
E[g(X)]=x∈Val(X)∑1{x∈A}PX(x)dx=x∈A∑PX(x)dx=P(x∈A)
连续情况:
E
[
g
(
X
)
]
=
∫
−
∞
∞
1
{
x
∈
A
}
f
X
(
x
)
d
x
=
∫
x
∈
A
f
X
(
x
)
d
x
=
P
(
x
∈
A
)
E[g(X)]=\int_{-\infty}^{\infty} 1\{x \in A\} f_{X}(x) d x=\int_{x \in A} f_{X}(x) d x=P(x \in A)
E[g(X)]=∫−∞∞1{x∈A}fX(x)dx=∫x∈AfX(x)dx=P(x∈A)
2.6 一些常见的随机变量
离散随机变量
-
伯努利分布:硬币掷出正面的概率为 p p p(其中: 0 ≤ p ≤ 1 0 \leq p \leq 1 0≤p≤1),如果正面发生,则为1,否则为0。
p ( x ) = { p if p = 1 1 − p if p = 0 p(x)=\left\{\begin{array}{ll}{p} & {\text { if } p=1} \\ {1-p} & {\text { if } p=0}\end{array}\right. p(x)={p1−p if p=1 if p=0 -
二项式分布:掷出正面概率为 p p p(其中: 0 ≤ p ≤ 1 0 \leq p \leq 1 0≤p≤1)的硬币 n n n次独立投掷中正面的数量。
p ( x ) = ( n x ) p x ( 1 − p ) n − x p(x)=\left(\begin{array}{l}{n} \\ {x}\end{array}\right) p^{x}(1-p)^{n-x} p(x)=(nx)px(1−p)n−x -
几何分布:掷出正面概率为 p p p(其中: p > 0 p >0 p>0)的硬币第一次掷出正面所需要的次数。
-
泊松分布:用于模拟罕见事件频率的非负整数的概率分布(其中: λ > 0 \lambda >0 λ>0)。
p ( x ) = e − λ λ x x ! p(x)=e^{-\lambda} \frac{\lambda^{x}}{x !} p(x)=e−λx!λx
连续随机变量
-
均匀分布:在 a a a和 b b b之间每个点概率密度相等的分布(其中: a < b a<b a<b)。
f ( x ) = { 1 b − a if a ≤ x ≤ b 0 otherwise f(x)=\left\{\begin{array}{ll}{\frac{1}{b-a}} & {\text { if } a \leq x \leq b} \\ {0} & {\text { otherwise }}\end{array}\right. f(x)={b−a10 if a≤x≤b otherwise -
指数分布:在非负实数上有衰减的概率密度(其中: λ > 0 \lambda >0 λ>0)。
f ( x ) = { λ e − λ x if x ≥ 0 0 otherwise f(x)=\left\{\begin{array}{ll}{\lambda e^{-\lambda x}} & {\text { if } x \geq 0} \\ {0} & {\text { otherwise }}\end{array}\right. f(x)={λe−λx0 if x≥0 otherwise -
正态分布:又被称为高斯分布。
f ( x ) = 1 2 π σ e − 1 2 σ 2 ( x − μ ) 2 f(x)=\frac{1}{\sqrt{2 \pi} \sigma} e^{-\frac{1}{2 \sigma^{2}}(x-\mu)^{2}} f(x)=2πσ1e−2σ21(x−μ)2
一些随机变量的概率密度函数和累积分布函数的形状如图2所示。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-tyIgawhu-1603691202076)(images/b958c16cfdce9e6bd2b810b10d71416e.png)]
分布 | 概率密度函数(PDF)或者概率质量函数(PMF) | 均值 | 方差 |
---|---|---|---|
B e r n o u l l i ( p ) Bernoulli(p) Bernoulli(p)(伯努利分布) | { p if x = 1 1 − p if x = 0 \left\{\begin{array}{ll}{p} & {\text { if } x=1} \\ {1-p} & {\text { if } x=0}\end{array}\right. {p1−p if x=1 if x=0 | p p p | p ( 1 − p ) p(1-p) p(1−p) |
B i n o m i a l ( n , p ) Binomial(n,p) Binomial(n,p)(二项式分布) | ( n k ) p k ( 1 − p ) n − k \left(\begin{array}{l}{n} \\ {k}\end{array}\right) p^{k}(1-p)^{n-k} (nk)pk(1−p)n−k 其中: 0 ≤ k ≤ n 0 \leq k \leq n 0≤k≤n | n p np np | n p q npq npq |
G e o m e t r i c ( p ) Geometric(p) Geometric(p)(几何分布) | p ( 1 − p ) k − 1 p(1-p)^{k-1} p(1−p)k−1 其中: k = 1 , 2 , ⋯ k=1,2,\cdots k=1,2,⋯ | 1 p \frac{1}{p} p1 | 1 − p p 2 \frac {1-p}{p^2} p21−p |
P o i s s o n ( λ ) Poisson(\lambda) Poisson(λ)(泊松分布) | e − λ λ x / x ! e^{-\lambda} \lambda^{x} / x ! e−λλx/x! 其中: k = 1 , 2 , ⋯ k=1,2,\cdots k=1,2,⋯ | λ \lambda λ | λ \lambda λ |
U n i f o r m ( a , b ) Uniform(a,b) Uniform(a,b)(均匀分布) | 1 b − a \frac{1}{b-a} b−a1 存在 x ∈ ( a , b ) x \in (a,b) x∈(a,b) | a + b 2 \frac{a+b}{2} 2a+b | ( b − a ) 2 12 \frac{(b-a)^2}{12} 12(b−a)2 |
G a u s s i a n ( μ , σ 2 ) Gaussian(\mu,\sigma^2) Gaussian(μ,σ2)(高斯分布) | 1 2 π σ e − 1 2 σ 2 ( x − μ ) 2 \frac{1}{\sqrt{2 \pi} \sigma} e^{-\frac{1}{2 \sigma^{2}}(x-\mu)^{2}} 2πσ1e−2σ21(x−μ)2 | μ \mu μ | σ 2 \sigma^2 σ2 |
E x p o n e n t i a l ( λ ) Exponential(\lambda) Exponential(λ)(指数分布) | λ e − λ x \lambda e^{-\lambda x} λe−λx x ≥ 0 , λ > 0 x\geq0,\lambda>0 x≥0,λ>0 | 1 λ \frac{1}{\lambda} λ1 | 1 λ 2 \frac{1}{\lambda^2} λ21 |
3. 两个随机变量
到目前为止,我们已经考虑了单个随机变量。然而,在许多情况下,在随机实验中,我们可能有不止一个感兴趣的量。例如,在一个我们掷硬币十次的实验中,我们可能既关心 X ( ω ) = X(\omega) = X(ω)=出现的正面数量,也关心 Y ( ω ) = Y (\omega) = Y(ω)=连续最长出现正面的长度。在本节中,我们考虑两个随机变量的设置。
3.1 联合分布和边缘分布
假设我们有两个随机变量,一个方法是分别考虑它们。如果我们这样做,我们只需要
F
X
(
x
)
F_X (x)
FX(x)和
F
Y
(
y
)
F_Y (y)
FY(y)。但是如果我们想知道在随机实验的结果中,
X
X
X和
Y
Y
Y同时假设的值,我们需要一个更复杂的结构,称为
X
X
X和
Y
Y
Y的联合累积分布函数,定义如下:
F
X
Y
(
x
,
y
)
=
P
(
X
≤
x
,
Y
≤
y
)
F_{XY}(x,y)=P(X \leq x,Y \leq y)
FXY(x,y)=P(X≤x,Y≤y)
可以证明,通过了解联合累积分布函数,可以计算出任何涉及到 X X X和 Y Y Y的事件的概率。
联合CDF:
F
X
Y
(
x
,
y
)
F_{XY }(x,y)
FXY(x,y)和每个变量的联合分布函数
F
X
(
x
)
F_X(x)
FX(x)和
F
Y
(
y
)
F_Y (y)
FY(y)分别由下式关联:
F
X
(
x
)
=
lim
y
→
∞
F
X
Y
(
x
,
y
)
d
y
F_{X}(x)=\lim _{y \rightarrow \infty} F_{X Y}(x, y) d y
FX(x)=y→∞limFXY(x,y)dy
F
Y
(
y
)
=
lim
y
→
∞
F
X
Y
(
x
,
y
)
d
x
F_{Y}(y)=\lim _{y \rightarrow \infty} F_{X Y}(x, y) dx
FY(y)=y→∞limFXY(x,y)dx
这里我们称
F
X
(
x
)
F_X(x)
FX(x)和
F
Y
(
y
)
F_Y (y)
FY(y)为
F
X
Y
(
x
,
y
)
F_{XY }(x,y)
FXY(x,y)的边缘累积概率分布函数。
性质:
- 0 ≤ F X Y ( x , y ) ≤ 1 0 \leq F_{XY }(x,y) \leq 1 0≤FXY(x,y)≤1
- lim x , y → ∞ F X Y ( x , y ) = 1 \lim _{x, y \rightarrow \infty} F_{X Y}(x, y)=1 limx,y→∞FXY(x,y)=1
- lim x , y → − ∞ F X Y ( x , y ) = 0 \lim _{x, y \rightarrow -\infty} F_{X Y}(x, y)=0 limx,y→−∞FXY(x,y)=0
- F X ( x ) = lim y → ∞ F X Y ( x , y ) F_{X}(x)=\lim _{y \rightarrow \infty} F_{X Y}(x, y) FX(x)=limy→∞FXY(x,y)
3.2 联合概率和边缘概率质量函数
如果 X X X和 Y Y Y是离散随机变量,那么联合概率质量函数 p X Y : R × R → [ 0 , 1 ] p_{X Y} : \mathbb{R} \times \mathbb{R} \rightarrow [0,1] pXY:R×R→[0,1]由下式定义:
p X Y ( x , y ) = P ( X = x , Y = y ) p_{X Y}(x,y)=P(X=x,Y=y) pXY(x,y)=P(X=x,Y=y)
这里, 对于任意 x x x, y y y, 0 ≤ P X Y ( x , y ) ≤ 1 0 \leq P_{XY} (x,y) \leq 1 0≤PXY(x,y)≤1, 并且 ∑ x ∈ V a l ( X ) ∑ y ∈ V a l ( Y ) P X Y ( x , y ) = 1 \sum_{x \in V a l(X)} \sum_{y \in V a l(Y)} P_{X Y}(x, y)=1 ∑x∈Val(X)∑y∈Val(Y)PXY(x,y)=1
两个变量上的联合PMF分别与每个变量的概率质量函数有什么关系?事实上:
p
X
(
x
)
=
∑
y
p
X
Y
(
x
,
y
)
p_{X}(x)=\sum_{y} p_{X Y}(x, y)
pX(x)=y∑pXY(x,y)
对于 p Y ( y ) p_Y (y) pY(y)类似。在这种情况下,我们称 p X ( x ) p_X(x) pX(x)为 X X X的边际概率质量函数。在统计学中,将一个变量相加形成另一个变量的边缘分布的过程通常称为“边缘化”。
3.3 联合概率和边缘概率密度函数
假设
X
X
X和
Y
Y
Y是两个连续的随机变量,具有联合分布函数
F
X
Y
F_{XY}
FXY。在
F
X
Y
(
x
,
y
)
F_{XY}(x,y)
FXY(x,y)在
x
x
x和
y
y
y中处处可微的情况下,我们可以定义联合概率密度函数:
f
X
Y
(
x
,
y
)
=
∂
2
F
X
Y
(
x
,
y
)
∂
x
∂
y
f_{X Y}(x, y)=\frac{\partial^{2} F_{X Y}(x, y)}{\partial x \partial y}
fXY(x,y)=∂x∂y∂2FXY(x,y)
如同在一维情况下,
f
X
Y
(
x
,
y
)
≠
P
(
X
=
x
,
Y
=
y
)
f_{XY}(x,y)\not= P(X = x,Y = y)
fXY(x,y)=P(X=x,Y=y),而是:
∬
x
∈
A
f
X
Y
(
x
,
y
)
d
x
d
y
=
P
(
(
X
,
Y
)
∈
A
)
\iint_{x \in A} f_{X Y}(x, y) d x d y=P((X, Y) \in A)
∬x∈AfXY(x,y)dxdy=P((X,Y)∈A)
请注意,概率密度函数 f X Y ( x , y ) f_{XY}(x,y) fXY(x,y)的值总是非负的,但它们可能大于1。尽管如此,可以肯定的是 ∫ − ∞ ∞ ∫ − ∞ ∞ f X Y ( x , y ) = 1 \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} f_{X Y}(x, y)=1 ∫−∞∞∫−∞∞fXY(x,y)=1
与离散情况相似,我们定义:
f
X
(
x
)
=
∫
−
∞
∞
f
X
Y
(
x
,
y
)
d
y
f_{X}(x)=\int_{-\infty}^{\infty} f_{X Y}(x, y) d y
fX(x)=∫−∞∞fXY(x,y)dy
作为
X
X
X的边际概率密度函数(或边际密度),对于
f
Y
(
y
)
f_Y (y)
fY(y)也类似。
3.4 条件概率分布
条件分布试图回答这样一个问题,当我们知道
X
X
X必须取某个值
x
x
x时,
Y
Y
Y上的概率分布是什么?在离散情况下,给定
Y
Y
Y的条件概率质量函数是简单的:
p
Y
∣
X
(
y
∣
x
)
=
p
X
Y
(
x
,
y
)
p
X
(
x
)
p_{Y | X}(y | x)=\frac{p_{X Y}(x, y)}{p_{X}(x)}
pY∣X(y∣x)=pX(x)pXY(x,y)
假设分母不等于0。
在连续的情况下,在技术上要复杂一点,因为连续随机变量的概率等于零。忽略这一技术点,我们通过类比离散情况,简单地定义给定
X
=
x
X = x
X=x的条件概率密度为:
f
Y
∣
X
(
y
∣
x
)
=
f
X
Y
(
x
,
y
)
f
X
(
x
)
f_{Y | X}(y | x)=\frac{f_{X Y}(x, y)}{f_{X}(x)}
fY∣X(y∣x)=fX(x)fXY(x,y)
假设分母不等于0。
3.5 贝叶斯定理
当试图推导一个变量给定另一个变量的条件概率表达式时,经常出现的一个有用公式是贝叶斯定理。
对于离散随机变量
X
X
X和
Y
Y
Y:
P
Y
∣
X
(
y
∣
x
)
=
P
X
Y
(
x
,
y
)
P
X
(
x
)
=
P
X
∣
Y
(
x
∣
y
)
P
Y
(
y
)
∑
y
′
∈
V
a
l
(
Y
)
P
X
∣
Y
(
x
∣
y
′
)
P
Y
(
y
′
)
P_{Y | X}(y | x)=\frac{{P_{XY}}(x, y)}{P_{X}(x)}=\frac{P_{X | Y}(x | y) P_{Y}(y)}{\sum_{y^{\prime} \in V a l(Y)} P_{X | Y}\left(x | y^{\prime}\right) P_{Y}\left(y^{\prime}\right)}
PY∣X(y∣x)=PX(x)PXY(x,y)=∑y′∈Val(Y)PX∣Y(x∣y′)PY(y′)PX∣Y(x∣y)PY(y)
对于连续随机变量 X X X和 Y Y Y:
f Y ∣ X ( y ∣ x ) = f X Y ( x , y ) f X ( x ) = f X ∣ Y ( x ∣ y ) f Y ( y ) ∫ − ∞ ∞ f X ∣ Y ( x ∣ y ′ ) f Y ( y ′ ) d y ′ f_{Y | X}(y | x)=\frac{f_{X Y}(x, y)}{f_{X}(x)}=\frac{f_{X | Y}(x | y) f_{Y}(y)}{\int_{-\infty}^{\infty} f_{X | Y}\left(x | y^{\prime}\right) f_{Y}\left(y^{\prime}\right) d y^{\prime}} fY∣X(y∣x)=fX(x)fXY(x,y)=∫−∞∞fX∣Y(x∣y′)fY(y′)dy′fX∣Y(x∣y)fY(y)
3.6 独立性
如果对于 X X X和 Y Y Y的所有值, F X Y ( x , y ) = F X ( x ) F Y ( y ) F_{XY}(x,y) = F_X(x)F_Y(y) FXY(x,y)=FX(x)FY(y),则两个随机变量 X X X和 Y Y Y是独立的。等价地,
- 对于离散随机变量, 对于任意 x ∈ V a l ( X ) x \in Val(X) x∈Val(X), y ∈ V a l ( Y ) y \in Val(Y) y∈Val(Y) , p X Y ( x , y ) = p X ( x ) p Y ( y ) p_{XY}(x,y) = p_X (x)p_Y (y) pXY(x,y)=pX(x)pY(y)。
- 对于离散随机变量, p Y ∣ X ( y ∣ x ) = p Y ( y ) p_Y |X (y|x) = p_Y (y) pY∣X(y∣x)=pY(y)当对于任意 y ∈ V a l ( Y ) y \in Val(Y) y∈Val(Y)且 p X ( x ) ≠ 0 p_X (x) \not= 0 pX(x)=0。
- 对于连续随机变量, f X Y ( x , y ) = f X ( x ) f Y ( y ) f_{XY}(x,y) = f_X (x)f_Y(y) fXY(x,y)=fX(x)fY(y) 对于任意 x , y ∈ R x,y \in \mathbb{R} x,y∈R。
- 对于连续随机变量, f Y ∣ X ( y ∣ x ) = f Y ( y ) f_{Y |X} (y|x) = f_Y (y) fY∣X(y∣x)=fY(y) ,当 f X ( x ) ≠ 0 f_X (x)\not = 0 fX(x)=0对于任意 y ∈ R y \in \mathbb{R} y∈R。
非正式地说,如果“知道”一个变量的值永远不会对另一个变量的条件概率分布有任何影响,那么两个随机变量 X X X和 Y Y Y是独立的,也就是说,你只要知道 f ( x ) f(x) f(x)和 f ( y ) f(y) f(y)就知道关于这对变量 ( X , Y ) (X,Y) (X,Y)的所有信息。以下引理将这一观察形式化:
引理3.1
如果
X
X
X和
Y
Y
Y是独立的,那么对于任何
A
,
B
⊆
R
A,B⊆ \mathbb{R}
A,B⊆R,我们有:
P
(
X
∈
A
,
Y
∈
B
)
=
P
(
X
∈
A
)
P
(
Y
∈
B
)
P(X \in A, Y \in B)=P(X \in A) P(Y \in B)
P(X∈A,Y∈B)=P(X∈A)P(Y∈B)
利用上述引理,我们可以证明如果
X
X
X与
Y
Y
Y无关,那么
X
X
X的任何函数都与
Y
Y
Y的任何函数无关。
3.7 期望和协方差
假设我们有两个离散的随机变量
X
X
X,
Y
Y
Y并且
g
:
R
2
⟶
R
g : \mathbf{R}^{2} \longrightarrow \mathbf{R}
g:R2⟶R是这两个随机变量的函数。那么
g
g
g的期望值以如下方式定义:
E
[
g
(
X
,
Y
)
]
≜
∑
x
∈
V
a
l
(
X
)
∑
y
∈
V
a
l
(
Y
)
g
(
x
,
y
)
p
X
Y
(
x
,
y
)
E[g(X, Y)] \triangleq \sum_{x \in V a l(X)} \sum_{y \in V a l(Y)} g(x, y) p_{X Y}(x, y)
E[g(X,Y)]≜x∈Val(X)∑y∈Val(Y)∑g(x,y)pXY(x,y)
对于连续随机变量
X
X
X,
Y
Y
Y,类似的表达式是:
E
[
g
(
X
,
Y
)
]
=
∫
−
∞
∞
∫
−
∞
∞
g
(
x
,
y
)
f
X
Y
(
x
,
y
)
d
x
d
y
E[g(X, Y)]=\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} g(x, y) f_{X Y}(x, y) d x d y
E[g(X,Y)]=∫−∞∞∫−∞∞g(x,y)fXY(x,y)dxdy
我们可以用期望的概念来研究两个随机变量之间的关系。特别地,两个随机变量的协方差定义为:
C
o
v
[
X
,
Y
]
≜
E
[
(
X
−
E
[
X
]
)
(
Y
−
E
[
Y
]
)
]
{Cov}[X, Y] \triangleq E[(X-E[X])(Y-E[Y])]
Cov[X,Y]≜E[(X−E[X])(Y−E[Y])]
使用类似于方差的推导,我们可以将它重写为:
C
o
v
[
X
,
Y
]
=
E
[
(
X
−
E
[
X
]
)
(
Y
−
E
[
Y
]
)
]
=
E
[
X
Y
−
X
E
[
Y
]
−
Y
E
[
X
]
+
E
[
X
]
E
[
Y
]
]
=
E
[
X
Y
]
−
E
[
X
]
E
[
Y
]
−
E
[
Y
]
E
[
X
]
+
E
[
X
]
E
[
Y
]
]
=
E
[
X
Y
]
−
E
[
X
]
E
[
Y
]
\begin{aligned} {Cov}[X, Y] &=E[(X-E[X])(Y-E[Y])] \\ &=E[X Y-X E[Y]-Y E[X]+E[X] E[Y]] \\ &=E[X Y]-E[X] E[Y]-E[Y] E[X]+E[X] E[Y]] \\ &=E[X Y]-E[X] E[Y] \end{aligned}
Cov[X,Y]=E[(X−E[X])(Y−E[Y])]=E[XY−XE[Y]−YE[X]+E[X]E[Y]]=E[XY]−E[X]E[Y]−E[Y]E[X]+E[X]E[Y]]=E[XY]−E[X]E[Y]
在这里,说明两种协方差形式相等的关键步骤是第三个等号,在这里我们使用了这样一个事实,即 E [ X ] E[X] E[X]和 E [ Y ] E[Y] E[Y]实际上是常数,可以被提出来。当 c o v [ X , Y ] = 0 cov[X,Y] = 0 cov[X,Y]=0时,我们说 X X X和 Y Y Y不相关。
性质:
- (期望线性) E [ f ( X , Y ) + g ( X , Y ) ] = E [ f ( X , Y ) ] + E [ g ( X , Y ) ] E[f(X,Y ) + g(X,Y)] = E[f(X,Y )] + E[g(X,Y)] E[f(X,Y)+g(X,Y)]=E[f(X,Y)]+E[g(X,Y)]
- V a r [ X + Y ] = V a r [ X ] + V a r [ Y ] + 2 C o v [ X , Y ] V ar[X + Y ] = V ar[X] + V ar[Y ] + 2Cov[X,Y] Var[X+Y]=Var[X]+Var[Y]+2Cov[X,Y]
- 如果 X X X和 Y Y Y相互独立, 那么 C o v [ X , Y ] = 0 Cov[X,Y ] = 0 Cov[X,Y]=0
- 如果 X X X和 Y Y Y相互独立, 那么 E [ f ( X ) g ( Y ) ] = E [ f ( X ) ] E [ g ( Y ) ] E[f(X)g(Y )] = E[f(X)]E[g(Y)] E[f(X)g(Y)]=E[f(X)]E[g(Y)].
4. 多个随机变量
上一节介绍的概念和想法可以推广到两个以上的随机变量。特别是,假设我们有 n n n个连续随机变量, X 1 ( ω ) , X 2 ( ω ) , ⋯ X n ( ω ) X _1 (\omega),X_2 (\omega),\cdots X_n (\omega) X1(ω),X2(ω),⋯Xn(ω)。在本节中,为了表示简单,我们只关注连续的情况,对离散随机变量的推广工作类似。
4.1 基本性质
我们可以定义
X
1
,
X
2
,
⋯
,
X
n
X_1,X_2,\cdots,X_n
X1,X2,⋯,Xn的联合累积分布函数、联合概率密度函数,以及给定
X
2
,
⋯
,
X
n
X_2,\cdots,X_n
X2,⋯,Xn时
X
1
X_1
X1的边缘概率密度函数为:
F
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
=
P
(
X
1
≤
x
1
,
X
2
≤
x
2
,
…
,
X
n
≤
x
n
)
F_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right)=P\left(X_{1} \leq x_{1}, X_{2} \leq x_{2}, \ldots, X_{n} \leq x_{n}\right)
FX1,X2,…,Xn(x1,x2,…xn)=P(X1≤x1,X2≤x2,…,Xn≤xn)
f
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
=
∂
n
F
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
∂
x
1
…
∂
x
n
f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right)=\frac{\partial^{n} F_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right)}{\partial x_{1} \ldots \partial x_{n}}
fX1,X2,…,Xn(x1,x2,…xn)=∂x1…∂xn∂nFX1,X2,…,Xn(x1,x2,…xn)
f
X
1
(
X
1
)
=
∫
−
∞
∞
⋯
∫
−
∞
∞
f
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
d
x
2
…
d
x
n
f_{X_{1}}\left(X_{1}\right)=\int_{-\infty}^{\infty} \cdots \int_{-\infty}^{\infty} f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right) d x_{2} \ldots d x_{n}
fX1(X1)=∫−∞∞⋯∫−∞∞fX1,X2,…,Xn(x1,x2,…xn)dx2…dxn
f X 1 ∣ X 2 , … , X n ( x 1 ∣ x 2 , … x n ) = f X 1 , X 2 , … , X n ( x 1 , x 2 , … x n ) f X 2 , … , X n ( x 1 , x 2 , … x n ) f_{X_{1} | X_{2}, \ldots, X_{n}}\left(x_{1} | x_{2}, \dots x_{n}\right)=\frac{f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \dots x_{n}\right)}{f_{X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right)} fX1∣X2,…,Xn(x1∣x2,…xn)=fX2,…,Xn(x1,x2,…xn)fX1,X2,…,Xn(x1,x2,…xn)
为了计算事件
A
⊆
R
n
A \subseteq \mathbb{R}^{n}
A⊆Rn的概率,我们有:
P
(
(
x
1
,
x
2
,
…
x
n
)
∈
A
)
=
∫
(
x
1
,
x
2
,
…
x
n
)
∈
A
f
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
d
x
1
d
x
2
…
d
x
n
P\left(\left(x_{1}, x_{2}, \ldots x_{n}\right) \in A\right)=\int_{\left(x_{1}, x_{2}, \ldots x_{n}\right) \in A} f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right) d x_{1} d x_{2} \ldots d x_{n}
P((x1,x2,…xn)∈A)=∫(x1,x2,…xn)∈AfX1,X2,…,Xn(x1,x2,…xn)dx1dx2…dxn
链式法则:
从多个随机变量的条件概率的定义中,可以看出:
f
(
x
1
,
x
2
,
…
,
x
n
)
=
f
(
x
n
∣
x
1
,
x
2
…
,
x
n
−
1
)
f
(
x
1
,
x
2
…
,
x
n
−
1
)
=
f
(
x
n
∣
x
1
,
x
2
…
,
x
n
−
1
)
f
(
x
n
−
1
∣
x
1
,
x
2
…
,
x
n
−
2
)
f
(
x
1
,
x
2
…
,
x
n
−
2
)
=
⋯
=
f
(
x
1
)
∏
i
=
2
n
f
(
x
i
∣
x
1
,
…
,
x
i
−
1
)
\begin{aligned} f\left(x_{1}, x_{2}, \ldots, x_{n}\right) &=f\left(x_{n} | x_{1}, x_{2} \ldots, x_{n-1}\right) f\left(x_{1}, x_{2} \ldots, x_{n-1}\right) \\ &=f\left(x_{n} | x_{1}, x_{2} \ldots, x_{n-1}\right) f\left(x_{n-1} | x_{1}, x_{2} \ldots, x_{n-2}\right) f\left(x_{1}, x_{2} \ldots, x_{n-2}\right) \\ &=\cdots=f\left(x_{1}\right) \prod_{i=2}^{n} f\left(x_{i} | x_{1}, \ldots, x_{i-1}\right) \end{aligned}
f(x1,x2,…,xn)=f(xn∣x1,x2…,xn−1)f(x1,x2…,xn−1)=f(xn∣x1,x2…,xn−1)f(xn−1∣x1,x2…,xn−2)f(x1,x2…,xn−2)=⋯=f(x1)i=2∏nf(xi∣x1,…,xi−1)
独立性:对于多个事件,
A
1
,
⋯
,
A
k
A_1,\cdots ,A_k
A1,⋯,Ak,我们说
A
1
,
⋯
,
A
k
A_1,\cdots ,A_k
A1,⋯,Ak 是相互独立的,当对于任何子集
S
⊆
{
1
,
2
,
⋯
,
k
}
S \subseteq \{1,2,\cdots,k\}
S⊆{1,2,⋯,k},我们有:
P
(
∩
i
∈
S
A
i
)
=
∏
i
∈
S
P
(
A
i
)
P\left(\cap_{i \in S} A_{i}\right)=\prod_{i \in S} P\left(A_{i}\right)
P(∩i∈SAi)=i∈S∏P(Ai)
同样,我们说随机变量 X 1 , X 2 , ⋯ , X n X_1,X_2,\cdots,X_n X1,X2,⋯,Xn是独立的,如果:
f
(
x
1
,
⋯
,
x
n
)
=
f
(
x
1
)
f
(
x
2
)
⋯
f
(
x
n
)
f(x_1,\cdots,x_n)=f(x_1)f(x_2)\cdots f(x_n)
f(x1,⋯,xn)=f(x1)f(x2)⋯f(xn)
这里,相互独立性的定义只是两个随机变量独立性到多个随机变量的自然推广。
独立随机变量经常出现在机器学习算法中,其中我们假设属于训练集的训练样本代表来自某个未知概率分布的独立样本。为了明确独立性的重要性,考虑一个“坏的”训练集,我们首先从某个未知分布中抽取一个训练样本
(
x
(
1
)
,
y
(
1
)
)
(x^{ (1)},y^{(1)})
(x(1),y(1)),然后将完全相同的训练样本的
m
−
1
m-1
m−1个副本添加到训练集中。在这种情况下,我们有:
P
(
(
x
(
1
)
,
y
(
1
)
)
,
…
.
(
x
(
m
)
,
y
(
m
)
)
)
≠
∏
i
=
1
m
P
(
x
(
i
)
,
y
(
i
)
)
P\left(\left(x^{(1)}, y^{(1)}\right), \ldots .\left(x^{(m)}, y^{(m)}\right)\right) \neq \prod_{i=1}^{m} P\left(x^{(i)}, y^{(i)}\right)
P((x(1),y(1)),….(x(m),y(m)))=i=1∏mP(x(i),y(i))
尽管训练集的大小为 m m m,但这些例子并不独立!虽然这里描述的过程显然不是为机器学习算法建立训练集的明智方法,但是事实证明,在实践中,样本的不独立性确实经常出现,并且它具有减小训练集的“有效大小”的效果。
4.2 随机向量
假设我们有n个随机变量。当把所有这些随机变量放在一起工作时,我们经常会发现把它们放在一个向量中是很方便的…我们称结果向量为随机向量(更正式地说,随机向量是从 Ω \Omega Ω到 R n \mathbb{R}^n Rn的映射)。应该清楚的是,随机向量只是处理 n n n个随机变量的一种替代符号,因此联合概率密度函数和综合密度函数的概念也将适用于随机向量。
期望:
考虑
g
:
R
n
→
R
g : \mathbb{R}^n \rightarrow \mathbb{R}
g:Rn→R中的任意函数。这个函数的期望值 被定义为
E
[
g
(
X
)
]
=
∫
R
n
g
(
x
1
,
x
2
,
…
,
x
n
)
f
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
d
x
1
d
x
2
…
d
x
n
E
[
g
(
X
)
]
=
∫
R
n
g
(
x
1
,
x
2
,
…
,
x
n
)
f
X
1
,
X
2
,
…
,
X
n
(
x
1
,
x
2
,
…
x
n
)
d
x
1
d
x
2
…
d
x
n
E[g(X)]=\int_{\mathbb{R}^{n}} g\left(x_{1}, x_{2}, \ldots, x_{n}\right) f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right) d x_{1} d x_{2} \ldots d x_{n}E[g(X)]\\=\int_{\mathbb{R}^{n}} g\left(x_{1}, x_{2}, \ldots, x_{n}\right) f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots x_{n}\right) d x_{1} d x_{2} \ldots d x_{n}
E[g(X)]=∫Rng(x1,x2,…,xn)fX1,X2,…,Xn(x1,x2,…xn)dx1dx2…dxnE[g(X)]=∫Rng(x1,x2,…,xn)fX1,X2,…,Xn(x1,x2,…xn)dx1dx2…dxn
其中,
∫
R
n
\int_{\mathbb{R}^{n}}
∫Rn是从
−
∞
-\infty
−∞到
∞
\infty
∞的
n
n
n个连续积分。如果
g
g
g是从
R
n
\mathbb{R}^n
Rn到
R
m
\mathbb{R}^m
Rm的函数,那么
g
g
g的期望值是输出向量的元素期望值,即,如果
g
g
g是:
g
(
x
)
=
[
g
1
(
x
)
g
2
(
x
)
⋮
g
m
(
x
)
]
g(x)=\left[\begin{array}{c}{g_{1}(x)} \\ {g_{2}(x)} \\ {\vdots} \\ {g_{m}(x)}\end{array}\right]
g(x)=⎣⎢⎢⎢⎡g1(x)g2(x)⋮gm(x)⎦⎥⎥⎥⎤
那么,
E
[
g
(
X
)
]
=
[
E
[
g
1
(
X
)
]
E
[
g
2
(
X
)
]
⋮
E
[
g
m
(
X
)
]
]
E[g(X)]=\left[\begin{array}{c}{E\left[g_{1}(X)\right]} \\ {E\left[g_{2}(X)\right]} \\ {\vdots} \\ {E\left[g_{m}(X)\right]}\end{array}\right]
E[g(X)]=⎣⎢⎢⎢⎡E[g1(X)]E[g2(X)]⋮E[gm(X)]⎦⎥⎥⎥⎤
协方差矩阵:对于给定的随机向量
X
:
Ω
→
R
n
X:\Omega\rightarrow \mathbb{R}^n
X:Ω→Rn,其协方差矩阵
Σ
\Sigma
Σ是
n
×
n
n \times n
n×n平方矩阵,其输入由
Σ
i
j
=
C
o
v
[
X
i
,
X
j
]
\Sigma_{i j}={Cov}\left[X_{i}, X_{j}\right]
Σij=Cov[Xi,Xj]给出。从协方差的定义来看,我们有:
$$
\begin{aligned}
\begin{equation}
\Sigma=\left[\begin{array}{ccc}{{Cov}\left[X_{1}, X_{1}\right]} & {\cdots} & {{Cov}\left[X_{1}, X_{n}\right]} \ {\vdots} & {\ddots} & {\vdots} \ {{Cov}\left[X_{n}, X_{1}\right]} & {\cdots} & {{Cov}\left[X_{n}, X_{n}\right]}\end{array}\right]\
=\left[\begin{array}{ccc}{E\left[X_{1}^{2}\right]-E\left[X_{1}\right] E\left[X_{1}\right]} & {\cdots} & {E\left[X_{1} X_{n}\right]-E\left[X_{1}\right] E\left[X_{n}\right]} \ {\vdots} & {\ddots} & {\vdots} \ {E\left[X_{n} X_{1}\right]-E\left[X_{n}\right] E\left[X_{1}\right]} & {\cdots} & {E\left[X_{n}^{2}\right]-E\left[X_{n}\right] E\left[X_{n}\right]}\end{array}\right]\
=\left[\begin{array}{ccc}{E\left[X_{1}^{2}\right]} & {\cdots} & {E\left[X_{1} X_{n}\right]} \ {\vdots} & {\ddots} & {\vdots} \ {E\left[X_{n} X_{1}\right]} & {\cdots} & {E\left[X_{n}^{2}\right]}\end{array}\right]-\left[\begin{array}{ccc}{E\left[X_{1}\right] E\left[X_{1}\right]} & {\cdots} & {E\left[X_{1}\right] E\left[X_{n}\right]} \ {\vdots} & {\ddots} & {\vdots} \ {E\left[X_{n}\right] E\left[X_{1}\right]} & {\cdots} & {E\left[X_{n}\right] E\left[X_{n}\right]}\end{array}\right]\
=E\left[X X^{T}\right]-E[X] E[X]{T}=\ldots=E\left[(X-E[X])(X-E[X]){T}\right]
\end{equation}
\end{aligned}
$$
其中矩阵期望以明显的方式定义。
协方差矩阵有许多有用的属性:
- Σ ⪰ 0 \Sigma \succeq 0 Σ⪰0;也就是说, Σ \Sigma Σ是正半定的。
- Σ = Σ T \Sigma=\Sigma^T Σ=ΣT;也就是说, Σ \Sigma Σ是对称的。
4.3 多元高斯分布
随机向量上概率分布的一个特别重要的例子叫做多元高斯或多元正态分布。随机向量 X ∈ R n X\in \mathbb{R}^n X∈Rn被认为具有多元正态(或高斯)分布,当其具有均值 μ ∈ R n \mu \in \mathbb{R}^n μ∈Rn和协方差矩阵 Σ ∈ S + + n \Sigma \in \mathbb{S}_{++}^{n} Σ∈S++n(其中$ \mathbb{S}_{++}^{n} 指 对 称 正 定 指对称正定 指对称正定n \times n$矩阵的空间)
f X 1 , X 2 , … , X n ( x 1 , x 2 , … , x n ; μ , Σ ) = 1 ( 2 π ) n / 2 ∣ Σ ∣ 1 / 2 exp ( − 1 2 ( x − μ ) T Σ − 1 ( x − μ ) ) f_{X_{1}, X_{2}, \ldots, X_{n}}\left(x_{1}, x_{2}, \ldots, x_{n} ; \mu, \Sigma\right)=\frac{1}{(2 \pi)^{n / 2}|\Sigma|^{1 / 2}} \exp \left(-\frac{1}{2}(x-\mu)^{T} \Sigma^{-1}(x-\mu)\right) fX1,X2,…,Xn(x1,x2,…,xn;μ,Σ)=(2π)n/2∣Σ∣1/21exp(−21(x−μ)TΣ−1(x−μ))
我们把它写成 X ∼ N ( μ , Σ ) X \sim \mathcal{N}(\mu, \Sigma) X∼N(μ,Σ)。请注意,在 n = 1 n = 1 n=1的情况下,它降维成普通正态分布,其中均值参数为 μ 1 \mu_1 μ1,方差为 Σ 11 \Sigma_{11} Σ11。
一般来说,高斯随机变量在机器学习和统计中非常有用,主要有两个原因:
首先,在统计算法中对“噪声”建模时,它们非常常见。通常,噪声可以被认为是影响测量过程的大量小的独立随机扰动的累积;根据中心极限定理,独立随机变量的总和将趋向于“看起来像高斯”。
其次,高斯随机变量便于许多分析操作,因为实际中出现的许多涉及高斯分布的积分都有简单的封闭形式解。我们将在本课程稍后遇到这种情况。
5. 其他资源
一本关于CS229所需概率水平的好教科书是谢尔顿·罗斯的《概率第一课》(A First Course on Probability by Sheldon Ross)。