计量经济学 第六章 自相关.

Slides:



Advertisements
Similar presentations
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
Advertisements

第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二节 换元积分法 一、第一类换元积分 法(凑微分法) 二、第二类换元积分法. 问题 解决方法 利用复合函数,设置中间变量. 过程令 一、第一类换元积分法(凑微分法)
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
计量经济学 第五章 异 方 差 性.
3.4 空间直线的方程.
第三章 异方差和自相关.
圆的一般方程 (x-a)2 +(y-b)2=r2 x2+y2+Dx+Ey+F=0 Ax2+Bxy+Cy2+Dx+Ey+ F=0.
第五章 二次型. 第五章 二次型 知识点1---二次型及其矩阵表示 二次型的基本概念 1. 线性变换与合同矩阵 2.
一、能线性化的多元非线性回归 二、多元多项式回归(线性化)
第三章 函数逼近 — 最佳平方逼近.
6.6 单侧置信限 1、问题的引入 2、基本概念 3、典型例题 4、小结.
《高等数学》(理学) 常数项级数的概念 袁安锋
第二章 经典单方程计量经济学模型: 一元线性回归模型
量化视角下的豆粕投资机会分析 格林期货研发培训中心 郭坤龙.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
第一章 商品 第一节 价值创造 第二节 价值量 第三节 价值函数及其性质 第四节 商品经济的基本矛盾与利己利他经济人假设.
初中数学八年级下册 (苏科版) 10.4 探索三角形 相似的条件(2).
第2章 Z变换 Z变换的定义与收敛域 Z反变换 系统的稳定性和H(z) 系统函数.
第三章 多维随机变量及其分布 §2 边缘分布 边缘分布函数 边缘分布律 边缘概率密度.
§3.3 多元线性回归模型的统计检验 一、拟合优度检验 二、方程的显著性检验(F检验) 三、变量的显著性检验(t检验) 四、参数的置信区间.
引子: 国内生产总值增加会减少财政收入吗?
第六章 正态条件下回归的推论.
第二章 回归模型 法、参数的普通最小二乘估计式及相关性质、对模型的经济意 义检验和统计检验,能应用Eviews软件进行最小二乘估计与统
一元线性回归模型 § 1 回归分析概述 § 2 一元线性回归模型的参数估计 § 3 一元线性回归模型的统计检验
计量软件实验:EViews软件应用 基本回归模型.
第2章 一元线性回归 2 .1 一元线性回归模型 2 .2 参数 的估计 2 .3 最小二乘估计的性质 2 .4 回归方程的显著性检验
计算机数学基础 主讲老师: 邓辉文.
第十章 方差分析.
宁波市高校慕课联盟课程 与 进行交互 Linux 系统管理.
第七章 参数估计 7.3 参数的区间估计.
若2002年我国国民生产总值为 亿元,如果 ,那么经过多少年国民生产总值 每年平均增长 是2002年时的2倍? 解:设经过 年国民生产总值为2002年时的2倍, 根据题意有 , 即.
习题 一、概率论 1.已知随机事件A,B,C满足 在下列三种情况下,计算 (1)A,B,C相互独立 (2)A,B独立,A,C互不相容
抽样和抽样分布 基本计算 Sampling & Sampling distribution
3.8.1 代数法计算终点误差 终点误差公式和终点误差图及其应用 3.8 酸碱滴定的终点误差
概 率 统 计 主讲教师 叶宏 山东大学数学院.
线 性 代 数 厦门大学线性代数教学组 2019年4月24日6时8分 / 45.
5.2 常用统计分布 一、常见分布 二、概率分布的分位数 三、小结.
第六章 自相关.
§4.2 序列相关性 Serial Correlation.
第五章 异方差.
定理21.9(可满足性定理)设A是P(Y)的协调子集,则存在P(Y)的解释域U和项解释,使得赋值函数v(A){1}。
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
相关与回归 非确定关系 在宏观上存在关系,但并未精确到可以用函数关系来表达。青少年身高与年龄,体重与体表面积 非确定关系:
第五节 缓冲溶液pH值的计算 两种物质的性质 浓度 pH值 共轭酸碱对间的质子传递平衡 可用通式表示如下: HB+H2O ⇌ H3O++B-
第4课时 绝对值.
第三章 两变量线性回归.
第四章 多元线性回归分析.
第一部分:概率 产生随机样本:对分布采样 均匀分布 其他分布 伪随机数 很多统计软件包中都有此工具 如在Matlab中:rand
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
§5.2 抽样分布   确定统计量的分布——抽样分布,是数理统计的基本问题之一.采用求随机向量的函数的分布的方法可得到抽样分布.由于样本容量一般不止2或 3(甚至还可能是随机的),故计算往往很复杂,有时还需要特殊技巧或特殊工具.   由于正态总体是最常见的总体,故本节介绍的几个抽样分布均对正态总体而言.
导 言 经济学的基本问题 经济学的基本研究方法 需求和供给.
概率论与数理统计B.
§2 方阵的特征值与特征向量.
第二节 简单线性回归模型的最小二乘估计 用样本去估计总体回归函数,总要使用特定的方法,而任何估 计参数的方法都需要有一定的前提条件——假定条件 一、简单线性回归的基本假定 为什么要作基本假定? ●只有具备一定的假定条件,所作出的估计才具有良好的统计性质。 ●模型中有随机扰动项,估计的参数是随机变量,显然参数估计值的分布与扰动项的分布有关,只有对随机扰动的分布作出假定,才能比较方便地确定所估计参数的分布性质,也才可能进行假设检验和区间估计等统计推断。
多元线性回归分析.
第三章 从概率分布函数的抽样 (Sampling from Probability Distribution Functions)
本底对汞原子第一激发能测量的影响 钱振宇
第三节 数量积 向量积 混合积 一、向量的数量积 二、向量的向量积 三、向量的混合积 四、小结 思考题.
数学模型实验课(二) 最小二乘法与直线拟合.
一元一次方程的解法(-).
Presentation transcript:

