TensorFlow训练Logistic回归

Logistic回归

在用线性模型进行回归训练时，有时需要根据这个线性模型进行分类，则要找到一个单调可微的用于分类的函数将线性回归模型的预测值关联起来。这时就要用到逻辑回归，之前看吴军博士的《数学之美》中说腾讯和谷歌广告都有使用logistics回归算法。

如下图，可以清晰看到线性回归和逻辑回归的关系，一个线性方程被逻辑方程归一化后就成了逻辑回归。.

Logistic模型

对于二分类，输出$y \in \{0,1\}$，假如线性回归模型为$z = \theta^Tx$，则要将z转成y，即$y=g(z)$。于是最直接的方式是用单位阶跃函数来表示，即

$$y = \left\{\begin{matrix} 0, & z<0;\\ 0.5, & z=0;\\ 1, & z>0; \end{matrix}\right.$$
如图，

但阶跃函数不连续，于是用sigmoid函数替代之，为

$$y = \frac{1}{1+e^{-z}}$$
如图，

则有，
$$y = \frac{1}{1+e^{-(\theta^Tx)}}$$
即logistics函数，可化为，
$$\ln \frac{y}{1 - y} = \theta_o + \theta_1 x_1 +\theta_2 x_2 +...+\theta_mx_m$$
此即为对数几率回归模型，其中y看成是样本x正例的概率，1-y则为样本x负例的概率，则
$$\ln \frac{p(y=1|x)}{1-p(y=1|x)}=\theta^T x$$

现在要解决的问题是如何求得$\theta$。对于给定样本集$\{(x_i,y_i)\}_{i=1}^{m}$，每个样本出现的概率为，

$$p(y_i,x_i)=p(y_i=1|x_i)^{y_i} (1-p(y_i=1|x_i))^{1-y_i}$$
其中$y_i$为1或0。则样本集出现的似然函数为
$$L(\theta) = \prod_{i=1}^{m}p(y_i,x_i)=\prod_{i=1}^{m}p(y_i=1|x_i)^{y_i} (1-p(y_i=1|x_i))^{1-y_i}$$

对数似然为：

$$l(\theta) = \sum_{i=1}^{m} \ln p(y_i,x_i)$$
$$= \sum_{i=1}^{m} y_i \ln p(y_i=1|x_i) + (1-y_i) \ln (1-p(y_i=1|x_i))$$
$$=\sum_{i=1}^{m} y_i \ln \frac {p(y_i=1|x_i)}{1-p(y_i=1|x_i)} + \sum_{i=1}^{m} \ln (1-p(y_i=1|x_i))$$
$$= \sum_{i=1}^{m} y_i \theta^T x_i - \sum_{i=1}^{m} \ln (1+e^{\theta^T x_i})$$

求对数似然最大化的$\theta$。其中通过求导没办法求得解，所以肯定要用迭代去逼近最优解，可以用梯度下降法或者牛顿法求的解。

实现代码

import tensorflow as tf
from numpy import *

x_train = [[1.0, 2.0], [2.0, 1.0], [2.0, 3.0], [3.0, 5.0], [1.0, 3.0], [4.0, 2.0], [7.0, 3.0], [4.0, 5.0], [11.0, 3.0],
           [8.0, 7.0]]
y_train = [1, 1, 0, 1, 0, 1, 0, 1, 0, 1]
y_train = mat(y_train)

theta = tf.Variable(tf.zeros([2, 1]))
theta0 = tf.Variable(tf.zeros([1, 1]))
y = 1 / (1 + tf.exp(-tf.matmul(x_train, theta) + theta0))

loss = tf.reduce_mean(- y_train.reshape(-1, 1) * tf.log(y) - (1 - y_train.reshape(-1, 1)) * tf.log(1 - y))
train = tf.train.GradientDescentOptimizer(0.01).minimize(loss)

init = tf.initialize_all_variables()

sess = tf.Session()
sess.run(init)
for step in range(1000):
    sess.run(train)
print(step, sess.run(theta).flatten(), sess.run(theta0).flatten())

---

欢迎关注：