2009-7-30 1
回归分析实验目的实验内容
2、掌握用数学软件求解回归分析问题。
1、直观了解回归分析基本内容。
1,回归分析的基本理论 。
3,实验作业。
2,用数学软件求解回归分析问题。
2009-7-30 3
一元线性回归 多元线性回归回归分析数学模型及定义
*
模型参数估计
*
检验
、
预测与控制可线性化的一元非线性回归
（
曲线回归
）
数学模型及定义
*
模型参数估计
*
多元线性回归中的检验与预测逐步回归分析
2009-7-30 4
一、数学模型例 1 测 16名成年女子的身高与腿长所得数据如下：
身高 143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164
腿长 88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102
以身高 x为横坐标，以腿长 y为纵坐标将这些数据点（ xI,yi）
在平面直角坐标系上标出,
140 145 150 155 160 165
84
86
88
90
92
94
96
98
100
102
散点图
xy 10
解答
2009-7-30 5
一般地，称由 xy
10
确定的模型为 一元线性回归模型,
记为
2
10
,0
DE
xy
固定的未知参数 0?,1? 称为回归系数，自变量 x 也称为回归变量,
一元线性回归分析的 主要任务 是：
1,用试验值 （样本值）对
0
,
1
和? 作点估计；
2,对回归系数
0
,
1
作假设检验；
3,在 x= 0x 处对 y 作预测，对 y 作区间估计,
xY 10,称为 y 对 x 的回归直线方程,
返回
2009-7-30 6
二、模型参数估计
1、回归系数的最小二乘估计有 n 组独立观测值,（ x
1
,y
1
）,（ x
2
,y
2
），?,（ x
n
,y
n
）
设
相互独立且,
nii
ii
DE
nixy
...,,0
,...,2,1,
21
2
10
记
n
i
ii
n
i
i
xyQQ
1
2
10
1
2
10
),(
最小二乘法 就是选择
0
和 1
的估计
0
，
1
使得
),(m i n)
,
(
10
,
10
10
QQ?
2009-7-30 7
221
10
xx
yxxy
xy
解得其中
n
i
i
n
i
i ynyxnx
11
1,1,
n
i
ii
n
i
i yxnxyxnx
11
22 1,1,
（经验）回归方程为,)( 110 xxyxy
或
n
i
i
n
i
ii
xx
yyxx
1
2
1
1
2009-7-30 8
2,2? 的无偏估计记
n
i
n
i
iiiie
yyxyQQ
1 1
2
2
1010
)?(
)
,
(
称 Q
e
为 残差平方和 或 剩余平方和,
2
的无偏估计 为 )2(?
2
nQ
ee
称
2
e
为 剩余方差 （残差的方差）,
2
e
分别与 0
、
1
独立 。
e
称为 剩余标准差,
返回
2009-7-30 9
三、检验、预测与控制
1、回归方程的显著性检验对回归方程 xY
10
的显著性检验，归结为对假设
0:;0:
1110
HH
进行检验,
假设 0,10H 被拒绝，则回归显著，认为 y 与 x 存在线性关系，所求的线性回归方程有意义；否则回归不显著,y 与 x 的关系不能用一元线性回归模型来描述，所得的回归方程也无意义,
2009-7-30 10
（ Ⅰ ） F检验法当
0
H 成立时,
)2/(?
nQ
U
F
e
~F （ 1,n - 2 ）
其中
n
i
i
yyU
1
2
（ 回归平方和）
故 F> )2,1(1 nF?,拒绝 0H,否则就接受 0H,
（ Ⅱ ） t检验法
n
i
i
n
i
ixx xnxxxL
1
22
1
2)(其中当 0H 成立时，
e
xxL
T
1
~t （ n - 2 ）
故 )2(
2
1
ntT?,拒绝 0H,否则就接受 0H,
2009-7-30 11
（ Ⅲ ） r检验法当 | r | > r 1- α 时，拒绝 H 0 ；否则就接受 H 0,
记
n
i
n
i
ii
n
i
ii
yyxx
yyxx
r
1 1
22
1
)()(
))((
其中2,121 1
1
1
nFnr
2009-7-30 12
2、回归系数的置信区间
0
和
1
置信水平为 1- α 的置信区间分别为
xx
e
xx
e
L
x
n
nt
L
x
n
nt
2
2
1
0
2
2
1
0
1
)2(
,
1
)2(
和?
xxexxe
LntLnt /?)2(
,/?)2(
2
1
1
2
1
1
2
的置信水平为 1-? 的置信区间为
)2(
,
)2(
2
2
2
2
1
n
Q
n
Q
ee
2009-7-30 13
3、预测与控制
（ 1）预测用 y 0 的回归值 0100 xy 作为 y 0 的预测值,
0y 的置信水平为1 的 预测区间 为
)(?),(? 0000 xyxy
其中
xx
e L
xx
nntx
2
0
21
0
11)2(?)(
特别，当 n 很大且 x 0 在 x 附近取值时,
y 的置信水平为1 的 预测区间近似为
2
1
2
1
, uyuy ee
2009-7-30 14
（ 2）控制要求, xy 10 的值以1 的概率落在指定区间yy,
只要控制 x 满足以下两个不等式
yxyyxy )(?,)(
要求 )(2 xyy,若 yxyyxy )(?,)( 分别有解 x?
和 x,即 yxyyxy )(?,)(,
则xx,就是所求的 x 的控制区间,
返回
2009-7-30 15
四、可线性化的一元非线性回归
（曲线回归）
例 2 出钢时所用的盛钢水的钢包，由于钢水对耐火材料的侵蚀，
容积不断增大,我们希望知道使用次数与增大的容积之间的关系,对一钢包作试验，测得的数据列于下表：
使用次数 增大容积 使用次数 增大容积
2
3
4
5
6
7
8
9
6,4 2
8,2 0
9,5 8
9,5 0
9,7 0
1 0,0 0
9,9 3
9,9 9
10
11
12
13
14
15
16
1 0,4 9
1 0,5 9
1 0,6 0
1 0,8 0
1 0,6 0
1 0,9 0
1 0,7 6
解答
2009-7-30 16
2 4 6 8 10 12 14 16
6
6,5
7
7,5
8
8,5
9
9,5
10
1 0,5
11
散点图此即 非线性回归 或 曲线回归 问题 （ 需要配曲线 ）
配曲线的一般方法是：
先对两个变量 x 和 y 作 n 次试验观察得 niyx ii,...,2,1),,(? 画出散点图，
根据散点图确定须配曲线的类型,然后由 n 对试验数据确定每一类曲线的未知参数 a 和 b,采用的方法是通过变量代换把非线性回归化成线性回归，即采用非线性回归线性化的方法,
2009-7-30 17
通常选择的六类曲线如下：
（ 1 ） 双曲线 xbay1
（ 2 ） 幂函数曲线 y = a bx,其中 x > 0,a > 0
（ 3 ） 指数曲线 y = a bxe 其中参数 a > 0,
（ 4 ） 倒指数曲线 y= a xbe / 其中 a > 0,
（ 5 ） 对数曲线 y = a + b l o g x,x > 0
（ 6 ） S 型曲线 xbeay 1
返回解例 2,由散点图我们选配倒指数曲线 y =a
xb
e
/
根据线性化方法，算得 4 5 8 7.2?,1 1 0 7.1 Ab
由此 6 7 8 9.11?
A
ea
最后得
xey
1 1 0 7.1
6 7 8 9.11
2009-7-30 18
一、数学模型及定义一般称
n
IC O VE
XY
2
),(,0)(
为高斯 — 马尔柯夫线性模型 ( k 元线性回归模型 ),并简记为 ),,(
2
n
IXY
ny
y
Y
...
...
1
，
nknn
k
k
xxx
xxx
xxx
X
...1
...............
...1
...1
21
22221
11211
，
k?
.,,
1
0
，
n?
.,,
2
1
kk xxy,..110 称为 回归平面方程,
返回线性模型 ),,(
2
nIXY 考虑的主要问题是,
（ 1 ） 用试验值（样本值）对未知参数? 和
2?
作点估计和假设检验，从而建立 y 与
k
xxx,...,,
21
之间的数量关系；
（ 2 ）在,,.,,,,0022011 kk xxxxxx 处对 y 的值作预测与控制，即对 y 作区间估计,
2009-7-30 19
二、模型参数估计
1,对 i? 和 2? 作估计用最小二乘法求
k
,.,,,
0
的估计量：作离差平方和
n
i
ikkii
xxyQ
1
2
110
.,,
选择 k,.,,,0 使 Q 达到最小。
解得估计值YXXX TT 1
得到的 i 代入回归平面方程得：
kk xxy.,, 110
称为 经验回归平面方程,i 称为 经验回归系数,
注意, 服从 p+ 1 维正态分布，且为? 的无偏估计，协方差阵为 C
2
,
C = L -1 =( c ij ),L = X ’ X
2009-7-30 20
2,多项式回归设变量 x,Y 的回归模型为
p
p
xxxY,,,
2
210
其中 p 是已知的,),,2,1( pi
i
是未知参数,? 服从正态分布
),0(
2
N,
令 ii xx?,i =1,2,…,k 多项式回归模型变为多元线性回归模型,
返回
kk xxxY,,,2210
称为 回归多项式,上面的回归模型称为 多项式回归,
2009-7-30 21
三、多元线性回归中的检验与预测
1,线性模型和回归系数的检验假设 0...,100 kH
（ Ⅰ ） F检验法
（ Ⅱ ） r检验法定义
eyy QU
U
L
U
R
为 y 与 x 1,x 2,.,,,x k 的 多元相关系数 或 复相关系数 。
由于 2
2
1
1
R
R
k
kn
F
,故用 F 和用 R 检验是等效的。
当 H 0 成立时,)1,(~
)1/(
/
knkF
knQ
kU
F
e
如果 F > F 1- α （ k,n - k - 1 ），则拒绝 H 0,认为 y 与 x 1,…,x k 之间显著地有线性关系；否则就接受 H 0,认为 y 与 x 1,…,x k 之间线性关系不显著,
其中
n
i
i yyU
1
2? （ 回归平方和）?
n
i
iie yyQ
1
2)?((残差平方和）
2009-7-30 22
2、预测
（ 1）点预测求出回归方程 kk xxy... 110,对于给定自变量的值?kxx,.,,,*1,用 **110*?... kk xxy 来预测
**110,,,kk xxy,称 *?y 为 *y 的点预测,
（ 2）区间预测
y 的1 的预测区间 （置信）区间为 )?,?(
21
yy,其中
)1(1
)1(1
2/1
0 0
2
2/1
0 0
1
kntxxcyy
kntxxcyy
k
i
k
j
jiije
k
i
k
j
jiije
C = L
-1
=( c
ij
),L = X ’ X
1 kn
Q e
e?
返回
2009-7-30 23
四、逐步回归分析
（ 4）“有进有出”的逐步回归分析。
（ 1）从所有可能的因子（变量）组合的回归方程中选择最优者；
（ 2）从包含全部变量的回归方程中逐次剔除不显著因子；
（ 3）从一个变量开始，把变量逐个引入方程；
选择“最优”的回归方程有以下几种方法：
“最优,的回归方程 就是包含所有对 Y有影响的变量,而不包含对 Y影响不显著的变量回归方程 。
以第四种方法，即 逐步回归分析法 在筛选变量方面较为理想,
2009-7-30 24
这个过程反复进行，直至既无不显著的变量从回归方程中剔除，又无显著变量可引入回归方程时为止。
逐步回归分析法 的思想：
从一个自变量开始，视自变量 Y作用的显著程度，从大到地依次逐个引入回归方程。
当引入的自变量由于后面变量的引入而变得不显著时，
要将其剔除掉。
引入一个自变量或从回归方程中剔除一个自变量，为逐步回归的一步。
对于每一步都要进行 Y值检验，以确保每次引入新的显著性变量前回归方程中只包含对 Y作用显著的变量。
返回
2009-7-30 25
1、考察温度 x对产量 y的影响，测得下列 10组数据：
温度 （ ℃ ） 20 25 30 35 40 45 50 55 60 65
产量 （ kg ） 1 3,2 1 5,1 1 6,4 1 7,1 1 7,9 1 8,7 1 9,6 2 1,2 2 2,5 2 4,3
求 y关于 x的线性回归方程，检验回归效果是否显著，并预测
x=42℃ 时产量的估值及预测区间（置信度 95%）,
2、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标 xi处测得纵坐标 yi共 11对数据如下：
x i 0 2 4 6 8 10 12 14 16 18 20
y i 0,6 2,0 4,4 7,5 1 1,8 1 7,1 2 3,3 3 1,2 3 9,6 4 9,7 6 1,7
求这段曲线的纵坐标 y关于横坐标 x的二次多项式回归方程,
2009-7-30 26
3,在研究化学动力学反应过程中，建立了一个反应速度和反应物含量的数学模型，形式为
342312
5
3
21
1 xxx
x
x
y
其中
51
,, 是未知参数，
321
,,xxx 是三种反应物 （氢,n 戊烷，
异构戊烷）的含量,y 是反应速度,今测得一组数据如表 4,试由此确定参数
51
,,,并给出置信区间,
51
,, 的参考值为
（ 1,0,05,0,02,0,1,2 ）,
序号 反应速度 y 氢 x
1
n 戊烷 x
2
异构戊烷 x
3
1 8,55 470 300 10
2 3,79 285 80 10
3 4,82 470 300 120
4 0,02 470 80 120
5 2,75 470 80 10
6 14,39 100 190 10
7 2,54 100 80 65
8 4,35 470 190 65
9 13,00 100 300 54
10 8,50 100 300 120
11 0,05 100 80 120
12 1 1,32 285 300 10
13 3,13 285 190 120
2009-7-30 27
4、混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成 12个试块，记录了养护日期 x（日）及抗压强度 y
（ kg/cm2）的数据：
养护时间 x 2 3 4 5 7 9 12 14 17 21 28 56
抗压强度 y 35 42 47 53 59 65 68 73 76 82 86 99
试求 xbay ln 型回归方程,
回归分析实验目的实验内容
2、掌握用数学软件求解回归分析问题。
1、直观了解回归分析基本内容。
1,回归分析的基本理论 。
3,实验作业。
2,用数学软件求解回归分析问题。
2009-7-30 3
一元线性回归 多元线性回归回归分析数学模型及定义
*
模型参数估计
*
检验
、
预测与控制可线性化的一元非线性回归
（
曲线回归
）
数学模型及定义
*
模型参数估计
*
多元线性回归中的检验与预测逐步回归分析
2009-7-30 4
一、数学模型例 1 测 16名成年女子的身高与腿长所得数据如下：
身高 143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164
腿长 88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102
以身高 x为横坐标，以腿长 y为纵坐标将这些数据点（ xI,yi）
在平面直角坐标系上标出,
140 145 150 155 160 165
84
86
88
90
92
94
96
98
100
102
散点图
xy 10
解答
2009-7-30 5
一般地，称由 xy
10
确定的模型为 一元线性回归模型,
记为
2
10
,0
DE
xy
固定的未知参数 0?,1? 称为回归系数，自变量 x 也称为回归变量,
一元线性回归分析的 主要任务 是：
1,用试验值 （样本值）对
0
,
1
和? 作点估计；
2,对回归系数
0
,
1
作假设检验；
3,在 x= 0x 处对 y 作预测，对 y 作区间估计,
xY 10,称为 y 对 x 的回归直线方程,
返回
2009-7-30 6
二、模型参数估计
1、回归系数的最小二乘估计有 n 组独立观测值,（ x
1
,y
1
）,（ x
2
,y
2
），?,（ x
n
,y
n
）
设
相互独立且,
nii
ii
DE
nixy
...,,0
,...,2,1,
21
2
10
记
n
i
ii
n
i
i
xyQQ
1
2
10
1
2
10
),(
最小二乘法 就是选择
0
和 1
的估计
0
，
1
使得
),(m i n)
,
(
10
,
10
10
QQ?
2009-7-30 7
221
10
xx
yxxy
xy
解得其中
n
i
i
n
i
i ynyxnx
11
1,1,
n
i
ii
n
i
i yxnxyxnx
11
22 1,1,
（经验）回归方程为,)( 110 xxyxy
或
n
i
i
n
i
ii
xx
yyxx
1
2
1
1
2009-7-30 8
2,2? 的无偏估计记
n
i
n
i
iiiie
yyxyQQ
1 1
2
2
1010
)?(
)
,
(
称 Q
e
为 残差平方和 或 剩余平方和,
2
的无偏估计 为 )2(?
2
nQ
ee
称
2
e
为 剩余方差 （残差的方差）,
2
e
分别与 0
、
1
独立 。
e
称为 剩余标准差,
返回
2009-7-30 9
三、检验、预测与控制
1、回归方程的显著性检验对回归方程 xY
10
的显著性检验，归结为对假设
0:;0:
1110
HH
进行检验,
假设 0,10H 被拒绝，则回归显著，认为 y 与 x 存在线性关系，所求的线性回归方程有意义；否则回归不显著,y 与 x 的关系不能用一元线性回归模型来描述，所得的回归方程也无意义,
2009-7-30 10
（ Ⅰ ） F检验法当
0
H 成立时,
)2/(?
nQ
U
F
e
~F （ 1,n - 2 ）
其中
n
i
i
yyU
1
2
（ 回归平方和）
故 F> )2,1(1 nF?,拒绝 0H,否则就接受 0H,
（ Ⅱ ） t检验法
n
i
i
n
i
ixx xnxxxL
1
22
1
2)(其中当 0H 成立时，
e
xxL
T
1
~t （ n - 2 ）
故 )2(
2
1
ntT?,拒绝 0H,否则就接受 0H,
2009-7-30 11
（ Ⅲ ） r检验法当 | r | > r 1- α 时，拒绝 H 0 ；否则就接受 H 0,
记
n
i
n
i
ii
n
i
ii
yyxx
yyxx
r
1 1
22
1
)()(
))((
其中2,121 1
1
1
nFnr
2009-7-30 12
2、回归系数的置信区间
0
和
1
置信水平为 1- α 的置信区间分别为
xx
e
xx
e
L
x
n
nt
L
x
n
nt
2
2
1
0
2
2
1
0
1
)2(
,
1
)2(
和?
xxexxe
LntLnt /?)2(
,/?)2(
2
1
1
2
1
1
2
的置信水平为 1-? 的置信区间为
)2(
,
)2(
2
2
2
2
1
n
Q
n
Q
ee
2009-7-30 13
3、预测与控制
（ 1）预测用 y 0 的回归值 0100 xy 作为 y 0 的预测值,
0y 的置信水平为1 的 预测区间 为
)(?),(? 0000 xyxy
其中
xx
e L
xx
nntx
2
0
21
0
11)2(?)(
特别，当 n 很大且 x 0 在 x 附近取值时,
y 的置信水平为1 的 预测区间近似为
2
1
2
1
, uyuy ee
2009-7-30 14
（ 2）控制要求, xy 10 的值以1 的概率落在指定区间yy,
只要控制 x 满足以下两个不等式
yxyyxy )(?,)(
要求 )(2 xyy,若 yxyyxy )(?,)( 分别有解 x?
和 x,即 yxyyxy )(?,)(,
则xx,就是所求的 x 的控制区间,
返回
2009-7-30 15
四、可线性化的一元非线性回归
（曲线回归）
例 2 出钢时所用的盛钢水的钢包，由于钢水对耐火材料的侵蚀，
容积不断增大,我们希望知道使用次数与增大的容积之间的关系,对一钢包作试验，测得的数据列于下表：
使用次数 增大容积 使用次数 增大容积
2
3
4
5
6
7
8
9
6,4 2
8,2 0
9,5 8
9,5 0
9,7 0
1 0,0 0
9,9 3
9,9 9
10
11
12
13
14
15
16
1 0,4 9
1 0,5 9
1 0,6 0
1 0,8 0
1 0,6 0
1 0,9 0
1 0,7 6
解答
2009-7-30 16
2 4 6 8 10 12 14 16
6
6,5
7
7,5
8
8,5
9
9,5
10
1 0,5
11
散点图此即 非线性回归 或 曲线回归 问题 （ 需要配曲线 ）
配曲线的一般方法是：
先对两个变量 x 和 y 作 n 次试验观察得 niyx ii,...,2,1),,(? 画出散点图，
根据散点图确定须配曲线的类型,然后由 n 对试验数据确定每一类曲线的未知参数 a 和 b,采用的方法是通过变量代换把非线性回归化成线性回归，即采用非线性回归线性化的方法,
2009-7-30 17
通常选择的六类曲线如下：
（ 1 ） 双曲线 xbay1
（ 2 ） 幂函数曲线 y = a bx,其中 x > 0,a > 0
（ 3 ） 指数曲线 y = a bxe 其中参数 a > 0,
（ 4 ） 倒指数曲线 y= a xbe / 其中 a > 0,
（ 5 ） 对数曲线 y = a + b l o g x,x > 0
（ 6 ） S 型曲线 xbeay 1
返回解例 2,由散点图我们选配倒指数曲线 y =a
xb
e
/
根据线性化方法，算得 4 5 8 7.2?,1 1 0 7.1 Ab
由此 6 7 8 9.11?
A
ea
最后得
xey
1 1 0 7.1
6 7 8 9.11
2009-7-30 18
一、数学模型及定义一般称
n
IC O VE
XY
2
),(,0)(
为高斯 — 马尔柯夫线性模型 ( k 元线性回归模型 ),并简记为 ),,(
2
n
IXY
ny
y
Y
...
...
1
，
nknn
k
k
xxx
xxx
xxx
X
...1
...............
...1
...1
21
22221
11211
，
k?
.,,
1
0
，
n?
.,,
2
1
kk xxy,..110 称为 回归平面方程,
返回线性模型 ),,(
2
nIXY 考虑的主要问题是,
（ 1 ） 用试验值（样本值）对未知参数? 和
2?
作点估计和假设检验，从而建立 y 与
k
xxx,...,,
21
之间的数量关系；
（ 2 ）在,,.,,,,0022011 kk xxxxxx 处对 y 的值作预测与控制，即对 y 作区间估计,
2009-7-30 19
二、模型参数估计
1,对 i? 和 2? 作估计用最小二乘法求
k
,.,,,
0
的估计量：作离差平方和
n
i
ikkii
xxyQ
1
2
110
.,,
选择 k,.,,,0 使 Q 达到最小。
解得估计值YXXX TT 1
得到的 i 代入回归平面方程得：
kk xxy.,, 110
称为 经验回归平面方程,i 称为 经验回归系数,
注意, 服从 p+ 1 维正态分布，且为? 的无偏估计，协方差阵为 C
2
,
C = L -1 =( c ij ),L = X ’ X
2009-7-30 20
2,多项式回归设变量 x,Y 的回归模型为
p
p
xxxY,,,
2
210
其中 p 是已知的,),,2,1( pi
i
是未知参数,? 服从正态分布
),0(
2
N,
令 ii xx?,i =1,2,…,k 多项式回归模型变为多元线性回归模型,
返回
kk xxxY,,,2210
称为 回归多项式,上面的回归模型称为 多项式回归,
2009-7-30 21
三、多元线性回归中的检验与预测
1,线性模型和回归系数的检验假设 0...,100 kH
（ Ⅰ ） F检验法
（ Ⅱ ） r检验法定义
eyy QU
U
L
U
R
为 y 与 x 1,x 2,.,,,x k 的 多元相关系数 或 复相关系数 。
由于 2
2
1
1
R
R
k
kn
F
,故用 F 和用 R 检验是等效的。
当 H 0 成立时,)1,(~
)1/(
/
knkF
knQ
kU
F
e
如果 F > F 1- α （ k,n - k - 1 ），则拒绝 H 0,认为 y 与 x 1,…,x k 之间显著地有线性关系；否则就接受 H 0,认为 y 与 x 1,…,x k 之间线性关系不显著,
其中
n
i
i yyU
1
2? （ 回归平方和）?
n
i
iie yyQ
1
2)?((残差平方和）
2009-7-30 22
2、预测
（ 1）点预测求出回归方程 kk xxy... 110,对于给定自变量的值?kxx,.,,,*1,用 **110*?... kk xxy 来预测
**110,,,kk xxy,称 *?y 为 *y 的点预测,
（ 2）区间预测
y 的1 的预测区间 （置信）区间为 )?,?(
21
yy,其中
)1(1
)1(1
2/1
0 0
2
2/1
0 0
1
kntxxcyy
kntxxcyy
k
i
k
j
jiije
k
i
k
j
jiije
C = L
-1
=( c
ij
),L = X ’ X
1 kn
Q e
e?
返回
2009-7-30 23
四、逐步回归分析
（ 4）“有进有出”的逐步回归分析。
（ 1）从所有可能的因子（变量）组合的回归方程中选择最优者；
（ 2）从包含全部变量的回归方程中逐次剔除不显著因子；
（ 3）从一个变量开始，把变量逐个引入方程；
选择“最优”的回归方程有以下几种方法：
“最优,的回归方程 就是包含所有对 Y有影响的变量,而不包含对 Y影响不显著的变量回归方程 。
以第四种方法，即 逐步回归分析法 在筛选变量方面较为理想,
2009-7-30 24
这个过程反复进行，直至既无不显著的变量从回归方程中剔除，又无显著变量可引入回归方程时为止。
逐步回归分析法 的思想：
从一个自变量开始，视自变量 Y作用的显著程度，从大到地依次逐个引入回归方程。
当引入的自变量由于后面变量的引入而变得不显著时，
要将其剔除掉。
引入一个自变量或从回归方程中剔除一个自变量，为逐步回归的一步。
对于每一步都要进行 Y值检验，以确保每次引入新的显著性变量前回归方程中只包含对 Y作用显著的变量。
返回
2009-7-30 25
1、考察温度 x对产量 y的影响，测得下列 10组数据：
温度 （ ℃ ） 20 25 30 35 40 45 50 55 60 65
产量 （ kg ） 1 3,2 1 5,1 1 6,4 1 7,1 1 7,9 1 8,7 1 9,6 2 1,2 2 2,5 2 4,3
求 y关于 x的线性回归方程，检验回归效果是否显著，并预测
x=42℃ 时产量的估值及预测区间（置信度 95%）,
2、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标 xi处测得纵坐标 yi共 11对数据如下：
x i 0 2 4 6 8 10 12 14 16 18 20
y i 0,6 2,0 4,4 7,5 1 1,8 1 7,1 2 3,3 3 1,2 3 9,6 4 9,7 6 1,7
求这段曲线的纵坐标 y关于横坐标 x的二次多项式回归方程,
2009-7-30 26
3,在研究化学动力学反应过程中，建立了一个反应速度和反应物含量的数学模型，形式为
342312
5
3
21
1 xxx
x
x
y
其中
51
,, 是未知参数，
321
,,xxx 是三种反应物 （氢,n 戊烷，
异构戊烷）的含量,y 是反应速度,今测得一组数据如表 4,试由此确定参数
51
,,,并给出置信区间,
51
,, 的参考值为
（ 1,0,05,0,02,0,1,2 ）,
序号 反应速度 y 氢 x
1
n 戊烷 x
2
异构戊烷 x
3
1 8,55 470 300 10
2 3,79 285 80 10
3 4,82 470 300 120
4 0,02 470 80 120
5 2,75 470 80 10
6 14,39 100 190 10
7 2,54 100 80 65
8 4,35 470 190 65
9 13,00 100 300 54
10 8,50 100 300 120
11 0,05 100 80 120
12 1 1,32 285 300 10
13 3,13 285 190 120
2009-7-30 27
4、混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成 12个试块，记录了养护日期 x（日）及抗压强度 y
（ kg/cm2）的数据：
养护时间 x 2 3 4 5 7 9 12 14 17 21 28 56
抗压强度 y 35 42 47 53 59 65 68 73 76 82 86 99
试求 xbay ln 型回归方程,
