(精选)ECM误差修正模型
- 格式:docx
- 大小:44.23 KB
- 文档页数:6
Error Correction Model 用EVIEWS怎么做一、利用EG两步法做协整检验。
在两个变量情况下(设为Y、X),包括两序列单整检验、两变量最小二乘法回归并得到残差序列并命名为e、对e作单位根检验。
二、在证明Y、X两序列间存在协整后,才可以建立ECM。
其中,误差修正项ecm的值就是之前的回归模型的残差序列e。
三、直接输入以下命令:ls y c y(-1) x x(-1)得到的估计结果在实际预测时比较方便,不过需要计算得到ecm项的系数。
四、也可以直接输入以下命令:ls y c x e(-1)其中,e(-1)项的系数就是ecm项的系数。
这个模型的优点是直观,但是不便于预测。
五、两种估计是等价的。
六、建议参考阅读易丹辉:《数据分析与EViews应用》,中国统计出版社2002年版。
(也许有新版也不一定)对于误差修正模型,需要先建立一个模型,然后进行回归分析,分析它的短期均衡关系。
操作:举个例子说,比如试图建立y对y(-1)和x的误差修正模型。
STEP1 建立长期关系ls y c y(-1) xSTEP2 对残差进行单位根检验来检验协整关系ecm=residuroot(10,h) ecmSTEP3 建立误差修正模型ls d(y) c d(y(-1)) d(x) ecm(-1)教程:案例1上面的分析可以证明序列lconsume、lincome及lconsme(-1)之间存在协整关系,故可以建立ecm(误差修正模型)。
先分别对序列lconsume、lincome及lconsme(-1)进行一阶差分,然后对误差修正模型进行估计。
在主窗口命令行中输入:ls d(lconsume) c d(lincome) d(lconsume(-1)) ecm(-1)此时的常数项系数不明显,我们去掉常数项后再进行回归,结果如下图8.6所示图8.6从上式可以看出上式中的T检验值均显著,误差修正项的系数为-0.252,这说明长期均衡对短期波动的影响不大。
实验报告(二)——误差修正模型(ECM)的建立与分析一、单位根检验:1、绘制cons与GDP的时间序列图:从时间序列图中可以看出,cons与GDP随时间增加都呈上升趋势,表现出非平稳性。
2、对cons进行单位根检验:先选择对原序列(level)进行单位根检验,根据cons与GDP的时间序列图的走势,选择trend and intercept的检验方法,在maximum lags中填写ADF 检验方法的滞后期为0,从上表中可以看出,P值为0.9888,大于0.05的显著性水平,说明原序列是非平稳的。
选择cons的一阶差分(1st)和trend and intercept,从上表中可以看出,经过一阶差分后,P值(=0.5099)仍然没有通过0.05的置信水平检验,说明是不平稳的,需要继续改进。
再试用ADF检验,在滞后期(maximum lags)中填入8,选择一阶差分和trend and intercept,得出上表,可以看出P值=0.0801,大于0.05,没有通过0.05的置信水平检验,说明是不平稳的,需要继续改进。
再试用ADF检验,在滞后期(maximum lags)中填入6,选择二阶差分和trend and intercept,得出上表,可以看出P值=0.0137,小于0.05,通过0.05的置信水平检验,说明是平稳的。
3、对GDP进行单位根检验:先选择对原序列(level)进行单位根检验,根据cons与GDP的时间序列图的走势,选择trend and intercept的检验方法,在maximum lags中填写ADF 检验方法的滞后期为0,从上表中可以看出,P值为1.0000,大于0.05的显著性水平,说明原序列是非平稳的。
选择GDP的一阶差分(1st)和trend and intercept,从上表中可以看出,经过一阶差分后,P值(=0.5574)仍然没有通过0.05的置信水平检验,说明是不平稳的,需要继续改进。
实验八:协整关系检验与误差修正模型(ECM)new实验八:协整关系检验与误差修正模型(ECM)一、实验目的通过上机实验,使学生加深对时间序列之间协整关系的理解,能够运用Eviews 软件检验时间序列数据之间的协整关系并以此估计误差修正模型(ECM)。
二、预备知识(1)用EViews估计线性回归模型的基本操作;(2)时间序列数据的协整关系及其检验方法;(3)误差修正模型的结构及估计方法。
三、实验内容(1)用EViews检验两个时间序列数据的协整关系;(2)用EViews估计误差修正模型;四、实验步骤(一)、建立工作文件sy8.wf1及导入数据打开sy8.xls文件,运用前面学过的方法,在EViews新建一个工作文件sy8.wf1,把sy8.xls的数据导入到EViews,并根据得到人均消费(consp)和人均GDP(gdpp)两个序列,分别计算对应的自然对数,即lnc=log(consp)、lngdp=log(gdpp)。
(二)、分别检验序列lnc和lngdp的单整阶数。
运用图示法观察序列的时间路径图,如图8-1所示。
可见,lnc和lngdp都随时间不断上升,表明两者都是非平稳的。
(再运用自相关函数法,判断lnc 的平稳性。
打开lnc 序列的窗口,点击view\Correlogram ,设定滞后阶数为12,可得样本自相关系数图,操作和结果分别如图8-2和图8-3所示。
可见,lnc 是非平稳的。
再分析lnc 的一阶差分是否平稳。
在自相关函数图中,设定显示序列的一阶差分(1st differenc )后,再观察其样本自相关函数图,设定和结果如图8-4和图8-5所示。
可见,lnc 取一阶差分后就达到平稳,因此,lnc 是一阶单整序列,即I(1)序列。
如果采用单位根检验,结果相同。
同理,也可检验得到lngdp 序列是I(1)序列。
(三)运用Engle-Granger 方法(即EG 检验)检验consp 与gdpp 的协整关系。
误差修正模型(Error Correction Model)误差修正模型的产生原因对于非稳定时间序列,可通过差分的方法将其化为稳定序列,然后才可建立经典的回归分析模型。
如:建立人均消费水平(Y)与人均可支配收入(X)之间的回归模型:Y t = α0 + α1X t + μt如果Y与X具有共同的向上或向下的变化趋势,进行差分,X,Y成为平稳序列,建立差分回归模型得:ΔY t = α1ΔX t + v t式中,v t = μt−μt− 1然而,这种做法会引起两个问题:(1)如果X与Y间存在着长期稳定的均衡关系Y t = α0 + α1X t + μt且误差项μt不存在序列相关,则差分式ΔY t = α1ΔX t + v t中的v t是一个一阶移动平均时间序列,因而是序列相关的;(2)如果采用差分形式进行估计,则关于变量水平值的重要信息将被忽略,这时模型只表达了X与Y间的短期关系,而没有揭示它们间的长期关系。
因为,从长期均衡的观点看,Y在第t期的变化不仅取决于X本身的变化,还取决于X 与Y在t-1期末的状态,尤其是X与Y在t-1期的不平衡程度。
另外,使用差分变量也往往会得出不能令人满意回归方程。
例如,使用ΔY1 = ΔX t + v t回归时,很少出现截距项显著为零的情况,即我们常常会得到如下形式的方程:式中,(*)在X保持不变时,如果模型存在静态均衡(static equilibrium),Y也会保持它的长期均衡值不变。
但如果使用(*)式,即使X保持不变,Y也会处于长期上升或下降的过程中,这意味着X与Y间不存在静态均衡。
这与大多数具有静态均衡的经济理论假说不相符。
可见,简单差分不一定能解决非平稳时间序列所遇到的全部问题,因此,误差修正模型便应运而生。
误差修正模型的概述误差修正模型(Error Correction Model,简记为ECM)是一种具有特定形式的计量经济学模型,它的主要形式是由Davidson、Hendry、Srba和Yeo于1978年提出的,称为DHSY 模型。
r语言ecm模型R语言ECM模型:时间序列分析中的必备工具引言时间序列分析是经济学、金融学、统计学等领域中常用的一种分析方法,通过对一系列按时间顺序排列的数据进行分析,揭示出数据中的趋势、周期性和随机性等重要特征。
在时间序列分析中,误差修正模型(Error Correction Model,ECM)是一种非常常用且有效的模型。
本文将介绍R 语言中的ECM模型,从安装R软件到具体应用,逐步介绍ECM模型的原理和使用方法。
一、安装R软件R语言是一种开源的统计计算和图形化编程语言,它提供了丰富的数据分析和统计建模功能,被广泛应用于学术研究和工业界。
R语言的优点在于它具有强大的社区支持和丰富的扩展包,可以满足不同领域、不同需求的数据分析任务。
要使用R语言进行ECM模型的分析,首先需要安装R软件。
在官网(二、载入数据在使用ECM模型前,首先需要准备好要分析的时间序列数据。
R语言提供了丰富的数据获取和导入方法,可以从本地文件、网络数据源等获取数据,并将数据导入到R的工作环境中进行分析和建模。
以一个简单的例子为例,假设我们想分析某公司每月销售额的变化。
我们已经准备好了一份名为“sales.csv”的数据文件,其中包含两列数据,分别是日期和销售额。
首先,我们需要将这个数据文件导入到R中进行分析。
使用如下代码将数据文件读取到R的工作环境中:Rsales <- read.csv("sales.csv")销售数据现在存储在一个名为“sales”的数据框中,我们可以使用该数据框进行后续的分析。
三、数据预处理在进行ECM模型的建模前,我们通常需要对原始数据进行预处理。
这包括数据清洗、缺失值处理、时序特性检验等步骤,以确保数据的准确性和可靠性。
1. 数据清洗在进行数据清洗时,我们需要检查数据是否存在异常值、重复值等问题,并进行相应的处理或修正。
R语言提供了强大的数据操作和处理函数,可以快速、高效地对数据进行清洗。
误差修正模型(Error Correction Model,ECM)第二节 误差修正模型(Error Correction Model ,ECM )一、误差修正模型的构造对于y t 的(1,1)阶自回归分布滞后模型:t t t t t y x x y εβββα++++=--12110在模型两端同时减y t-1,在模型右端10-±t x β,得:tt t t t t t t tt t t t x y x x y x y x x y εααγβεββββαββεββββα+--+∆=+---+--+∆=+-+++∆+=∆------)(])1()1()[1()1()(1101012120120121100其中,12-=βγ,)1/()(200ββαα-+=,)1/(211ββα-=。
记 11011-----=t t t x y ecm αα (5-5) 则 t t t t ecm x y εγβ++∆=∆-10 (5-6) 称模型(5-6)为“误差修正模型”,简称ECM 。
二、误差修正模型的含义如果y t ~ I(1),x t ~ I(1),则模型(5-6)左端)0(~I y t ∆,右端)0(~I x t ∆,所以只有当y t 和x t 协整、即y t 和x t 之间存在长期均衡关系时,式(5-5)中的ecm~I(0),模型(5-6)两端的平稳性才会相同。
当y t 和x t 协整时,设协整回归方程为:t t t x y εαα++=10它反映了y t 与x t 的长期均衡关系,所以称式(5-5)中的ecm t -1是前一期的“非均衡误差”,称误差修正模型(5-6)中的1-t ecm γ是误差修正项,12-=βγ是修正系数,由于通常1||2<β,这样0<γ;当ecm t -1 >0时(即出现正误差),误差修正项1-t ecm γ< 0,而ecm t -1 < 0时(即出现负误差),1-t ecm γ> 0,两者的方向恰好相反,所以,误差修正是一个反向调整过程(负反馈机制)。
协整与误差修正模型在处理时间序列数据时,我们还得考虑序列的平稳性。
如果一个时间序列的均值或自协方差函数随时间而改变,那么该序列就是非平稳的。
对于非平稳的数据,采用传统的估计方法,可能会导致错误的推断,即伪回归。
若非平稳序列经过一阶差分变为平稳序列,那么该序列就为一阶单整序列。
对一组非平稳但具有同阶的序列而言,若它们的线性组合为平稳序列,则称该组合序列具有协整关系。
对具有协整关系的序列,我们算出误差修正项,并将误差修正项的滞后一期看做一个解释变量,连同其他反映短期波动关系的变量一起。
建立误差修正模型。
建立误差修正模型的步骤如下:首先,对单个序列进行单根检验,进行单根检验有两种:ADF( Augument Dickey-Fuller )和DF(Dickey-Fuller) 检验法。
若序列都是同阶单整,我们就可以对其进行协整分析。
在此我们只介绍单个方程的检验方法。
对于多向量的检验参见Johensen 协整检验。
我们可以先求出误差项,再建立误差修正模型,也可以先求出向量误差修正模型,然后算出误差修正项。
补充一点的是,误差修正模型反映的是变量短期的相互关系,而误差修正项反映出变量长期的关系。
下面我们给出案例分析。
案例分析在此,我们考虑从1978 年到2002 年城镇居民的人均可支配收入income 与人均消费水平consume 的关系,数据来自于《中国统计年鉴》,如表8.1 所示。
根据相对收入假设理论,在一定时期,人们的当期的消费水平不仅与当期的可支配收入、而且受前期的消费水平的影响,具有一定的消费惯性,这就是消费的棘轮效应。
从这个理论出发,我们可以建立如下 ( 8.1 ) 式的模型。
同时根据生命周期假设理论,消费者的消费不仅与当期收入有关,同时也受过去各项的收入以及对将来预期收入的限制和影响。
从我们下面的数据分析中,我们可以把相对收入假设理论与生命周期假设理论联系起来,推出如下的结果:当期的消费水平不仅与当期的可支配收入有关,而且还与前期的可支配收入、前两期的消费水平有关。
在此先对人均可支配收入和人均消费水平取对数,同时给出如下的模型Iconsume t 0 1lconsume t 1 2lincome t t=1,2,…,n (8.1 ) 如果当期的人均消费水平与当期的人均可支配收入及前期的人均消费水平均为一阶单整序列,而它们的线性组合为平稳序列,那么我们可以求出误差修正序列,并建立误差修正模型,如下:Iconsume0 1 Iincome t 2 Iconsum t e13ecm t 1 4t t=1,2,…,nt(8.2)ecm t =Iconsume t 0 1Iconsume t 12Iincome t 1t=1,2, … ,n(8.3)从( 8.2 )式我们可以推出如下的方程:Iconsume t(1 1 3)Iconsume t 1 ( 1 3 1 )Iconsume t 22Iincome t2 t(8.4 ))Iincome t 1 0 3 0 4t( 2 32在(8.2 )中Iconsume lincome分别为变量对数滞后一期的值,ecm( 1)为误差修正项,如( 8.3 )式所示。
(8.2 )式为含有常数项和趋势项的形式,我们省略了只含趋势项或常数项及二项均无的形式。
分析步骤:1、单位根检验。
我们先介绍ADF检验。
在检验过程中,若ADF检验值的绝对值大于临界值的绝对值,则认为被检验的序列为平稳序列。
在此我们先以对lincome的检验为例,在主菜单中选择Quick/Series Statistics/Unit Root Test,屏幕提示用户输入待检验序列名,输入lincome,单击OK进入单位根检验定义的对话框,如图8.1。
Tesl Ty p _________________ I rclude in test equationmini HHI ■■■■ VIIBI mmaawiiBimi■■ ■■■■■■ -r:q* jAug inn ent ed 0 icIcey-Fuller;Iniercept11 IIU■ in IAJI ■ HE iMiaiL&BiiiBJuiBiiuii■■JIIL mi Jir? Hhillips^l-'error* and inter csptTesthriritroot in:〉tlone*L BYC I1st difference Lagged diflerences:2nd difference图8.1对话框由三部分构成。
检验类型(Test Type )中默认项是ADF检验。
Test for unit root In 中可选择的是对原序列、一阶差分序列或是二阶差序列做单位根检验,在此我们保持默认的level,即原序列。
右上方的In clude in test equation 中,有三个选项,依次为含常数项,含常数项和趋势项,没有常数项且没有趋势。
在右下方的空格里默认为2,但我们一般根据AIC最小来确定滞后期数,本文选定为滞后一期。
检验的顺序为:先选含趋势项和常数项的检验,如果趋势项的T统计量不明显,就再选只含常图8.3数项的,如果常数项的T 统计量不明显,就选择常数项和趋势项均不包括的一项。
们选含趋势项和常数项的检验时,会出现下面的结果,如图8.2所示。
AiigmeniCifrl」Ljllei U'nll lluut 】uE uri I IF4CT% Chine al Value* 5% CMicnl Value 1[J%CnhLSi Wiiw *IMs_ krnnon cntiiL al values |口r rejecti sn ul ll^'pulheBi^ ol 勺 jrirt raol. AjjgrnB-led Dickey-尸 ul er Tesl Equa B ion De pt* nilen I! Var lab I B D(LlMC) Mnt Kfid L 前S<3UlinDate O2/13)Qd Time: 22:D1 Sample adjusledJ I9EO ^JUZIrtCludiJd uL'atruaiiuni! 23 Jfttr bi^u^l in<jj u r d^Oi Hl图8.2在检验的结果输出窗口中, 左上方为ADF 检验值,右上方为1%5囁口1 0%勺显著水平下的临界值,从图8.1中可以看出ADF 统计的检验值为-3.117,其绝对值小于10%勺显著水平的临界值-3.2856的绝对值。
同时趋势值的T 统计来看,在5% 的水平下显著。
注意,这里的T 统计量不同于我们在做最小二乘时用的T 统计值。
这些T 统计检验的临界值在F uller (1976)中给出.从上面的分析我们可以认为 该序列为非平稳的序列,且该序列有趋势项和常数项。
在下文中我们会进行一步介 绍只含常数项的和常数项与趋势项均不包括的ADF 检验的过程。
在上面分析的基础上,我们回到图 8.1的窗口,检验lincome 差分一阶的平稳 性。
在图8.1中的Test for unit root In 中选差分一阶,同时在 Include in test equation 中选取含趋势项和常数项这一项,我们同样根据 AIC 和SC 最小来选择滞 后两期。
此时会出现如下图8.3的结果:AiJipnufiieil Dlcldiy 」till^i ll»H Rom lisjjr 4111 DCLINCOMH) ADF Tesl Statiistic-3.3330231 % 亡mi 匚si ValuH*1-4.4E91 5% CihlcAl Value -3.6454ID% CiW»l Vteiia-3 2CT12*M^cKinr on critical </alu )e3 rar rejaction of hypothecs af m uni1 rootAuiiment^d Oicfc&y-Fuller Teal Equ'iiiDnOependem X^roble- ID(LlNCOME,2)i Mel hod: Le :3st Squiar-ssDate: 02^1 印04 T mt 22: ce S&itifi le( MjiNtsd) 1962 2002Inc fajd «rl mtar R r*atin nr ?1 rflRi ^riuFlinq R ntrlpciinisVdrlabls CoefTisieirt Std IzrtQi t-Statistic Arob. D(JJNCOME(1)) U3.-333D2-3Qg?D(UMCOME(r1]^)□ 31392S □,3J2T51 1.553530 □.14QD□J531J0 □.1SB35Z DElil 犯" UL03E3cU 'I02S75J □ [DGOJB 2 07OBII5 aonni^TF?EMD(197S) ・U 0DJQ J 3 □.□□1717 皿日旧d □.5Sd3R-s-qu-are J0 434dZ4 h/eai dep e nde nt o.ooaaaaAdjd£t&d iR-iquai-ed 0 299DE5 S.D. /arO.OSG405 S.E. of r»gr»ttion 0 0-RdS2 info ■J.0C22ZI2Sunn 二qu 就注时 D 03-FTE6 匚lifFimnLoq likelihQiid3T D48E5 Fstaiistic3.QT2ES1Cjui^in-WstMin stst1 S524I3Proh(F-<tsrtislie)Q,tM6&44当我ADF Test 3抽血亡 ・3.11 -4.416? ■3.6219 -3.2474从上图中可以看出ADF的绝对值小于5%水平下的临界值的绝对值,大于10%勺检验值的绝对值。
但此时趋势项的T检验值不明显。
所以我们回到图8.1的窗口,在In elude in test equation中选取含常数项这一项。
其结果如下图8.4所示,结果显示ADF的绝对值为3.4546大于5%水平下的临界值的绝对值,此时常数项的T检验值为3.34572,大于在显著水平为5%水平下的T临界值为 2.61,所以常数项T检验值很明显。
我们认为lincome序列差分一阶后为平稳的。
值得注意的是,我们在此选择10%为临界值来判断非平稳的情况,而选择5%勺临界值来判断平稳的情况,也就是,当ADF检验值的绝对值大于5%水平下的临界的绝对值。