回归方法

1 逻辑斯蒂回归

因变量是定性变量(假设只取两个值),自变量可以定量、定性(哑变量)。输出结果在0到1之间,用逻辑斯蒂函数

\[ p(X)=\frac{e^{\beta_{0}+\beta_{1} X}}{1+e^{\beta_{0}+\beta_{1} X}} \]

逻辑斯蒂函数可以等价地写成 \[ \frac{p(X)}{1-p(X)}=e^{\beta_{0}+\beta_{1} X} \] 称左边为“发生比", 它的取值范围为 \((0, \infty)\) 取对数,得: \[ \log \left(\frac{p(X)}{1-p(X)}\right)=\beta_{0}+\beta_{1} X \text { . } \] 极大似然估计参数,似然函数: \[ L\left(\beta_{0}, \beta_{1}\right)=\prod_{i: y_{i}=1} p\left(x_{i}\right) \prod_{i^{\prime}: y_{i^{\prime}}=0}\left(1-p\left(x_{i^{\prime}}\right)\right) \] 梯度下降法或牛顿迭代法求系数估计。

2 判别分类

贝叶斯定理: \[ p_{k}(x):=P(Y=k \mid X=x)=\frac{\pi_{k} f_{k}(x)}{\sum_{l=1}^{K} \pi_{l} f_{l}(x)} \] 估计\(\pi_k\):取随机样本,计算属于第\(k\)类的样本占总样本的比例。

将一个待判别的\(x\)分类到使得\(p_k(x)\)达到最大的那个类,这种方法被称为贝叶斯分类器。

2.1 线性判别(LDA)

假设预测变量只有1个。此外,假设 \(f_{k}(x)\) 是正态的: \[ f_{k}(x)=\frac{1}{\sqrt{2 \pi} \sigma_{k}} \exp \left\{-\frac{1}{2 \sigma_{k}^{2}}\left(x-\mu_{k}\right)^{2}\right\}, \quad x \in R \] 再假设 \(\sigma_{1}^{2}=\cdots=\sigma_{K}^{2}\), 简记为 \(\sigma^{2}\), 那么: \[ p_{k}(x)=\frac{\pi_{k} \frac{1}{\sqrt{2 \pi} \sigma} \exp \left\{-\frac{1}{2 \sigma^{2}}\left(x-\mu_{k}\right)^{2}\right\}}{\sum_{l=1}^{K} \pi_{l} \frac{1}{\sqrt{2 \pi} \sigma} \exp \left\{-\frac{1}{2 \sigma^{2}}\left(x-\mu_{l}\right)^{2}\right\}} \]\(p_{k}(x)\) 两边取对数,可知贝叶斯分类器其实是将观测x分到使得 \[ \delta_{k}(x)=x \frac{\mu_{k}}{\sigma^{2}}-\frac{\mu_{k}^{2}}{2 \sigma^{2}}+\log \pi_{k} \] 达到最大的那一类。 假设\(K =2,\ \pi_{1}=\pi_{2}\), 则当 \(2x\left(\mu_{1}-\mu_{2}\right)>\mu_{1}{ }^{2}-\mu_{2}{ }^{2}\) 时,贝叶斯分类器将 观测分入类1, 否则分入类2. 贝叶斯决策边界: \[ x=\frac{\mu_{1}^{2}-\mu_{2}^{2}}{2\left(\mu_{1}-\mu_{2}\right)}=\frac{\mu_{1}+\mu_{2}}{2} \text { . } \]

2.2 二次判别分析(QDA)

假设来自第k类的随机观测服从\(N\left(\mu_{k}, \sigma_{k}^{2}\right)\), 在这种假设下,贝叶斯 分类器把观测分入使得 \[ \begin{aligned} \delta_{k}(x) &=-\frac{1}{2 \sigma_{k}^{2}}\left(x-\mu_{k}\right)^{2}+\log \pi_{k} \\ &=-\frac{1}{2 \sigma_{k}^{2}} x^{2}+\frac{x \mu_{k}}{\sigma_{k}^{2}}-\frac{\mu_{k}^{2}}{2 \sigma_{k}^{2}}+\log \pi_{k} \end{aligned} \] 达到最大的那一类。

