如果你在运营一个2C的平台，那么你肯定关心用户流失的问题。腾讯有个产品叫信鸽Pro，它能够通过对用户往期行为的挖掘，预测用户潜在的流失（付费）行为，进而实现精准营销。据说，腾讯自己的手游就是用这个系统做用户分析的。

信鸽Pro获取大量用户数据，提取用户特征，然后通过算法建模，评估出用户可能的行为。算法建模中最基础的一步就是对用户进行分类。这里就介绍一种常用的分类算法 － 逻辑回归。

模型

用户数据比较复杂，这里用平面上的点举例。假设平面上有一些点，如图所示：

整个平面上只有两种图形，一种是三角形，另一种是圆形。可以把它们想象为两种不同的用户，比如活跃用户／非活跃用户。

问题：如果随意在这个平面新增加一个点, 比如点P(5,19)，那怎知把它归到哪一组更合适？可以想象为对新用户的预测。

思路

我们发现，三角形大都位于左上方，而圆形大都位于右下方。我们可以用尺子在图上画一条直线，该直线尽可能的将三角形和圆形分到两边。然后观察新点位于哪一侧。若与三角形在同一侧，则它应该属于三角形；若位于圆形一侧，则应属于圆形。在本例中，坐标P应该属于三角形更合适。

这个问题似乎很简单。但是，如果三维空间存在类似的问题，答案就没有那么显而易见了。那4维空间呢? 1024维空间呢?

不过别担心！借助计算机算法，N维空间分类的问题已经很容易解决，逻辑回归就是常用的一种。

逻辑回归

逻辑回归的核心思想就是通过现有数据，对分类边界线建立回归公式，以此进行分类。

在介绍算法之前，需要先介绍一个函数: Sigmoid函数。

Sigmoid函数

Sigmoid函数的表达式为：

在坐标系中的图形为：

x>0时，x越大y越接近于1；x<0时，x越小，y越接近于0。如果把坐标拉长，曲线中间就会很“陡”。直观上x的“轻微”变化，都会导致y接近于0或1。

Sigmoid函数的作用是将任意实数转换成0～1的数，而0和1刚好可以用做分类，比如，用1表示三角形，用0表示圆形。小于0.5的可以划分为0类，大于0.5的划分为1类。（注：Sigmoid是单调增长函数，因而多个数字通过Sigmoid转换后相对位置不变，这是选择该函数的重要原因。）

分析步骤

 

简化模型

为便于分析，把模型中的坐标简化一些。下面的六个坐标点和一条分割线：

其中红色三角形坐标分别是(1,2)、(1.5,7)和(2,6)。绿色圆点坐标分别是（1,0）、(2,3)和(2.5,6)。分割线的函数为y=4x-3. 它的形式还可以转换成：3-4x+y=0 。

我们设表达式f(x,y) = 3-4x+y

把六个点的坐标代到这个方程式里,有

<表1>

(注：标识1表示三角形；标识0表示圆形)

f(x,y)>0的点在分割线上方，是三角形；f(x,y)<0的点在分割线下方，是圆形.

如果有个三角形的坐标是(2,4.5),那这个点的f(x,y)值等于-0.5，这个点就被分割线错误划分了。

现在的问题是，我们只有一些坐标以及这些坐标的分类信息，如何找到一条最优的分割线,使得尽可能少的点被错误划分？

 

损失函数

损失函数 (Loss Function) 的作用是判断直线错误划分数据的程度。一种方法是计算被错误划分的点的个数，错误点越少，直线越好。但，这种方法很难优化。另一种方法是计算点到直线的距离。

如果是一个平面来划分三维空间的点，那距离公式为

一般的，n维空间上一个点到超平面的距离为

w是超平面的参数向量,

x是超平面的自变量,

b是截距

超平面函数：

表示x向量的第i个元素（特性）；后面会用到，表示空间中的第i个点。

为了方便计算，一般在x中增加一个元素1，w中增加一个元素w0=b

于是超平面函数变为：

距离公式变为：

