PCA无监督,LDA有监督
在PCA一文中,我们简单地提到了
如果是二维空间中的样本点,那么我们就是求解出一条直线使得样本投影到该直线上的方差最大。从回归的角度来看其实就是求解出一个线性函数来拟合样本点集合。
所以我们可以从维度降低的角度来考察线性分类器。考虑二分类的情形,假设我们有一个D维输入向量
x,然后我们使用
y=wTx投影到一维。我们设置一个阈值,有
N1个点大于它是
C1类,有
N2个点小于是
C2类,这就得到了一个标准的线性分类器。两类的均值向量分别为
m1=N11n∈C1∑xn
m2=N21n∈C2∑xn投影之后最简单的度量类别之间分开程度的方法就是类别均值投影之后的距离,我们要选择一个合适的
w使得距离最大即
warg max wT(m2−m1)由于只考虑了均值即偏差并没有考虑方差(通常需要权衡利弊两者的关系),而本身又是降维的方式可以想象得到这种方法会存在一个问题即高维中可分,降维后两个类别可能会挤在一起。如果类概率分布的协方差矩阵与对角化矩阵差距较大时,这个问题尤为突出。
Fisher提出的思想是最大化一个函数,这个函数能够让类均值的投影分开得尽可能大,同时能够让每个类别内部的方差较小,从而最小化了类别的重叠。
经过变换后,类内方差可以表示为
sk2=n∈Ck∑(wTxn−wTmk)2Fisher准则根据类间方差与类内方差的比值定义
J(w)=s12+s22(wT(m2−m1))2=wTSWwwTSBw
SB=(m2−m1)(m2−m1)T
SW=n∈C1∑(xn−m1)(xn−m1)T+n∈C2∑(xn−m2)(xn−m2)T对J求w的导数
∂w∂J(w)=0
(wTSBw)SWw=(wTSWw)SBw忽略标量因子
(wTSBw)和
(wTSWw),
SBw是在
(m2−m1)方向上。所以**
w∝SW−1(m2−m1)**,这就是Fisher线性判别函数。
对于二分类问题,Fisher准则可以看成最小平方的一个特例
我们假设
C1类的目标值为
N1N,
C2类的目标值为
−N2N,平方和误差函数
E=21n=1∑N(wTxn+wo−tn)2求偏导
∂w∂E=n=1∑N(wTxn+wo−tn)xn=0
∂wo∂E=n=1∑N(wTxn+wo−tn)=0先化解第二个式子,
Nwo+n=1∑NwTxn=n=1∑Ntn=N1N1N−N2N2N=0
wo=−wTmm=N1n=1∑Nxn=N1(N1m1+N2m2)然后把
wo带入到
∂w∂E
n=1∑N(wTxn−wTm−tn)xn=0
n∈C1∑(xnxnT−mxn)w−n∈C1∑tnxn+n∈C2∑(xnxnT−mxn)w−n∈C2∑tnxn=0
(n∈C1∑xnxnT−N1m1mT)w−N1NN1m1+(n∈C2∑xnxnT−N2m2mT)w+N2NN2m2=0
{(n∈C1∑xnxnT+n∈C2∑xnxnT)−(N1m1+N2m2)mT}w=N(m1−m2)我们已经推导出了
SW=n∈C1∑(xn−m1)(xn−m1)T+n∈C2∑(xn−m2)(xn−m2)T
SW=(n∈C1∑xnxnT+n∈C2∑xnxnT)−(n∈C1∑m1xnT+xnm1T)−(n∈C2∑m2xnT+xnm2T)+(n∈C1∑m1m1T+n∈C2∑m2m2T)
SW=(n∈C1∑xnxnT+n∈C2∑xnxnT)−(N1m1m1T+N1m1m1T)−(N2m2m2T+N2m2m2T)+(N1m1m1T+N2m2m2T)
(n∈C1∑xnxnT+n∈C2∑xnxnT)=SW+N1m1m1T+N2m2m2T
把这个代入
{(n∈C1∑xnxnT+n∈C2∑xnxnT)−(N1m1+N2m2)mT}w=N(m1−m2)可得
{SW+N1m1m1T+N2m2m2T−(N1m1+N2m2)N1(N1m1+N2m2)T}w=N(m1−m2)
{SW+NN1N2[N2Nm1m1T+N1Nm2m2T−N2N1m1m1T−m1m2T−m2m1T−N1N2m2m2T]}w=N(m1−m2)
{SW+NN1N2[m1m1T+m2m2T−m1m2T−m2m1T]}w=N(m1−m2)
{SW+NN1N2SB}w=N(m1−m2)所以从这个式子中还是能看出
w∝SW−1(m2−m1)