【深度学习实践】02. Softmax 回归

虽然称作 Softmax 回归,但本质上 Softmax 处理的是分类问题。分类,顾名思义,即对样本数据集进行“哪一个”的解答,通常类别的范围在分类前是已知的,例如一封电子邮件,我们需要将其分类为垃圾邮件或者非垃圾邮件。通常,机器学习实践者用分类这个词来描述以下两个有微妙差别的问题:

  1. 只对样本的“硬性”类别感兴趣,即想要知道其属于哪个类别;
  2. “软性”的类别,即需要得到属于每个类别的概率(置信度)。

这两者的界限往往比较模糊。即使我们只关心所属的类别,我们仍然会使用类别概率的模型。以下展示分类问题的简单模型,每个 o 代表每个类别的置信度:


分类模型

通常我们通过 独热编码(one-hot encoding)对标签类别进行编码。独热编码是一个向量,它的分量和类别一样多。类别对应的分量设置为1,其他所有分量设置为0。在我们的例子中,标签y将是一个三维向量,其中(1, 0, 0)对应于“猫”、(0, 1, 0)对应于“鸡”、(0, 0, 1)对应于“狗”:
y \in \{(1, 0, 0), (0, 1, 0), (0, 0, 1)\}.

为什么不适用 1,2,3 代表猫、鸡,狗的类别,转化为回归问题呢?主要是因为三种类别互相之间没有顺序,而1和2的距离比1和3的距离更近,在回归问题中先验不同。但如果我们试图预测\{\text{婴儿}, \text{儿童},\text{青少年}, \text{青年人}, \text{中年人}, \text{老年人}\},那么将这个问题转变为回归问题,并且保留这种格式是有意义的。

softmax

既然我们需要根据输出类别的置信度来判断预测类别,并且通过独热编码来计算损失函数,那么接下来我们需要将输出的各类别的实数域上的值缩放到“概率分布”中,也就是确保各类别的输出和为 1,相对关系保持不变(原先相对较大的值在转化后仍然较大,负数要转化为正值)。softmax 函数可以帮助我们完成次转化:
\hat{\mathbf{y}} = \mathrm{softmax}(\mathbf{o}) \quad \text{其中}\quad \hat{y}_j = \frac{\exp(o_j)}{\sum_k \exp(o_k)}
这里,对于所有的j总有0 \leq \hat{y}_j \leq 1。因此,\hat{\mathbf{y}}可以视为一个正确的概率分布。softmax运算不会改变未规范化的预测\mathbf{o}之间的顺序,只会确定分配给每个类别的概率。因此,在预测过程中,我们仍然可以用下式来选择最有可能的类别。
\operatorname*{argmax}_j \hat y_j = \operatorname*{argmax}_j o_j.

尽管softmax是一个非线性函数,但softmax回归的输出仍然由输入特征的仿射变换决定。因此,softmax回归是一个线性模型

损失函数

我们的标签使用了独热编码,模型输出通过softmax转化为概率分布,此时我们要计算损失函数,即需要衡量两者概率之间的区别。这里我们使用最大似然估计,根据最大似然估计,对于任何标签\mathbf{y}和模型预测\hat{\mathbf{y}},损失函数为:
l (\mathbf{y}, \hat{\mathbf{y}}) = - \sum_{j=1}^q y_j \log \hat{y}_j

通常此又被成为交叉熵损失。由于\mathbf{y}是一个长度为q的独热编码向量,所以除了类别项以外的所有项j都消失了。由于所有\hat{y}_j都是预测的概率,所以它们的对数永远不会大于0。因此,如果正确地预测实际标签,即如果实际标签P(\mathbf{y} \mid \mathbf{x})=1,则损失函数不能进一步最小化。注意,这往往是不可能的。例如,数据集中可能存在标签噪声(比如某些样本可能被误标),或输入特征没有足够的信息来完美地对每一个样本分类(简单来说,模型输出不可能与独热编码的标签完全相同)。

对于 损失函数的导数:
\partial_{o_j} l(\mathbf{y}, \hat{\mathbf{y}}) =\mathrm{softmax}(\mathbf{o})_j - y_j.
换句话说,导数是我们softmax模型分配的概率与独热标签之间的差异。这使梯度计算在实践中变得容易很多。

Pytorch 实现

import torch
from torch import nn
from d2l import torch as d2l

batch_size = 256
train_iter, test_iter = d2l.load_data_fashion_mnist(batch_size)

# 我们在线性层前定义了展平层(flatten),来调整网络输入的形状
net = nn.Sequential(nn.Flatten(), nn.Linear(784, 10))

def init_weights(m):
    if type(m) == nn.Linear:
        nn.init.normal_(m.weight, std=0.01)

net.apply(init_weights)

loss = nn.CrossEntropyLoss()
trainer = torch.optim.SGD(net.parameters(), lr=0.1)

num_epochs = 10
d2l.train_ch3(net, train_iter, test_iter, loss, num_epochs, trainer)
输出结果
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 220,192评论 6 511
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,858评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 166,517评论 0 357
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,148评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,162评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,905评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,537评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,439评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,956评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,083评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,218评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,899评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,565评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,093评论 0 23
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,201评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,539评论 3 375
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,215评论 2 358

推荐阅读更多精彩内容