超平面上方的点f(x)>0, 下方的点f(x)<0,因此点到超平面的距离（分正负）:

d是一个负无穷到正无穷的数。

通过sigmoid函数，将d变成一个0～1的值，设h = sigmoid(d)。若d为正且越大，h越接近于1，也就越应该属于三角形（分类1）；若d为负，且绝对值越大，h越接近于0，该点也就越应该属于圆形(分类0)。因此，h越接近于分类标识，划分的准确性越高。

设第j个点的分类表示为，那么下面的公式就表示点j被错误划分的概率。

我们把损失函数设定为所有点被错误划分的平均概率

平方是为了保证概率为正，前面的1/2是为了求导数后消除参数。

那么，问题转化成：找到w的一个值，使得损失函数的值最小。

 

用梯度下降法求w

所谓梯度，就是函数在某个点增长最快的方向，有时称为斜度。如果函数是一个曲线，某个点的梯度就是该点的斜率，或导数。

如果是曲面，梯度是在两个方向上求偏导。

梯度下降法的核心思想是：欲求函数的最小值，先从某一点出发，沿着函数的梯度的反方向探寻，直到找到最小值。设每次探寻Delta(w)，步长为alpha，则梯度下降的算法的公式为：

 

求导

用梯度下降法需要先对损失函数求导，我们的损失函数被分成三部分：

•  ——— (1)

•  ——— (2)

•  ———- (3)

可以通过复合函数求导法对损失函数求偏导：

梯度公式重点关注的是导数的符号，所以这里可以简化一下。函数（2）是单调递增函数，所以导数是正数，不影响整体导数的符号，可以去除。 公式（3）的分母是正数，也不影响导数的符号，也可以去掉。最后得：

代入梯度下降算法公式得：

1/m为正数，也可用去掉。

 

代码

loadData()函数返回坐标值和分类标识。第一个返回值取前三列 x0,x1,x2；第二个返回值取第四列，即label

from numpy import *def loadData():  data = [[1, 1, 2, 1],          [1, 2, 6, 1],          [1, 1.5, 7, 1],          [1, 1, 0, 0],          [1, 2, 3, 0],          [1, 2.5, 6, 0]]  ds = [e[0:3] for e in data]  label = [e[-1] for e in data]  return ds, label

Sigmoid函数

def sigmoid(x):    return 1.0/(1+exp(-x))

梯度下降算法

input ds: 坐标数据; label: 标签return w: 系数向量, nx1的矩阵

def reduce(ds, label):  #转换成矩阵  dmat = mat(ds)  lmat = mat(label).T  #mxn的矩阵的行数和列数  m,n = shape(dmat)  #步长  alpha = 0.1  #循环次数  loops = 200  #初始化w为[1,1,1],即分割线为 1+x+y=0  w = ones((n,1))  for i in range(loops):     h = sigmoid(dmat*w)    err = (h - lmat)    w = w - alpha * dmat.T* err   return w.A[:,0]

测试

def test():  ds,l = loadData()  print reduce(ds,l)

运行结果

>>> import lr>>> lr.test()[ 3.1007773  -5.54393712  1.60563033]

 

也就是说w=(3.1, -5.5, 1.6), 即w0=3.1, w1=-5.5, w2 = 1.6

分割线的表达式为：w0+w1x+w2y=0， 代入w后得 3.1-5.5x+1.6y=0, 即y=3.44x-1.9 。 见下图，该直线正确地将图形划分开。

执行过程

w的初始值为(1,1,1),也就是0号线。

每次循环都会调整分割线的位置，执行到第200次的时候，分割线已经能够很好对坐标分组了。 <分割线调整图> 线编号n表示第n次调整(循环)之后的位置

应用

把上面的x,y转换成用户特征，比如登录时间，登录频率等等。把三角形和圆形转换成付费用户和免费用户，就得到了付费用户预测模型；把三角形和圆形转换成流失用户和有效用户，就得到了流失用户预测模型。

当然，这只是个理论模型，实际应用要比这复杂的多的多。

 

相关文章