计量经济学 第六章 自相关

引子 T检验和F检验一定就可靠吗? 研究居民储蓄存款Y与居民收入X的关系: 用普通最小二乘法估计其参数,结果为 (1.8690) (0.0055) t= (14.9343) (64.2069) F=4122.531 检验结果:回归系数的标准误差非常小,t统计量较大,说明居 民收入X对居民储蓄存款Y的影响非常显著。同时可决系数也非 常高,F统计量=4122.531,也表明模型异常的显著。 但此估计结果可能是虚假的,t统计量和F统计量都被虚假 地夸大,因此所得结果是不可信的。为什么?

经济系统中的自相关 一、自相关的概念 二、自相关产生的背景与原因 三、自相关性的后果 四、自相关性的检验 五、自相关问题的处理方法

第一节 什么是自相关 一、自相关的概念 自相关(auto correlation),又称序列相关(serial correlation)是指总体回归模型的随机误差项之间存在相关关系。即不同观测点上的误差项彼此相关。

一阶自相关系数 自相关系数ρ的定义与普通相关系的公式形式相同,ρ的取值范围为 -1≤ ρ ≤ 1。 式(6.1)中ut-1是ut滞后一期的随机误差项,因此,将式(6.1)计算的自相关系数ρ称为一阶自相关系数。 式(6.1)

二、自相关产生的原因 1、经济系统的惯性 自 相 关 产 2、经济活动的滞后效应 生 的 原 因 3、数据处理造成的相关 4、蛛网现象 5、模型设定偏误

自相关现象大多出现在时间序列数据中,而经济系统的经济行为都具有时间上的惯性。 原因1-经济系统的惯性 自相关现象大多出现在时间序列数据中,而经济系统的经济行为都具有时间上的惯性。 如GDP、价格、就业等经济指标都会随经济系统的周期而波动。例如,在经济高涨时期,较高的经济增长率会持续一段时间,而在经济衰退期,较高的失业率也会持续一段时间,这种现象就会表现为经济指标的自相关现象。

滞后效应是指某一指标对另一指标的影响不仅限于当期而是延续若干期。由此带来变量的自相关。 原因2- 经济活动的滞后效应 滞后效应是指某一指标对另一指标的影响不仅限于当期而是延续若干期。由此带来变量的自相关。 例如,居民当期可支配收入的增加,不会使居民的消费水平在当期就达到应有水平,而是要经过若干期才能达到。因为人的消费观念的改变客观上存在自适应期。

