统计学习是关于计算机 基于数据构建概率统计模型并运用模型对数据进行预测与分析的一门学科。总的来说,统计学习也称为统计机器学习。
以计算机以及网络为平台;
以数据为对象,是数据驱动的学科;
目的是对数据进行预测和分析;
以方法为中心,构建模型并应用模型进行预测和分析;
是概率论、统计学、信息论、计算理论、最优化理论以及计算机科学等多个领域的交叉学科,具有独立的理论体系和方法论。
统计学习的对象是数据(Data).从数据出发,提取数据特征,抽象出数学模型,发现数学中的知识,又回到对数据分析和预测中去。数据是多样的,包括存在计算机上的一些数字、图片、音频等。
监督学习:监督学习的任务就是学习一个模型,使得模型能够对任意给定的输入,对其相应的输出做出一个好的预测。
非监督学习
半监督学习
强化学习
模型——Model
策略——Strategy
算法——Algorithm
统计学习首要考虑的问题就是学习什么样的模型。在监督学习过程中,模型就是所要学习的条件概率分布或者决策函数。模型的假设空间包含所有可能条件概率分布或者决策函数。
1.损失函数和风险函数
监督学习问题实在假设空间中选取模型作为决策函数,对于给定的输入X,由函数给出相应的输出Y,这个输出的值可能和预测值不一样。用一个损失函数或者代价函数来度量预测的错误的程度。统计学习常用的损失函数有如下几种:
0-1损失函数
L ( Y , f ( X ) ) = { 1 , Y ≠ f ( X ) 0 , Y = f ( X ) L(Y,f(X)) = \left\{ \begin{aligned} 1 & , & Y \neq f(X) \\ 0 & , & Y = f(X) \end{aligned} \right. L(Y,f(X))={10,,Y=f(X)Y=f(X)
平方损失函数
L ( Y , f ( X ) ) = ( Y − f ( X ) ) 2 L(Y,f(X)) =(Y-f(X))^2 L(Y,f(X))=(Y−f(X))2
绝对损失函数
L ( Y , f ( X ) ) = ∣ Y − f ( X ) ∣ L(Y,f(X)) = | Y - f(X)| L(Y,f(X))=∣Y−f(X)∣
对数损失函数或者对数似然损失函数
L ( Y , f ( X ) ) = − lg P ( Y ∣ X ) L(Y,f(X)) = -\lg P(Y|X) L(Y,f(X))=−lgP(Y∣X)
损失函数值越小,说明模型越好。由于模型的输入、输出(X,Y)是随机变量,遵循联合分布P(X,Y),所以损失函数的期望是:
R e x p ( f ) = E p [ L ( Y , f ( X ) ) ] = ∫ x × y L ( y , f ( x ) ) P ( x , y ) d x d y R_{exp}(f) = E_p[L(Y,f(X))]=\int_{x\times y} L(y,f(x))P(x,y)dxdy Rexp(f)=Ep[L(Y,f(X))]=∫x×yL(y,f(x))P(x,y)dxdy
这称为风险函数或者期望损失(Expected Loss)。但是在实际情况下,由于联合分布P(X,Y)是未知的,期望损失还是难以计算的。此时需要计算经验风险,也即是模型关于选练级的平均损失,当样本容量趋于无穷大时,经验风险时趋近于期望风险。实际上,这也是统计学习方法中的两种策略:经验风险最小化和结构风险最小化。
经验风险最小化(empirical risk minimization,ERM)策略表示经验风险最小的模型是最优的模型。即:
min f ϵ F 1 N ∑ i = 1 N L ( y i , f ( x i ) ) \min \limits_{f \epsilon F} \frac{1}{N} \sum^N_{i=1} L(y_i,f(x_i)) fϵFminN1i=1∑NL(yi,f(xi))
当样本容量较大时,经验风险最小化能提供一个很好的效果,但是当样本容量过小时,模型容易出现“过拟合”的现象。此时需要考虑结构风险最小化(structural risk minimization,SRM),它是为了防止过拟合而提出来的策略。结构风险最小化等于正则化。结构风险需要在经验风险上再加上正则化项或者惩罚项,可定义如下:
经 验 风 险 : R e m p ( f ) = 1 N ∑ i = 1 N L ( y i , f ( x i ) ) 结 构 风 险 : R s r m ( f ) = 1 N ∑ i = 1 N L ( y i , f ( x i ) ) + λ J ( f ) 经验风险: R_{emp}(f)= \frac{1}{N} \sum^N_{i=1} L(y_i,f(x_i)) \\ 结构风险: R_{srm}(f) = \frac{1}{N} \sum^N_{i=1} L(y_i,f(x_i)) + \lambda J(f) 经验风险:Remp(f)=N1i=1∑NL(yi,f(xi))结构风险:Rsrm(f)=N1i=1∑NL(yi,f(xi))+λJ(f)
其中J(f)表示模型的复杂度。复杂度表示了对复杂模型的惩罚。如贝叶斯估计就是最好的结构最小化策略的例子。
最终的问题编程求解最优解问题:
min f ϵ F 1 N ∑ i = 1 N L ( y i , f ( x i ) ) + λ J ( f ) \min \limits_{f \epsilon F} \frac{1}{N} \sum^N_{i=1} L(y_i,f(x_i)) + \lambda J(f) fϵFminN1i=1∑NL(yi,f(xi))+λJ(f)
算法是指学习模型的具体计算的方法。
训练误差的大小反应的是判断给定的问题是不是一个容易学习的问题是有意义的,测试误差反映了学习方法对未知的测试数据集的预测能力。测试误差小的方法具有更好的预测的能力,通常将学习方法对未知数据的预测能力称为泛化能力。
当所选的模型的复杂度比真模型更高,这种现象称为过拟合。过拟合表示学习时选择的模型包含的参数过多,以至于出现这一模型对已知数据预测得更好,但对未知测试数据得预测能力差。此时采取的策略可以是使用正则化等。
模型选择得典型方法是正则化。正则化是结构风险最小化策略得实现。
交叉验证
简单交叉验证
S折交叉验证
留一交叉验证:S折交叉验证的特殊情形。
如果学到的模型是f,那么这个模型对未知数据预测的误差是:
R e x p ( f ) = E p [ L ( Y , f ( X ) ) ] = ∫ x × y L ( y , f ( x ) ) P ( x , y ) d x d y R_{exp}(f) = E_p[L(Y,f(X))]=\int_{x\times y} L(y,f(x))P(x,y)dxdy Rexp(f)=Ep[L(Y,f(X))]=∫x×yL(y,f(x))P(x,y)dxdy
泛化误差反映了学习方法的泛化能力。实际上,泛化误差也是期望误差。
学习方法的泛化能力往往是通过分析研究泛化误差的概率上界进行的,简称为泛化误差上界。
定理:对于而分类问题,当假设空间是有限个函数的集合F={f1,f2,…fd}时,对任意一个函数f属于F,至少满足以下不等式:
R ( f ) ≤ R ‾ ( f ) + ε ( d , N , δ ) 其 中 , ε ( d , N , δ ) = 1 2 N ( lg d + lg 1 δ ) R(f) \leq \overline {R}(f) + \varepsilon(d,N,\delta) \\ 其中,\varepsilon(d,N,\delta)= \sqrt{\frac{1}{2N}(\lg d+ \lg \frac{1}{\delta})} R(f)≤R(f)+ε(d,N,δ)其中,ε(d,N,δ)=2N1(lgd+lgδ1)
右边式子中第一项是训练误差,训练误差越小,泛化误差也就越小。第2项是单调递减函数,当N趋于无穷时趋于0.
证明过程如下:
提 示 : 在 证 明 过 程 中 需 要 用 到 H o e f f d i n g 不 等 式 。 H o e f f d i n g 不 等 式 描 述 如 下 : 设 S n = ∑ i = 1 X i 是 随 机 变 量 X 1 、 X 2 、 . . . X n 之 和 , X i ϵ [ a i , b i ] , 则 对 任 意 t > 0 , 以 下 不 等 式 成 立 : P ( S n − E S n ≥ t ) ≤ e x p ( − 2 t 2 ∑ i = 1 ( b i − a i ) ) P ( E S n − S n ≥ t ) ≤ e x p ( − 2 t 2 ∑ i = 1 ( b i − a i ) ) 对 于 任 意 函 数 f ϵ F , R ‾ ( f ) 是 N 个 独 立 随 机 变 量 L ( Y , f ( X ) ) 的 样 本 均 值 , R ( f ) 是 随 机 变 量 值 的 期 望 值 。 对 于 样 本 空 间 中 的 所 有 随 机 变 量 , ∵ S n = X 1 + X 2 + . . . X n , X ‾ = 1 N S n , E S n = E ( ∑ i = 1 X i ) ∴ E ( X ‾ ) = 1 N E S n ∴ P ( X n ‾ − E ( X n ‾ ) ≥ t ) = P ( S n − E S n ≥ n t ) ≤ e x p ( − 2 t 2 ∑ i = 1 ( b i − a i ) ) 即 : P ( X n ‾ − E ( X n ‾ ) ≥ t ) = P ( S n − E S n ≥ n t ) ≤ e x p ( − 2 n 2 t 2 ∑ i = 1 ( b i − a i ) ) 对 于 二 分 类 任 务 , 在 样 本 空 间 中 : R ( f ) = E [ L ( Y , f ( X ) ) ] 代 表 着 期 望 风 险 , 表 示 备 选 模 型 在 测 试 集 上 的 表 现 R ‾ ( f ) = 1 N ∑ i = 1 n L ( y i , f ( x i ) ) , 代 表 经 验 风 险 , 表 示 着 备 选 模 型 在 训 练 集 上 的 表 现 ∵ 损 失 函 数 取 值 于 [ 0 , 1 ] , 即 对 所 有 i , [ a i , b i ] = [ 0 , 1 ] 。 ∴ P ( X n ‾ − E ( X n ‾ ) ≥ t ) = P ( S n − E S n ≥ n t ) ≤ e x p ( − 2 n 2 t 2 ∑ i = 1 ( b i − a i ) ) = e x p ( − 2 n 2 t 2 n ) = e x p ( − 2 n t 2 ) ∴ 对 于 ∀ ε > 0 , 以 下 满 足 : P ( R ( f ) − R ‾ ( f ) ≥ ε ) ≤ e x p ( − 2 n ε 2 ) ∵ 命 题 是 对 任 意 一 个 函 数 f 属 于 F , 至 少 都 有 概 率 1 − δ 满 足 不 等 式 成 立 ∴ 对 于 d 个 备 选 模 型 , 有 P ( ∃ f ∈ F , P ( R ( f ) − R ‾ ( f ) ≥ ε ) = P ( R ( f 1 ) − R ‾ ( f 1 ) ≥ ε ⋃ R ( f 2 ) − R ‾ ( f 2 ) ≥ ε . . . ⋃ R ( f n ) − R ‾ ( f n ) ≥ ε ) ≤ ∑ i = 1 d P ( R ( f i ) − R ‾ ( f i ) ≥ ε ) ≤ d e x p ( − 2 n ε 2 ) ∴ 原 命 题 满 足 : P ( ∀ f ∈ F , P ( R ( f ) − R ‾ ( f ) ≤ ε ) ≥ 1 − d e x p ( − 2 n ε 2 ) 令 1 − d e x p ( − 2 n ε 2 ) = δ , 则 ε = ε ( d , n , δ ) = 1 2 n ( lg d + lg 1 δ ) , 其 中 ε 是 一 个 三 元 组 , d 是 备 选 模 型 个 数 , n 是 样 本 数 即 : 上 述 不 等 式 P ( ∀ f ∈ F , P ( R ( f ) − R ‾ ( f ) ≤ ε ) ≥ 1 − δ 成 立 提示:在证明过程中需要用到Hoeffding不等式。\\ Hoeffding不等式描述如下:\\ 设S_n=\sum_{i=1}Xi 是随机变量X_1、X_2、...X_n之和,X_i \epsilon[a_i,b_i],则对任意t>0,以下不等式成立:\\ P(S_n-ES_n \geq t)\leq exp(\frac{-2t^2}{\sum_{i=1}(b_i - a_i)})\\ P(ES_n-S_n \geq t)\leq exp(\frac{-2t^2}{\sum_{i=1}(b_i - a_i)})\\ 对于任意函数f \epsilon F,\overline {R}(f)是N个独立随机变量L(Y,f(X))的样本均值,R(f)是随机变量值的期望值。\\ 对于样本空间中的所有随机变量,\\ \because S_n = X_1 + X_2+...X_n,\overline{X}=\frac{1}{N}S_n,ES_n=E(\sum_{i=1} Xi) \\ \therefore E(\overline{X})= \frac{1}{N}ES_n\\ \therefore P(\overline{X_n}-E(\overline{X_n})\geq t)=P(S_n-ES_n \geq nt) \leq exp(\frac{-2t^2}{\sum_{i=1}(b_i - a_i)})\\ 即:P(\overline{X_n}-E(\overline{X_n})\geq t)=P(S_n-ES_n \geq nt) \leq exp(\frac{-2n^2t^2}{\sum_{i=1}(b_i - a_i)})\\ 对于二分类任务,在样本空间中:\\ R(f)= E[L(Y,f(X))]代表着期望风险,表示备选模型在测试集上的表现\\ \overline{R}(f)=\frac{1}{N}\sum^n_{i=1}L(y_i,f(x_i)),代表经验风险,表示着备选模型在训练集上的表现\\ \because 损失函数取值于[0,1],即对所有i,[a_i,b_i]=[0,1]。\\ \therefore P(\overline{X_n}-E(\overline{X_n})\geq t)\\ =P(S_n-ES_n \geq nt)\\ \leq exp(\frac{-2n^2t^2}{\sum_{i=1}(b_i - a_i)}) = exp(\frac{-2n^2t^2}{n})=exp(-2nt^2)\\ \therefore 对于 \forall \varepsilon > 0,以下满足:\\ P(R(f)-\overline{R}(f) \geq \varepsilon) \leq exp(-2n\varepsilon ^2)\\ \because 命题是对任意一个函数f属于F,至少都有概率1-\delta满足不等式成立\\ \therefore 对于d个备选模型,有P(\exists f \in F,P(R(f)-\overline{R}(f) \geq \varepsilon)\\ = P(R(f_1)-\overline{R}(f_1) \geq \varepsilon \bigcup R(f_2)-\overline{R}(f_2) \geq \varepsilon ... \bigcup R(f_n)-\overline{R}(f_n) \geq \varepsilon) \\ \leq \sum_{i=1}^d P(R(f_i)-\overline{R}(f_i) \geq \varepsilon)\\ \leq d exp(-2n\varepsilon ^2)\\ \therefore 原命题满足:P(\forall f \in F,P(R(f)-\overline{R}(f) \leq \varepsilon) \geq 1-d exp(-2n\varepsilon ^2)\\ 令1-d exp(-2n\varepsilon ^2)=\delta,则\varepsilon=\varepsilon(d,n,\delta)= \sqrt{\frac{1}{2n}(\lg d+ \lg \frac{1}{\delta})},\\ 其中\varepsilon是一个三元组,d是备选模型个数,n是样本数\\ 即:上述不等式P(\forall f \in F,P(R(f)-\overline{R}(f) \leq \varepsilon) \geq 1- \delta成立 提示:在证明过程中需要用到Hoeffding不等式。Hoeffding不等式描述如下:设Sn=i=1∑Xi是随机变量X1、X2、...Xn之和,Xiϵ[ai,bi],则对任意t>0,以下不等式成立:P(Sn−ESn≥t)≤exp(∑i=1(bi−ai)−2t2)P(ESn−Sn≥t)≤exp(∑i=1(bi−ai)−2t2)对于任意函数fϵF,R(f)是N个独立随机变量L(Y,f(X))的样本均值,R(f)是随机变量值的期望值。对于样本空间中的所有随机变量,∵Sn=X1+X2+...Xn,X=N1Sn,ESn=E(i=1∑Xi)∴E(X)=N1ESn∴P(Xn−E(Xn)≥t)=P(Sn−ESn≥nt)≤exp(∑i=1(bi−ai)−2t2)即:P(Xn−E(Xn)≥t)=P(Sn−ESn≥nt)≤exp(∑i=1(bi−ai)−2n2t2)对于二分类任务,在样本空间中:R(f)=E[L(Y,f(X))]代表着期望风险,表示备选模型在测试集上的表现R(f)=N1i=1∑nL(yi,f(xi)),代表经验风险,表示着备选模型在训练集上的表现∵损失函数取值于[0,1],即对所有i,[ai,bi]=[0,1]。∴P(Xn−E(Xn)≥t)=P(Sn−ESn≥nt)≤exp(∑i=1(bi−ai)−2n2t2)=exp(n−2n2t2)=exp(−2nt2)∴对于∀ε>0,以下满足:P(R(f)−R(f)≥ε)≤exp(−2nε2)∵命题是对任意一个函数f属于F,至少都有概率1−δ满足不等式成立∴对于d个备选模型,有P(∃f∈F,P(R(f)−R(f)≥ε)=P(R(f1)−R(f1)≥ε⋃R(f2)−R(f2)≥ε...⋃R(fn)−R(fn)≥ε)≤i=1∑dP(R(fi)−R(fi)≥ε)≤dexp(−2nε2)∴原命题满足:P(∀f∈F,P(R(f)−R(f)≤ε)≥1−dexp(−2nε2)令1−dexp(−2nε2)=δ,则ε=ε(d,n,δ)=2n1(lgd+lgδ1),其中ε是一个三元组,d是备选模型个数,n是样本数即:上述不等式P(∀f∈F,P(R(f)−R(f)≤ε)≥1−δ成立
监督学习包括生成模型和判别模型。
生成模型:由数据学习联合概率分布P(X,Y),然后求出条件概率分布P(Y|X)作为预测的模型。由给定输入的X产生输出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)
判别模型:由数据直接学习方法决策函数F(X)或者条件概率分布P(Y|X)作为预测模型。判别方法关心给定的输入X,应该如何预测如何的模型。常见的判别模型:K近邻方法、感知机、决策树、逻辑回归、最大熵模型、支持向量机、条件随机场等。
生成方法:
判别方法: