线性回归推导整理

记录一下线性回归推导。以后多写写博客,多记录

线性回归公式 y=w_0x_0+w_1x_1+...+w_nx_n+b{} 其中,w0为参数,x0 为样本值,b 为偏执项

可以记为y = w^Tx+b{}   (1)其中,wT为转置矩阵。

预测样本和真实值之间存在误差 \left | y-\hat{y} \right | = \varepsilon  其中\varepsilon 为误差

对于每个样本都存在误差 \varepsilon_i =\hat{y}_i -w^Tx_i    (2)

假设误差\varepsilon 是服从独立分布的,并且服从高斯分布,则有

p(\varepsilon_i )=\frac{1}{\sqrt{2\pi }\delta }exp(-\frac{\varepsilon_i^2 }{2\delta^2})    (3)

将(2)代入(3)则有条件概率

p\left ( Y=y_i|X=x_i;W=w_i \right )=\frac{1}{\sqrt{2\pi }\delta }exp(-\frac{(\hat{y}_i -w^Tx_i )^2}{2\delta^2})  (4)

在已知条件概率的情况下,可以使用最大似然函数来估计参数,也就说在知道y 和x 的情况下,可以估计w是真实样本的最大概率是多少。

L(w)=\prod p\left ( Y=y_i|X=x_i;W=w_i \right )=\prod \frac{1}{\sqrt{2\pi }\delta }exp(-\frac{(\hat{y}_i -w^Tx_i )^2}{2\delta^2}) (5)

因为求累乘很麻烦,因此转为求对数

logL(w)=log\prod \frac{1}{\sqrt{2\pi }\delta }exp(-\frac{(\hat{y}_i -w^Tx_i )^2}{2\delta^2})(6)

左边转化:

log\sum \frac{1}{\sqrt{2\pi }\delta }exp(-\frac{(\hat{y}_i -w^Tx_i )^2}{2\delta^2})(7)

继续化简

\sum log\frac{1}{\sqrt{2\pi }\delta }+\sum log exp(-\frac{(\hat{y}_i -w^Tx_i )^2}{2\delta^2})(8)

下一步:

m*log\frac{1}{\sqrt{2\pi }\delta }+\sum(-\frac{(\hat{y}_i -w^Tx_i )^2}{2\delta^2})(9)

最终化简为:

logL(w)=m*log\frac{1}{\sqrt{2\pi }\delta }-\frac{1}{2} \frac{1}{\delta^2}\sum(\hat{y}_i -w^Tx_i )^2(10)

因为L(w)是概率值,所以应该要求这个概率值越大越好,公式右边第一项是常数,没有影响,所以就要求第二项越小越好

因此得到目标函数:

J(w) = \frac{1}{2}\sum(\hat{y}_i -w^Tx_i )^2(11)

为了让公式(11)越小越好,因此应该就是求函数的极小值。因为公式(11)是凸函数,因此求极小值就是求导数为0的点

将公式(11)展开

(公式打不上去,就去别去扒了一张图来)

对其求偏导

\frac{\partial J(w)}{\partial w} = \frac{1}{2}\frac{\partial}{\partial w} (y-Xw)^T(y-Xw)(12)

\frac{\partial J(w)}{\partial w} = \frac{1}{2}\frac{\partial}{\partial w} (y^T-w^TX^T)(y-Xw)(13)

对公式(13)化简整理:

\frac{\partial J(w)}{\partial w} = \frac{1}{2}\frac{\partial}{\partial w} (y^Ty-w^TX^Ty-y^TXw+w^TX^TXw)(14)

对矩阵求偏导,最终得到

\frac{\partial J(w)}{\partial w} =X^TXw-Xy(15)

求偏导等于0的点因此为

X^TXw=Xy(16)

w=(X^TX)^{-1}Xy

到此求出参数w 与样本之间的关系。

但这里

(X^TX)^{-1}不一定可逆,因此就涉及到了使用梯度下降的方式来求解。之后写关于梯度下降的博客

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值