因为某些原因对数据进行了修整和内插处理,在这样的数据序列中就会有自相关。 原因3-数据处理造成的相关 因为某些原因对数据进行了修整和内插处理,在这样的数据序列中就会有自相关。 例如,将月度数据调整为季度数据,由于采用了加合处理,修匀了月度数据的波动,使季度数据具有平滑性,这种平滑性产生自相关。对缺失的历史资料,采用特定统计方法进行内插处理,使得数据前后期相关,产生了自相关。

供给t=β1+β0价格t-1+ut 原因4-蛛网现象 蛛网现象是微观经济学中的一个概念。它表示某种商品的供给量受前一期价格影响而表现出来的某种规律性,即呈蛛网状收敛或发散于供需的均衡点。 许多农产品的供给呈现为蛛网现象,供给对价格的反应要滞后一段时间,因为供给需要经过一定的时间才能实现。如果时期t的价格Pt低于上一期的价格Pt-1,农民就会减少时期t+1的生产量。如此则形成蛛网现象,此时的供给模型为: 供给t=β1+β0价格t-1+ut (6.2)

原因5-模型设定偏误 如果模型中省略了某些重要的解释变量或者模型函数形式不正确,都会产生系统误差,这种误差存在于随机误差项中,从而带来了自相关。由于该现象是由于设定失误造成的自相关,因此,也称其为虚假自相关。 例如,应该用两个解释变量解释Y,即: (6.3) 而建立模型时,模型设定为: (6.4) 则X3t对Yt的影响在式(6.4)中便归入随机误差项ut中,由于X3t在不同观测点上是相关的,这就造成了ut在不同观测点是相关的,呈现出系统模式,此时ut是自相关的。

模型形式设定偏误也会导致自相关现象。如将U形成本曲线设定为线性成本曲线,则必定会导致自相关。由设定偏误产生的自相关是一种虚假自相关,可通过改变模型设定予以消除。 自相关关系主要存在于时间序列数据中,但是在横截面数据中,也可能会出现自相关,通常称其为空间自相关(Spatial auto correlation)。

例如,在消费行为中,一个家庭、一个地区的消费行为可能会影响另外一些家庭和另外一些地区,就是说不同观测点的随机误差项可能是相关的。 多数经济时间序列在较长时间内都表现为上升或下降的超势,因此大多表现为正自相关。但就自相关本身而言是可以为正相关也可以为负相关。

三、自相关的表现形式 自相关的性质可以用自相关系数ρ的符号判断自相关的方向,即ρ <0为负相关, ρ >0为正相 关。当| ρ |接近1时,表示相关的程度很高。但自相关是u1,u2,…,un序列自身的相关,因n个随机误差项的关联形式不同而具有不同的自相关形式。自相关大多出现在时间序列数据中。

ut =ρut-1 + vt ; -1<ρ< 1 对于样本观测期为n的时间序列数据,可得到总体回归模型(PRF)的随机误差项为u1,u2,…,un,如果自相关形式为 ut =ρut-1 + vt ; -1<ρ< 1 (6.5) 其中ρ 为自相关系数,vt 为经典误差项,即E(vt) = 0,Var (vt) =σ2,Cov (vt,vt+s) = 0,s≠0。则式(6.5)称为一阶自回归模式,记为AR (1)。因为模型(6.5)中ut-1是ut 滞后一期的值,因此称为一阶。式(6.5)中的ρ也称为一阶自相关系数。

如果式(6.5)中的随机误差项vt不是经典误差项,即vt中包含有ut的成份,如包含有ut-2则需将ut-2显含在回归模型中,其为 ut = ρ1ut-1+ρ2ut-2 +vt′ (6.6) 其中,ρ1为一阶自相关系数, ρ2为二阶自相关系数, vt′是经典误差项。式(6.6)称为二阶自回归模式,记为AR(2)。一般地,如果u1,u2,…,ut之间的关系为

(6.7) 其中,vt 为经典误差项。则称式(6.7)为m阶自回归模式,记为AR(m)。 在经济计量分析中,通常采用一阶自回归形式,即假定自回归形式为一阶自回归AR(1)。