当有\(p\)个自变量时, LDA的协方差矩阵有\(p(p+1)/2\)个参数, 而QDA的\(K\)个协方差矩阵有\(Kp(p+1)/2\)个参数。所以LDA没有QDA分类器光滑,LDA拥有更小的方差和更大的预测偏差。样本量小,用LDA;样本量多,用QDA。

3 K最近邻分类

对新的观测,根据其\(k\)个最近邻的训练数据的类别,通过多数表决等方式进行类别预判。因此,\(k\)最近邻方法不具有显式学习过程。

\(k\)最近邻法的三个基本要素:

涵盖最邻近\(k\)个点的 \(x\) 的邻域记作 \(N_{k}(x)\) ,在 \(N_{k}(x)\) 中根据分类决策规则(如多数表决)决定 \(x\) 的类別 \(y\) \[ y=\arg \max _{c_{j}} \sum_{x_{i} \in N_{k}(x)} I\left\{y_{i}=c_{j}\right\}, \quad i=1, \cdots, N ; j=1, \cdots, L ; \]\(k=1\)时,最近邻分类器偏差较小但方差很大,决策边界很不规则。 当\(k\)变大时,方差较低但偏差却增大,将得到一个接近线性的决策边界。 在实际中,可用交叉验证的方法选择\(k\)的大小。

4 决策树

分而治之

因变量数值型 - 回归问题 - 回归树

因变量类别型 - 分类问题 - 分类树

4.1 回归树

出于简化模型和增加模型的可解释性的考虑,通常将自变量空间划为高维矩形,或称为盒子。划分区域的目标是找到使模型的残差平方和RSS最小的矩形区域 $R_{1}, , R_{J} $. RSS的定义是: \[ \sum_{j=1}^{J} \sum_{i \in R_{j}}\left(y_{i}-\hat{y}\left(R_{j}\right)\right)^{2} \] 其中, \(\hat{y}\left(R_{j}\right)\) 是第\(j\)个矩形区域中训练观测的平均响应值。

将自变量空间划分为\(J\)个矩形区域,一般采用自上而下、贪婪的方法:递归二叉分裂。

在执行递归二叉分裂时,先选择自变量\(X_{j}\) 和分割点\(s\),将自变量空间分为两个区域: \[ \begin{array}{l} \left\{X \mid X_{j}<s\right\} \\ \left\{X \mid X_{j} \geq s\right\} \end{array} \] 选择出来的两个区域要能够使RSS尽可能地减少。

重复至某个停止标准,如所有区域包含的观测个数都不大于5.

成本复杂性剪枝

不是考虑每一棵可能的子树,而是考虑以非负调节参数 \(\alpha\) 标记的一列子树。每一个 \(\alpha\) 的取值对应一棵子树 \(T \subset T_{0}\) 。当 \(\alpha\) 值给定时,其对应的子树需使下式最小 \[ \sum_{m=1}^{|T|} \sum_{i: \boldsymbol{x}_{i} \in R_{m}}\left(y_{i}-\hat{y}\left(R_{m}\right)\right)^{2}+\alpha|T| \] 这里的\(|T|\)表示树的叶节点个数, \(R_m\) 是第\(m\)个叶节点对应的盒子, \(\hat{y}\left(R_{m}\right)\) 是相应的响应预测值。

交叉验证挑选\(\alpha\),使均方预测误差达到最小,从而确定最优子树。

4.2 分类树

递归二叉分裂。RSS的替代指标:

\[ E_{m}=1-\max _{k} \hat{p}_{m k} \] 其中, \(\hat{p}_{\mathrm{mk}}\) 表示第\(m\)个区域的训练观测中第\(k\)类所占的比例, \(\max _{k} \hat{p}_{m k}\) 是分类正确率。

但分类错误率在分类树的构建中不够敏感,常用基尼指数和互熵衡量节点纯度。

\[ G_{m}=\sum_{k=1}^{K} \hat{p}_{m k}\left(1-\hat{p}_{m k}\right) . \]

\[ D_{m}=-\sum_{k=1}^{K} \hat{p}_{m k} \log \hat{p}_{m k} \]

5 集成学习

两大类集成树的产生方法:

5.1 Bagging

