机器学习八股文
Machine Leaning
ML 基础概念
-
Overfitting / Underfitting
过拟合指模型与数据的匹配程度过高,将训练数据一些正常的起伏、波动、异常值也当作是数据的特征,导致模型对新数据的泛化能力变差。具体的表现为在训练集上表现非常优秀,而在验证集/测试集中表现非常差。
- 解决过拟合的方法一般有:1) 适量减少特征的数量;2) 添加正则化项(Regularization)。正则化,顾名思义,目的是为了降低特征对于预测结果的影响能力。常见的正则化项有L1正则项和L2正则项。详见正则化。
欠拟合与过拟合相反,指的是模型缺乏足够的泛化能力。
- 解决欠拟合的方法有:1) 增加训练轮数;2) 增加模型特征;3) 减少正则项。
-
Bias / Variance trade-off
偏差(Bias)指模型预测结果与真实值的差异程度,描述了模型的拟合能力;方差(Varience)指模型面对不同数据集时的差异程度,描述了数据扰动对模型的影响。
一般来说,越简单模型的偏差越高,方差越低;越复杂模型的偏差越低,方差越高。这同样也对应着模型的过拟合与欠拟合。权衡偏差与方差的常见方法有交叉认证(Cross-Validation)。K折交叉验证的基本方法为:将训练集平均分为 k k k份,每次训练取其中一份作为验证集,剩下 k − 1 k-1 k−1份作为训练集,重复 k k k次,直到每一份小数据集都被作为过验证集。最终的损失为 k k k次训练的损失取平均。
正则化 Regularization
-
L1 vs L2
-
L1正则化,又称LASSO、L1范数,是所有参数的绝对值之和。
∥ x ∥ 1 = ∑ i = 1 m ∣ x i ∣ \lVert x \lVert_1=\sum_{i=1}^m\lvert x_i \lvert ∥x∥1=i=1∑m∣xi∣ -
L2正则化,又称Ridge,岭回归,是所有参数的平方和的平方根。
∥ x ∥ 2 = ∑ i = 1 m x i 2 \lVert x \lVert_2=\sqrt{\sum_{i=1}^m x_i^2} ∥x∥2=i=1∑mxi2
-
两种范数都有助于降低过拟合风险。L1范数可以用于特征选择,但不能直接求导,因此不能使用常规的梯度下降法/牛顿法等进行优化(常见方法有坐标轴下降法和 Lasso 回归法);L2范数方便求导。
-
-
L1范数的稀疏性 / 为何L1正则化可以用于特征选择?
L1范数相比于L2范数,更容易得到稀疏解,即L1范数可以将不重要的特征参数优化至0.
- 如何理解?
假设损失函数 L L L 与某个参数 x x x 的关系如下图所示:此时最优点位于绿色点处, x < 0 x<0 x<0.
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-g2CYlP6P-1631367951014)(imgs/l1vsl2_01.jpg)]
这时施加 L2 正则化,新的损失函数 ( L + C x 2 ) (L+Cx^2) (L+Cx2) 如下图蓝线所示,可以看到最优的 x x x 在黄点处, x x x 的绝对值减小了,但依然非零。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-SpLZ1mwd-1631367951016)(imgs/l1vsl2_02.jpg)]
而如果施加 L1 正则化,则新的损失函数 ( L + C ∣ x ∣ ) (L+C\lvert x \lvert) (L+C∣x∣) 如下图粉线所示,最优的 x x x 就变成了 0。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-2dgHoqZv-1631367951017)(imgs/l1vsl2_03.jpg)]
略加推导可以得到,当施加 L2 正则化时,当且仅当损失函数原本的导数为 0 时,损失函数才会在 x = 0 x=0 x=0 处极小;而施加 L1 正则化时,参数 C C C 与损失函数的导数仅需满足 C > ∣ L ∣ C>\lvert L \lvert C>∣L∣ 的关系, x = 0 x=0 x=0 便会成为损失函数的一个极小值点。
上面只分析了一个参数 x x x。事实上 L1 正则化会使得许多参数的最优值变成 0,使得模型变得稀疏。利用这样的特性,我们便可以使用L1正则化来帮助筛选特征。
机器学习中的评估指标 Metrics
-
Precision / Recall / F 1 F_1 F1 Score
对于二分类问题,我们常常使用精确率(Precision)、召回率(Recall)以及 F 1 F_1 F1 Score来评估二分类模型的性能。对于一个二分类器,在数据集上的预测情况可以分为以下4种:
- TP(True Positive),将正类正确预测为正类;
- TN(True Negative),将负类正确预测为负类;
- FP(False Positive),将负类错误预测为正类;
- FN(False Negative),将正类错误预测为负类;
有了以上概念,我们可以给出以下评估指标的定义:
- 精确率定义为:
P = T P T P + F P P=\frac{TP}{TP+FP} P=TP+FPTP
即在模型预测为正类的样本中,预测正确的比例。可以看到,精确率更加关注于模型认为是正类样本的结果。 - 召回率定义为:
R = T P T P + F N R=\frac{TP}{TP+FN} R=TP+FNTP
即在正类的样本中,模型预测正确的比例。相比之下,召回率更加关注于那些真实值为正类的样本。 - 此外,
F
1
F_1
F1 值定义为精确率与召回率的调和均值,即
2 F 1 = 1 P + 1 R \frac{2}{F_1}=\frac{1}{P}+\frac{1}{R} F12=P1+R1
F 1 = 2 × P × R P + R = 2 T P 2 T P + F P + F N F_1 = \frac{2 \times P \times R}{P + R} = \frac{2TP}{2TP+FP+FN} F1=P+R2×P×R=2TP+FP+FN2TP
当精确率和召回率都高时, F 1 F_1 F1 值也会高。
-
混淆矩阵 Confusion Matrix
分类结果的混淆矩阵如下表所示。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-jmzKbNu1-1631367951020)(imgs/ConfusionMatrix.jpeg)]
-
macro- F 1 F_1 F1 vs micro- F 1 F_1 F1
很多时候我们有多个二分类混淆矩阵(例如多次训练与测试 / 多个数据集 / 多分类任务中每两两类别的组合等),这是我们希望在 n n n 个二分类混淆矩阵上综合考察模型性能。
-
macro- F 1 F_1 F1
一种直接的做法是直接计算各个混淆矩阵的精确率和召回率,再计算平均值,分别得到 macro- P P P、macro- R R R和对应的macro- F 1 F_1 F1.
macro- P = 1 n ∑ i = 1 n P i , macro- R = 1 n ∑ i = 1 n R i , \text{macro-}P = \frac{1}{n}\sum_{i=1}^n P_i, \qquad \text{macro-}R = \frac{1}{n}\sum_{i=1}^n R_i, macro-P=n1i=1∑nPi,macro-R=n1i=1∑nRi,
macro- F 1 = 2 × macro- P × macro- R macro- P + macro- R \text{macro-}F_1 = \frac{2 \times \text{macro-}P \times \text{macro-}R}{\text{macro-}P + \text{macro-}R} macro-F1=macro-P+macro-R2×macro-P×macro-R -
micro- F 1 F_1 F1
另一种做法是先将各个混淆矩阵的对应元素进行平均,得到 T P ‾ \overline{TP} TP、 T N ‾ \overline{TN} TN、 F P ‾ \overline{FP} FP和 F N ‾ \overline{FN} FN,再基于这些值计算出micro- P P P、micro- R R R和对应的micro- F 1 F_1 F1.
micro- P = T P ‾ T P ‾ + F P ‾ , micro- R = T P ‾ T P ‾ + F N ‾ , \text{micro-}P = \frac{\overline{TP}}{\overline{TP}+\overline{FP}}, \qquad \text{micro-}R = \frac{\overline{TP}}{\overline{TP}+\overline{FN}}, micro-P=TP+FPTP,micro-R=TP+FNTP,
micro- F 1 = 2 × micro- P × micro- R micro- P + micro- R \text{micro-}F_1 = \frac{2 \times \text{micro-}P \times \text{micro-}R}{\text{micro-}P + \text{micro-}R} micro-F1=micro-P+micro-R2×micro-P×micro-R
-
-
ROC 曲线 / AUC 面积
ROC 曲线(Receiver Operating Characteristic)与 AUC (Area Under ROC Curve)是面对不平衡分类问题时最常用的评估指标。要了解 ROC 是什么,首先我们根据混淆矩阵再定义两个指标:True Positive Rate(TPR) 以及 False Positive Rate(FPR).
T P R = R = T P T P + F N , F P R = F P T N + F P , TPR = R = \frac{TP}{TP+FN}, \qquad FPR = \frac{FP}{TN+FP}, TPR=R=TP+FNTP,FPR=TN+FPFP,
正常来说,一个好的模型应该满足高 TPR 和低 FPR。对于任意一个训练好的模型,在给定测试数据上我们都能计算出它的 TPR 和 FPR。以 FPR 为横坐标,TPR 为纵坐标,我们可以将任意模型的一对 (FPR, TPR) 画在该坐标图中,如下图所示。同时我们将由该坐标轴构成的空间称为 ROC 空间。图1中假设有 A、B、C、D、E 共计五个模型。在 ROC 空间中,模型越靠近左上角,表明模型效果越好。[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-RJBQh1xv-1631367951022)(imgs/ROC_01.png)]
在二分类问题的大多数情况中(尤其是神经网络中),我们判定一个样本是正类和负类的依据是设置一个阈值,超过该阈值的样本被标记为正类,反之则为负类。一般而言,这个阈值被设置为0.5。那么如果我们尝试使用不同的阈值来划分正负类,我们就能得到多组 (FPR, TPR)。我们可以根据这些坐标近似地在 ROC 空间中画出一条曲线,即 ROC 曲线。只要 (FPR, TPR) 点足够多,我们可以计算出曲线下的面积,即 AUC面积,如下图所示。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-A1l3w7l1-1631367951024)(imgs/ROC_02.png)]
Loss与优化
-
凸优化问题
对于一个优化问题,如果其目标函数是凸函数,且可行域是凸集(集合中任意两点连线上的任意点都在集合内),那么它就是一个凸优化问题。
定义域 D \mathbb{D} D 是一个凸集的函数 f f f 是凸函数,当且仅当对于任意的 x , y ∈ D x,y \in \mathbb{D} x,y∈D 和 θ ∈ [ 0 , 1 ] \theta \in [0,1] θ∈[0,1],都有:
f ( θ x + ( 1 − θ ) y ) ≤ θ f ( x ) + ( 1 − θ ) f ( y ) f(\theta x+(1-\theta)y) \le \theta f(x)+(1-\theta) f(y) f(θx+(1−θ)y)≤θf(x)+(1−θ)f(y)
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-BnweXJgW-1631367951027)(imgs/convex_func.jpg)]
数学中强调凸优化问题的重要性,在于凸优化问题的局部最优解必然也是其全局最优解。这个特性使得我们可以使用贪心算法、梯度下降法、牛顿法等方法来求解凸优化问题。事实上,我们求解许多非凸优化问题,也是通过将其拆解为若干个凸优化问题,再分别进行求解。
-
MSE / MSELoss
均方误差 (Mean Square Error, MSE),是回归任务中最常见的度量指标。
E ( f ; D ) = ∑ i = 1 m ( f ( x i ) − y i ) 2 E(f;D)=\sum_{i=1}^m (f(x_i) - y_i)^2 E(f;D)=i=1∑m(f(xi)−yi)2
-
以 MSELoss 为损失函数的逻辑回归是凸优化问题吗?
不是。逻辑回归将线性模型通过 sigmoid 非线性函数映射为分类问题,其 MSE 是一个非凸函数,优化时可能得到局部最优解而得不到全局最优解,所以以 MSELoss 为损失函数的逻辑回归不是凸优化问题。
-
线性回归,最小二乘法与最大似然估计的关系?
求解线性回归常用的方法有最小二乘法 (OLS)和最大似然估计 (MLE)。
-
最小二乘法以预测值和真实值的平方和作为损失函数 (MSELoss)。
J ( w ) = ∑ i = 1 m ( h w ( x i ) − y i ) 2 J(w)=\sum_{i=1}^m (h_w(x_i) - y_i)^2 J(w)=i=1∑m(hw(xi)−yi)2
-
最大似然估计在已知 x x x 与 y y y 的情况下,以概率最大的角度,估计模型可能性最大的参数 h w h_w hw。设误差 ϵ i = y i − h w ( x i ) \epsilon_i = y_i - h_w(x_i) ϵi=yi−hw(xi), 由于 ϵ i \epsilon_i ϵi 符合高斯分布,可得概率密度函数:
p ( ϵ i ) = 1 σ 2 π e − ( ϵ i ) 2 2 σ 2 p(\epsilon_i) = \frac{1}{\sigma \sqrt{2\pi}} e^{-\frac{(\epsilon_i)^2}{2\sigma^2}} p(ϵi)=σ2π1e−2σ2(ϵi)2
将 ϵ i = y i − h w ( x i ) \epsilon_i = y_i - h_w(x_i) ϵi=yi−hw(xi) 代入,可得:
p ( y i ∣ h w ( x i ) ) = 1 σ 2 π e − ( y i − h w ( x i ) ) 2 2 σ 2 p(y_i | h_w(x_i)) = \frac{1}{\sigma \sqrt{2\pi}} e^{-\frac{(y_i - h_w(x_i))^2}{2\sigma^2}} p(yi∣hw(xi))=σ2π1e−2σ2(yi−hw(xi))2
则似然函数公式如下:
L ( h w ( x i ) ) = ∏ i = 1 m p ( y i ∣ h w ( x i ) ) = ∏ i = 1 m 1 σ 2 π e − ( y i − h w ( x i ) ) 2 2 σ 2 \begin{aligned} L(h_w(x_i)) &= \prod_{i=1}^m p(y_i | h_w(x_i))\\ &= \prod_{i=1}^m \frac{1}{\sigma \sqrt{2\pi}} e^{-\frac{(y_i - h_w(x_i))^2}{2\sigma^2}} \\ \end{aligned} L(hw(xi))=i=1∏mp(yi∣hw(xi))=i=1∏mσ2π1e−2σ2(yi−hw(xi))2
等号两边取对数,不影响函数的极值点。
log L ( h w ( x i ) ) = ∑ i = 1 m log 1 σ 2 π e − ( y i − h w ( x i ) ) 2 2 σ 2 = m log 1 σ 2 π − 1 2 σ 2 ∑ i = 1 m ( y i − h w ( x i ) ) 2 \begin{aligned} \log L(h_w(x_i)) &= \sum_{i=1}^m \log \frac{1}{\sigma \sqrt{2\pi}} e^{-\frac{(y_i - h_w(x_i))^2}{2\sigma^2}}\\ &= m \log \frac{1}{\sigma \sqrt{2\pi}} - \frac{1}{2\sigma^2} \sum_{i=1}^m (y_i - h_w(x_i))^2 \end{aligned} logL(hw(xi))=i=1∑mlogσ2π1e−2σ2(yi−hw(xi))2=mlogσ2π1−2σ21i=1∑m(yi−hw(xi))2
我们知道 h w ( x ) h_w(x) hw(x) 是关于权重 w w w 的函数,不妨设为 l ( w ) l(w) l(w)。因此有:
l ( w ) = m log 1 σ 2 π − 1 2 σ 2 ∑ i = 1 m ( y i − h w ( x i ) ) 2 \begin{aligned} l(w) = m \log \frac{1}{\sigma \sqrt{2\pi}} - \frac{1}{2\sigma^2} \sum_{i=1}^m (y_i - h_w(x_i))^2 \end{aligned} l(w)=mlogσ2π1−2σ21i=1∑m(yi−hw(xi))2
去除前面的常数项和常数系数,可以看到与最小二乘法的的公式一致,之后求解过程同最小二乘法。因此得出结论,最小二乘法与最大似然估计从两个不同的角度出发,得到一致的结果。
-
-
相对熵与交叉熵 Ralative-Entropy / Cross-Entropy
我们常用信息量来量化数据中的信息。设事件 x 0 x_0 x0 发生的概率为 p ( x 0 ) p(x_0) p(x0),则其信息量为:
I ( x 0 ) = − log p ( x 0 ) I(x_0) = -\log p(x_0) I(x0)=−logp(x0)
熵 (Entropy) 被用来度量一个系统的混乱程度,代表一个系统中所有事件信息量的期望。熵越大,该系统的不确定性也越大。
H ( X ) = − ∑ x ∈ X p ( x i ) log p ( x i ) H(X) = -\sum_{x \in X}p(x_i) \log p(x_i) H(X)=−x∈X∑p(xi)logp(xi)
相对熵 (Ralative Entropy),又称 KL 散度 (Kullback-Leibler Divergence),是两个随机分布 p p p 与 q q q 之间的对数差值的期望。
D K L ( p ∣ ∣ q ) = ∑ x ∈ X p ( x ) log p ( x ) q ( x ) = − ∑ x ∈ X p ( x ) [ log q ( x ) − log p ( x ) ] D_{KL}(p||q)=\sum_{x\in X} p(x)\log\frac{p(x)}{q(x)}=-\sum_{x\in X} p(x)[\log q(x) - \log p(x)] DKL(p∣∣q)=x∈X∑p(x)logq(x)p(x)=−x∈X∑p(x)[logq(x)−logp(x)]
交叉熵 (Cross-Entropy),与 KL 散度类似,是两个随机分布 p p p 与 q q q 之间距离的另一种度量。
C E H ( p , q ) = − ∑ x ∈ X p ( x ) l o g q ( x ) CEH(p,q)=−\sum_{x \in X}p(x)logq(x) CEH(p,q)=−x∈X∑p(x)logq(x)
为何在机器学习中常使用交叉熵而不是 KL 散度作为损失函数?
可以看到,相对熵、交叉熵之间存在以下关系:
D K L ( p ∣ ∣ q ) = C E H ( p , q ) − H ( p ) D_{KL}(p||q) = CEH(p,q) - H(p) DKL(p∣∣q)=CEH(p,q)−H(p)
在机器学习中,可以将 p p p 看作真实分布, q q q 为预测分布。则当 p p p 的分布已知时, H ( p ) H(p) H(p) 为常数,交叉熵与 KL 散度等价。
朴素贝叶斯 Naive Bayes
-
概率相关公式及贝叶斯定理
- 条件概率:事件 A 在另外一个事件 B 已经发生条件下的发生概率,即 p ( A ∣ B ) p(A|B) p(A∣B);
- 联合概率:事件 A 和事件 B 同时发生的概率,即 p ( A , B ) = p ( A ∣ B ) ∗ p ( B ) p(A, B) = p(A|B) * p(B) p(A,B)=p(A∣B)∗p(B);
- 全概率:若事件 B1, B2, …, Bn 构成一个完备事件组,即他们两两不相容,且和为全集,则对于任意事件 A 有:
p ( A ) = ∑ i = 1 n [ p ( A ∣ B i ) ∗ p ( B i ) ] p(A)=\sum^n_{i=1} [p(A|B_i)*p(B_i)] p(A)=∑i=1n[p(A∣Bi)∗p(Bi)] - 贝叶斯概率:在日常生活中,有时候我们难以直接求出
p
(
A
i
∣
B
)
p(A_i|B)
p(Ai∣B),但我们若已知
p
(
B
∣
A
i
)
p(B|A_i)
p(B∣Ai),
p
(
A
i
)
p(A_i)
p(Ai) 和
p
(
B
)
p(B)
p(B),则有:
p ( A i ∣ B ) = p ( B ∣ A i ) ∗ p ( A i ) p ( B ) = p ( B ∣ A i ) ∗ p ( A i ) ∑ j = 1 n p ( B ∣ A j ) ∗ p ( A j ) p(A_i|B) = \frac{p(B|A_i)*p(A_i)}{p(B)} = \frac{p(B|A_i)*p(A_i)}{\sum^n_{j=1}p(B|A_j)*p(A_j)} p(Ai∣B)=p(B)p(B∣Ai)∗p(Ai)=∑j=1np(B∣Aj)∗p(Aj)p(B∣Ai)∗p(Ai)
其中, p ( A i ∣ B ) p(A_i|B) p(Ai∣B) 被称为后验概率, p ( A i ) p(A_i) p(Ai) 被称为先验概率。
-
朴素贝叶斯分类器
- 设 x = { a 1 , a 2 , . . . , a m } x=\{a_1,a_2,...,a_m\} x={a1,a2,...,am} 为一个待分类项,其中 a i a_i ai 是 x x x 的特征属性。
- 有类别集合 C = { y 1 , y 2 , . . . , y n } C=\{y_1,y_2,...,y_n\} C={y1,y2,...,yn}。
- 对每一个类别 y i y_i yi,统计各个特征属性的条件概率,即 p ( a 1 ∣ y 1 ) p(a_1|y_1) p(a1∣y1), p ( a 2 ∣ y 1 ) p(a_2|y_1) p(a2∣y1), …, p ( a m ∣ y 1 ) p(a_m|y_1) p(am∣y1)。
- 根据贝叶斯公式,求得 p ( y i ∣ x ) = p ( x ∣ y i ) ∗ p ( y i ) p ( x ) p(y_i|x)=\frac{p(x|y_i)*p(y_i)}{p(x)} p(yi∣x)=p(x)p(x∣yi)∗p(yi)。
- 计算所有类别的 p ( y i ∣ x ) p(y_i|x) p(yi∣x),概率最大的 y k y_k yk 即为预测的类别。
- 该分类器之所以被称为“朴素”贝叶斯,是因为模型假设待分类项 x x x 的所有特征都是独立的事件。
- 常见的朴素贝叶斯分类器有:
- GaussianNB,该分类器用高斯分布来假设类别的先验概率分布,一般用于连续型数据。
- MultinomialNB,该分类器用多项式分布来假设类别的先验概率分布,用于多项式数据。
- BernoulliNB,该分类器用伯努利分布来假设类别的先验概率分布,用于二项分布数据。
-
朴素贝叶斯分类器的优缺点
朴素贝叶斯分类器的主要优点有:
- 模型发源于古典数学概论,算法比较简单,且有稳定的分类效率;
- 对小规模的数据表现好,适合多分类任务;
- 对缺失数据不敏感,例如文本分类等任务;
- 不存在过拟合的说法。
朴素贝叶斯分类器的缺点有:
- 朴素贝叶斯假设特征之间相互独立,但在现实中这个假设往往不成立;
- 朴素贝叶斯需要先估计先验概率,如果估计不准确容易影响分类结果;
- 基于概率的分类有可能会不准确。
-
Generative Model vs Discriminative Model 生成模型 / 判别模型
生成模型通过学习联合概率 P ( X , Y ) P(X,Y) P(X,Y),即特征 x x x 与类别 y y y 同时出现的概率,再对每一个类别求条件概率,取概率最大的类别作为预测结果,即 P ( Y ∣ X ) = P ( X , Y ) P ( X ) P(Y|X) = \frac{P(X,Y)}{P(X)} P(Y∣X)=P(X)P(X,Y)。
- 生成模型能学习到更多信息,如每个特征的边缘分布 p ( x ) p(x) p(x);
- 生成模型收敛速度快,且对小规模数据或稀疏数据表现较好;
- 生成模型不容易出现过拟合现象;
- 生成模型的效果一般没有判别模型好。
判别模型则是通过学习条件概率,即直接预测特征 x x x 下类别 y y y 的概率。
- 判别模型的分类边界更加灵活,能够拟合更加复杂的边界;
- 只用学习分类的信息,问题得到简化;
- 准确率普遍较生成模型较高。
支持向量机 Support Vector Machine, SVM
-
SVM 的核心思想?
SVM 是定义在特征空间中的线性分类器,目标是使得分隔超平面两边的数据点离超平面的间隔最大。当数据线性可分时,可以通过硬间隔最大化来学习一个线性的超平面;当数据近似线性可分时,则可以通过加上一个松弛变量,通过软间隔最大化来学习超平面;当线性不可分时,也可以使用核函数,将输入空间映射到高维的特征空间。
-
硬间隔最大化
首先定义超平面 ( w , b ) (w,b) (w,b) 关于样本点 ( x i , y i ) (x_i,y_i) (xi,yi) 的函数间隔为 γ ^ i = y i ( w x i + b ) \hat{\gamma}_i=y_i(wx_i+b) γ^i=yi(wxi+b) 但是,如果成比例地改变 w w w 和 b b b 的值,函数间隔会发生改变,但其超平面并没有变。为了解决这个问题,引入几何间隔为超平面与样本点的直线距离 γ i = y i ( w ∣ ∣ w ∣ ∣ x i + b ∣ ∣ w ∣ ∣ ) \gamma_i=y_i(\frac{w}{||w||}x_i+\frac{b}{||w||}) γi=yi(∣∣w∣∣wxi+∣∣w∣∣b)。硬间隔最大化的优化目标就是找到合适的超平面 ( w , b ) (w,b) (w,b),使得点到超平面的间隔最大,即:
max w , b γ , s.t., y i ( w ∣ ∣ w ∣ ∣ x i + b ∣ ∣ w ∣ ∣ ) ≥ γ , i = 1 , 2 , . . . , N \max_{w,b} \gamma, \text{s.t.,} y_i(\frac{w}{||w||}x_i+\frac{b}{||w||}) \geq \gamma, i=1,2,...,N w,bmaxγ,s.t.,yi(∣∣w∣∣wxi+∣∣w∣∣b)≥γ,i=1,2,...,N
我们可以将这个优化问题转化为:
min w , b 1 2 ∣ ∣ w ∣ ∣ 2 , s.t., y i ( w x i + b ) − 1 ≥ 0 , i = 1 , 2 , . . . , N \min_{w,b} \frac{1}{2}||w||^2, \text{s.t.,} y_i(wx_i+b)-1 \geq 0,i=1,2,...,N w,bmin21∣∣w∣∣2,s.t.,yi(wxi+b)−1≥0,i=1,2,...,N
可以使用拉格朗日乘子将该问题转换为对偶问题,这样可以便于求解:
min α 1 2 ∑ i N ∑ j N α i α j y i y j ( x i ⋅ x j ) − ∑ i N α i , s.t., ∑ i N α i y i = 0 , α i ≥ 0 , i = 1 , 2 , . . . , N \min_\alpha \frac{1}{2}\sum_i^N \sum_j^N \alpha_i \alpha_j y_iy_j(x_i\cdot x_j)-\sum_i^N \alpha_i, \text{s.t.,} \sum_i^N \alpha_i y_i=0, \alpha_i \geq 0,i=1,2,...,N αmin21i∑Nj∑Nαiαjyiyj(xi⋅xj)−i∑Nαi,s.t.,i∑Nαiyi=0,αi≥0,i=1,2,...,N
转换为对偶问题还可以便于引入核技巧来解决非线性问题,即将内积 ( x i ⋅ x j ) (x_i\cdot x_j) (xi⋅xj) 用核函数 K ( x i , x j ) = ϕ ( x i ) ⋅ ϕ ( x j ) K(x_i,x_j)=\phi(x_i)\cdot\phi(x_j) K(xi,xj)=ϕ(xi)⋅ϕ(xj) 来代替。 -
软间隔最大化
在现实中数据并不可能总是线性可分(可能存在噪音点、离群点等脏数据),线性不可分意味着有些样本点 ( x i , y i ) (x_i,y_i) (xi,yi) 不能满足函数间隔大于 1 的约束条件。为了使得模型更加鲁棒,我们需要将硬间隔转换为软间隔。对于每一个样本点,我们引进一个松弛变量 ξ i \xi_i ξi,使得函数间隔加上松弛变量大于 1,即 y i ( w x i + b ) ≥ 1 − ξ i y_i(wx_i+b)\geq1-\xi_i yi(wxi+b)≥1−ξi。
同时,添加一个惩罚参数 C C C,目标函数变为 1 2 ∣ ∣ w ∣ ∣ 2 + C ∑ i N ξ i \frac{1}{2}||w||^2 + C\sum_i^N \xi_i 21∣∣w∣∣2+C∑iNξi。惩罚参数的添加是使得函数间隔尽可能大,而误分类点尽可能少,惩罚参数是调和二者的系数。当 ξ i = 0 \xi_i = 0 ξi=0 时,样本点在间隔之外;当 0 < ξ i < 1 0 <\xi_i < 1 0<ξi<1 时,分类正确,样本点在超平面与间隔边界之间;当 ξ i > 1 \xi_i > 1 ξi>1 时,样本点被误分。
-
Hinge Loss 合页损失函数
Hinge Loss 的图像如图所示,横轴是函数间隔 γ ^ i = y i ( w x i + b ) \hat{\gamma}_i=y_i(wx_i+b) γ^i=yi(wxi+b),当间隔 γ ^ i ≥ 1 \hat{\gamma}_i\geq 1 γ^i≥1 时,表示正确分类,损失值为 0;当间隔 γ ^ i < 0 \hat{\gamma}_i < 0 γ^i<0 时分类错误;当间隔 0 < γ ^ i < 1 0 < \hat{\gamma}_i < 1 0<γ^i<1 时,分类正确,但样本点在间隔边界与超平面之间。Hinge Loss 对这样的样本点也会进行惩罚,提高了样本点的利用率,这使得 SVM 对训练样本数的依赖大大减少。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-bPdD0B3j-1631367951028)(imgs/hinge-loss.jpg)]
逻辑回归 Logistic Regression, LR
-
逻辑回归的核心思想?
逻辑回归主要用于分类问题,对于所给数据集,认为可以用一条直线将数据线性分类。逻辑回归与线性回归的最主要区别在于逻辑回归主要关注于分类概率与输入向量的关系,即 P ( Y = 1 ) P(Y=1) P(Y=1) 与 x x x 的直接关系,然后通过概率值来判断是否属于某一类别。
逻辑回归主要考虑二分类问题。给定数据集:
D = ( x 1 , y 1 ) , ( x 2 , y 2 ) , . . . , ( x N , y N ) , x i ∈ R n D=(x_1,y_1), (x_2,y_2),...,(x_N,y_N), x_i\in R^n D=(x1,y1),(x2,y2),...,(xN,yN),xi∈Rn
由于 w T x + b w^T x+b wTx+b取值连续,因此可以用它来拟合条件概率 p ( Y = 1 ∣ x ) p(Y=1|x) p(Y=1∣x). 最理想的函数其实是:
p ( Y = 1 ∣ x ) = { 0 , z < 0 0.5 , z = 0 1 , z > 0 , z = w T x + b p(Y=1|x)= \begin{cases} 0, & z < 0\\ 0.5, & z = 0 \\ 1, & z > 0 \end{cases} , z = w^T x + b p(Y=1∣x)=⎩⎪⎨⎪⎧0,0.5,1,z<0z=0z>0,z=wTx+b
但是,这个函数不可微,因此,我们使用 sigmoid 函数来拟合概率:
y = 1 1 + e z , z = w T x + b y=\frac{1}{1+e^z}, z=w^T x + b y=1+ez1,z=wTx+b
将 y y y 视为类后验概率估计,则重写公式有:
P ( Y = 1 ∣ x ) = 1 1 + e z , z = w T x + b , P(Y=1|x)=\frac{1}{1+e^z}, z=w^T x + b,\\ P(Y=1∣x)=1+ez1,z=wTx+b,
z = ln P ( Y = 1 ∣ x ) 1 − P ( Y = 1 ∣ x ) z = \ln\frac{P(Y=1|x)}{1-P(Y=1|x)} z=ln1−P(Y=1∣x)P(Y=1∣x)
因此,逻辑回归实际上是使用线性回归模型的预测值来逼近分类任务的对数几率,其优点有:
- 不仅能够预测出类别,还能预测出属于该类别的概率,对于一些需要预测概率的任务很适用;
- 对数几率函数在任意阶都是连续可导的凸函数,因此可以使用许多优化算法求解。
-
逻辑回归的损失函数与梯度
设: P ( Y = 1 ∣ x ) = p ( x ) , P ( Y = 0 ∣ x ) = 1 − p ( x ) P(Y=1|x) = p(x), P(Y=0|x) = 1 - p(x) P(Y=1∣x)=p(x),P(Y=0∣x)=1−p(x),则似然函数可以写为:
L ( w ) = ∏ [ p ( x i ) ] y i [ 1 − p ( x i ) ] 1 − y i L(w) = \prod[p(x_i)]^{y_i}[1-p(x_i)]^{1-y_i} L(w)=∏[p(xi)]yi[1−p(xi)]1−yi
为了方便求解,一般取对数似然函数:
l ( w ) = ln L ( w ) = ∑ [ y i ln p ( x i ) + ( 1 − y i ) ln ( 1 − p ( x i ) ) ] = ∑ [ y i ln p ( x i ) 1 − p ( x i ) + ln ( 1 − p ( x i ) ) ] = ∑ [ y i z i − ln ( 1 + e z i ) ] \begin{aligned} l(w) = \ln L(w) & = \sum[y_i\ln p(x_i) + (1-y_i)\ln (1-p(x_i))] \\ & = \sum[y_i\ln\frac{p(x_i)}{1-p(x_i)} + \ln(1-p(x_i))] \\ & = \sum[y_i z_i - \ln(1 + e^{z_i})] \end{aligned} l(w)=lnL(w)=∑[yilnp(xi)+(1−yi)ln(1−p(xi))]=∑[yiln1−p(xi)p(xi)+ln(1−p(xi))]=∑[yizi−ln(1+ezi)]
使用梯度下降法求解时,取似然函数的相反值进行优化,求其梯度为:
∂ J ( w ) ∂ w j = − ∑ i ∂ [ y i z i − ln ( 1 + e z i ) ] ∂ z i ⋅ ∂ z i ∂ w j = − ∑ i ( y i − p ( x i ) ) ⋅ x j \begin{aligned} \frac{\partial J(w)}{\partial w_j} & = -\sum_i \frac{\partial [y_i z_i - \ln(1 + e^{z_i})]}{\partial z_i} \cdot \frac{\partial z_i}{\partial w_j} \\ & = -\sum_i (y_i - p(x_i)) \cdot x_j \end{aligned} ∂wj∂J(w)=−i∑∂zi∂[yizi−ln(1+ezi)]⋅∂wj∂zi=−i∑(yi−p(xi))⋅xj
权重更新为:
w j : = w j + η ( ∑ i ( y i − p ( x i ) ) ⋅ x j ) , for i in range ( n ) w_j := w_j + \eta(\sum_i (y_i - p(x_i)) \cdot x_j),\text{ for }i\text{ in range}(n) wj:=wj+η(i∑(yi−p(xi))⋅xj), for i in range(n)
-
与 SVM 的区别?
- 同样是更加关注于分类边界附近的数据点,不同之处在于,LR 通过非线性映射减小离分类平面远的数据点的权重,提升离分类平面近的数据点的权重;而 SVM 则是完全只关注于分类错误的点以及靠近分类平面的数据点,而不关注那些离分类平面较远的普通样本。
- LR 是参数模型,而 SVM 是非参数模型。参数模型的前提是假设数据符合某种分布(如 LR 假设数据服从伯努利分布),该分布由一些参数决定。依赖数据分布会使得 LR 模型在样本不均衡时影响分类效果;而非参数模型不直接依赖于分布。
- LR 能够产生概率,而 SVM 不能直接产生概率。
决策树 Decision Tree
-
什么是决策树?
决策树以树状的结构对数据特征进行分类。建立一棵树,树上的每一个节点都是一个决策的判断,记为分支。将数据输入决策树,当在某节点数据不再分裂,则形成一个叶子节点,完成最终的分类。
-
如何确定分类的指标?
简要来说,我们在每次选择划分特征时,都会选择对分类最有帮助的特征来作为接下来的划分特征。具体来说,我们通过计算特征的熵来评价特征的重要程度。
熵 (Entropy) 被用来度量一个系统的混乱程度,代表一个系统中所有事件信息量的期望。
H ( X ) = − ∑ x ∈ X p ( x i ) log p ( x i ) H(X) = -\sum_{x \in X}p(x_i) \log p(x_i) H(X)=−x∈X∑p(xi)logp(xi)
H ( X ∣ A ) = − ∑ i = 1 d p ( A = a i ) H ( X ∣ A = a i ) H(X|A) = -\sum_{i=1}^d p(A=a_i)H(X|A=a_i) H(X∣A)=−i=1∑dp(A=ai)H(X∣A=ai)
熵越大,该系统的不确定性也越大。构造树的基本思想时随着树深度的增加,节点的熵迅速降低。熵降低的速度越快越好,以便生成一个高度尽可能矮的决策树。首先,我们根据分类的 label,计算不选择任何特征时,系统的熵值 H ( D ) H(D) H(D)。随后,我们分别计算考虑了各个特征之后,系统的熵值。我们计算每个特征的信息增益 G a i n ( D ∣ A ) = H ( D ) − H ( D ∣ A ) Gain(D|A) = H(D) - H(D|A) Gain(D∣A)=H(D)−H(D∣A),并选择信息增益最大的特征作为新节点的划分特征。
-
决策树的分类
- ID3。使用信息增益作为特征选择的评估方式。缺点主要是:当遇到稀疏特征时(如用户 ID),由于每项特征的样本比较少,可能会出现信息增益特别大的情况,导致决策树错误选择该稀疏特征作为划分特征。通常来说,ID3决策树偏向选择那些取值较多的特征属性。
- C4.5。使用信息增益率作为特征选择的评估方式。信息增益率计算方式如下:
G a i n R a t i o ( D ∣ A ) = G a i n ( D ∣ A ) H ( A ) = H ( D ) − H ( D ∣ A ) H ( A ) GainRatio(D|A) = \frac{Gain(D|A)}{H(A)} = \frac{H(D)-H(D|A)}{H(A)} GainRatio(D∣A)=H(A)Gain(D∣A)=H(A)H(D)−H(D∣A)
然而,信息增益率偏向选择取值较少的特征。因此,C4.5决策树会选择信息增益大于平均水平的特征来进行信息增益率的计算。此外,C4.5相对于ID3决策树,增加了对连续值和缺失值的处理过程。 - CART。对于分类问题,CART 决策树使用 Gini 系数作为评估标准;而对于回归问题,CART 决策树使用 MSE 作为评估标准。公式如下:
G i n i ( X ) = ∑ x ∈ X p ( x i ) ( 1 − p ( x i ) ) Gini(X) = \sum_{x\in X} p(x_i)(1-p(x_i)) Gini(X)=x∈X∑p(xi)(1−p(xi))
G i n i ( X ∣ A ) = ∑ i = 1 d p ( A = a i ) G i n i ( X ∣ A = a i ) Gini(X|A) = \sum_{i=1}^d p(A=a_i) Gini(X|A=a_i) Gini(X∣A)=i=1∑dp(A=ai)Gini(X∣A=ai)
-
随机森林 Random Forest, RF
使用 Bagging 算法,即有放回采样地选取 n n n 个样本,建立 m m m 个决策树分类器。多个分类器采用投票机制来产生最终的分类结果。
随机森林的“随机”有两重含义:- 样本随机选择,即 n < ∥ D ∥ n < \lVert D\rVert n<∥D∥,在每个决策树中选择一部分样本,有利于增加模型对离群点、噪音点的鲁棒性。
- 特征随机选择,每个决策树中只挑选一部分特征进行分类,有利于筛选不重要的特征或无效的特征。
XGBoost
-
Boosting 思想
当每个弱分类器的分类效果不理想时,我们可以尝试将多个不同的弱分类器组合起来,形成效果更好的强分类器。Boosting 着重于减小模型的 bias。Boosting 与 Bagging 的区别如下:
- 样本选择上,Bagging 在每次训练中是通过 Bootstrap 有放回选取不同的数据,每轮训练的数据集是独立的;Boosting 中每轮的训练集不变,只是训练集中每个样本的权重在训练中发生变化,样本的权重根据上一轮的分类结果进行调整。
- 模型预测上,Bagging 中每个模型的权重相等,通过投票或取均值的方式进行结合;Boosting 中每个弱分类器都有对应的权重,对于分类误差小的分类器分配更大的权重。
- 计算方式上,Bagging 的各个模型可以并行训练;Boosting 中每个模型只能串行训练,因为后一个模型的样本权重需要通过上一轮的训练结果来调整。
- Bagging 着重于减小模型的 Variance;Boosting 着重于减小模型的 bias。
为何说 Bagging 是减小 Variance,Bagging 是减小 bias?
在 Bagging 中,对样本重采样,每个子样本训练一个强(甚至过强)模型,最后取平均。由于子样本之间的相似程度高,多个模型的类型、结构与分布也非常相似,故多个模型有着近似的效果。而经过 Bagging 后,可以有效减小子模型的过拟合情况,因此可以显著降低模型的 Variance。
在 Boosting 中,把许多个弱分类器组合成一个强的分类器,着重于提高弱模型的性能,即着重降低模型的 bias。而 Variance 则不在 Boosting 的考虑范围之内。在 Boosting 中,多个模型有着同样的优化目标,并且通过在每一轮中的不断优化,能够达到降低 bias 的目的。 -
XGBoost 基本原理
在第 t t t 轮训练中,在保留前 t − 1 t-1 t−1 轮训练结果的前提下,加入一棵树 f t f_t ft,使得目标函数尽可能地降低。用公式表达如下:
O b j t = ∑ i = 1 n l ( y i , y ^ i t ) = ∑ i = 1 n l ( y i , y ^ i t − 1 + f t ( x i ) ) \begin{aligned} Obj_t & = \sum_{i=1}^n l(y_i, \hat{y}_i^t) \\ & = \sum_{i=1}^n l(y_i, \hat{y}_i^{t-1} + f_t(x_i)) \\ \end{aligned} Objt=i=1∑nl(yi,y^it)=i=1∑nl(yi,y^it−1+ft(xi))
设损失函数为 MSE,则原目标函数写为:
O b j t = ∑ i = 1 n ( y i − ( y ^ i t − 1 + f t ( x i ) ) ) 2 = ∑ i = 1 n [ 2 ( y ^ i t − 1 − y i ) f t ( x i ) + f t ( x i ) 2 ] + ∑ i = 1 n ( y i − y ^ i t − 1 ) 2 \begin{aligned} Obj_t &= \sum_{i=1}^n (y_i - (\hat{y}_i^{t-1} + f_t(x_i)))^2 \\ & = \sum_{i=1}^n[2(\hat{y}_i^{t-1} - y_i)f_t(x_i)+f_t(x_i)^2] + \sum_{i=1}^n ({y_i - \hat{y}_i^{t-1}})^2 \end{aligned} Objt=i=1∑n(yi−(y^it−1+ft(xi)))2=i=1∑n[2(y^it−1−yi)ft(xi)+ft(xi)2]+i=1∑n(yi−y^it−1)2
其中, ∑ i = 1 n ( y i − y ^ i t − 1 ) 2 \sum_{i=1}^n ({y_i - \hat{y}_i^{t-1}})^2 ∑i=1n(yi−y^it−1)2 与本轮无关,可以视为常数, ( y ^ i t − 1 − y i ) (\hat{y}_i^{t-1} - y_i) (y^it−1−yi) 一般被叫做残差,表示了上一轮预测值与真实值之间的差异,也是 XGBoost 算法在每一轮中预测的主要目标。即,将上一轮的训练结果看作一个整体,而新的一轮则对残差值进行预测。[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-rJT1oaOG-1631367951029)(imgs/XGBoost.png)]
此外,XGBoost 在每个叶子节点上都增加了正则化项 Ω ( f t ) = γ T + λ 1 2 ∑ j = 1 T w j 2 \Omega(f_t) = \gamma T + \lambda\frac{1}{2}\sum^T_{j=1} w_j^2 Ω(ft)=γT+λ21∑j=1Twj2,其中, T T T 代表叶子节点数量, λ 1 2 ∑ j = 1 T w j 2 \lambda\frac{1}{2}\sum^T_{j=1} w_j^2 λ21∑j=1Twj2 为 L2 正则化项。
Deep Learning
DL 基础概念
-
为什么神经网络需要偏置项?
对于神经网络中的每一个神经元,都有 y i = W T X i + b y_i = W^TX_i + b yi=WTXi+b。这个式子本质上就是要用这个函数在空间中划分决策面。而如果没有偏置项,那么划分的超平面就只能经过原点。偏置项的加入使得神经网络的拟合更加灵活,如果没有偏置项,训练可能难以收敛或出现其他 bug。
在所有场合都可以使用偏置项吗?
不是。例如在卷积层之后,如果要添加 Batch Normalization 层,最好不添加偏置项,因为不起作用,且会占用显卡内存。
在 BN 中,有一步关键操作为:
x i ^ = x i − μ B σ B 2 + ϵ \hat{x_i} = \frac{x_i - \mu_\mathcal{B}}{\sqrt{\sigma^2_{\mathcal{B}} + \epsilon}} xi^=σB2+ϵxi−μB
其中, μ B \mu_\mathcal{B} μB 为均值, σ B 2 \sigma^2_{\mathcal{B}} σB2 为方差。在该操作中,偏置项在计算中会被抵消掉,故偏置项不起作用。
-
Back Propagation
BP 神经网络是由一个输入层、一个输出层和若干个隐藏层构成的。输入信号从输入层进入,经过隐藏层计算,并由输出层输出。将输出的结果和真实值进行比对得到训练的误差。该误差沿着输出层,经过隐藏层,最终传播到输入层的权值参数。由于误差传播方向和训练方向相反,故称“反向传播”。
反向传播是为了解决神经网络无法直接应用梯度下降法的问题。由于梯度下降法只能用于“能够通过得到误差”的情况,例如逻辑回归。但隐藏层并不存在所谓“误差”,因此只能通过先将误差反向传播到隐藏层,应用链式法则得到求导函数,再使用梯度下降法进行优化。反向传播算法可以看作是梯度下降法在链式法则(Chain Rule)中的应用。
-
梯度消失和梯度爆炸问题
在反向传播的梯度更新中,若更新的梯度一直小于 0,就可能触发连乘效应,在之后的传播中越传越小,导致靠近输入层的权值几乎不更新,训练收敛速度变慢,这便是梯度消失。与之相反,若梯度过大则会触发梯度爆炸,以致于溢出,出现梯度为 NaN 的问题。
当激活函数为 Sigmoid 时,容易触发梯度消失问题。因为 Sigmoid 函数的导数最大值只有 0.25,如图所示。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-J8IOPS58-1631367951031)(imgs/sigmoid.jpg)]
常见的缓解梯度消失 / 梯度爆炸的方法有:
- 使用其他激活函数,如ReLU等;
- 用更合理的权值初始化方式,如 Xavier 初始化,He 初始化。这两种初始化方法都能保证在传播时权值的方差不变。
- Batch Normalization. 梯度的更新与 x x x 的值也有关系,因此用 BN 限制 x x x 的分布也有利于缓解梯度消失 / 梯度爆炸问题;
- 对权重进行正则化(L1、L2);
- 使用 ResNet 网络。ResNet 通过添加 Shortcut Connections,使得层与层之间可以跨层连接,减少了梯度消失 / 梯度爆炸的问题。
- 通过梯度截断(Gradient Truncation)手动防止梯度爆炸。
-
能不能将神经网络的所有权值都初始化为 0?
不能。事实上,不能将神经网络的所有权值都设置为同一值。否则,在神经网络的更新中,两权值的更新将一模一样。多个相同的神经元相当于只有一个神经元,会使得神经网络无法拟合。
因此,一般我们选择随机初始化,或是使用其他初始化方法,如 Xavier 初始化,He 初始化。这两种初始化方法都能保证在传播时权值的方差不变。
-
在深度学习中缓解过拟合问题
深度学习中防止过拟合常见的方法有:
- 获取更多、质量更高的数据
- 采集新的数据
- 数据增强(图片镜像、翻转等)
- 利用对抗网络生成数据
- 正则化(L1、L2)
- Dropout
- Early Stopping
- 集成学习,如 Bagging、Boosting 等。
- 获取更多、质量更高的数据
-
Dropout 是什么?
Dropout 是在每次训练过程中都随机舍弃一些神经元之间的连接。这样做可以降低对部分上层神经元的依赖关系,迫使模型去学习一些更具有鲁棒性的特征,使得模型泛化能力更强。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-WvRrQBTL-1631367951034)(imgs/dropout.png)]
-
Dropout 和 Batch Normalization 在训练和预测中的区别?
Dropout 在训练时采用,是为了减少神经元对部分上层神经元的依赖,减少过拟合的风险。而在预测中,应该用训练完成的模型,不需要 Dropout。
对于 BN,在训练时使用每一批数据的均值和方差进行计算,对每一批数据单独进行归一化。但在测试时,可能不存在 batch 的概念(例如预测单条数据)。因此在测试时,一般使用所有训练数据的均值和方差进行计算。
为什么不在训练的时候使用所有训练数据的均值和方差?
因为在训练中使用所有数据的均值和方差容易出现过拟合现象。
BN 的原理就是将每一批数据都归一到相同的分布。而每一批数据的均值和方差都不相同,这个差别能够增加模型的鲁棒性,在一定程度上减少模型的过拟合现象。
也正是因此,当应用 BN 时,一般要求将训练集完全打乱,并用一个较大的 batch size,否则,一批数据可能无法较好得代表训练集的分布,会影响模型训练的效果。
-
常见的 Non-Linear Activation Function 及其优缺点
非线性激活函数是神经网络与感知机网络最大的区别,即将非线性特性引入到网络中。如果不用非线性激活函数,则每一层都是上一层的线性变换,无论网络有多少层,输出都是输入的线性组合。而加入非线性层后,神经网络便拥有了学习非线性关系的能力,这使得神经网络可以逼近任意形状的函数。
函数名 函数表达式 优点 缺点 Sigmoid f ( z ) = 1 1 + e − z \displaystyle f(z)=\frac{1}{1+e^{-z}} f(z)=1+e−z1 1. 将输入转换为 (0, 1) 的区间。 1. 在神经网络反向传播中可能出现梯度消失问题;
2. 函数均值不为 0,使得权值总往同一方向更新,
收敛速度慢。tanh f ( z ) = e z − e − z e z + e − z \displaystyle f(z)=\frac{e^z-e^{-z}}{e^z+e^{-z}} f(z)=ez+e−zez−e−z 1. 将输入转换为 (0, 1) 的区间;
2. 解决 Sigmoid 均值非 0 问题。1. 依然存在梯度消失问题;
2. 幂函数计算复杂,训练时间大。ReLU f ( z ) = m a x ( 0 , x ) f(z)=max(0,x) f(z)=max(0,x) 1. 解决了梯度消失的问题;
2. 计算与收敛速度快。1. 函数均值不为 0;
2. Dead ReLU Problem. 当 x < 0 x<0 x<0 时梯度为 0,且
该神经元之后的神经元梯度永远为 0,即神经元直接
失效。通过合理的初始化,或是降低学习率来解决。Leaky ReLU f ( z ) = m a x ( α ⋅ x , x ) f(z)=max(\alpha \cdot x,x) f(z)=max(α⋅x,x) 1. 使得神经元在负数区域偏向于
激活而不是直接失效。1. 函数均值不为 0;
2. α \alpha α 的值需要确定且比较敏感(通常是一个非常小
的值,如 0.01)。 -
Gradient Descent vs Stochastic Gradient Descent vs Mini-Batch Gradient Descent
方法 优点 缺点 Gradient Descent 1. 参数梯度更新方向大致确定;
2. 适合并行化计算。1. 训练时收敛速度慢;
2. 当数据量大时,需要大量显存。Stochastic Gradient
Descent (SGD)1. 每次只随机抽取一条数据进行梯度更新,
花费代价小;
2. 适合大量数据的训练。1. 需要更多迭代次数;
2. 参数更新的过程震荡很大,参数更新方向有很大的波动;
3. 不适合并行化计算。Mini-Batch Gradient
Descent (MBGD)结合了前两种方法的优势:
1. 比 GD 收敛速度快,比 SGD 更加稳定;
2. 能利用高度优化的矩阵运算,适合并行化;1. 难以选择合适的学习率。太小的学习率会导致收敛缓慢;
太大的学习率会导致波动过大,可能跳出局部最优解。
可以采用动态调整学习率的方法(learning rate decay)。 -
常见的优化器及其对比
方法 特点 GD / SGD / MBGD 1. 难以选择合适的学习率。学习率太小会导致收敛缓慢;学习率太大会导致波动过大,可能跳出局部最优解。
2. 每个参数的学习率都是相同的。如果数据是稀疏的,且不同特征的出现频率相差较大,则会出现部分特征
学习不足的问题;
3. 在训练中容易陷入鞍点,即局部最优点,在这些点的梯度为 0,无法继续训练。Momentum 1. 借鉴了物理中的动量概念,模拟物体运动时的惯性,即更新的时候在一定程度上保留之前更新的方向,而不
是像 GD 算法一样完全按照新的梯度方向更新。这样可以增加稳定性,并且有一定的摆脱局部最优解的能力。
2. Momentum 算法会观察上一步的梯度,若当前梯度方向与历史梯度一致,则增强该方向的梯度,否则则削
弱该方向的梯度。AdaGrad 1. 针对 GD 算法中对于每个参数都保持同一学习率的问题,AdaGrad 算法能在训练中自动对不同参数的学习率
进行调整。对于出现频率比较低的特征,加大更新的学习率;对于出现频率高的,减小学习率。
2. 由于这个特性,AdaGrad 非常适合用于处理稀疏的数据。RMSprop 1. Root Mean Square prop. 对 AdaGrad 算法的改进,把 AdaGrad 的将历史梯度相加变成对历史梯度求均值;
2. 这种方法可以缓解 AdaGrad 算法学习率下降较快的问题。Adam 1. Adam 算法结合了 AdaGrad 和 RMSprop 的优点,即动态更新参数的学习率。不同于 RMSprop 只参考了参
数的历史平均值,Adam 同时参考了梯度的平均值和方差。
2. 在各大机器学习库中,两次估计的衰减率默认值 β 1 \beta_1 β1 和 β 2 \beta_2 β2 分别为 0.9 和 0.999.AdamW 1. 针对 Adam 算法中先进行梯度衰减再进行正则化,使得梯度大的参数无法正常被正则化的问题,在 AdamW
中将梯度衰减的步骤移到正则化后,解决了这一问题。 -
如何正确使用迁移学习?
通过使用之前在大数据集上经过训练的预训练模型,我们可以直接使用相应的结构和权重,将它们应用到我们正在面对的问题上。
-
场景一:现有数据集***大***,数据与原数据相似度***高***
这是最理想的情况,采用预训练模型会变得非常高效。最好的运用方式是保持模型原有的结构和初始权重不变,随后在新数据集的基础上重新训练 / 微调。
-
场景二:现有数据集***小***,数据与原数据集相似度***高***
在这种情况下,由于数据和原数据集相似度高,我们不需要重新训练模型,只需要将输出层改为新问题的结构即可。
-
场景三:现有数据集***大***,数据与原数据集相似度***低***
因为实际数据与预训练模型的训练数据之间存在很大差异,采用预训练模型将不会是一种高效的方式。因此最好的方法还是只沿用预训练模型的结构。将预处理模型中的权重全都初始化后,在新数据集上重新开始训练。
-
场景四:现有数据集***小***,数据与原数据集相似度***低***
这是最糟糕的一种情况。为了防止过拟合,我们不能从头开始训练。我们可以利用预训练模型较低的层进行特征提取,弥补数据集大小不足的缺陷,再利用较高的层进行训练(一般而言,神经网络较高的层具有较高的区分度,更适合用来训练数据本身)。因此,我们冻结预训练模型前 k k k 层的权重,用于提取数据的特征,然后训练后 n − k n-k n−k 层,并将原输出层改为新问题的结构。
-
Natural Language Processing
NLP 基本概念
-
常见的文本相似度计算方法
- 欧式距离,用于计算两等长句子向量的相似度。 distance = ( A − B ) ∗ ( A − B ) T \text{distance} = \sqrt{(A-B)*(A-B)^T} distance=(A−B)∗(A−B)T;
- 余弦距离,用于计算两等长句子向量的相似度。 distance = A ∗ B T ∣ A ∣ ∗ ∣ B ∣ \text{distance} = \frac{A*B^T}{|A|*|B|} distance=∣A∣∗∣B∣A∗BT;
- Jaccard 相似度。将句子看作单词的集合。则 A 与 B 的 Jaccard 相似度为: similarity = ∣ A ∩ B ∣ ∣ A ∪ B ∣ \text{similarity} = \frac{|A\cap B|}{|A\cup B|} similarity=∣A∪B∣∣A∩B∣;
- TF-IDF。TF 是词频 (Term Frequency),表示在一个文章中某个单词出现的频率;IDF 是逆文本频率指数 (Inverse Document Frequency),表示含有该关键词的文档占所有文档的比例。TF-IDF 建立在以下假设上:对区别文档最有意义的词语应该是那些在文档中出现频率高,而在整个文档集合的其他文档中出现频率少的词语;
- 最小编辑距离。一种经典的距离计算方法,用来度量字符串之间的差异。将字符串 A 不断修改(增删改),直至成为字符串 B,所需要的修改次数代表了字符串 A 和 B 的差异大小。常使用动态规划来计算最小编辑距离。
-
word2vec 模型
在 NLP 中,我们希望用一个数学形式表示不同的单词,于是便有了词向量。最初的词向量是 one-hot 词向量,但这种向量维度过大,非常稀疏,且不能反映词与词之间的关系。于是便有了分布式词向量,即固定 embedding 的维度,embedding 中的每一个值都是通过计算不同单词的贡献得到的。
训练 word2vec 模型主要有两种方式:CBOW 和 Skip-Gram。
- CBOW 是让模型根据某个词前面的 C 个词和之后的 C 个词,预测这个词出现的概率。如图,训练过程其实就是学习这两个矩阵
W
W
W 和
W
′
W'
W′,其中,
W
W
W 矩阵又被叫做 lookup table,即所有词嵌入向量的词表。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Rph14TiM-1631367951035)(imgs/word2vec-CBOW.jpg)] - Skip-Gram 和 CBOW 相反,是根据某一个词来预测它的前 C 个词和后 C 个词。同样训练两个矩阵
W
W
W 和
W
′
W'
W′,其中,
W
W
W 矩阵是 lookup table。一般来说,Skip-Gram 的训练时间比 CBOW 要慢。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-9vp8ZOKQ-1631367951036)(imgs/word2vec-skip-gram.jpg)]
为了加快训练速度,word2vec 采用了两种优化方式。
- Hierarchical Softmax,用霍夫曼树代替神经网络,本质上是将 n 分类变成 log(n) 次二分类。
- Negative Sampling,由于霍夫曼树中高频词离根结点较近,但是如果中心词是较生僻的词,那么就要消耗很长时间。简要来说就是从负样本中选取一部分来更新,而不是更新全部的权重。
- CBOW 是让模型根据某个词前面的 C 个词和之后的 C 个词,预测这个词出现的概率。如图,训练过程其实就是学习这两个矩阵
W
W
W 和
W
′
W'
W′,其中,
W
W
W 矩阵又被叫做 lookup table,即所有词嵌入向量的词表。
-
GloVe 模型
GloVe 模型利用了词语的共现频率来计算相关性。首先引入词语的共现矩阵 X X X,其中 X i j X_{ij} Xij 是在 word i 的上下文中 word j 的出现次数, X i = ∑ k X i k X_i = \sum_k X_{ik} Xi=∑kXik 是出现在 word i 的上下文中所有词的出现次数,则共现概率为 P i j = P ( j ∣ i ) = X i j X i P_ij = P(j|i) = \frac{X_{ij}}{X_i} Pij=P(j∣i)=XiXij,是word j 出现在 word i 上下文的概率。可以发现,共现概率的比例可以反映两个词的相关度。
HMM / CRF
-
隐马尔可夫模型 Hidden Markov Model, HMM
- 马尔可夫链,即一个状态序列,满足在任意时刻 t t t 的状态仅与其前一时刻的状态有关。隐马尔可夫链,则是无法直接观测到某一时刻的状态,而是要通过其他的观测状态才能预测隐藏的状态。
- 隐马尔可夫模型的两个基本假设:
- 齐次性假设:即隐藏的马尔可夫状态在任意时刻的状态只依赖于前一时刻的状态,与其他时刻的状态无关;
- 观测独立性假设:任意时刻的观测状态只取决与当前状态的隐藏状态,和其他时刻的观测状态或隐藏状态无关。
- 隐马尔可夫模型的五个要素:
- 隐藏状态集 Q Q Q = { q 1 q_1 q1, q 2 q_2 q2, …, q N q_N qN},即隐藏节点只能取值于该集合中的元素。
- 观测状态集 V V V = { v 1 v_1 v1, v 2 v_2 v2, …, v M v_M vM},即观测节点的状态也只能取值于该集合中的元素。
- 隐藏状态转移矩阵 A A A = [ a i j ] N × N [a_{ij}]_{N\times N} [aij]N×N,表示从一种隐藏状态到另一种隐藏状态的转移概率。
- 观测概率矩阵 B B B = [ b i j ] N × M [b_{ij}]_{N\times M} [bij]N×M,表示对于某一种隐藏状态,其观测状态的分布概率。
- 初始隐藏状态概率 π \pi π = [ p 1 , p 2 , . . . , p n ] [p_1, p_2, ..., p_n] [p1,p2,...,pn],表示初始时刻处于各个隐藏状态的概率。
- 隐马尔可夫模型要求解的基本问题:
- 概率计算问题。对于已知模型 λ \lambda λ = ( A , B , π ) (A, B, \pi) (A,B,π),和已知观测序列 O O O = { o 1 o_1 o1, o 2 o_2 o2, …, o M o_M oM},求产生这种观测序列的概率是多少,即求 p ( O ∣ λ ) p(O|\lambda) p(O∣λ)。
- 学习问题。对于已知观测序列 O O O = { o 1 o_1 o1, o 2 o_2 o2, …, o M o_M oM},求解模型 λ \lambda λ = ( A , B , π ) (A, B, \pi) (A,B,π) 的参数,使得产生这种观测序列的概率 p ( O ∣ λ ) p(O|\lambda) p(O∣λ) 最大,即用最大似然估计方法估计模型的参数。
- 解码问题。同样对于已知模型 λ \lambda λ = ( A , B , π ) (A, B, \pi) (A,B,π),和已知观测序列 O O O = { o 1 o_1 o1, o 2 o_2 o2, …, o M o_M oM},求解最优的隐藏状态序列 I I I = { i 1 i_1 i1, i 2 i_2 i2, …, i N i_N iN},使得 p ( I ∣ O ) p(I|O) p(I∣O) 最大。
- 对于基本问题的解法:
-
对第一个问题的解法:
- 暴力解法:时间复杂度为 O ( T N T ) O(TN^T) O(TNT);
- 前向算法:利用动态规划思想,将前面时刻计算过的概率保存下来。
- 对于第一个时刻的状态: a 1 ( i ) = π i b i ( o 1 ) a_1(i) = \pi_ib_i(o_1) a1(i)=πibi(o1), i ∈ [ 1 , N ] i\in [1,N] i∈[1,N];
- 对于第 t t t 个时刻的状态: a t ( i ) = [ ∑ j = 1 N a t − 1 ( j ) a j i ] b i ( o t ) a_t(i) = [\sum_{j=1}^N a_{t-1}(j)a_{ji}]b_i(o_t) at(i)=[∑j=1Nat−1(j)aji]bi(ot)。
-
对第二个问题的解法:
Baum-Welch 算法:与 EM 算法相似,在 E-step 中,计算联合分布 P ( O , I ∣ λ ) P(O,I|\lambda) P(O,I∣λ) 和条件分布 P ( I ∣ O , λ ˉ ) P(I|O,\bar{\lambda}) P(I∣O,λˉ),根据联合分布和条件分布计算期望表达式 L ( λ , λ ˉ ) L(\lambda,\bar{\lambda}) L(λ,λˉ);在 M-step 中最大化 λ \lambda λ 的值,使得 λ ˉ = arg max λ L ( λ , λ ˉ ) \bar{\lambda} = \argmax_\lambda L(\lambda,\bar{\lambda}) λˉ=λargmaxL(λ,λˉ)。
-
对第三个问题的解法:
Viterbi 维特比算法:可以看作一个求最长路径的动态规划算法。
初始化两个状态变量: δ t ( i ) \delta_t(i) δt(i) 表示在 t t t 时刻隐藏状态为 i i i 的所有状态转移路径中概率最大值,初始化 δ 1 ( i ) = π i b i ( o 1 ) \delta_1(i) = \pi_i b_i(o_1) δ1(i)=πibi(o1)。 ψ t ( i ) \psi_t(i) ψt(i) 则是在 t t t 时刻使得隐藏状态为 i i i 的转移路径中概率最大的前一时刻的隐藏状态,初始化为 0。则两状态变量的递推表达式为:
{ δ t ( i ) = max 1 ≤ j ≤ N [ δ t − 1 ( j ) a j i ] b i ( o t ) ψ t ( i ) = arg max 1 ≤ j ≤ N [ δ t − 1 ( j ) a j i ] \begin{cases} \delta_t(i) = \max_{1\leq j\leq N}[\delta_{t-1}(j)a_{ji}]b_i(o_t) \\ \psi_t(i) = \argmax_{1\leq j\leq N}[\delta_{t-1}(j)a_{ji}] \end{cases} ⎩⎨⎧δt(i)=max1≤j≤N[δt−1(j)aji]bi(ot)ψt(i)=1≤j≤Nargmax[δt−1(j)aji]
在第 T T T 步,取 δ T ( i ) \delta_T(i) δT(i) 最大值即为最可能隐藏序列出现的概率,此时最大的 ψ T ( i ) \psi_T(i) ψT(i) 即为第 T T T 的状态。
随后,从第 T T T 步开始回溯,即 i t − 1 ∗ i^*_{t-1} it−1∗ = ψ t ( i t ∗ ) \psi_t(i^*_t) ψt(it∗),得到完整的隐藏序列 I = ( i 1 ∗ , i 2 ∗ , . . . , i T ∗ ) I=(i^*_1, i^*_2, ..., i^*_T) I=(i1∗,i2∗,...,iT∗)。
-
-
条件随机场 Conditional Random Field, CRF
首先介绍随机场。一组随机变量,分布在同一个样本空间,那么它们就组成了一个随机场。我们希望利用这些随机变量之间的关系来解决实际问题。
马尔可夫随机场是满足马尔可夫独立性的随机场,即每个节点仅与其相邻的节点有关系。并不像贝叶斯网络(有向无环图)一样,通过变量之间的条件分布建模(节点与节点之间的依赖关系),马尔可夫随机场是根据变量之间的联合分布来建模的。当知道了变量之间的联合分布,则它们之间的条件分布也就能求解出来。因此,马尔可夫随机场是生成模型。
条件随机场则是对条件概率建模。即已经观测到部分点的前提下,求解整张图的分布。
HMM 是通过对可观测序列和隐藏序列的联合概率建模,估计模型的隐含变量的分布,再计算概率最大的隐藏序列,是一个生成模型。CRF 则是直接对条件概率 P ( I ∣ O ) P(I|O) P(I∣O) 建模,通过可观测序列直接判别出隐藏序列,是一个判别模型。
比较常见的条件随机场是线性链条件随机场。设 X = ( X 1 , X 2 , . . . , X n ) X = (X_1, X_2, ..., X_n) X=(X1,X2,...,Xn), Y = ( Y 1 , Y 2 , . . . , Y n ) Y = (Y_1, Y_2, ..., Y_n) Y=(Y1,Y2,...,Yn) 均为线性链表示的随机变量序列。在给定 X X X 的情况下,随机变量序列 Y Y Y 的条件概率分布构成线性链条件随机场,即 Y Y Y 的分布只与其相邻的节点有关。
与 HMM 相似,条件随机场主要求解的也是三个问题:
- 概率计算问题:给定条件随机场 P ( Y ∣ X ) P(Y|X) P(Y∣X), 观测序列 x x x 和隐藏序列 y y y,计算条件概率 P ( Y i = y i ∣ x ) P(Y_i=y_i|x) P(Yi=yi∣x),可以通过前向后向解法求解
- 学习问题:已知观测序列和隐藏序列,通过极大似然估计法来学习模型的最大概率参数。
- 预测问题:给定条件随机场 Y = ( Y ∣ X ) Y=(Y|X) Y=(Y∣X) 和观测序列 x x x,求条件概率最大的隐藏序列 y ∗ y^* y∗,即对观测序列进行标注。预测问题的常用算法是维特比算法。
CRF 是一个序列化标注算法,接收一个输入序列 X = ( x 1 , x 2 , . . . , x n ) X = (x_1, x_2, ..., x_n) X=(x1,x2,...,xn),输出目标序列 Y = ( y 1 , y 2 , . . . , y n ) Y = (y_1, y_2, ..., y_n) Y=(y1,y2,...,yn),可以被看作是一个 Seq2Seq 模型。在词性标注任务中,输入为文本序列,输出则为对应的词性序列。
相比于 HMM 需要对状态转移矩阵和观测概率矩阵建模,CRF 属于判别模型,其直接对 P ( I ∣ O ) P(I|O) P(I∣O) 建模:
P ( I ∣ O ) = 1 Z ( O ) e ∑ i T ∑ k M λ k f k ( O , I i − 1 , I i , i ) P(I|O) = \frac{1}{Z(O)}e^{\sum_i^T \sum_k^M \lambda_k f_k(O, I_{i-1}, I_i, i)} P(I∣O)=Z(O)1e∑iT∑kMλkfk(O,Ii−1,Ii,i)
其中,下标 i 表示当前所在的节点(token)位置,下标 k 表示第 k 个特征函数,并且每个特征函数都附属一个权重 λ k \lambda_k λk, 1 Z ( O ) \frac{1}{Z(O)} Z(O)1 是归一化系数。
RNN / LSTM
-
为什么需要 RNN?
循环神经网络(Recurrent Neural Network, RNN)。当给定的数据是序列型数据,如文本、音频等数据时,我们往往希望模型能够学习到给定数据的上下文环境。例如,对于一个句子,序列模型试图从同一个句子前面的单词推导出关系。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-glizhUB6-1631367951037)(imgs/rnn.png)]
在循环神经网络的每一个时间步骤(time step)中,我们取一个输入 x i x_i xi 和上一个节点的权值 h i − 1 h_{i-1} hi−1 作为输入,并产生一个输出 y i y_i yi 和权值 h i h_i hi,这个权值又被传递到下一个时间步骤,直到输入序列被读取完毕。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-cEkIJBuo-1631367951039)(imgs/multi-tasks.jpeg)]
普通的 RNN(Vanilla RNN)常使用 BP 算法来训练权值,但由于梯度消失 / 梯度爆炸问题,RNN 会丧失学习远距离信息的能力。为了解决远距离依赖问题,提出了 LSTM(Long Short-Term Memory)。
-
LSTM 网络
LSTM(Long Short-Term Memory)相对于普通 RNN 网络,能够显著的缓解长期依赖关系丢失的问题。LSTM 的主要思想是利用门结构来去除或添加单元之间信息传递的能力。LSTM 拥有三个门,来保护和控制单元状态,分别为遗忘门、输入门和输出门。
-
遗忘门
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-mxip4YKk-1631367951040)(imgs/forget-gate.png)]
第一步是决定从上一个单元中保留多少消息。将上一单元的状态 h t − 1 h_{t-1} ht−1 和这一层的输入 x i x_i xi 经过 sigmoid 层,输出一个 0-1 的值,代表要从上一层的单元状态保留多少信息。 -
输入门
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-bQrE3D6E-1631367951041)(imgs/input-gate.png)]
这一步是决定在这一层的单元状态中保留多少信息。将上一单元的状态 h t − 1 h_{t-1} ht−1 和这一层的输入 x i x_i xi 分别经过 sigmoid 层和 tanh 层,得到一个候选的单元状态 C ~ t \tilde{C}_t C~t。[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-HqxiKGdk-1631367951042)(imgs/update-cell.png)]
随后,根据遗忘门得到的遗忘比例 f t f_t ft 和这一层要输入的单元状态 C ~ t \tilde{C}_t C~t,得到这一层的最终单元状态 C t = f t ∗ C t − 1 + i t ∗ C ~ t C_t = f_t*C_{t-1} + i_t*\tilde{C}_t Ct=ft∗Ct−1+it∗C~t。 -
输出门
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-zrEEAcyQ-1631367951043)(imgs/output-gate.png)]
最终,我们需要决定这一层的单元的输出状态。将上一单元的状态 h t − 1 h_{t-1} ht−1 和这一层的输入 x i x_i xi 经过 sigmoid 层,确定要输出的部分 o t o_t ot,再将这一层的单元状态 C t C_t Ct 经过 tanh 层,再与 o t o_t ot 结合,得到最终的输出状态 h t h_t ht。
-
-
GRU 网络
与 LSTM 对比,GRU 网络更加简单,训练更加高效。GRU 去除了单元状态,将 LSTM 的 3 个门减少到 2 个,分别为更新门和重置门,分别决定了应该让多少信息通过单元,以及应该丢弃多少信息。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-5RlanzpX-1631367951044)(imgs/gru.jpg)]
-
如何计算 LSTM 和 GRU 的参数量?
一个单元内一共有四个非线性门 ( W [ h t − 1 , x t ] + b W[h_{t-1},x_t] + b W[ht−1,xt]+b),每一个门内的可训练变量包括一个矩阵 W W W 和一个置项 b b b。
因此,一个 LSTM 非线性门的参数即为 (embed_size + hidden_size) * hidden_size +hidden_size,LSTM 四个门的总参数量为 ((embed_size + hidden_size) * hidden_size +hidden_size) * 4。
同理,一个 GRU 单元的参数量为 ((embed_size + hidden_size) * hidden_size + hidden_size) * 3。
-
RNN / LSTM 的局限性
- 对于 RNN 来说,在访问一个单元前需要遍历之前所有的单元,使得在长文本下极度容易出现梯度消失问题。
- LSTM 利用门机制稍微缓解了 RNN 的梯度消失问题,但在超长文本前仍然存在该问题。
- 在 LSTM 中,在每一个单元中都有 4 个 MLP 层,需要消耗大量的计算资源,且模型本身不利于并行化。
TextCNN
-
如何卷积?
输入一个长度为 s s s 的句子,将其分词后映射到词向量,假设词向量的维度为 d d d,那么该句子可以表示为一个 s × d s\times d s×d 的矩阵。将该矩阵看作一张图像,用卷积神经网络提取特征。
文本卷积和图像卷积的区别在与文本序列只在垂直方向做卷积,即卷积核的宽度固定为词向量的维度 d d d。
-
TextCNN 的优缺点?
- 优点:网络结构简单,训练速度快,适合进行并行化,对短文本效果好;使用 Max-Pooling,便于提取最关键信息,因此适用于文本分类等任务。
- 缺点:全局 Max-Pooling 丢失了结构信息,很难发现文本中的依赖关系;只能学习到关键词是什么,无法学习到关键词的频率和顺序。
Attention Mechanism / Transformer
-
Seq2Seq 中的 Attention 机制
在 Seq2Seq 中,我们使用 encoder 将输入文本转化为一个定长向量,再用 decoder 将该向量转换为输出文本。但是,在面对长文本时,我们很难在定长向量中保留完整的输入文本信息,因此在 decode 时会存在信息丢失的问题。为了缓解这个问题,我们引入了 Attention 机制。
以机器翻译任务举例,在翻译到某一个单词时,我们希望能够注意到这个单词所对应的上下文,并结合之前已翻译的部分作出相应的翻译。这样,我们在 decoder 中就可以注意到输入文本的全部信息,而不只局限于那个定长的向量。
Attention 的计算过程如下:
- 得到 encoder 中的 hidden state h e → = ( h 1 , h 2 , . . . , h n ) \overrightarrow{h_e} = (h_1, h_2, ..., h_n) he=(h1,h2,...,hn)。
- 假设当前翻译到的 decoder state 为
s
t
−
1
→
\overrightarrow{s_{t-1}}
st−1,则可以计算该状态与输入的每一个单元
h
j
h_j
hj 状态的关联性
e
t
j
=
a
(
s
t
−
1
,
h
j
)
e_{tj} = a(s_{t-1},h_j)
etj=a(st−1,hj),写成向量形式则为
e
t
→
=
a
(
s
t
−
1
→
,
h
→
)
\overrightarrow{e_t} = a(\overrightarrow{s_{t-1}}, \overrightarrow{h})
et=a(st−1,h),其中,
a
a
a 是相关性的计算,常见的计算方式有:
- 直接点乘 a ( s t − 1 , h j ) = s t − 1 → T ⋅ h → a(s_{t-1},h_j)=\overrightarrow{s_{t-1}}^T\cdot\overrightarrow{h} a(st−1,hj)=st−1T⋅h;
- 加权点乘 a ( s t − 1 , h j ) = s t − 1 → T ⋅ W ⋅ h → a(s_{t-1},h_j)=\overrightarrow{s_{t-1}}^T\cdot W \cdot\overrightarrow{h} a(st−1,hj)=st−1T⋅W⋅h,其中, W W W 是可训练矩阵;
- 多层感知机 a ( s t − 1 , h j ) = V ⋅ tanh ( W 1 ⋅ s t − 1 → + W 2 ⋅ h → ) a(s_{t-1},h_j)=V \cdot \tanh(W_1 \cdot \overrightarrow{s_{t-1}} + W_2 \cdot \overrightarrow{h}) a(st−1,hj)=V⋅tanh(W1⋅st−1+W2⋅h),其中, V V V、 W 1 W_1 W1、 W 2 W_2 W2 都是可训练矩阵;
- 缩放的点乘 a ( s t − 1 , h j ) = s t − 1 → T ⋅ h → ∣ h → ∣ a(s_{t-1},h_j)=\frac{\overrightarrow{s_{t-1}}^T\cdot\overrightarrow{h}}{\sqrt{|\overrightarrow{h}|}} a(st−1,hj)=∣h∣st−1T⋅h。Softmax 函数对非常大的输入很敏感。这会使得梯度的传播出现问题,并且会导致学习的速度下降,甚至会导致学习的停止。那如果我们使用 ∣ h → ∣ \sqrt{|\overrightarrow{h}|} ∣h∣ 来对输入的向量做缩放,就能够防止进入到 softmax 函数的饱和区,使梯度过小。
- 对 e t → \overrightarrow{e_t} et 进行 softmax 操作后得到 Attention 分布 α t → = s o f t m a x ( e t → ) \overrightarrow{\alpha_t} = softmax(\overrightarrow{e_t}) αt=softmax(et),其中, α t j = exp ( e t j ) ∑ i = 1 n exp ( e t i ) \alpha_{tj} = \frac{\exp(e_{tj})}{\sum_{i=1}^n \exp(e_{ti})} αtj=∑i=1nexp(eti)exp(etj)。
- 计算得到上下文表示 c t → = ∑ j = 1 n α t j ⋅ h j \overrightarrow{c_t}=\sum_{j=1}^n \alpha_{tj}\cdot h_j ct=∑j=1nαtj⋅hj。
- 我们可以将该上下文表示利用到下一个时间步的状态生成 s t = f ( s t − 1 , y t − 1 , c t ) s_t = f(s_{t-1}, y_{t-1}, c_t) st=f(st−1,yt−1,ct)。
-
Q(Query), K(Key), V(Value)
在 Attention 中,Q(Query) 指的是被查询的向量,即根据什么来关注其他的单词;K(Key) 指的是查询的向量,即被关注的向量的关键词;V(Value) 则是的被关注的信息本身。
使用 Q 和 K 计算了相似度之后得到相似度评分,之后有了相似度评分,就可以把内容 V 加权回去了。
-
Transformer
既然我们知道 Attention 机制本身就可以获取上下文信息,那么我们可不可以将原本的 RNN 结构完全去掉,仅仅依赖于 Attention 模型呢?这样我们可以使得训练并行化,并且可以拥有全局的信息。根据这个思想,产生了 Transformer 模型。其模型结构如下:
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-2jhqiRcJ-1631367951045)(imgs/transformer.jpg)]
-
Self-Attention 机制
Seq2Seq 中的 Attention 机制是在 decode 过程中,逐步计算对应的上下文表示,仿照这个思想,Self-Attention 就是在 encode 阶段,便考虑到每个输入单词与其他单词的关联性,从而得到具有上下文信息的 input embedding 信息。因此,对于 Self-Attention,其 Q, K, V 都来自于同一个输入矩阵,即 Q=K=V。
Self-Attention 的计算过程如下:
- 输入序列 x → \overrightarrow{x} x;
- 将 x → \overrightarrow{x} x 分别与对应 Q, K, V 的三个可训练矩阵 W q W_q Wq, W k W_k Wk, W v W_v Wv 点乘,得到 Q = x → ⋅ W q Q=\overrightarrow{x}\cdot W_q Q=x⋅Wq, K = x → ⋅ W k K=\overrightarrow{x}\cdot W_k K=x⋅Wk, V = x → ⋅ W v V=\overrightarrow{x}\cdot W_v V=x⋅Wv;
- 计算 A t t e n t i o n ( Q , K , V ) = s o f t m a x ( Q ⋅ K T d K ) ⋅ V Attention(Q,K,V)=softmax(\frac{Q\cdot K^T}{\sqrt{d_K}})\cdot V Attention(Q,K,V)=softmax(dKQ⋅KT)⋅V,其中, d K d_K dK 为 K K K 的维度。
-
Multi-Head Attention
为了使模型能够从不同角度获取输入序列的上下文信息,同时引入多组 ( W q i W_{qi} Wqi, W k i W_{ki} Wki, W v i W_{vi} Wvi) 矩阵,分别得到多个 ( Q i Q_i Qi, K i K_i Ki, V i V_i Vi),再将它们按列拼接,之后经过一个联合矩阵 W o W_o Wo,得到最终的 Attention 表示。过程如图所示:
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-C14bChnH-1631367951046)(imgs/multi-head.jpg)]
注意,在 Transformer 的模型中,有多个 Multi-Head Attention 步骤。其中,encoder 中的 Attention 和 decoder 中的第一步 Attention 的步骤都仅以前一级的输出作为输入,而在 decoder 中的第二步 Attention 则不仅接受来自前一级的输出,还要接收 encoder 的输出。
即,在第一种 Multi-Head Attention 中,有 Q = K = V Q = K = V Q=K=V,在第二种 Multi-Head Attention 中,则 Q ≠ K = V Q \neq K = V Q=K=V。
-
Positional Encoding
由于 Transformer 模型没有循环结构或卷积结构,为了使模型能够学习到输入序列的顺序,我们需要插入一些关于 tokens 位置的信息。因此提出了 Positional Encoding 的概念,其与 input embedding 具有相同的维度,便于相加。
但是,如果直接使用计数的方式来进行 encoding,即 p o s = 1 , 2 , . . . , n − 1 pos = 1, 2, ..., n - 1 pos=1,2,...,n−1,那么最后一个 token 的encoding 将会比第一个 token 大很多,与原 embedding 相加后会造成数据不平衡的现象。原论文作者们的方法是使用了不同频率的正弦和余弦函数来作为位置编码:
P E ( p o s , 2 i ) = s i n ( p o s / 1000 0 2 i / d m o d e l ) P E ( p o s , 2 i + 1 ) = c o s ( p o s / 1000 0 2 i / d m o d e l ) \begin{aligned} PE_{(pos,2i)} & = sin(pos/10000^{2i/d_{model}}) \\ PE_{(pos,2i+1)} & = cos(pos/10000^{2i/d_{model}}) \\ \end{aligned} PE(pos,2i)PE(pos,2i+1)=sin(pos/100002i/dmodel)=cos(pos/100002i/dmodel) -
Add & Norm 层
- Add 指的是 Residual Connection,与 ResNet 的原理相似,是将上一层的信息直接传到下一层,可以帮助解决多层神经网络训练困难的问题。
- Norm 指的是 Layer Normalization,在层与层之间对每一行数据进行缩放。这样可以缓解梯度消失的状况,同时使模型更快收敛。
Batch Normalization 和 Layer Normalization 的区别?
在 BN 中,我们将每一个 batch 中的数据按列进行缩放。而在 NLP 任务中,由于输入序列的长度是不确定的,且不同行同一位置的单词直接并没有直接联系,直接做缩放可能会影响原语义表达。因此,在 NLP 等序列型任务中,我们一般采用 Layer Normalization,即对每一行数据进行缩放。
-
-
BERT: Bi-directional Encoder Representation from Transformers
-
双向表示
区别于 Bi-LSTM 的双向表示,分别正序和反序得到表示再进行拼接,BERT 中的双向指的是根据前文和后文来预测被 masked 的单词。
-
Embedding
BERT 中的 embedding 由三个部分组成:Token Embedding,Segment Embedding,Position Embedding。
- Token Embedding 是词向量,其中,第一个词为 [CLS] 标记,可以用于之后的下游任务。
- Segment Embedding 用于区分 BERT 输入的两个句子,之后的 pre-training 将会用到。
- Position Embedding 由学习得到,而不是普通 Transformer 中的三角函数。
-
Pre-training Tasks
-
Masked LM
在训练过程中,将 15% 的单词用 [mask] 代替,让模型去预测被遮挡的单词,最终的损失函数只计算被遮盖的 token。
但是如果一直用 [mask] 表示(实际预测时并不会遇到 [mask] 标记)会影响模型,因此作者设置了一下规则:80% 的时间用 [mask] 来代替被遮盖的单词,10% 的时间随机用另一个单词代替,剩下 10% 的时间保留原单词。
值得注意的是,模型并不知道哪些单词被遮盖了,这使得模型能够关注到每一个单词,依赖上下文信息预测单词,赋予了模型一定的纠错能力。
-
Next Sentence Prediction
对于输入的两个句子 A 和 B,让模型预测 B 是否应该是 A 的后一句。该任务的目的是让模型理解两个句子直接的关系。
-
-
为什么BERT在第一句前会加一个 [CLS] 标志?
为了获得整个句子的语义表示,用于其他任务。一个没有明显语义的 [CLS] 标记会更加公平地融合句子中每个单词的语义,从而获得更加完整的句子表示。
-
BERT 的优缺点?
优点是建立在 Transformer 上,相对rnn更加高效,具有强大的信息提取能力,能捕捉更长距离的依赖。且双向模型比单向的 Transformer 效果更好;
缺点则是该模型几乎无法修改,只能拿来直接用。由于只能预测 15% 的词,模型收敛较慢,需要强大算力支撑。
-
使用BERT预训练模型为什么最多只能输入 512 个词,最多只能两个句子合成一句?
这是由于在预训练的时候,在参数中设置了 position embedding 的大小和 segment embedding 的大小,分别为 512 和 2。在这之外的单词和句子会没有与之对应的 embedding。
-
BERT 的输入和输出分别是什么?
输入是 token embedding,segment embedding 和 position embedding,输出是文本中各个字 / 词融合了全文语义信息后的向量表示。
-
计算 BERT 模型的参数数量?
- 词向量参数:vocab_size=30522, hidden_size=768, max_position_embedding=512, token_type_embedding=2,因此参数量为 (30522 + 512 + 2) * 768。
- Multi-head Attention:len = hidden_size = 768, d k d_k dk = d q d_q dq = d v d_v dv = d m o d e l / n h e a d d_{model}/n_{head} dmodel/nhead = 768 / 12 = 64,将12个头进行拼接后还要进行线性变换,因此参数量为 768 * 64 * 12 * 3 + 768 * 768。
- 前馈网络参数: FFN ( x ) = max ( 0 , x W 1 + b 1 ) W 2 + b 2 \text{FFN}(x)=\max(0, xW_1+b_1)W_2 + b_2 FFN(x)=max(0,xW1+b1)W2+b2,W_1 和 W_2 的参数量均为 768 * (768 * 4),总参数量为 768 * 768 * 4 * 2。
总参数量 = 词向量参数 + 12 (层数) * (Multi-head + 前馈网络) = 110M
-
-
ALBERT
-
Factorized Embedding Parameterization
在 BERT 中,模型直接将词表对应到 word embedding 中,embedding 的维度大小和隐藏层 H 的维度大小相等。这是没有必要的,因为当维度大小 d H d_H dH 增加时,word embedding 维度的增加没有意义。因此引入多一层转换矩阵 E,让词表 V 先通过转换矩阵,再转换为隐藏层的维度大小。这样可以明显减小参数量,由之前的 ( d V ∗ d H ) (d_V * d_H) (dV∗dH) 减少为 ( d V ∗ d E + d E ∗ d H ) (d_V * d_E + d_E * d_H) (dV∗dE+dE∗dH)。
-
Cross-Layer Parameter Sharing
BERT 框架中的参数主要包括 Attention 层的参数和 Feed Forward 网络的参数,ALBERT 将这些参数都共享,大大减小了参数量,为了弥补性能的损失,ALBERT 加大了隐藏层的维度大小,由“窄而深”变成“宽而浅”。
-
Sentence Order Prediction
针对 BERT 的第二个训练任务,即判断 A 是否是 B 的下一句话,过于简单的问题,ALBERT 增加了预训练的难度,即将负样本换成了两个句子的逆序排列。
[NSP 任务] 正样本:同一个文档的两个连续句子;负样本:两个连续句子交换顺序
[SOP 任务] 正样本:同一个文档的两个连续句子;负样本:不同文档的句子
-
-
XLNet
由于 BERT 在预训练过程中需要加入 [mask],而在下游任务及预测过程中都没有这样的标记,因此会造成性能损失。XLNet 则通过自回归语言模型的思想来解决,即从左到右依次生成。为了保持模型仍然是双向的,能够同时从前文和后文获取信息,XLNet 引入了 Attention Mask 机制。
假设模型在预训练过程中需要预测第 k k k 个词,那么首先先将序列随机打乱,再取前 k − 1 k-1 k−1 个词进行预测,这样既可以读到前后文的信息,又可以省去 [mask] 标记。
这样的预训练模式天然符合下游任务序列生成的任务,因此可以预计 XLNet 在文本摘要,机器翻译,信息检索等领域具有优势。
-
TinyBERT
由于 BERT 模型过于庞大,很难实际应用落地。因此提出了一种蒸馏 BERT 的方法 TinyBERT,它的大小不到 BERT 的 1/7,但速度提高了 9 倍。
知识蒸馏的基本思想是使用一个大的训练好的模型来知道小模型更好的训练。TinyBERT 的基本思想是减少 Transformer 的层数以及降低 hidden_size 的大小。模型结构如下:
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-F1CSYTkj-1631367951047)(imgs/TinyBERT.jpg)]
TinyBERT 的 loss 分为三部分:
-
Embedding Layer Distillation
TinyBERT 的 embedding 大小比教师模型更小,因此需要通过一个维度变换矩阵来把学生模型的 embedding 映射到教师模型所在空间,再通过 MSE 来计算 loss:
L e m b d = MSE ( E S W e , E T ) \mathcal{L}_{embd}= \text{MSE}(E^SW_e, E^T) Lembd=MSE(ESWe,ET) -
Transformer Layer Distillation
TinyBERT 的知识蒸馏采取每隔 k 层蒸馏的方式。设 Teacher BERT 有 12 层,TinyBERT 有 4 层,则学生模型每隔 3 层就与教师模型计算一次 loss,其中,loss 又分为 Attention Loss 和 Hidden Loss:
L a t t n = 1 h ∑ i = 1 h MSE ( A i S , A i T ) \mathcal{L}_{attn} = \frac{1}{h}\sum_{i=1}^h \text{MSE}(A_i^S, A_i^T) Lattn=h1i=1∑hMSE(AiS,AiT)
其中,h 为 Attention 头数, A i ∈ { A q , A k , A v } A_i\in \{A_q,A_k,A_v\} Ai∈{Aq,Ak,Av}。L h i d n = MSE ( H S W h , H T ) \mathcal{L}_{hidn} = \text{MSE}(H^SW_h, H^T) Lhidn=MSE(HSWh,HT)
-
Prediction Layer Distillation
在预测层的 loss 计算取决于不同的具体任务,但都需要结合教师模型和学生模型的 loss。
-
-
RoBERTa
- 去除 NSP 任务
- 动态掩码。RoBERTa的做法是将训练数据复制多份,每份采用不同的随机挑选 token 进行掩码。这样有助于模型适应不同的掩码策略,学习不同的语言表征。
- 文本编码。使用了更大的词汇表来训练。