第二节 自相关的后果 自相关的后果 对参数估计的影响 对模型检验的影响 对模型预测的影响

一、对参数估计的影响 在有自相关的条件下,仍然使用普通最小二乘法将低估估计量 的方差 。 将低估真实的 。

(6.9) 对于一元线性回归模型: (6.8) 假定随机误差项u存在一阶自相关: 其中, 为现期随机误差, 为前期随机误差。 其中, 为现期随机误差, 为前期随机误差。 是经典误差项,满足零均值 ,同方差 ,无自相关 的假定。

将随机误差项 的各期滞后值: 逐次代入式(6.9)可得: (6.10)

式(6.10)表明随机误差项 可表示为独立同分布的随机误差序列 , , ,… 的加权和,权数分别为1, , ,…,当0<ρ<1时,这些权数是随时间推移而呈几何衰减的;而当 -1<ρ <0时,这些权数是随时间推移而交错振荡衰减的。

由式(6.10)可推得: (6.11) (6.12) 式(6.11)、式(6.12)表明,在 为一阶自回归的相关形式时,随机误差 依然是零均值、同方差的误差项。

(6.13) 由于现期的随机误差项 并不影响回归模型中随机误差项 的以前各期值 (k>0),所以 与 不 相关,即有 。因此,由式(6.9)可得随机误差项 与其以前各期 的协方差分别为: (6.13)

(6.14)

以此类推,可得 : (6.15) 这些协方差分别称为随机误差项 的一阶自协方差、二阶自协方差和 k 阶自协方差。 对于一元线性回归模型(6.8),当 u 为经典误差项时,普通最小二乘估计量 的方差为: (6.16)

当式(6.8)的随机误差项u有自相关时, 依然是无偏的,即 ,这一点在普通最小二乘法无偏性证明中可以看到。因为,无偏性证明并不需要 u 满足无自相关的假定。那么,最小二乘估计量 是否是有效的,即在所有线性无偏估计量中方差最小呢?下面我们将以说明。

式(6.8)中的最小二乘估计量为 (6.17) 因此: (6.18)

当存在自相关时,普通最小二乘估计量不再是最佳线性无估计量,即它在线性无偏估计量中不是方差最小的。在实际经济系统中,通常存在正的自相关,即ρ>0,同时X 序列自身也呈正相关,因此式(6.18)右边括号内的值通常大于0。因此,在有自相关的条件下,仍然使用普通最小二乘法将低估估计量 的方差 。 将低估真实的 。

二、对模型检验的影响 考虑自相关时的检验 对模型检验的影响 忽视自相关时的检验

(一)考虑自相关时的检验 由于 并不是所有线性无偏估计量中最小的,使用t 检验判断回归系数的显著性时就可能得到错误的结论。

(二)忽视自相关时的检验 如果我们忽视自相关问题依然假设经典假定成立,使用 ,将会导致错误结果。 如果我们忽视自相关问题依然假设经典假定成立,使用 ,将会导致错误结果。 当ρ >0,即有正相关时,对所有的j 都有ρ >0。另外回归模型中的解释变量在不同时期通常是正相关的,对于xt和xt+j来说 是大于0的。

因此,由式(6.18)可以看出,普通最小二乘法的方差 通常会低估 的真实方差。当ρ 较大和Xt 有较强的正自相关时,普通最小二乘估计量的方差会有很大偏差,我们会夸大估计量的估计精度,即得到较小的标准误。 因为 的标准误是 的标准差的估计值,在有自相关时,普通最小二乘估计 的标准误就不可靠了。

一个被低估了的标准误意味着一个较大的t 统计量。因此,当ρ>0时,通常t 统计量都很 大。这种有偏的t 统计量是不能用来判断回归系数的显著性的。 类似地,由于自相关的存在,参数的最小二乘估计量是无效的,使得F 检验和 检验也是不可靠的。

三、对模型预测的影响 模型预测的精度决定于抽样误差和总体误差项的方差 。抽样误差来自于对 的估计,在自相关情形下, 的方差的最小二乘估计变得不可靠,由此必定加大抽样误差。同时,在自相关情形下,对 的估计 也会不可靠。由此可看出,影响预测精度的两大因素都会因自相关的存在而加大不确定性,使预测的置信区间不可靠,从而降低预测的精度。

