文章目录
1.过拟合(over-fitting)问题
到现在为止,我们已经学习了几种不同的学习算法,包括线性回归和逻辑回归,它们能够有效地解决许多问题,但是当将它们应用到某些特定的机器学习应用时,会遇到过度拟合(over-fitting)的问题,可能会导致它们效果很差。
下面通过两个例子来解释什么过度拟合问题:
例1:这是一个回归问题的例子。
- 第一个模型是一个线性模型,欠拟合(under-fitting)(或者说“高偏差”(“High bias”)),不能很好地适应我们的训练集;
- 第三个模型是一个四次方的模型,过于强调拟合原始数据,而丢失了算法的本质:预测新数据。我们可以看出,若给出一个新的值使之预测,它将表现的很差,是过拟合(或者说“高方差”(“High variance”)),虽然能非常好地适应我们的训练集,但在新输入变量进行预测时可能会效果不好;
- 中间的模型似乎最合适;
例2:这是一个分类问题的例子。
就以多项式理解,x 的次数越高,拟合的越好,但相应的预测的能力就可能变差。
我们来总结一下什么是过拟合问题:如果我们有非常多的特征(特征少的时候可以通过现有特征构造新特征),我们通过学习得到的假设可能能够非常好地适应训练集(代价函数可能几乎为 0( J ( θ ) ≈ 0 J(θ) ≈0 J(θ)≈0)),但是可能无法泛化到新样本。(泛化:指一个假设模型应用到新样本的能力。)
问题是,如果我们发现了过拟合问题,应该如何处理?
解决过拟合问题方法:
1.丢弃一些不能帮助我们正确预测的特征:
- 可以手工选择保留哪些特征
- 用模型选择算法(Model selection algorithm)帮忙(如PCA)
2.正则化(Regularization):
保留所有的特征,但减少梯度(magnitude)或模型参数 θ j θ_j θj的大小
2.正则化线性回归(Regularized linear regression)
2.1 正则化线性回归的代价函数
上面的回归问题中如果我们的模型是:
h
θ
(
x
)
=
θ
0
+
θ
1
x
1
+
θ
2
x
2
2
+
θ
3
x
3
3
+
θ
4
x
4
4
h_θ(x) = θ_0+θ_1x_1+θ_2x_2^2+θ_3x_3^3+θ_4x_4^4
hθ(x)=θ0+θ1x1+θ2x22+θ3x33+θ4x44我们可以从之前的事例中看出,正是那些高次项导致了过拟合的产生,所以如果我们能让这些高次项的系数接近于 0 的话,我们就能很好的拟合了。
所以我们要做的就是在一定程度上减小这些参数 θ 的值,这就是正则化的基本方法。
假如我们有非常多的特征,我们并不知道其中哪些特征我们要惩罚,我们将对所有的特征进行惩罚,并且让代价函数最优化的软件来选择这些惩罚的程度。这样的结果是最终得到了一个较为简单的能防止过拟合问题的假设函数(这里是正则化线性回归的代价函数):
J
(
θ
)
=
1
2
m
[
∑
i
=
1
m
(
h
θ
(
x
(
i
)
)
−
y
(
i
)
)
2
+
λ
∑
∗
∗
j
=
1
∗
∗
n
θ
j
2
]
J(θ)=\frac{1}{2m}[\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})^2+\lambda\sum_{**j=1**}^{n}θ_j^2]
J(θ)=2m1[i=1∑m(hθ(x(i))−y(i))2+λ∗∗j=1∗∗∑nθj2]
其中,
λ
∑
j
=
1
n
θ
j
2
\lambda\sum_{j=1}^{n}θ_j^2
λ∑j=1nθj2称为正则化项,
λ
\lambda
λ称为正则化参数(regularization parameter)
总的来说,加入正则化项,可以使我们最终得到一个相对简单的假设模型(“Simpler” hypothesis)。
还是以假设函数 h θ ( x ) = θ 0 + θ 1 x 1 + θ 2 x 2 2 + θ 3 x 3 3 + θ 4 x 4 4 h_θ(x) = θ_0+θ_1x_1+θ_2x_2^2+θ_3x_3^3+θ_4x_4^4 hθ(x)=θ0+θ1x1+θ2x22+θ3x33+θ4x44为例,若 λ \lambda λ过大,则 θ = [ θ 1 θ 2 ⋮ θ n ] → 0 θ=\begin{bmatrix} θ_1 \\ θ_2 \\ \vdots \\ θ_n \end{bmatrix} \to 0 θ=⎣⎢⎢⎢⎡θ1θ2⋮θn⎦⎥⎥⎥⎤→0,进而 h θ ( x ) → θ 0 h_θ(x)\toθ_0 hθ(x)→θ0,模型欠拟合。
故为了让正则化起到应有的效果,要选择一个合适的 λ \lambda λ。
2.2 正则化线性回归的梯度下降法
梯度下降法的更新公式:
Repeat until convergence{
θ 0 = θ 0 − α 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x 0 ( i ) \qquad\qquad\qquad\qquad\quadθ_0 = θ_0-α\frac{1}{m}\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})x_0^{(i)} θ0=θ0−αm1∑i=1m(hθ(x(i))−y(i))x0(i)
θ
j
=
θ
j
−
α
[
1
m
∑
i
=
1
m
(
h
θ
(
x
(
i
)
)
−
y
(
i
)
)
x
j
(
i
)
+
λ
m
θ
j
]
\qquad\qquad\qquad\qquad\quadθ_j = θ_j-α[\frac{1}{m}\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})x_j^{(i)}+\frac{\lambda}{m}θ_j]
θj=θj−α[m1∑i=1m(hθ(x(i))−y(i))xj(i)+mλθj],(j=1,2,…n)
}
对更新公式进行合并整理之后,有:
Repeat until convergence{
θ 0 = θ 0 − α 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x 0 ( i ) \qquad\qquad\qquad\qquad\quadθ_0 = θ_0-α\frac{1}{m}\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})x_0^{(i)} θ0=θ0−αm1∑i=1m(hθ(x(i))−y(i))x0(i)
θ
j
=
θ
j
(
1
−
α
λ
m
)
−
α
1
m
∑
i
=
1
m
(
h
θ
(
x
(
i
)
)
−
y
(
i
)
)
x
j
(
i
)
\qquad\qquad\qquad\qquad\quadθ_j = θ_j(1-α\frac{\lambda}{m})-α\frac{1}{m}\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})x_j^{(i)}
θj=θj(1−αmλ)−αm1∑i=1m(hθ(x(i))−y(i))xj(i),(j=1,2,…n)
}
其中, ( 1 − α λ m ) (1-α\frac{\lambda}{m}) (1−αmλ)是一个比1略小一点点的数,比如0.99。
2.3 正则化线性回归的正规方程
2.3.1 正则化线性回归的正规方程表示
优化目标:
m
i
n
θ
J
(
θ
)
min_θJ(θ)
minθJ(θ)
其中的矩阵是一个(n+1)*(n+1)的矩阵。
2.3.2 正规方程中的不可逆(Non-invertibility)问题(选学)
在多变量线性回归中讲过,假如 m ≤ n m\leq n m≤n, θ = ( X T X ) − 1 X T y θ=(X^TX)^{-1}X^Ty θ=(XTX)−1XTy 中 X T X X^TX XTX不可逆。
现在在J(θ)中加入正则化项之后,如果
λ
>
0
\lambda >0
λ>0,中要求逆的矩阵一定是可逆的。
3.正则化逻辑回归(Regularized logistic regression)
3.1 正则化逻辑回归的代价函数
J ( θ ) = − 1 m ∑ i = 1 m [ y ( i ) l o g ( h θ ( x ( i ) ) ) + ( 1 − y ( i ) ) l o g ( 1 − h θ ( x ( i ) ) ) ] + λ 2 m ∑ ∗ ∗ j = 1 ∗ ∗ n θ j 2 J(θ)=-\frac{1}{m}\sum_{i=1}^{m}[y^{(i)}log(h_θ(x^{(i)}))+(1-y^{(i)})log(1-h_θ(x^{(i)}))]+\frac{\lambda}{2m}\sum_{**j=1**}^{n}θ_j^2 J(θ)=−m1i=1∑m[y(i)log(hθ(x(i)))+(1−y(i))log(1−hθ(x(i)))]+2mλ∗∗j=1∗∗∑nθj2
3.2 正则化逻辑回归的梯度下降法
梯度下降法的更新规则:
Repeat until convergence{
θ 0 = θ 0 − α 1 m ∑ i = 1 m ( h θ ( x ( i ) ) − y ( i ) ) x 0 ( i ) \qquad\qquad\qquad\qquad\quadθ_0 = θ_0-α\frac{1}{m}\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})x_0^{(i)} θ0=θ0−αm1∑i=1m(hθ(x(i))−y(i))x0(i)
θ
j
=
θ
j
−
α
[
(
1
m
∑
i
=
1
m
(
h
θ
(
x
(
i
)
)
−
y
(
i
)
)
x
j
(
i
)
)
+
λ
m
θ
j
]
\qquad\qquad\qquad\qquad\quadθ_j = θ_j-α[(\frac{1}{m}\sum_{i=1}^{m}(h_θ(x^{(i)})-y^{(i)})x_j^{(i)})+\frac{\lambda}{m}θ_j]
θj=θj−α[(m1∑i=1m(hθ(x(i))−y(i))xj(i))+mλθj],(j=1,2,…n)
}
需要注意的三点:
- 虽然这里的的梯度下降和正则化的线性回归中的表达式表面上“看起来”一样,但由于两者的假设函数 h θ ( x ) h_θ(x) hθ(x)不同,所以还是有很大差别。 这里的假设函数 h θ ( x ) = 1 1 + e − θ T x h_θ(x)=\frac{1}{1+e^{-θ^Tx}} hθ(x)=1+e−θTx1
- 上述更新规则中,方括号里面的项为 ∂ ∂ θ j J ( θ ) \frac{\partial}{\partialθ_j}J(θ) ∂θj∂J(θ),J(θ)为正则化代价函数
- θ 0 θ_0 θ0不参与其中的任何一个正则化
3.3 如何在高级优化算法中使用正则化
在逻辑回归模型的文章中我们提到了,使用优化算法时,那么我们需要做的是编写代码,当输入参数 θ 时,它们会计算出两样东西:
- J(θ)
- ∂ ∂ θ j J ( θ ) ( f o r j = 0 , 1 , 2 , . . . , n + 1 ) \frac{\partial}{\partial θ_j}J(θ)(for \ j = 0,1,2,...,n+1) ∂θj∂J(θ)(for j=0,1,2,...,n+1)
完成上述编码之后,就可以用梯度下降法或者共轭梯度法、BFGS、L-BFGS等高级优化算法来拟合模型参数θ了。
下图为在Octave中使用fminunc(这一算法的作用是求函数在无约束条件下的最小值)这一高级优化算法来拟合模型参数θ的伪代码: