最大似然估计是机器学习中最常用的参数估计方法之一。整个建模过程需要一个似然函数来描述在不同模型参数下真实数据发生的概率,似然函数是关于模型参数的函数。最大似然估计就是寻找最优参数,使得观测数据发生的概率最大、统计模型与真实数据最相似。
本专栏之前的文章介绍了线性回归以及最小二乘法的数学推导过程。对于一组训练数据,使用线性回归建模,可以有不同的模型参数来描述数据,这时候可以用最小二乘法来选择最优参数来拟合训练数据,即使用误差的平方作为损失函数。机器学习求解参数的过程被称为参数估计,机器学习问题也变成求使损失函数最小的最优化问题。最小二乘法比较直观,很容易解释,但不具有普遍意义,对于更多其他机器学习问题,比如二分类和多分类问题,最小二乘法就难以派上用场了。本文将给大家介绍一个具有普遍意义的参数估计方法:最大似然估计。
我们以一个赌博的例子来模拟机器学习的概率推理过程。假如你参与了一场赌博,你会被告知一个硬币抛掷10次的正反情况,接下来由你下注,而你只有一次机会,猜对硬币下一次正反情况则赢得100元,猜错则损失100元。这时,你会如何决策?
概率和似然
一般地,硬币有正反两面,如果硬币正反两面是均匀的,即每次抛掷后硬币为正的概率是0.5。使用这个硬币,很可能抛10次,有5次是正面。但是假如有人对硬币做了手脚,比如提前对硬币做了修改,硬币每次都会正面朝上,现在抛10次,10次都是正面,那么下次你绝对不会猜它是反面,因为前面的10次结果摆在那里,直觉上你不会相信这是一个普通的硬币。现在有一人抛了10次硬币,得到6正4反的结果,如何估算下次硬币为正的概率呢?
因为硬币并不是我们制作的,我们不了解硬币是否是完全均匀的,只能根据现在的观察结果来反推硬币的情况。假设硬币上有个参数θ \theta θ ,它决定了硬币的正反均匀程度,θ = 0.5 \theta = 0.5 θ = 0 . 5 表示正反均匀,每次抛硬币为正的概率为0.5,θ = 1.0 \theta = 1.0 θ = 1 . 0 表示硬币只有正面,每次抛硬币为正的概率为1。那么,从观察到的正反结果,反推硬币的构造参数θ \theta θ 的过程,就是一个参数估计的过程。
概率
抛掷10次硬币可能出现不同的情况,可以是“5正5反”、“4正6反”,“10正0反”等。假如我们知道硬币是如何构造的,即已知硬币的参数θ \theta θ ,那么出现“6正4反”的概率为:
P ( 6 正 4 反 ∣ θ = 0.5 ) = C 10 6 × 0. 5 6 × ( 1 − 0.5 ) 4 ≈ 0.2051 P ( 6 正 4 反 ∣ θ = 0.6 ) = C 10 6 × 0. 6 6 × ( 1 − 0.6 ) 4 ≈ 0.2508 P ( 6 正 4 反 ∣ θ = 0.9 ) = C 10 6 × 0. 9 6 × ( 1 − 0.9 ) 4 ≈ 0.0112
P(6正4反 \ |\ \theta=0.5)=C_{10}^{6}\times 0.5^6 \times (1-0.5)^4 \approx 0.2051 \\
P(6正4反 \ |\ \theta=0.6)=C_{10}^{6}\times 0.6^6 \times (1-0.6)^4 \approx 0.2508 \\
P(6正4反 \ |\ \theta=0.9)=C_{10}^{6}\times 0.9^6 \times (1-0.9)^4 \approx 0.0112
P ( 6 正 4 反 ∣ θ = 0 . 5 ) = C 1 0 6 × 0 . 5 6 × ( 1 − 0 . 5 ) 4 ≈ 0 . 2 0 5 1 P ( 6 正 4 反 ∣ θ = 0 . 6 ) = C 1 0 6 × 0 . 6 6 × ( 1 − 0 . 6 ) 4 ≈ 0 . 2 5 0 8 P ( 6 正 4 反 ∣ θ = 0 . 9 ) = C 1 0 6 × 0 . 9 6 × ( 1 − 0 . 9 ) 4 ≈ 0 . 0 1 1 2
上面这个公式是概率函数,表示已知参数θ \theta θ ,事实“6正4反”发生的概率。参数θ \theta θ 取不同的值时,事情发生的概率不同。在数学上一般使用P P P 或P r Pr P r 表示概率(Probability)函数。
上述过程中,抛10次硬币,要选出6次正面,使用了排列组合。因为“6正4反”可能会出现正正正正正正反反反反
、正正正正正反正反反反
、正正正正反正正反反反
等共210种组合,要在10次中选出6次为正面。假如每次正面的概率是0.6,那么反面的概率就是(1-0.6)。每次抛掷硬币的动作是相互独立,互不影响的,“6正4反”发生的概率就是各次抛掷硬币的概率乘积,再乘以210种组合。
概率反映的是:已知背后原因,推测某个结果发生的概率 。
似然
与概率不同,似然反映的是:已知结果,反推原因 。具体而言,似然(Likelihood)函数表示的是基于观察的数据,取不同的参数θ \theta θ 时,统计模型以多大的可能性接近真实观察数据。这就很像开篇提到的赌局,已经给你了一系列硬币正反情况,但你并不知道硬币的构造,下次下注时你要根据已有事实,反推硬币的构造。例如,当观察到硬币“10正0反”的事实,猜测硬币极有可能每次都是正面;当观察到硬币“6正4反”的事实,猜测硬币有可能不是正反均匀的,每次出现正面的可能性是0.6。
似然函数与前面的概率函数的计算方式极其相似,与概率函数不同的是,似然函数是θ \theta θ 的函数,即θ \theta θ 是未知的。似然函数衡量的是在不同参数θ \theta θ 下,真实观察数据发生的可能性。似然函数通常是多个观测数据发生的概率的联合概率,即多个观测数据都发生的概率。在机器学习里可以这样理解,目标y y y 和特征x \boldsymbol{x} x 同时发生,这些数值被观测到的概率。单个观测数据发生的可能性为P ( θ ) P(\theta) P ( θ ) ,如果各个观测之间是相互独立的,那么多个观测数据都发生的概率可表示为各个样本发生的概率的乘积 。
这里稍微解释一下事件独立性与联合概率之间的关系。如果事件A和事件B相互独立,那么事件A和B同时发生的概率是P ( A ) × P ( B ) P(A) \times P(B) P ( A ) × P ( B ) 。例如,事件“下雨”与事件“地面湿”就不是相互独立的,“下雨”与”地面湿”是同时发生、高度相关的,这两个事件都发生的概率就不能用概率的乘积来表示。两次抛掷硬币相互之间不影响,因此硬币正面朝上的概率可以用各次概率的乘积来表示。
似然函数通常用L L L 表示,对应英文Likelihood。观察到抛硬币“6正4反”的事实,硬币参数θ \theta θ 取不同值时,似然函数表示为:
L ( θ ; 6 正 4 反 ) = C 10 6 × θ 6 × ( 1 − θ ) 4
L(\theta ; 6正4反)=C_{10}^{6}\times \theta^6 \times (1-\theta)^4 \\
L ( θ ; 6 正 4 反 ) = C 1 0 6 × θ 6 × ( 1 − θ ) 4
这个公式的图形如下图所示。从图中可以看出:参数θ \theta θ 为0.6时,似然函数最大,参数为其他值时,“6正4反”发生的概率都相对更小。在这个赌局中,我会猜测下次硬币为正,因为根据已有观察,硬币很可能以0.6的概率为正。
“6正4反”的似然函数
推广到更为一般的场景,似然函数的一般形式可以用下面公式来表示,也就是之前提到的,各个样本发生的概率的乘积。
L ( θ ; X ) = P 1 ( θ ; X 1 ) × P 2 ( θ ; X 2 ) . . . × P n ( θ ; X n ) = ∏ P i ( θ ; X i )
L(\theta ; \mathbf{X}) = P_1(\theta ; X_1) \times P_2(\theta ; X_2) ... \times P_n(\theta ; X_n) = \prod P_i(\theta; X_i)
L ( θ ; X ) = P 1 ( θ ; X 1 ) × P 2 ( θ ; X 2 ) . . . × P n ( θ ; X n ) = ∏ P i ( θ ; X i )
最大似然估计
理解了似然函数的含义,就很容易理解最大似然估计的机制。似然函数是关于模型参数的函数,是描述观察到的真实数据在不同参数下发生的概率。最大似然估计要寻找最优参数,让似然函数最大化。或者说,使用最优参数时观测数据发生的概率最大。
线性回归的最大似然估计
之前的文章提到,线性回归的误差项ε ε ε 是预测值与真实值之间的差异,如下面公式所示。它可能是一些随机噪音,也可能是线性回归模型没考虑到的一些其他影响因素。
y ( i ) = ϵ ( i ) + ∑ j = 1 n w j x j ( i ) = ϵ ( i ) + w ⊤ x ( i )
y^{(i)} = \epsilon^{(i)} + \sum_{j=1}^n w_jx_{j}^{(i)} =\epsilon^{(i)} + \boldsymbol{w}^\top \boldsymbol{x}^{(i)}
y ( i ) = ϵ ( i ) + j = 1 ∑ n w j x j ( i ) = ϵ ( i ) + w ⊤ x ( i )
线性回归的一大假设是:误差服从均值为0的正态分布,且多个观测数据之间互不影响,相互独立。正态分布(高斯分布)的概率密度公式如下面公式,根据正态分布的公式,可以得到ϵ \epsilon ϵ 的概率密度。
假设x x x 服从正态分布,它的均值为μ \mu μ ,方差为σ \sigma σ ,它的概率密度公式如下。公式左侧的P ( x ; μ , σ ) P(x ; \mu, \sigma) P ( x ; μ , σ ) 表示x x x 是随机变量,; ; ; 分号强调μ \mu μ 和σ \sigma σ 不是随机变量,而是这个概率密度函数的参数。条件概率函数中使用的∣ | ∣ 竖线有明确的意义,P ( y ∣ x ) P(y|x) P ( y ∣ x ) 表示给定x x x (Given x x x ),y y y 发生的概率(Probability of y y y )。
P ( x ; μ , σ ) = 1 2 π σ 2 exp ( − ( x − μ ) 2 2 σ 2 )
P(x ; \mu, \sigma) = \frac{1}{\sqrt{2 \pi \sigma^2}}\exp{({-\frac{(x - \mu)^2}{2 \sigma^2}})}
P ( x ; μ , σ ) = 2 π σ 2 1 exp ( − 2 σ 2 ( x − μ ) 2 )
既然误差项服从正态分布,那么:
P ( ϵ ( i ) ) = 1 2 π σ 2 exp ( − ( ϵ ( i ) ) 2 2 σ 2 )
P(\epsilon^{(i)}) = \frac{1}{\sqrt{2 \pi \sigma^2}}\exp{({-\frac{(\epsilon^{(i)})^2}{2 \sigma^2}})}
P ( ϵ ( i ) ) = 2 π σ 2 1 exp ( − 2 σ 2 ( ϵ ( i ) ) 2 )
由于ϵ ( i ) = y ( i ) − w ⊤ x ( i ) \epsilon^{(i)} = y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)} ϵ ( i ) = y ( i ) − w ⊤ x ( i ) ,并取均值μ \mu μ 为0,可得到:
P ( y ( i ) ∣ x ( i ) ; w ) = 1 2 π σ 2 exp ( − ( y ( i ) − w ⊤ x ( i ) ) 2 2 σ 2 )
P(y^{(i)}|\boldsymbol{x}^{(i)}; \boldsymbol{w}) = \frac{1}{\sqrt{2 \pi \sigma^2}}\exp{({-\frac{(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2}{2 \sigma^2}})}
P ( y ( i ) ∣ x ( i ) ; w ) = 2 π σ 2 1 exp ( − 2 σ 2 ( y ( i ) − w ⊤ x ( i ) ) 2 )
上式表示给定x ( i ) \boldsymbol{x}^{(i)} x ( i ) ,y ( i ) y^{(i)} y ( i ) 的概率分布。w \boldsymbol{w} w 并不是随机变量,而是一个参数,所以用; ; ; 分号隔开。或者说w \boldsymbol{w} w 和x ( i ) \boldsymbol{x}^{(i)} x ( i ) 不是同一类变量,需要分开单独理解。P ( y ( i ) ∣ x ( i ) , w ) P(y^{(i)}|\boldsymbol{x}^{(i)}, \boldsymbol{w}) P ( y ( i ) ∣ x ( i ) , w ) 则有完全不同的意义,表示x ( i ) \boldsymbol{x}^{(i)} x ( i ) 和w \boldsymbol{w} w 同时发生时,y ( i ) y^{(i)} y ( i ) 的概率分布。
前文提到,似然函数是所观察到的各个样本发生的概率的乘积。一组样本有m m m 个观测数据,其中单个观测数据发生的概率为刚刚得到的公式,m m m 个观测数据的乘积如下所示。
L ( w ) = L ( w ; X , y ) = ∏ i = 1 m P ( y ( i ) ∣ x ( i ) ; w )
L(\boldsymbol{w}) = L(\boldsymbol{w}; \boldsymbol{X}, \boldsymbol{y}) = \prod_{i=1}^{m}P(y^{(i)}|\boldsymbol{x}^{(i)}; \boldsymbol{w})
L ( w ) = L ( w ; X , y ) = i = 1 ∏ m P ( y ( i ) ∣ x ( i ) ; w )
最终,似然函数可以表示成:
L ( w ) = ∏ i = 1 m 1 2 π σ 2 exp ( − ( y ( i ) − w ⊤ x ( i ) ) 2 2 σ 2 )
L(\boldsymbol{w}) = \prod_{i=1}^{m} \frac{1}{\sqrt{2 \pi \sigma^2}}\exp{({-\frac{(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2}{2 \sigma^2}})}
L ( w ) = i = 1 ∏ m 2 π σ 2 1 exp ( − 2 σ 2 ( y ( i ) − w ⊤ x ( i ) ) 2 )
其中,x ( i ) \boldsymbol{x}^{(i)} x ( i ) 和y ( i ) y^{(i)} y ( i ) 都是观测到的真实数据,是已知的,w \boldsymbol{w} w 是需要去求解的模型参数。
给定一组观测数据X \boldsymbol{X} X 和y \boldsymbol{y} y ,如何选择参数w \boldsymbol{w} w 来使模型达到最优的效果?最大似然估计法告诉我们应该选择一个w \boldsymbol{w} w ,使得似然函数L L L 最大。L L L 中的乘积符号和exp \exp exp 运算看起来就非常复杂,直接用L L L 来计算十分不太方便,于是统计学家在原来的似然函数基础上,取了log \log log 对数。log \log log 的一些性质能大大化简计算复杂程度,且对原来的似然函数增加log \log log 对数并不影响参数w w w 的最优值。通常使用花体的ℓ \ell ℓ 来表示损失函数的对数似然函数。
ℓ ( w ) = log L ( w ) = log ∏ i = 1 m 1 2 π σ 2 exp ( − ( y ( i ) − w ⊤ x ( i ) ) 2 2 σ 2 ) = ∑ i = 1 m log [ 1 2 π σ 2 ⋅ exp ( − ( y ( i ) − w ⊤ x ( i ) ) 2 2 σ 2 ) ] = ∑ i = 1 m log [ 1 2 π σ 2 ] + ∑ i = 1 m log [ exp ( − ( y ( i ) − w ⊤ x ( i ) ) 2 2 σ 2 ) ] = m log 1 2 π σ 2 − 1 2 σ 2 ∑ i = 1 m ( y ( i ) − w ⊤ x ( i ) ) 2
\begin{aligned}
\ell(\boldsymbol{w}) &= \log\ L(\boldsymbol{w}) \\
& = \log \ \prod_{i=1}^{m} \frac{1}{\sqrt{2 \pi \sigma^2}}\exp{({-\frac{(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2}{2 \sigma^2}})} \\
&= \sum_{i=1}^{m}\log[\frac{1}{\sqrt{2 \pi \sigma^2}}\cdot\exp{({-\frac{(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2}{2 \sigma^2}})}] \\
&= \sum_{i=1}^{m}\log[\frac{1}{\sqrt{2 \pi \sigma^2}}] + \sum_{i=1}^{m}\log[\exp{({-\frac{(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2}{2 \sigma^2}})}] \\
&= m\log{\frac{1}{\sqrt{2 \pi \sigma^2}}} - \frac{1}{2 \sigma^2}\sum_{i=1}^{m}(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x^{(i)}})^2
\end{aligned}
ℓ ( w ) = log L ( w ) = log i = 1 ∏ m 2 π σ 2 1 exp ( − 2 σ 2 ( y ( i ) − w ⊤ x ( i ) ) 2 ) = i = 1 ∑ m log [ 2 π σ 2 1 ⋅ exp ( − 2 σ 2 ( y ( i ) − w ⊤ x ( i ) ) 2 ) ] = i = 1 ∑ m log [ 2 π σ 2 1 ] + i = 1 ∑ m log [ exp ( − 2 σ 2 ( y ( i ) − w ⊤ x ( i ) ) 2 ) ] = m log 2 π σ 2 1 − 2 σ 2 1 i = 1 ∑ m ( y ( i ) − w ⊤ x ( i ) ) 2
上面的推导过程主要利用了下面两个公式:
log ( a b ) = log ( a ) + log ( b ) log e ( exp ( a ) ) = a
\log (ab) = \log(a) + \log(b) \\
\log_e (\exp(a)) = a
log ( a b ) = log ( a ) + log ( b ) log e ( exp ( a ) ) = a
由于log \log log 对数可以把乘法转换为加法,似然函数中的乘积项变成了求和项。又因为log \log log 对数可以消去幂,最终可以得到上述结果。
由于我们只关心参数w \boldsymbol{w} w 取何值时,似然函数最大,标准差σ \sigma σ 并不会影响w \boldsymbol{w} w 取何值时似然函数最大,所以可以忽略掉带有标准差σ \sigma σ 的项m log 1 2 π σ 2 m\log{\frac{1}{\sqrt{2 \pi \sigma^2}}} m log 2 π σ 2 1 。再在− 1 2 σ 2 ∑ i = 1 m ( y ( i ) − w ⊤ x ( i ) ) 2 -\frac{1}{2 \sigma^2}\sum_{i=1}^{m}(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2 − 2 σ 2 1 ∑ i = 1 m ( y ( i ) − w ⊤ x ( i ) ) 2 加个负号,负负得正,原来似然函数ℓ \ell ℓ 最大化问题就变成了最小化问题,其实最后还是最小化:
∑ i = 1 m ( y ( i ) − w ⊤ x ( i ) ) 2
\sum_{i=1}^{m}(y^{(i)} - \boldsymbol{w}^\top \boldsymbol{x}^{(i)})^2
i = 1 ∑ m ( y ( i ) − w ⊤ x ( i ) ) 2
这与最小二乘法所优化的损失函数几乎一样,都是“真实值 - 预测值”的平方和,可以说是殊途同归。
接下来对公式参数求解,线性回归的求解 一文已经详细推导过:
对目标函数求导,让导数为0,得到一个矩阵方程,矩阵方程的解即为模型的最优解
使用梯度下降法,求模型的最优解。
最小二乘与最大似然
前面的推导中发现,最小二乘与最大似然的公式几乎一样。直观上来说,最小二乘法是在寻找观测数据与回归超平面之间的误差距离最小的参数。最大似然估计是最大化观测数据发生的概率。当我们假设误差是正态分布的,所有误差项越接近均值0,概率越大。正态分布是在均值两侧对称的,误差项接近均值的过程等同于距离最小化的过程。
总结
最大似然估计是机器学习中最常用的参数估计方法之一,逻辑回归、深度神经网络等模型都会使用最大似然估计。我们需要一个似然函数来描述真实数据在不同模型参数下发生的概率,似然函数是关于模型参数的函数。最大似然估计就是寻找最优参数,使得观测数据发生的概率最大、统计模型与真实数据最相似。
参考资料
如何通俗地理解概率论中的「极大似然估计法」? - 马同学的回答 - 知乎 https://www.zhihu.com/question/24124998/answer/242682386
Andrew Ng:CS229 Lecture Notes
留下评论