第三节 自相关的检验 当随机误差项存在自相关时会给普通最小二乘法的应用带来严重的后果。因此,如何诊断回归模型是否存在自相关就致关重要。

自相关的检验 图示检验法 DW检验法

一、图示检验法 图示法是一种直观的诊断方法,它是把给定的回归模型直接用普通最小二乘法估计参数,求出残差项et ,et作为ut随机项的真实估计值,再描绘et 的散点图,根据散点图来判断et的相关性。残差et的散点图通常有两种绘制方式 。

et et-1 图 6.1 et与et-1的关系 绘制et-1 ,et 的散点图。用(et-1 ,et )(t = 1,2,…n)作为散布点绘图,如果大部分点落在第Ⅰ、Ⅲ象限,表明随机误差项ut存在着正自相关,如图6.1所示。

如果大部分点落在第Ⅱ、Ⅳ象限,那么随机误差项ut 存在着负自相关,如图6.2所示。 et-1 et 图 6.2 et与et-1的关系 如果大部分点落在第Ⅱ、Ⅳ象限,那么随机误差项ut 存在着负自相关,如图6.2所示。

t et 图 6.3 et的分布 按照时间顺序绘制回归残差项 et 的图形。如果et (t=1,2,…,n)随着t 的变化逐次有规律地变化,呈现锯齿形或循环形状的变化,就可断言 et存在相关,表明存在着自相关;如果et随着t的变化逐次变化并不断地改变符号,那么随机误差项ut存在负自相关;如图6.3所示。

图 6.4 et的分布 t et 如果et随着t的变化逐次变化并不频繁地改变符号,而是几个正的et后面跟着几个负的,则表明随机误差项存ut 在正自相关,如图6.4所示。

二、DW检验法 DW检验是J.Durbin(杜宾)和G.S.Watson(沃特森)于1951年提出的一种适用于小样本的检验方法。DW检验只能用于检验随机误差项具有一阶自回归形式的自相关问题。这种检验方法是建立经济计量模型中最常用的方法,一般的计算机软件都可以计算出DW值。

随机误差项的一阶自回归形式为: (6.20) 为了检验序列的相关性,构造的原假设是: (6.21) 为了检验上述假设,构造DW 统计量首先要求出回归估计 式的残差et 定义DW统计量为 : (6.22)

DW值的取值范围 (6.23) 则: (6.24) 由于

4 (2,4) 2 (0,2) -1 (-1,0) (0,1) 1 DW 因此, 。所以,DW 值与 的对 应关系如表6.1所示。 表 6.1 4 (2,4) 2 (0,2) -1 (-1,0) (0,1) 1 DW

由上述讨论可知DW的取值范围为:0≤DW≤4 根据样本容量 n 和解释变量的数目k′(不包括常数项)查DW分布表,得临界值 dL 和 dU ,然后依下列准则考察计算得到的DW值,以决定模型的自相关状态。

表6.2 DW 检验决策规则 误差项u1,u2,…,un间存在负相关 4-dL≤DW≤4 不能判定是否有自相关 4-dU≤DW<4-dL 误差项u1,u2,…,un间无自相关 dU<DW<4-dU dL<DW≤dU 误差项u1,u2,…,un间存在正相关 0≤DW≤dL

表6.2可以用坐标图更加直观地表示出来: f (DW) 不能确定 正自相关 无自相关 负自相关 DW 4 4-dL 4- dU 2 dU

① DW检验有两个不能确定的区域,一旦DW值落在这两个区域,就无法判断。这时,只有增大样本容量或选取其他方法;

② DW统计量的上、下界表要求n≥15,这是因为样本如果再小,利用残差就很难对自相关的存在性做出比较正确的诊断; ④ 只适用于有常数项的回归模型并且解释变量中不能含滞后的被解释变量。

