第三章 回归分析
- 格式:doc
- 大小:548.50 KB
- 文档页数:16
22 第三章 回归分析
§1 一元线性回归
一、回归模型
设随机变量y与自变量x之间存在线性关系,它们的第i次观测数据是:(xi,yi)(i=1,2,…,n)那么这组数据可以假设具有如下的数学结构式:iiixy0(i=1,…,n),其中β0, β为待估参数,),0(~2Ni,且n,,,21相互独立,这就是一元线性回归的数学模型。
二、参数估计
1.回归系数
设b0和b分别是参数β0, β的最小二乘估计,于是一元线性回归方程为:
iibxby0ˆ (i=1,2,…,n)
b0,b叫做回归系数,它使偏差平方和niiiniiibxbyyyQ12012)()ˆ(取最小值。
由
0)(20)(210100niiiiniiixbxbybQbxbybQ
整理得正规方程组: 020()()()iiiiiinbxbyxbxbxy
解得 xxxySSbxbyb/,0
其中 222)(xnxxxSiixx
yxnyxyyxxSiiiixy))((
另外 ynyyySiiyy22)(
2.最小二乘估计b0,b的统计性质
(1)E(b)= β,E(b0)= β0
即b0,b分别是β0,β的无偏估计
(2)22()/()iDbxx
22201()[/()]iDbxxxn
即回归系数b0,b与σ2,x的波动大小有关,b0还与n有关,这就是说,x值越分散,数据越多,估计b0,b越精确。
三、假设检验 23 1.回归方程显著性检验
欲检验y与x之间是否有线性关系,即检验假设H0:β=0。
设(xi,yi)(i=1,2,…,n),iibxbyL0ˆ:,即)(ˆxxbyy
其中
xxxySSb,xbyb0
记 2)(yySi总
2)ˆ(iiyyS剩
2)ˆ(yySi回
则 剩回总SSS
且 yySS总,xybSS回 xyyybSSS剩
方差分析表
来源 平方和 自由度 均方 F F0
回归 S回 1 S回 2/1/nSS剩回 )2,1(nF
剩余 S剩 n-2 S剩/n-2
总和 S总 n-1
若 )2,1(nFF,则H0真
若 )2,1(nFF,则H0不真
四、预报和控制
对任一给定的x0由回归方程可得回归值0000ˆ,ˆybxby是x0处的观察值0000xy的一个估计值。
1.预报
在α下寻找正数δ,使得y0以1-α的概率落在]ˆ,ˆ[00yy内,即
)ˆˆ(000yyyP
由概率统计知:当N较大,2000ˆ~(0,),xxyyN时,且有
24 %99)ˆ3ˆˆ3ˆ(%95)ˆ2ˆˆ2ˆ(%68)ˆˆˆˆ(000000000yyyPyyyPyyyP
2.控制
欲使21ˆyyy,如何控制x?
2211)(ˆ,)(ˆyxyyxy
求出x1,x2的问题。
五、例题
为了获得小麦高产,在5块田上进行了对比试验,在同样肥力及管理水平下,取得了如下数据:
试验号(i) 播种量
(11月18日,斤/亩) 基本苗数xi
(12月19日,万/亩) 有效穗数y
(5月5日,万/亩)
1
2
3
4
5 25
30
35
40
45 15
25.8
30
36.6
44.4 39.4
42.9
41.0
43.1
49.2
试对基本苗数x与有效穗数y进行回归分析。
解
1.求回归直线方程
编号 x y x2 y2 xy
1
2
3
4
5 15.0
25.8
30.0
36.6
44.4 39.4
42.9
41.0
43.1
49.2 225.00
665.64
900.00
1339.56
1971.36 1552.36
1840.41
1681.00
1857.61
2420.64
591.00
1106.82
1230.00
1577.46
2184.48
∑ 151.8 215.8 5101.56 9352.02 6689.76
35.55)(115.144))((192.492)(12222ynySyxnxySxnxSyyxyxx
32.34,29.00xbybSSbxxxy
所求的回归方程为 xy29.032.34ˆ
2.回归方程显著性检验 25 55.13)ˆ(80.41)ˆ(3.55)(1)(2222212xyyyiixyxxiyyiibSSyySbSSbyySSynyyyS剩回总
方差分析表
来源 平方和 自由度 均方和 F F0
回归 41.80 1 41.80 9.25* F0.1(1,3)=5.53
剩余 13.55 3 4.52
总和 55.35 4
F检验结果表明,回归方程比较显著,其可信程度为90%,在一般情况下是有效的。
在例中通过取样的方法测得一块麦田的基本苗数x0=26万/亩,则可预报第二年成熟时的有效穗数y0近似的为80.41ˆ0y(万/亩)
§2 多元线性回归
一、回归模型
设随机变量y与m个自变量X1,X2,…,Xm之间存在线性关系,它们的第i次观测数据为:(Xi1,Xi2,…,Xim;yi)(i=1,2,…,n),那么,这一组数据可以假设具有如下的数学结构式:
),,2,1(22110niXXXyiimmiii,
其中),0((~2Ni,且n1相互独立,),,,(10m为待估参数,这就是多元线性回归模型,即1-多线性回归模型。
二、参数估计
设b0,b1,…,bm分别是参数m,,,10的最小二乘估计,则回归方程为:
),,2,1(ˆ110niXbXbbyimmii
b0,b1,…,bm称为偏回归系数,它使偏差平方和(残差或剩余)
2111012)()ˆ(niimmiiniiiXbXbbyyyQ
取最小值。
应用求极值的方法,由
),,2,1,0(,0mibQi
经整理化简后得正规方程组: 26 mymmmmmymmymmSbSbSbSSbSbSbSSbSbSbS22112222212111212111
其中 najiajainajajiaiijXXnXXXXXXS11))((
najaajnaajajjyyXnyXyyXXS11))((
naanaaiiynyXnX111,1
由方程组解出 ),2,1(mibi
再计算 mmXbXbyb110
记 1(),(,,)xxijmmxyymySSSSS
则正规方程组的矩阵形式为:
xyxxSbS*,即 xyxxSSb1*
其中 ),,,(21*mbbbb
三、回归方程显著性检验
欲检验y与x1,…,xm之间是否有线性关系,即检验假设
0,,0:10mH
其方法是对平方和进行方差分析,将y的总离平方和S总分解为
S总=S回+S剩
其中 naayyyySS12)(总
naayyUS12)ˆ(回
naaayyQS1)ˆ(剩
且 f总=n-1,f回=m, f剩=n-m-1
27 方差分析表
方差来源 平方和 自由度 均方 F 临界F0
回归 S回 m
S回/m F=[S回/m]/
[S剩/n-m-1]
剩余 S剩 n-m-1 S剩/n-m-1
总和 S总 n-1
若 )1,(mnmFFa,则拒绝H。
四、回归系数显著性检验
欲检验Xi对y的影响是否显著,即检验假设H0:i=0,(i=1,2,…,m)
采用F检验
)1,1(~)1/(/)(2mnFmnSCbFiiiii剩
在H0成立时:
)1/(/2mnSCbFiiii剩,
其中Cii为S-1xx=C=(Cij)主对角线上的元素。
若)1,1(mnFFai,则拒绝H0,即Xi对y有影响。
五、例题
为了估算猪的毛重,测得14头猪的体长X1(cm),胸围X2(cm),与体重y(kg)的数据如表,试建立y与X1,X2的回归方程,并进行显著性检验。
序号 X1 X2 y X21 X22 X1X2 X1y X2y y2
1 41 49 28 1681 2401 2009 1148 1372 784
2 45 58 39 2025 3364 2610 1755 2262 1521
3 51 62 41 2601 3844 3162 2091 2542 1681
4 52 71 44 2704 5041 3692 2288 3124 1936
5 59 62 43 3481 3844 3658 2537 2666 1849
6 62 74 50 3844 5476 4588 3100 3700 2500
7 69 71 51 4761 5041 4899 3519 3621 2601
8 72 74 57 5184 5476 5328 4104 4218 3249
9 78 79 63 6084 6241 6162 4914 4977 3969
10 80 81 66 6400 6561 6480 5280 5346 4356
11 90 85 70 8100 7225 7650 6300 5950 4900
12 92 94 76 8464 8836 8648 6992 7144 5776
13 98 91 80 9604 8281 8918 7840 7280 6400
14 103 95 84 10609 9025 9785 8652 7980 7056
992 1046 792 65993 80656 77589 60520 62182 48578
解 1270.8674.71,56.57XXy