学习笔记 TF010:softmax 分类

回答多选项问题，使用 softmax 函数，对数几率回归在多个可能不同值上的推广。函数返回值是 C 个分量的概率向量，每个分量对应一个输出类别概率。分量为概率，C 个分量和始终为 1。每个样本必须属于某个输出类别，所有可能样本均被覆盖。分量和小于 1,存在隐藏类别；分量和大于 1,每个样本可能同时属于多个类别。类别数量为 2,输出概率与对数几率回归模型输出相同。

变量初始化，需要 C 个不同权值组，每个组对应一个可能输出，使用权值矩阵。每行与输入特征对应，每列与输出类别对应。

鸢尾花数据集 Iris，包含 4 个数据特征、3 类可能输出，权值矩阵 4X3。

训练样本每个输出类别损失相加。训练样本期望类别为 1,其他为 0。只有一个损失值被计入，度量模型为真实类别预测的概率可信度。每个训练样本损失相加，得到训练集总损失值。TensorFlow 的 softmax 交叉熵函数，sparse_softmax_cross_entropy_with_logits 版本针对训练集每个样本只对应单个类别优化，softmax_cross_entropy_with_logits 版本可使用包含每个样本属于每个类别的概率信息的训练集。模型最终输出是单个类别值。

不需要每个类别都转换独立变量，需要把值转换为 0~2 整数(总类别数 3)。tf.stack 创建张量，tf.equal 把文件输入与每个可能值比较。tf.argmax 找到张量值为真的位置。

推断过程计算测试样本属于每个类别概率。tf. argmax 函数选择预测输出值最大概率类别。tf.equal 与期望类别比较。tf.reduce_meen 计算准确率。

import tensorflow as tf#导入 TensorFlow 库
import os#导入 OS 库
W = tf.Variable(tf.zeros([4, 3]), name="weights")#变量权值，矩阵，每个特征权值列对应一个输出类别
b = tf.Variable(tf.zeros([3], name="bias"))#模型偏置，每个偏置对应一个输出类别
def combine_inputs(X):#输入值合并
    print "function: combine_inputs"
    return tf.matmul(X, W) + b
def inference(X):#计算返回推断模型输出(数据 X)
    print "function: inference"
    return tf.nn.softmax(combine_inputs(X))#调用 softmax 分类函数
def loss(X, Y):#计算损失(训练数据 X 及期望输出 Y)
    print "function: loss"
    return tf.reduce_mean(tf.nn.sparse_softmax_cross_entropy_with_logits(logits=combine_inputs(X), labels=Y))#求平均值，针对每个样本只对应单个类别优化
def read_csv(batch_size, file_name, record_defaults):#从 csv 文件读取数据，加载解析，创建批次读取张量多行数据
    filename_queue = tf.train.string_input_producer([os.path.join(os.getcwd(), file_name)])
    reader = tf.TextLineReader(skip_header_lines=1)
    key, value = reader.read(filename_queue)
    decoded = tf.decode_csv(value, record_defaults=record_defaults)#字符串(文本行)转换到指定默认值张量列元组，为每列设置数据类型
    return tf.train.shuffle_batch(decoded, batch_size=batch_size, capacity=batch_size * 50, min_after_dequeue=batch_size)#读取文件，加载张量 batch_size 行
def inputs():#读取或生成训练数据 X 及期望输出 Y
    print "function: inputs"
    #数据来源： https://archive.ics.uci.edu/ml/datasets/Iris
    #iris.data 改为 iris.csv，增加 sepal_length, sepal_width, petal_length, petal_width, label 字段行首行
    sepal_length, sepal_width, petal_length, petal_width, label =\
        read_csv(100, "iris.csv", [[0.0], [0.0], [0.0], [0.0], [""]])
    #转换属性数据
    label_number = tf.to_int32(tf.argmax(tf.to_int32(tf.stack([
        tf.equal(label, ["Iris-setosa"]),
        tf.equal(label, ["Iris-versicolor"]),
        tf.equal(label, ["Iris-virginica"])
    ])), 0))#将类名称转抽象为从 0 开始的类别索引
    features = tf.transpose(tf.stack([sepal_length, sepal_width, petal_length, petal_width]))#特征装入矩阵，转置，每行一样本，每列一特征
    return features, label_number
def train(total_loss):#训练或调整模型参数(计算总损失)
    print "function: train"
    learning_rate = 0.01
    return tf.train.GradientDescentOptimizer(learning_rate).minimize(total_loss)
def evaluate(sess, X, Y):#评估训练模型
    print "function: evaluate"
    predicted = tf.cast(tf.arg_max(inference(X), 1), tf.int32)#选择预测输出值最大概率类别
    print sess.run(tf.reduce_mean(tf.cast(tf.equal(predicted, Y), tf.float32)))#统计所有正确预测样本数，除以批次样本总数，得到正确预测百分比
with tf.Session() as sess:#会话对象启动数据流图，搭建流程
    print "Session: start"
    tf.global_variables_initializer().run()
    X, Y = inputs()
    total_loss = loss(X, Y)
    train_op = train(total_loss)
    coord = tf.train.Coordinator()
    threads = tf.train.start_queue_runners(sess=sess, coord=coord)
    training_steps = 1000#实际训练迭代次数
    for step in range(training_steps):#实际训练闭环
        sess.run([train_op])
        if step % 10 == 0:#查看训练过程损失递减
            print str(step)+ " loss: ", sess.run([total_loss])
    print str(training_steps) + " final loss: ", sess.run([total_loss])
    evaluate(sess, X, Y)#模型评估
    coord.request_stop()
    coord.join(threads)
    sess.close()

参考资料： 《面向机器智能的 TensorFlow 实践》

欢迎加我微信交流：qingxingfengzi

我的微信公众号:qingxingfengzigz

我老婆张幸清的微信公众号：qingqingfeifangz