一、经典最小二乘法
①经典最小二乘法原理介绍
最小二乘法的原理实质比较简单,本质的目的是要确定如下建立的一元线性回归模型的两个回归参数
a1和
b1:
y=a1x+b1
若已知m组样本观测数据
(xi,yi)(i=1,2,3,4⋅⋅⋅m),经典的做法是根据离差平方和来达到一种最小的准则来进行确定的,即确定满足下面条件的
a1′和
b1′,它们使得下面函数取到最小值:
Q(a1′,b1′)=i=1∑m(yi−a1′xi−b1′)2
然后的问题就变得简单了,通过我们上述
Q函数的定义我们不难知道
Q(a1′,b1′)是一个非负函数,关于
a1′,b1′的导数是存在的,通过分别对
a1′,b1′来进行求导并令其为零会得到:
∂a1′∂Q(a1′,b1′)=i=1∑m−2xi(yi−a1′xi−b1′)=0⟶i=1∑mxiyi=i=1∑mb1′xi+a1′i=1∑mxi2⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅①
∂b1′∂Q(a1′,b1′)=i=1∑m−2(yi−a1′xi−b1′)=0⟶i=1∑myi=mb1′+a1′i=1∑mxi⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅②
联立
①,②分别对
m取平均我们很容易就求出来了
a1′和b1′:
a1′=x2−(x)2xy−xˉyˉ
b1′=y−a1′x
二、基于离差概率总和最小对最小二乘法的改进
①改进原理介绍
我们都知道最小二乘法并不能保证所有的样本测试数据点都在整个回归直线上面,而是比较“均匀”的分布在直线两边,我们可以考虑一种概率来重新对数据拟合误差的大小进行一种新规定:
若已知m组样本观测数据
(xi,yi),yi=0(i=1,2,3,4⋅⋅⋅m),我们新规定线性回归模型如下:
y=a2x+b2
下面定义概率
Pi如下:
Pi=yi∣yi−a2xi−b2∣
发现绝对值的存在显然不便于我们处理,我们不妨改进
Pi的定义如下:
Pi(a2,b2)=(yiyi−a2xi−b2)2
而我们所要求的即为满足如下的函数
R的最小值时的
a2和b2:
R(a2,b2)=i=1∑mPi
根据
R函数的定义我们知道
R(a2,b2)是一个非负函数,关于
a2,b2有偏导数令其偏导为零有:
∂a2∂R(a2,b2)=i=1∑m−2yixi(1−yixia2−yib2)=0⟶i=1∑myixi=a2i=1∑myi2xi2+b2i=1∑myi2xi⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅③
∂b2∂R(a2,b2)=i=1∑m−2yi1(1−yixia2−yib2)=0⟶i=1∑myi1=a2i=1∑myi2xi+b2i=1∑myi21⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅④
联立
③,④分别对
m取平均我们很容易就求出来了
a2和b2:
b2=(y2x)×(y2x)−(y21)×(y2x2)(yx)×(y2x)−(y1)×(y2x2)
a2=(y2x)(y1)−b2(y21)
之后通过利用编程来进行计算这两个参数即可,这种方法对我们的初始数据
(xi,yi)中的
yi有着特殊的要求。
三、基于垂线段总和最小对最小二乘法的改进
①改进原理介绍
我们仍假定所求回归直线为:
y=a3x+b3
经典最小二乘法仅仅是考虑的两点之间的直线距离,我们不妨这样考虑:可否将我们所有的样本点
(xi,yi)(i=1,2,3,4⋅⋅⋅m)向回归直线做垂线段,并将此平方总和作为度量指标来替代经典二乘法的距离呢?
据此想法,我们定义
S函数如下:
S(a3,b3)=i=1∑ma32+1(a3xi−yi+b3)2
而我们所要求的即为满足如下的函数
S的最小值时的
a3和b3:
根据
S函数的定义我们知道
S(a3,b3)是一个非负函数,关于
a3,b3有偏导数令其偏导为零有:
∂a3∂S(a3,b3)=i=1∑ma32+12xi(a3xi−yi+b3)(a32+1)−2a3(a3xi−yi+b3)2=0
∂b3∂S(a3,b3)=i=1∑ma32+12(a3xi−yi+b3)=0
即:
i=1∑ma3(xi2−yi2)+i=1∑m(a32−1)xiyi+i=1∑mb3(1−a32)xi+i=1∑m2a3b3yi=a3b32⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⑤
i=1∑m(a3xi−yi+b3)=0⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⑥
而由
⑥我们可知:
a3x−y+b3=0
b3=y−a3x⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⑦
联立
⑤,⑦分别对
m取平均我们很容易就求出来了
a3和b3:
a32(xy−xˉyˉ)+a3(x2−y2−xˉxˉ+yˉyˉ)+(xˉyˉ−xy)
a3=2(xy−xˉyˉ)−(x2−y2−xˉxˉ+yˉyˉ)±(x2−y2−xˉxˉ+yˉyˉ)2−4(xy−xˉyˉ)2
((xy−xˉyˉ)=0)
a3=(x2−y2−xˉxˉ+yˉyˉ)−(xˉyˉ−xy)((xy−xˉyˉ)=0)
再由
⑦可得
b3。
根据上面我们会发现有两个
a3供我们选择,当然我们选择哪一个是要根据题意而定的。并非都取,要结合实际数据来看。
上述为笔者考虑到的最小二乘法的一些改进方案,主要考虑了概率和垂线两种情况,未尝难免出现运算错误和稳定性考虑不周的情况,衷心恳请大家批评指正,笔者在此深表感谢!