第四节 自相关的补救 广义差分法 Cochrane - Orcutt迭代法 第四节 自相关的补救 对于虚假自相关,由于是模型设定偏误造成的,只能通过改变模型的设定去消除。对于设定正确的模型,如随机误差项有自相关,则为真实的自相关。对于真实的自相关可采用如下方法予以消除。 广义差分法 Cochrane - Orcutt迭代法 其它方法 :一阶差分法 德宾两步法

一、广义差分法 对于自相关的结构已知的情形可采用广义差分法解决。 由于随机误差项ut是不可观测的,通常我们假定ut为一阶自回归形式,即 ut =ρut-1 + vt (6.25) 其中,|ρ|<1,vt为经典误差项。 当自相关系数为已知时,使用广义差分法,自相关问题就可彻底解决。我们以一元线性回归模型为例说明广义差分法的应用。

对于一元线性回归模型 (6.26) 将模型(6.26)滞后一期可得 (6.27) 用ρ 乘式(6.27)两边,得 (6.28)

用式(6.26)减去式(6.28)可得 (6.29) 模型,随机误差项无序列相关。式(6.29)中,ut-ρut-1 = vt 是经典误差项。因此,模型(6.29)已经是经典线性回归 令

对模型(6.30)使用普通最小二乘估计就会得到参数估计的最佳线性无偏估计量。 则式(6.29) 可表示为: (6.30) 对模型(6.30)使用普通最小二乘估计就会得到参数估计的最佳线性无偏估计量。 式(6.29)称为广义差分方程,因为被解释变量与解释变量均为现期值减去前期值的一部分,由此而得名。

在进行广义差分时,解释变量X与被解释变量Y 均以差分形式出现,因而样本容量由n减少为 n-1,即丢失了第一个观测值。如果样本容量较大,减少一个观测值对估计结果影响不大。但是,如果样本容量较小,则对估计精度产生较大的影响。此时,可采用普莱斯-温斯滕(Prais-Winsten)变换,将第一个观测值变换为 补充到差分序列 中,再使用普通最小二乘法估计参数。

二、Cochrane - Orcutt 迭代法 在实际应用中,自相关系数 ρ 往往是未知的,ρ 必须通过一定的方法估计。最简单的方法是据DW统计量估计ρ 。由DW 与ρ 的关系可知 : 但是,式(6.31)得到的是一个粗略的结果, 是对ρ 精度不高的估计。其根本原因在于我们对有自相关的回归模型使用了普通最小二乘法。为了得到ρ的精确的估计值 ,人们通常采用科克伦—奥克特(Cochrane—Orcutt)迭代法。

该方法利用残差et 去估计未知的ρ。对于一元线性回归模型 (6.32) 假定ut 为一阶自回归形式,即 : (6.33)

科克伦-奥克特迭代法估计ρ 的步骤如下: 使用普遍最小二乘法估计模型 (6.32) 并获得残差: 1

利用残差 做如下的回归 (6.34) 2

3 利用式(6.34)得到的 ,对模型 (6.32)进行广义差分,即 对式(6.35)使用普通最小二乘法,可得样本回归函数为: 令 利用式(6.34)得到的 ,对模型 (6.32)进行广义差分,即 (6.35) (6.36) 令 对式(6.35)使用普通最小二乘法,可得样本回归函数为:

4 因为据式(6.34)得到 并不是对ρ 的最佳估计,必须进一步迭代,寻求最佳估计。由前一步估计的结果有: 因为据式(6.34)得到 并不是对ρ 的最佳估计,必须进一步迭代,寻求最佳估计。由前一步估计的结果有: 将 代入原回归方程(6.32),求得新的残差如下: (6.37)

5 利用残差 做如下的回归 (6.38) 这里得到的 就是 ρ 的第二轮估计值。

我们并不能确认 是否是ρ 的最佳估计值,我们还要继续估计ρ 的第三轮估计值 。这种方法是迭代法,当估计的 与 相差很小时,就找到了ρ 的最佳估计值。在模型实践中,通常采用科克伦—奥克特两步法。第一步,据式(6.34)估计ρ 。第二步,用估计值 做广义差分,估计广义差分方程的参数。这种两步法与上述迭代方法结果很相近。

三、其它方法简介 (一)一阶差分法 一阶差分法是消除序列相关的一种简单有效的方法。我们仍以一元线性回归模型来说明一阶差分法的应用。

