线性回归,加权回归

服了,简书又不支持latex公式了。
本文内容看博客:http://blog.csdn.net/longgb123/article/details/79079434
下面是markdown格式的原文:

## 一、普通线性回归(OLS)
损失函数:
$$ J(w)=\frac{1}{n}\sum_{i=1}^{n} (y_{i}-w*x_{i})^{2}=\frac{1}{n}||Y-X*w||^{2} $$
其中:$Y$、$w$、$x_{i}$为向量,$X$为矩阵。对该损失函数求解如下,即为对$J(w)$函数求$w$的偏导:
$$\begin{aligned} 
J(w)=\frac{1}{n}||Y-X*w||^{2}&=\frac{1}{n}(Y-Xw)^{T}(Y-Xw)
\\&=\frac{1}{n}(Y^{T}-w^{T}X^{T})(Y-Xw)
\\&=\frac{1}{n}(Y^{T}Y-w^{T}X^{T}Y-Y^{T}Xw+w^{T}X^{T}Xw)
\end{aligned}$$
需要用到的矩阵求导公式为:
$$\frac{dBA}{dA}=B^{T}$$
$$\frac{dA^{T}B}{dA}=B$$
$$\frac{dA^{T}BA}{dA}=2BA$$
所以,$J(w)$对$w$求导得到:
$$\begin{aligned} 
\frac{dJ(w)}{dw}&=\frac{1}{dw}\left(\frac{1}{n}(Y^{T}Y-w^{T}X^{T}Y-Y^{T}Xw+w^{T}X^{T}Xw)\right)
\\&=\frac{1}{n}(0-X^{T}Y-X^{T}Y+2X^{T}Xw)
\\&=\frac{1}{n}(-2X^{T}Y+2X^{T}Xw)=0
\end{aligned}$$
当$X^{T}X$为可逆矩阵的时候,有解:
$$w=(X^{T}X)^{-1}X^{T}Y$$
因为:
$$Xw=X(X^{T}X)^{-1}X^{T}Y=\hat{Y}=\hat{H}Y$$
所以,帽子矩阵为:
$$\hat{H}=X(X^{T}X)^{-1}X^{T}$$

## 二、加权回归
损失函数:
$$J(w)=\frac{1}{n}\sum_{i=1}^{n}\alpha_{i}(y_{i}-w*x_{i})^{2}=\frac{1}{n}\alpha||Y-X*w||^{2}$$
同理推导:
$$\begin{aligned} 
J(w)=\frac{1}{n}\alpha||Y-X*w||^{2}&=\frac{1}{n}(Y-Xw)^{T}\alpha(Y-Xw)
\\&=\frac{1}{n}(Y^{T}-w^{T}X^{T})\alpha(Y-Xw)
\\&=\frac{1}{n}(Y^{T}\alpha Y-w^{T}X^{T}\alpha Y-Y^{T}\alpha Xw+w^{T}X^{T}\alpha Xw)
\end{aligned}$$
其中,$\alpha$为是权重的对角矩阵。对$w$求导得到:
$$\begin{aligned} 
\frac{dJ(w)}{dw}&=\frac{1}{dw}\left(\frac{1}{n}(Y^{T}\alpha Y-w^{T}X^{T}\alpha Y-Y^{T}\alpha Xw+w^{T}X^{T}\alpha Xw)\right)
\\&=\frac{1}{n}(0-X^{T}\alpha Y-X^{T}\alpha Y+2X^{T}\alpha Xw)
\\&=\frac{1}{n}(-2X^{T}\alpha Y+2X^{T}\alpha Xw)=0
\end{aligned}$$
所以,得到:
$$w=(X^{T}\alpha X)^{-1}X^{T}\alpha Y$$
加权的帽子矩阵为:
$$\hat{H}=X(X^{T}\alpha X)^{-1}X^{T}\alpha$$
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • Android 自定义View的各种姿势1 Activity的显示之ViewRootImpl详解 Activity...
    passiontim阅读 176,621评论 25 709
  • 在此特此声明:一下所有链接均来自互联网,在此记录下我的查阅学习历程,感谢各位原创作者的无私奉献 ! 技术一点一点积...
    远航的移动开发历程阅读 13,898评论 12 197
  • 弥补黑暗的空虚 延伸白昼的希望 你是女神 你是救世主 让女人痛苦 让李太白着魔 你是祸水 你是扫帚星 给人团圆 给...
    年轻的风阅读 1,616评论 0 1
  • 楼中拟把同窗看,然已不见人些许。一届七十网新子。还可多少相记叙。 离歌本以尚且远,人却匆散忙终曲。驹隙一过十年期,...
    千年不陈皮阅读 1,648评论 0 2

友情链接更多精彩内容