R语言——回归(初级)

0. 回归的定义

回归分析(regression analysis)是确定两种或两种以上变量间相互依赖的定量关系的一种统计分析方法。运用十分广泛,回归分析按照涉及的变量的多少,分为一元回归和多元回归分析;按照自变量的多少,可分为简单回归分析和多重回归分析 [1] ;按照自变量因变量之间的关系类型,可分为线性回归分析和非线性回归分析。如果在回归分析中,只包括一个自变量和一个因变量,且二者的关系可用一条直线近似表示,这种回归分析称为一元线性回归分析。如果回归分析中包括两个或两个以上的自变量,且自变量之间存在线性相关,则称为多重线性回归分析。

R中做回归分析的函数已经超过200个了

1. OLS回归

OLS(普通最小二乘法),包括简单的线性回归、多项式回归和多元线性回归。OLS回归是目前最常见的统计分析方法。OLS回归是通过预测变量的加权和来预测量化的因变量,其中权重是通过数据估计而得的参数。

OLS回归你和模型的形式:
\hat{Y}_i=\hat{\beta}_0+\hat{\beta}_1X_{1i}+...+\hat{\beta}_kX_{ki}
其中i=1...n,n为观测的数目,k为预测变量的数目。

我们的目标是通过减少相应变量的真实值与预测值的差值来获得模型参数(截距项和斜率),具体而言就是使得残差平方和最小。
\sum_{i=1}^n(Y_i-\hat{Y}_i)=\sum_{i=1}^n(Y_i-\hat{\beta}_0-\hat{\beta}_1X_{1i}-...-\hat{\beta}_kX_{ki})^2=\sum_{i=1}^n\varepsilon^2_i
为了能够恰当的解释OLS模型的系数,数据必须满足以下的假设统计:

  • 正态性 对于固定的自变量的值,因变量值成正太分布。
  • 独立性 Y_i之间相互独立。
  • 线性 因变量和自变量之间为线性相关。
  • 同方差性 因变量的方差不随自变量的水平不同而发生变化。

1.1 使用lm()拟合回归模型

myfit <- lm(formula , data)

1.2 简单线性回归

R基础安装中的数据集women提供了15个年龄在30-39之间女性的身高和体重信息。我们通过身高来预测体重,获得一个灯饰可以帮助我们分辨出那些过重或者过轻的个体。

fit <- lm(weight ~ height, data=women)
summary(fit)
women$weight
fitted(fit)
residuals(fit)
plot(women$height,women$weight,
     main="Women Age 30-39", 
     xlab="Height (in inches)", 
     ylab="Weight (in pounds)")
# add the line of best fit
abline(fit)
简单线性回归反馈.png
  • summary()函数返回了回归拟合的基础汇总信息,我们可以从中得知截距、权重分别为多少Intercept=-87.51667,W(height)=3.45;在 Pr(>|t|)栏中可以看到回归系数(3.45)显著不为0(p<0.001),主要是后面***表示显著可信;Multiple R-squared参数为0.991,表明模型可以解释体重99.1%的方差。
  • fitted()函数返回的是回归拟合的预测值。
  • residuals()函数返回的是残差值。
    \hat{Weight}=-87.51667+3.45*Height

最后利用plot(0函数,对拟合进行可视化表述:


简单线性回归图例.jpeg

1.3 多项式回归

使用曲线往往比单纯的使用直线更为有效一些,可以通过增加一个平方项来提高预测的精度。

fit2 <- lm(weight ~ height + I(height^2), data=women)
summary(fit2)
plot(women$height,women$weight,
     main="Women Age 30-39",
     xlab="Height (in inches)",
     ylab="Weight (in lbs)")
lines(women$height,fitted(fit2))
多项式回归反馈.png

由summay()函数可以明确的看到Multiple R-squared: 0.9995,且在p<0.001水平下,回归系数均非常显著,综上可以得出包含二次项提高了模型的拟合程度。新的预测公式为:
\hat{Weight}=261.88-7.348×Height+0.083×Height^2

多项式回归.jpeg

1.4 多元线性回归

当预测变量不止一个时,简单线性回归见变成了多元线性回归。由于与多项式线性回归相似,此处就不再赘述。

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

推荐阅读更多精彩内容

  • 首页 资讯 文章 资源 小组 相亲 登录 注册 首页 最新文章 IT 职场 前端 后端 移动端 数据库 运维 其他...
    Helen_Cat阅读 3,837评论 1 10
  • R中的线性回归函数比较简单,就是lm(),比较复杂的是对线性模型的诊断和调整。这里结合Statistical Le...
    真依然很拉风阅读 65,141评论 1 64
  • 谷雨连春夏,景色最是佳。 无须添国色,一杯雨前茶。
    流云的天空阅读 449评论 2 3
  • 春天 春天是窗前景和惆怅 窗前景里有远山如黛 和一颗清风照料的种子 而惆怅却没有释义 追寻它的源头 可能是从你第一...
    郁结阅读 245评论 0 1
  • 我还在睡午觉。突然听到老张说:“外面下雨了”。老张急冲冲的跑到顶楼收起晾晒的被单。“还好我跑的快,不然就淋湿啦”。...
    言若NN阅读 561评论 3 2