版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/longgb123/article/details/79079434
一、普通线性回归(OLS)
损失函数:
J(w)=1n∑i=1n(yi−w∗xi)2=1n||Y−X∗w||2
其中:
Y
、
w
、
xi
为向量,
X
为矩阵。对该损失函数求解如下,即为对
J(w)
函数求
w
的偏导:
J(w)=1n||Y−X∗w||2=1n(Y−Xw)T(Y−Xw)=1n(YT−wTXT)(Y−Xw)=1n(YTY−wTXTY−YTXw+wTXTXw)
需要用到的矩阵求导公式为:
dBAdA=BT
dATBdA=B
dATBAdA=2BA
所以,
J(w)
对
w
求导得到:
dJ(w)dw=1dw(1n(YTY−wTXTY−YTXw+wTXTXw))=1n(0−XTY−XTY+2XTXw)=1n(−2XTY+2XTXw)=0
当
XTX
为可逆矩阵的时候,有解:
w=(XTX)−1XTY
因为:
Xw=X(XTX)−1XTY=Y^=H^Y
所以,帽子矩阵为:
H^=X(XTX)−1XT
二、加权回归
损失函数:
J(w)=1n∑i=1nαi(yi−w∗xi)2=1nα||Y−X∗w||2
同理推导:
J(w)=1nα||Y−X∗w||2=1n(Y−Xw)Tα(Y−Xw)=1n(YT−wTXT)α(Y−Xw)=1n(YTαY−wTXTαY−YTαXw+wTXTαXw)
其中,
α
为是权重的对角矩阵。对
w
求导得到:
dJ(w)dw=1dw(1n(YTαY−wTXTαY−YTαXw+wTXTαXw))=1n(0−XTαY−XTαY+2XTαXw)=1n(−2XTαY+2XTαXw)=0
所以,得到:
w=(XTαX)−1XTαY
加权的帽子矩阵为:
H^=X(XTαX)−1XTα