线性回归,加权回归,推导过程

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/longgb123/article/details/79079434

一、普通线性回归(OLS)

损失函数:

J(w)=1ni=1n(yiwxi)2=1n||YXw||2

其中: Y w xi 为向量, X 为矩阵。对该损失函数求解如下,即为对 J(w) 函数求 w 的偏导:
J(w)=1n||YXw||2=1n(YXw)T(YXw)=1n(YTwTXT)(YXw)=1n(YTYwTXTYYTXw+wTXTXw)

需要用到的矩阵求导公式为:
dBAdA=BT

dATBdA=B

dATBAdA=2BA

所以, J(w) w 求导得到:
dJ(w)dw=1dw(1n(YTYwTXTYYTXw+wTXTXw))=1n(0XTYXTY+2XTXw)=1n(2XTY+2XTXw)=0

XTX 为可逆矩阵的时候,有解:
w=(XTX)1XTY

因为:
Xw=X(XTX)1XTY=Y^=H^Y

所以,帽子矩阵为:
H^=X(XTX)1XT

二、加权回归

损失函数:

J(w)=1ni=1nαi(yiwxi)2=1nα||YXw||2

同理推导:
J(w)=1nα||YXw||2=1n(YXw)Tα(YXw)=1n(YTwTXT)α(YXw)=1n(YTαYwTXTαYYTαXw+wTXTαXw)

其中, α 为是权重的对角矩阵。对 w 求导得到:
dJ(w)dw=1dw(1n(YTαYwTXTαYYTαXw+wTXTαXw))=1n(0XTαYXTαY+2XTαXw)=1n(2XTαY+2XTαXw)=0

所以,得到:
w=(XTαX)1XTαY

加权的帽子矩阵为:
H^=X(XTαX)1XTα

猜你喜欢

转载自blog.csdn.net/longgb123/article/details/79079434