《回归分析》

  • 格式:doc
  • 大小:1.17 MB
  • 文档页数:22

第六章 回归分析 引 言 一、变量间的关系 1.确定关系(函数关系) 如∶正方形的面积y与边长x的关系、做匀速直线运动的质点的位移S与时间t的关系等等。 2.非确定关系(相关关系) 例1.作物的单位面积产量y与施肥量x的关系。 例2.人的体重y与身高x的关系。 例3.人的血压y与年龄x的关系。 例4.某种商品(如糖果)的消费量y与居民(按 人口计算)的平均收入x的关系。 共同点∶变量间虽然存在着密切的关系,但从一个变量的每一确定值,不能求出另一变量的确定的值。可是在大量的试验中,这种不确定的联系,具有统计规律性,表现为 yx|的一个确定的概率分布

定义. 设有两个随机变量,对其中某一个变量的每一个可能的值,另一个变量有一个确定的概率分布与之对应。则称这两个随机变量间存在着相关关系。 二、回归分析 1.回归关系 设y与x之间存在着某种相关关系。其中x是可以控制或可以精确观察的变量(如年龄、试验的温度等等),即,可以随意指定x的n个值nxxx,,,21

。因此视x为普通变量。此时称变量y

对x具有回归关系。 [问题]如何表达变量y对x的回归关系? 回归分析是研究相关(回归)关系的一种数学工具,它能帮助我们从一个变量的取值去估计另外一个变量所取得的值.在各个学科领域得到广泛应用。

2.回归分析要解决的一些问题 设通过试验或抽样调查得到了x和y的一组观察值 ),(,),,(),,(2211nnyxyxyx

利用这些数据研究y与x之间的关系,需解决如下几个问题∶ (1)如何判断y与x间是否存在相关关系(或者回归关系)? ————利用样本数据做y与x的散点图,由散点图的趋势初步判断两者是否有线性关系 (2)如何确定y与x之间的定量关系式(即回归函数)? ——————用最小二乘法 (3)如何对所确定的关系(回归模型)的可信程度作统计检验?——t分布检验法、相关系数r检验法 (4)如何利用所得到y与x的定量关系式(回归模型)来解决实际问题?————对变量进行预测和控制。 下面,将围绕这四个问题展开讨论 §6.1 一元线性回归 设通过试验或抽样调查得到了x和y的一组观察值 ),(,),,(),,(2211nnyxyxyx

利用这些数据研究y与x之间的关系。 一、散点图 利用x和y的一组观察值 ),(,),,(),,(2211nnyxyxyx

在直角坐标系上描点,得到散点图,根据散点图的趋势,来观察x和y是否想直线关系。如果是,可以考虑建立线性回归模型。

例如,下列散点图中,上面两个象直线关系,可考虑对变量x和y建立线性回归模型。

[例1]为研究全国人均国民收入x(元)与人均消费金额y(元)之间的关系,收集到1981—1990年10年的样本数据10,,2,1),,(iyxii如下表∶

年 份 1981 1982 1983 1984 1985 人均国民收入 393 419 461 544 668 人均消费金额 249 267 289 329 406 年 份 1986 1987 1988 1989 1990 人均国民收入 738 860 1069 1169 1251 人均消费金额 451 513 643 699 713

20040060080010001200X200400600Y 0 回归分析解法后面将介绍:

(n10, x757.2, 2xnS909739.6, y455.9,

niiy14559, 2ynS288428.9, niiiyx13963622,

,5623.06.9097394592.7573963622)(2111xniiniiinSyxyxb 126.302.7575623.09.45510xbyb 经验线性回归方程为 xy5623.0126.30ˆ )

二、模型 xy

10 (1) 其中,0、1为常数,x为普通变量,y、为随机变量,满足),0(~2N,即),(~210xNy。 x—自变量,y—因变量,—随机波动项(误差)

于是样本),(,),,(),,(2211nnyxyxyx满足∶

iiixy10,(ni,,2,1) (2) 假设n,,,21相互独立,从而nyyy,,,21

满足独立、正态、

等方差的前提条件。 “确定y与x之间的定量关系式(即回归函数)”的问题归结为“如何在上述假定下确定0、1和”的问题。 [注1]当0、1和已知时,由),(~210xNy可得

01(){||}1yxPU

对可靠性%951,96.1U 于是,y的点估计为xy10ˆ



误差限为96.1,y的置信区 间为∶ ]96.1,96.1[1010xx

虚线为y的95%的置信区间控制线

[注2]两类误差 实际问题中0、1

和未知的,要通过观测数据),(,),,(),,(2211nnyxyxyx进行估计。这样用x来预测y时,就面临两类误差∶ (i)所代表的误差—随机误差,由y与x的关系的紧密程度决定,与数据),(iiyx(ni,2,1)无关。 (ii)利用样本),(iiyx(ni,2,1)对0、1和的估计所产生的误差,与数据及数据量的大小有关。

三、最小二乘估计 线性回归模型: xy10

设 00ˆb,11ˆb,则iy的估计值为

iixbby10ˆ,),,2,1(ni

残差∶ )(ˆ10iiiiixbbyyye,),,2,1(ni 残差平方和∶

niiiniiiniiexbbyyyeSS12101212

)()ˆ(

它描述了用线性函数01ˆybbx来近似表示y产生的误差程度。 最小二乘原则∶选取0b和1b使eSS达到最小。按照这种方法确定

的00ˆb和11ˆb叫做0和1的最小二乘估计,这种方法也叫做最小二乘法。  0

和1的最小二乘估计确定

关于0b和1b二元函数eSS的极值在二阶偏导书为0的点取得:

0)1)((21100niiiexbbyb

SS

0))((21101niiiiexxbbyb

SS 



niiiniiniiniiniiyxbxbxybxbn1112011110)()()(

—正规方程

方程两边同除以n得 



niiiniiyxnbxnbxybxb11120101)1( 其解为





niiniiniiniiniiiniiniiixnxyxnyxxxyyxxbxbyb121211112111100)(1)()(1)())((ˆ

ˆ



——0和1的最小二乘估计 称: xbby10ˆ

为y对x经验线性回归方程。

四、2的无偏估计 (1) 2)2(][nSSEe (2) (重要!)记 )2/(2nSSSexy, ----y对x的回归剩余方差 ,显然22.[]yxES,即2yxS为2的无偏估计。 )2/(nSSSexy -----y对x的回归剩余标准差

[注]:简化最小二乘估计的计算步骤 记 niixxxnS122)(,niiyyynS122)( 注意: 2211()nxiiSxxn,2211()nyiiSyyn,表示样本二阶中心矩。不是样本方差,要注意区分。

xbybnSyxyxbxniiniii10021111ˆ,)(ˆ

niiiniiiniiexbbyyyeSS12101212

)()ˆ(

2212xynSbnS

, )2/(ˆnSSSexy

a.将nxxx,,,21输入计算器,求出x 和 2xnS

b.将nyyy,,,21输入计算器,求出y、niiy1 和 2ynS

c.将nnyxyxyx,,,2211 输入计算器,求出niiiyx1

d.计算xbybnSyxyxbxniiniii102111,)( e.写出经验线性回归方程 xbby10ˆ

f.计算 2212xyenSbnSSS, )2/(ˆnSSSexy 注:实际上,利用具有回归功能的计算器,输入相关数据后,就可