线性回归算法

线性回归算法

我的微信公众号: s406205391; 欢迎大家一起学习,一起进步!!!
线性回归是利用数理统计中回归分析,来确定两种或两种以上变量间相互依赖的定量关系的一种统计分析方法,运用十分广泛。

​ 我们看下面一组数据,这组数据通过工资和年龄两个特征来预测银行贷款的额度。

那么,工资和年龄对银行贷款给我们的额度的影响会有多大呢?这也就是我们所需要求的参数

工资 年龄 额度
4000 25 20000
8000 30 70000
5000 28 35000
7500 33 5000
12000 40 85000

通俗的讲:我们通过图片来展示我们的数据,图中X1和X2就是我们的两个特征(年龄,工资),Y是银行最终会借给我们多少钱。我们现在就需要找到最合适的一条线(想象一个高维)来最好的拟合我们的数据点。这时,我们假设θ1是年龄的参数,θ2是工资的参数。我们便可得到一个拟合平面:

h_\theta(x)=\theta_0 + \theta_1 x_1 + \theta_2 x20是偏置项)

​ 将上述公式用矩阵表达,即得:

h_\theta(x)=\sum\limits_{i=0}\limits^{n}\theta_{i}x_{i}=\theta^{T}x

​ 该式展开其实为:h_\theta(x)=\theta_{0}x_{0} + \theta_{1}x_{1} + \theta_{2}x_{2},相较于上式,多引入一个特征x0,这是为了后期方便矩阵的计算。所以我们可以把特征x0统一赋值为1。

​ 因为,真实值和预测值之间肯定是要存在误差的。我们引入ε来表示误差。所以对于每一个样本,即方程即为:y^{(i)}=\theta^{T}x^{(i)} + \epsilon^{(i)}

image

误差满足高斯分布

误差具有如下特点

  • 误差是独立并且具有相同的分布,并且服从均值为0,方差为θ2的高斯分布

  • 独立:误差在样本之间是相互独立的。

  • 同分布。

    既然,误差是服从高斯分布的,那么我们就可以用高斯分布的方程来表示: y^{(i)}=\theta^{T}x^{(i)} + \epsilon^{(i)}

求解目标函数

​ 我们首先将高斯函数,带入我们的方程,即可得:

预测值与误差y^{(i)}=\theta^{T}x^{(i)} + \epsilon^{(i)} (1)

由于误差服从高斯分布p(\epsilon^{(i)})=\frac{1}{\sqrt{2\pi\delta}}e^{-\frac{(\epsilon^{(i)})^2}{2\sigma^{2}}} (2)

将(1)式带入(2)式可得p(y^{(i)}|x^{(i)};\theta)=\frac{1}{\sqrt{2\pi\delta}}e^{-\frac{(y^{(i)}-\theta^{T}x^{(i)})^{2}}{2\sigma^{2}}}

​ 上式我们可以用最大函数求解

似然函数L(\theta)=\prod\limits_{i=1}\limits^{m} p(y^{(i)}|x^{(i)};\theta)= \prod\limits_{i=1}\limits^{m}\frac{1}{\sqrt{2\pi\delta}}e^{-\frac{(y^{(i)}-\theta^{T}x^{(i)})^{2}}{2\sigma^{2}}}

​ 为了方便求解,我们将上式取对数:

对数似然,并化简

logL(\theta)=log\prod\limits_{i=1}\limits^{m} p(y^{(i)}|x^{(i)};\theta)= log\prod\limits_{i=1}\limits^{m}\frac{1}{\sqrt{2\pi\delta}}e^{-\frac{(y^{(i)}-\theta^{T}x^{(i)})^{2}}{2\sigma^{2}}}

=\sum\limits_{i=1}\limits^{m}log\frac{1}{\sqrt{2\pi\delta}} e^{-\frac{(y^{(i)}-\theta^{T}x^{(i)})^{2}}{2\sigma^{2}}}

=mlog\frac{1}{\sqrt{2\pi\delta}} - \frac{1}{\delta^2}\frac{1}{2} \sum\limits_{i=1}\limits^{m}(y^{(i)}-\theta^Tx^{i})^2

目标:去除常数项之后,让目标函数越大越好

J(\theta)=\frac{1}{2}\sum\limits_{i=1}\limits^{m}(y^{(i)} - \theta^T x^{(i)})^2; 即让该式越小越好,可以用最小二乘解

应用最小二乘法求解上式

J(\theta)=\frac{1}{2}\sum\limits_{i=1}\limits^{m}(y^{(i)} - \theta^T x^{(i)})^2=\frac{1}{2}(X\theta-y)^T(X\theta-y)

​ 对该式求偏导即得:

\theta=(x^TX)^{-1}X^Ty

常用的评估方法,得到最佳目标函数的参数

​ 上述我们将目标函数直接求出了结果,但是在某些情况下,我们并不能直接求得结果,因为上述我们对X求其逆矩阵,但是逆矩阵并不是一定能求的。

​ 那么,我们如何评估目标函数拟合的效果呢?我们一般会用到R2项。R2的取值越接近于1,我们认为模型拟合的越好。

最常用的评估项R2: 1-\frac{\sum\limits_{i-1}\limits^{m}(\hat y_{i} - y_i)^2}{\sum\limits_{i=1}\limits^{m}(y_i-y\overline)^2}

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 200,527评论 5 470
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,314评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 147,535评论 0 332
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,006评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,961评论 5 360
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,220评论 1 277
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,664评论 3 392
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,351评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,481评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,397评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,443评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,123评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,713评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,801评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,010评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,494评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,075评论 2 341

推荐阅读更多精彩内容