式中,ut为一阶自回归AR(1)。将模型(6.39)变换为 : 如果原模型存在完全一阶正自相关,即ρ =1则 (6.41) 其中,vt 为经典误差项。则式(6.40)的随机误差项为经典误差项,无自相关问题。对式 (6.40)使用普通最小二乘法估计参数,可得到最佳线性无偏估计量。

(二)德宾两步法 当自相关系数未知时,也可采用德宾提出的两步法,消除自相关。将广义差分方程(6.29)表示为:

采用如下的两个步骤消除自相关。 第一步,把式(6.42)作为一个多元回归模型,使用普通最小二乘法估计参数。把Yt-1的回归系数 看作ρ 的一个估计值 。 第二步,求得 后,使用 进行广义差分,求得序列: 和 然后使用普通最小二乘法对广义差分方程估计参数,求得最佳线性无偏估计量。

第五节 案例分析 研究范围:中国农村居民收入-消费模型 (1985~2003) 第五节 案例分析 研究范围:中国农村居民收入-消费模型 (1985~2003) 研究目的:消费模型是研究居民消费行为的工具和手段。通过消费模型的分析可判断居民消费边际消费倾向,而边际消费倾向是宏观经济系统中的重要参数。 建立模型 Yt-居民消费,Xt-居民收入,ut-随机误差项。 数据收集:1985~2003年农村居民人均收入和消费 (见表6.3)

表6.3 1985-2003年农村居民人均收入和消费 单位:元 年份 全年人均纯收入 (现价) 全年人均消费 性支出 (现价) 消费价格 表6.3 1985-2003年农村居民人均收入和消费 单位:元 年份 全年人均纯收入 (现价) 全年人均消费 性支出 (现价) 消费价格 指数 (1985=100) 人均实际纯 收入 (1985可比价) 人均实际消费 性支出 1985 397.60 317.42 100.0 317.40 1986 423.80 357.00 106.1 399.43 336.48 1987 462.60 398.30 112.7 410.47 353.42 1988 544.90 476.70 132.4 411.56 360.05 1989 601.50 535.40 157.9 380.94 339.08 1990 686.30 584.63 165.1 415.69 354.11 1991 708.60 619.80 168.9 419.54 366.96 1992 784.00 659.80 176.8 443.44 373.19 1993 921.60 769.70 201.0 458.51 382.94

续 表 年份 全年人均纯 收入 (现价) 全年人均消费性支出 消费价格 指数 (1985=100) 人均实际纯 (1985可比价) 续 表 年份 全年人均纯 收入 (现价) 全年人均消费性支出 消费价格 指数 (1985=100) 人均实际纯 (1985可比价) 人均实际消费 性支出 1994 1221.00 1016.81 248.0 492.34 410.00 1995 1577.70 1310.36 291.4 541.42 449.69 1996 1923.10 1572.10 314.4 611.67 500.03 1997 2090.10 1617.15 322.3 648.50 501.77 1998 2162.00 1590.33 319.1 677.53 498.28 1999 2214.30 1577.42 314.3 704.52 501.75 2000 2253.40 1670.00 314.0 717.64 531.85 2001 2366.40 1741.00 316.5 747.68 550.08 2002 2475.60 1834.00 315.2 785.41 581.85 2003 2622.24 1943.30 320.2 818.86 606.81

据表6.3的数据使用普通最小二乘法估计消费模型得: (6.44) = 0.9788,F = 786.0548,df =17, DW = 0.7706,该回归方程可决系数较高,回归系数均显著。对样本量为19、一个解释变量的模型、5%显著水平,查DW 统计表可知,dL=1.18,dU= 1.40,模型中DW<dL,显然消费模型中有自相关。这一点也可从残差图中看出,点击EViews方程输出窗口的按钮Resids可得到残差图,如图6.6所示。

图6.6 残差图

自相关问题的处理: 使用科克伦-奥克特的两步法解决自相关问题据模型(6.44)可得残差序列et,在EViews中,每次回归的残差存放在resid序列中,为了对残差进行回归分析,需生成命名为e 的残差序列。在主菜单选择Quick/Generate Series或点击工作文件窗口工具栏中的Procs/Generate Series,在弹出的对话框中输入e = resid,点击OK得到残差序列et。使用et进行滞后一期的自回归,在EViews命今栏中输入ls e e (-1)可得回归方程:

