内容正文:
第九章 线性回归分析与方差分析
第一节 一元线性回归分析
第二节 可线性化的非线性回归
第三节 多元线性回归简介
第四节 方差分析
第二节 可线性化的非线性回归
在实际问题中,常常会遇到这样的情形:散点图上的几个样本数据点明显地不在一条直线附近,而在某曲线周围:
或者,用线性回归方程描述变量间的关系计算的结果与样本值误差较大,这表明变量之间不存在线性相关关系,而是一种非线性的相关关系.
下面举例说明对这类问题用线性化处理的方法.
例1 在彩色显像技术中,考虑析出银的光学密度x与形成染料光学密度Y之间的相关关系,其中11个样本数据如下所示:
xi 0.05 0.06 0.07 0.10 0.14 0.20 0.25 0.31 0.38 0.43 0.47
yi 0.10 0.14 0.23 0.37 0.59 0.79 1.00 1.12 1.19 1.25 1.29
解 根据这11个样本数据点
(xi,yi)作出散点图(图9-3).
图9-3
从散点图上看出,这些数据点在一条曲线L周围.
根据有关的专业知识,结合散点图,可以认为曲线L大致为:
对上式两边取对数:
令
即有:
0.25
0.22
0.17
0.11
0.00
-0.24
-0.53
-0.99
-1.47
-1.97
-2.30
2.13
2.33
2.63
3.23
4.00
5.00
7.14
10.00
14.29
16.67
20.00
于是数据( )相应地变换成( )
将变换后的数据点( )画出散点图(图9-4)
从散点图可以看出 与 具有线性相关关系,因此可用一元线性回归分析.
利用一元线性回归的方法可以计算出 与 的经验回归方程为
图9-4
可求得x与y之间相关关系的一个经验公式:
这里a=0.54,b= -0.15
所以
3′32.1″
2′13.9″
1′42.4″
43.86″
19.72″
9.95″
时间t(分′秒″)
1500
1000
800
400
200
100
距离x(m)
例2 赛跑是大家熟知的一种体育活动.下表给出了截至1997年底在6个不同的距离上中短跑成绩的世界记录:
试根据这些记录数据分析出运动员的赛跑成绩与所跑距离间的相关关系.
解 根据记录数据点(xi,ti)作出散点图 (图9-5)
图9-5
从散点图上看出,全部点(xi,ti)分布在一条曲线附近,因而x与t之间可以存在一种线性关系.
我们用一元线性回归分析,可计算出x与t间的线性回归模型为
t=-9.99+0.1455x
由此模型,当x=100,200,400,800,1000,1500(m)时,
t的理论值分别为:
4.56″, 19.10″,48.20″,1'46.4″,2'15.5″,3'28.2″
可以看出t的理论值与实际记录数据多数都比较接近。
仔细分析,可发现线性回归模型的一些不合理之处。
如:当赛跑距离小于68m时,所需时间为负值;
当赛跑距离为100m时所需时间只须4.56″.
再仔细分析,发现:短距离100m、200m及长距离1500米需要的时间实际值均高于线性模型的理论值,而中间的400m、800m、1000m需要的时间实际值均低于线性模型的理论值.
它告诉我们x与t的关系可能为一曲线,且曲线是下凸的。具有这种性质的最简单曲线当属幂函数:t=axb
它告诉我们x与t的关系可能为一曲线,且曲线是下凸的。
对上式二边取对数
lnt=lna+blnx
令
t′=lnt a′=lna x′=lnx
得
t′= a′+bx
为一线性关系
具有这种性质的最简单曲线当属幂函数:
t=axb
用一元线性回归分析估计a′,b,从而算出
最后可得t与x间的幂函数模型:
t=0.048x1.145
当x=100,200,400,800,1000,1500(m)时,
利用幂函数模型算出t的理论值分别为:
9.39″,20.78″,45.97″,1'41.68″, 2'11.29′,3'28.88″
比较计算结果可知:幂函数模型比线性回归模型更能确切地反映t与x间的关系.
第三节 多元线性回归简介
其中b0,b1,…,bp, 为与x1,…,xp无关的未知参数.
假定要考察p个自变量x1,x2,…,xp与因变量Y之间的相关关系.
设
这就是p元线性回归模型.
对变量x1,…,xp,Y作n次观测得到样本值:
( ) i=1,…,n
这里y1,…,yn独立、同分布,且