8.2 一元线性回归模型及其应用
散点大致排成一条斜线时,画一条「最贴近」它们的直线,就能用 x 估计 y。这一节讲:哪条线算最贴近、怎么算出来、用它预测时要小心什么。
本节知识点统计建模的检验与报告解释一元线性回归模型的含义最小二乘与经验回归方程残差及结构检查拟合效果、预测与外推变换后的非线性经验回归
第 1 步
从例子进入
小卖部记下了 5 天的最高气温和卖出的冰棍数。先点下面的按钮,看顺着这 5 天的散点画出趋势线,延伸到 30 ℃ 会指向哪里:
| 气温 / ℃ |
20 |
22 |
24 |
26 |
28 |
| 冰棍 / 根 |
40 |
46 |
50 |
55 |
59 |
想一想
看答案与解析
答案:C
C。气温每高 2 ℃,大约多卖 5 根,从 28 ℃ 的 59 根再往上走一步,约 64 根。
你其实在心里画了一条直线:五个点并不正好在一条线上(每一步多卖 6,4,5,4 根),但大致沿着一条斜线往上走。
这一节把这件事做准:在所有直线里,按一个明确的标准挑出最贴近这些点的那一条(最小二乘法),再用它估计。也要讲清楚:这样的估计是「平均来说大约多少」,不是「一定是多少」。
第 2 步
讲清楚
1. 先认识几个词
要用到:8.1 散点图与相关系数(回看 8.1);必修一 2.3 二次函数最值(回看 2.3)。线性回归是用一条最优直线去拟合散点的宏观规律。
2. 经验回归主公式与残差计算法则
想一想
自己试试:某校调查了 50 名高三男生的身高 x(cm)与体重 y(kg),建立经验回归方程 y^=0.8x−85。关于式中的 0.8,下列说法正确的是?
看答案与解析
答案:B
B。b^ 带帽子,是样本算出的估计值,换一个样本,估计值可能改变。A 把估计当成了总体参数;C 漏了「平均」:具体两个人之间,体重差多少还受 e 影响。
想一想
自己试试:一组数据 xˉ=4,yˉ=10,算得 b^=2。a^ 是?
看答案与解析
答案:A
A。a^=yˉ−b^xˉ=10−2×4=2。可以验:x=4 时 y^=2×4+2=10=yˉ,直线过样本中心。选 B 的同学把减号写成了加号;选 C 的同学忘了乘 b^。
想一想
自己试试:经验回归方程 y^=0.5x+3,数据点 (4,6) 的残差是?
看答案与解析
答案:A
A。y^=0.5×4+3=5,残差 =y−y^=6−5=1,点在直线上方。选 B 是减反了(残差是「真实减预测」);选 C 只算了预测值。
3. 三步破题思维链例题
4. 三大翻车现场红绿对比卡
第 3 步
动手
① 自己画一条线,和最小二乘比一比
拖动绿线两端的空心圆(或选中后按 ↑↓)改变直线;金色竖线是残差。黑点也能拖(选中后按方向键移动)。菱形是样本中心 (xˉ,yˉ)。试着让 Q 尽量小,再点「显示最小二乘直线」对答案。
② 进阶:残差图、R2 和外推
选一组数据。上图是散点和经验回归直线(虚线部分在数据范围之外),下图是残差图:横轴 x,纵轴残差。拖滑块改变要预测的 x0,看它跑出数据范围会怎样。
第 4 步
练一练
从易到难排好了。前两道直接用,后三道是最常见的错。
直接用
想一想
- 一组成对数据算得 xˉ=2,yˉ=5,∑(xi−xˉ)(yi−yˉ)=12,∑(xi−xˉ)2=4。经验回归方程是?
看答案与解析
答案:A
A。b^=412=3,a^=5−3×2=−1。验:x=2 时 y^=5,过样本中心。B 是 a^ 的减号写成了加号;C 是把分子分母颠倒了;D 是把 yˉ 直接当成了 a^。
想一想
- 某店广告费 x(万元,数据在 2~8 之间)与销售额 y(万元)的经验回归方程是 y^=2.5x+4。广告费 6 万元时,销售额估计约多少万元?
看答案与解析
答案:B
B。2.5×6+4=19,6 在数据范围里。意思是「广告费 6 万元时,销售额平均约 19 万元」,某一次具体是多少会上下浮动。A 漏加了 4;D 是把 6 当成了 y 去反解。
易错
想一想
看答案与解析
答案:C
C。最小二乘量的是竖直方向的偏差:同一个 xi,真实 yi 和直线上的 bxi+a 差多少,再平方相加。「点到直线的距离」是垂直于直线量的,是另一回事,A、B 都不对。
D 也不行:只要直线过样本中心,不管斜率多少,残差之和都是 0,挑不出最好的那条。在「动手」①里点「画出残差的平方」,看到的金色方块就是要压小的东西。
想一想
- 调查 7~12 岁的孩子,得到身高 y(cm)对年龄 x(岁)的经验回归方程 y^=6x+76。哪句话对?
看答案与解析
答案:C
C。经验回归方程给的是平均趋势的估计。A 错在「一定」:个体有随机误差 e,同岁孩子高矮不一。B 错在外推:25 岁远在 7~12 岁的数据范围外,那里身高早已不按每年 6 cm 涨了。「动手」②里把 x0 拖出数据范围,可以看到弯曲数据的直线预测差得多远。
想一想
- 三对数据 (1,1),(2,3),(3,2),y 对 x 的经验回归方程是 y^=0.5x+1。现在知道 y=3,想估计 x,应该用哪个结果?
看答案与解析
答案:B
B。最小二乘压的是 y 方向的偏差;要估计 x,该压的是 x 方向的偏差,得重新拟合:xˉ=yˉ=2,∑(xi−xˉ)(yi−yˉ)=1,∑(yi−yˉ)2=2,所以 x^=21(y−2)+2=0.5y+1,y=3 时 x^=2.5。
反解得 4,差得很远(数据都在 1~3 之间)。只有点全在一条直线上时,两种做法才一样。
第 5 步
高考真题
先独立作答,做完再对照分步解析。选择题点选项自动判对错;填空题写法不同时,看答案自己判。要用到后面章节的题收在最后,学到那里再做。
2016 全国Ⅲ卷(理科) 第 18 题 · 解答
我国2008—2014年生活垃圾无害化处理量 y(单位:亿吨)的折线图近似呈线性上升趋势,年份代码 t=1,2,…,7 分别对应2008—2014年。题给参考数据为 ∑i=17yi=9.32,∑i=17tiyi=40.17,∑i=17(yi−yˉ)2=0.55,7≈2.646。 建立y关于t的回归方程(系数精确到0.01),预测2016年我国生活垃圾无害化处理量。附:y^=a^+b^t,b^=∑(ti−tˉ)2∑(ti−tˉ)(yi−yˉ),a^=yˉ−b^tˉ。
答案与解析
答案:y^=0.92+0.10t;按该方程预测2016年为1.82亿吨。
思路
先计算最小二乘系数,再按照题目要求将系数保留到0.01,最后代入正确年份代码。
- 计算中心化统计量
tˉ=4,yˉ=9.32/7,Stt=28,Sty=40.17−4×9.32=2.89。
- 求斜率与截距
b^=2.89/28≈0.103214,a^=9.32/7−4b^≈0.918571。系数精确到0.01,得 y^=0.92+0.10t。
- 换算年份并预测
2016年对应 t=2016−2007=9,由上述回归方程得 y^=0.92+0.10×9=1.82 亿吨。
易错点
- 时间变量是年份代码,不是实际公历年份。
- 题目要求系数精确到0.01,应说明预测采用的是舍入后的方程,避免把精度差异误判为算错。
小结
2017 山东卷(理科) 第 5 题 · 单选
为了研究某班学生的脚长x(单位:厘米)和身高y(单位:厘米)的关系,从该班随机抽取10名学生,根据测量数据的散点图可以看出y与x之间有线性相关关系。设其回归直线方程为 y^=b^x+a^,已知 ∑i=110xi=225,∑i=110yi=1600,b^=4。该班某学生的脚长为24,据此估计其身高为
看答案与解析
答案:C
带截距的最小二乘回归直线经过样本中心,用均值求截距即可。
求样本中心
xˉ=22.5 cm,yˉ=160 cm。
利用回归直线的均值性质
yˉ=b^xˉ+a^,所以 a^=160−4×22.5=70。回归方程为 y^=4x+70。
代入待预测脚长
当 x=24 cm,y^=4×24+70=166 cm,选C。
已知斜率和样本均值时,可用样本中心迅速恢复回归直线。