(6.45) 由式(6.45)可知 =0.4960,对原模型进行 广义差分,得到广义差分方程: (6.39) 由式(6.45)可知 =0.4960,对原模型进行 广义差分,得到广义差分方程: (6.39) 对式(6.39)的广义差分方程进行回归,在EViews命令栏中输入ls Y-0.4960*Y (-1) c X-0.4960*X (-1)回车后可得方程输出结果如表6.4。

表6.4 广义差分方程输出结果 Dependent Variable: Y-0.496014*Y(-1) 表6.4 广义差分方程输出结果 Dependent Variable: Y-0.496014*Y(-1) Method: Least Squares Date: 03/26/05 Time: 12:32 Sample(adjusted): 1986 2003 Included observations: 18 after adjusting endpoints Variable Coefficient Std. Error t-Statistic Prob. C 60.44431 8.964957 6.742287 0.0000 X-0.496014*X(-1) 0.583287 0.029410 19.83325 R-squared 0.960914 Mean dependent var 231.9218 Adjusted R-squared 0.958472 S.D. dependent var 49.34525 S.E. of regression 10.05584 Akaike info criterion 7.558623 Sum squared resid 1617.919 Schwarz criterion 7.657554 Log likelihood -66.02761 F-statistic 393.3577 Durbin-Watson stat 1.397928 Prob(F-statistic) 0.000000

由表6.4可得回归方程为: (6.46) R2 = 0.9609 F = 393.3577 df = 16 DW = 1.3979 式中, 由于使用了广义差分数据,样本容量减少了1个,为18个。查5%显著水平的DW 统计表可知 dL= 1.16,dU = 1.39,模型中DW =1.3979>dU,说明广义差分模型中已无自相关。同时,可决系数R2、t、F 统计量均达到理想水平。

如果我们使用普莱斯-温斯腾变换补充第一个观测值,即 对比模型(6.44)和(6.46),很明显普通最小二乘法低估了回归系数的标准误。原模型中SE( )= 0.0214,广义差分模型中为SE ( )= 0.0294。 如果我们使用普莱斯-温斯腾变换补充第一个观测值,即 得到普莱斯-温斯腾变换的广义差分模型为:

对比模型(6.47)和(6.46)可发现两者的参数估计值和各检验统计量的差别很微小,说明在本例中使用普莱斯-温斯腾变换与直接使用科克伦-奥克特两步法的估计结果无显著差异,这是因为本例中的样本还不算太小。如果实际应用中样本较小,则两者的差异就会较大。通常对于小样本,应采用普莱斯-温斯腾变换补充第一个观测值。

由差分方程(6.46)可知: (6.48) 由此,我们得到最终的中国农村居民消费模型: 由模型(6.49)的中国农村居民消费模型可知,中国农村居民的边际消费倾向为0.5833,即中国农民每增加收入1元,将增加消费支出0.5833元。

本章小结 1.当总体回归模型的随机误差项在不同观测点上彼此相关时就产生了自相关问题。 2.自相关的出现有多种原因。时间序列的惯性、模型设定错误、数据的处理等等。 3.在出现自相关时,普通最小二乘估计量依然是无偏、一致的,但不再是有效的。通常的t 检验和F 检验都不能有效地使用。 4.为了研究问题的方便和考虑实际问题的代表意义,我们通常将自相关设定为一阶自相关即AR(1)模式。用一阶自相关系数 ρ 表示自相关的程度与方向。当然,实际问题也存在AR(m)模式或其它模式。

5.由于ut是不可观测的,通常我们使用ut的估计量et判断ut的特性。我们可通过et的图形判断自相关的存在,也可使用依据et计算的DW 统计量判断自相关的存在。 6.如果自相关系数 ρ 是已知的,我们可以使用广义差分法消除序列相关。 7.如果自相关系数是 ρ 未知的,我们可采用科克伦-奥克特迭代法求得 ρ 的估计值,然后用广义差分法消除序列相关。

第 六 章 结 束 了! THANKS