Bagging主要关注降低预测模型的方差。 给定\(n\)个独立随机变量\(Z_{1}, \cdots, Z_{n}\), 假设它们的方差都为 \(\sigma^{2}\), 那么样本 均值 \(\bar{Z}=\frac{1}{n} \sum_{i=1}^{n} Z_{i}\) 的方差为 \(\sigma^{2} / \mathrm{n}_{\text {。 }}\) 启发:从总体中抽取多个训练集,对每个训练集分别建立预测模型,再对由此得到的全部预测模型求平均,从而降低方差,得到一个集成模型。

即,可以用B个独立的训练集训练出B个模型: \(\hat{f}^{1}(x), \cdots, \hat{f}^{B}(x)\), 然 后求平均,得到一个低方差的模型: \[ \hat{f}_{\text {avg }}(\boldsymbol{x})=\frac{1}{B} \sum_{b=1}^{B} \hat{f}^{b}(\boldsymbol{x}) . \] 在实际中,不容易得到多个训练集。自助抽样法(Bootstrap)可以解决这个问题。

\(B\)足够大,可稳定预测准确率。

可以对某一自变量在一棵个体分类上因分裂导致的基尼指数减少量加总,再在所有\(B\)棵个体分类树上求平均值。平均值越大就说明这个自变量越重要。

5.2 随机森林

以决策树为基础构建Bagging分类树的基础上,进一步在决策树的训练过程中引入了自变量的随机选择,从而达到对树的去相关(decorrelating),实现对Bagging的改进。

在建立这些个体分类树时,每考虑树上的一个分裂点,都要从全部 的p个自变量中选出一个包含\(q\ (1 \leq q \leq p)\)个自变量的随机样本作为候选变量。这个分裂点所用的自变量只能从这\(q\)个变量中选择。在每个分裂点处都重新进行抽样,选出\(q\)个自变量。若\(q=p\),则随机森林就是Bagging. 通常取\(q\)\(p\)的平方根。

若个体分类树有高度相关性,对高度相关的变量求平均无法大幅度减少方差。而随机森林对不相关变量求平均,可大幅度降低方差。

样本多样性来自

5.3 Adaboost

Boosting:可将弱分类器提升为强分类器。根据分类器的表现对训练样本分布进行调整,使先前分类器错分的训练样本在后续得到更多的关注。

Adaboost是Boosting中的一种,算法:

输入:训练集 \(D=\left\{\left(\boldsymbol{x}_{1}, y_{1}\right), \cdots,\left(\boldsymbol{x}_{m}, y_{m}\right)\right\}\); 分类器算法 \(\mathfrak{L} ;\) 训练轮数\(T\);

过程

(a)\(\mathcal{D}_{1}(\boldsymbol{x})=1 / m\). (b) 对 \(t=1, \cdots, T\), 执行: (c) \(h_{t}=\mathfrak{L}\left(D, \mathcal{D}_{t}\right)\) (d) \(\epsilon_{t}=P_{\boldsymbol{x} \sim \mathcal{D}_{t}}\left(h_{t}(\boldsymbol{x}) \neq f(\boldsymbol{x})\right)\); (e) 如果 \(\epsilon_{t}>0.5\), 则停止;否则,继续执行; (f) \(\alpha_{t}=\frac{1}{2} \ln \left(\frac{1-\epsilon_{t}}{\epsilon_{t}}\right)\); (g) 令 \[ \begin{aligned} \mathcal{D}_{t+1}(\boldsymbol{x}) &=\frac{\mathcal{D}_{f}(\boldsymbol{x}) \exp \left(-\alpha_{t} f(\boldsymbol{x}) h_{t}(\boldsymbol{x})\right)}{Z_{t}} \\ &=\frac{\mathcal{D}_{t}(\boldsymbol{x})}{Z_{t}} \times\left\{\begin{array}{ll} \exp \left(-\alpha_{t}\right), & \text { 如果 } h_{t}(\boldsymbol{x})=f(\boldsymbol{x}) \\ \exp \left(\alpha_{t}\right), & \text { 如果 } h_{t}(\boldsymbol{x}) \neq f(\boldsymbol{x}) \end{array}\right. \end{aligned} \] 其中 \(Z_{t}\) 是归一化常数; (h) 循环结束.

输出\(H(\boldsymbol{x})=\operatorname{sign}\left(\sum_{t=1}^{T} \alpha_{t} h_{t}(\boldsymbol{x})\right)\).

从偏差-方差权衡角度,Adaboost更关注降低偏差。