计量经济学软件包Eviews使用说明
一、启动软件包
假定用户有Windows95/98的操作经验,我们通过一个实际问题的处理过程,使用户对EViews的应用有一些感性认识,达到速成的目的。
1、 Eviews的启动步骤:
进入Windows /双击Eviews快捷方式,进入EViews窗口;或点击开始 /程序/Econometric Views/Eviews,进入EViews窗口。
菜单栏 |
标题栏 |
命令窗口 |
控制按钮 |
信息栏 |
路径 |
状态栏 |
主显示窗口 |
(图一) |
2、EViews窗口介绍
标题栏:窗口的顶部是标题栏,标题栏的右端有三个按钮:最小化、最大化(或复原)和关闭,点击这三个按钮可以控制窗口的大小或关闭窗口。
菜单栏:标题栏下是主菜单栏。主菜单栏上共有7个选项: File,Edit,Objects,View,Procs,Quick,Options,Window,Help。用鼠标点击可打开下拉式菜单(或再下一级菜单,如果有的话),点击某个选项电脑就执行对应的操作响应(File,Edit的编辑功能与Word, Excel中的相应功能相似)。
命令窗口:主菜单栏下是命令窗口,窗口最左端一竖线是提示符,允许用户在提示符后通过键盘输入EViews(TSP风格)命令。如果熟悉MacroTSP(DOS)版的命令可以直接在此键入,如同DOS版一样地使用EViews。按F1键(或移动箭头),键入的历史命令将重新显示出来,供用户选用。
主显示窗口:命令窗口之下是Eviews的主显示窗口,以后操作产生的窗口(称为子窗口)均在此范围之内,不能移出主窗口之外。
状态栏:主窗口之下是状态栏,左端显示信息,中部显示当前路径,右下端显示当前状态,例如有无工作文件等。
二、创建工作文件
工作文件是用户与EViews对话期间保存在RAM之中的信息,包括对话期间输入和建立的全部命名对象,所以必须首先建立或打开一个工作文件用户才能与Eviews对话。工作文件好比你工作时的桌面一样,放置了许多进行处理的东西(对象),像结束工作时需要清理桌面一样,允许将工作文件保存到磁盘上。如果不对工作文件进行保存,工作文件中的任何东西,关闭机器时将被丢失。
进入EViews后的第一件工作应从创建新的或调入原有的工作文件开始。只有新建或调入原有工作文件, EViews才允许用户输入开始进行数据处理。
建立工作文件的方法:点击File/New/Workfile。选择数据类型和起止日期,并在出现的对话框中提供必要的信息:适当的时间频率(年、季度、月度、周、日);确定起止日期或最大处理个数(开始日期是项目中计划的最早的日期;结束日期是项目计划的最晚日期,非时间序列提供最大观察个数,以后还可以对这些设置进行更改)。
下面我们通过研究我国城镇居民消费与可支配收入的关系来学习Eviews的应用。数据如下:
表一
1998年我国城镇居民人均可支配收入与人均消费性支出 单位:元
地 区 | 可支配收 入(inc) | 消 费 性支 出(consum) | 地 区 | 可支配收 入(inc) | 消 费 性支 出(consum) |
北 京 | 8471.98 | 6970.83 | 河 南 | 4219.42 | 3415.65 |
天 津 | 7110.54 | 5471.01 | 湖 北 | 4826.36 | 4074.38 |
河 北 | 5084.64 | 3834.43 | 湖 南 | 5434.26 | 4370.95 |
山 西 | 4098.73 | 3267.70 | 广 东 | 8839.68 | 7054.09 |
内蒙古 | 4353.02 | 3105.74 | 广 西 | 5412.24 | 4381.09 |
辽 宁 | 4617.24 | 3890.74 | 海 南 | 4852.87 | 3832.44 |
吉 林 | 4206.64 | 3449.74 | 重 庆 | 5466.57 | 4977.26 |
黑龙江 | 4268.50 | 3303.15 | 四 川 | 5127.08 | 4382.59 |
上 海 | 8773.10 | 6866.41 | 贵 州 | 4565.39 | 3799.38 |
江 苏 | 6017.85 | 4889.43 | 云 南 | 6042.78 | 5032.67 |
浙 江 | 7836.76 | 6217.93 | 陕 西 | 4220.24 | 3538.52 |
安 徽 | 4770.47 | 3777.41 | 甘 肃 | 4009.61 | 3099.36 |
福 建 | 6485.63 | 5181.45 | 青 海 | 4240.13 | 3580.47 |
江 西 | 4251.42 | 3266.81 | 宁 夏 | 4112.41 | 3379.82 |
山 东 | 5380.08 | 4143.96 | 新 疆 | 5000.79 | 3714.10 |
(数据来源:中国统计年鉴-1999光盘J10、J11,中国统计出版社)
下面的图片说明了具体操作过程。
1、打开新建对象类型对话框,选择工作文件Workfile,见图二。
菜单栏 |
标题栏 |
命令窗口 |
控制按钮 |
信息栏 |
路径 |
状态栏 |
主显示窗口 |
(图一) |
(图二) |
(图三) |
2、打开工作文件时间频率和样本区间对话框,输入频率和样本区间,见图三。
信息栏兰 |
标题栏 |
工具条 |
控制框昂 |
对象 |
小图标 |
(图四) |
3、点击OK确认,得新建工作文件窗口,见图四。
工作文件窗口:工作文件窗口是EViews的子窗口。它有标题栏、控制按钮和工具条。标题栏指明窗口的类型workfile、工作文件名。标题栏下是工作文件窗口的工具条,工具条上有一些按钮。Views观察按钮、Procs过程按钮、Save(保存)工作文件、Sample(设置观察值的样本区间)、Gener(利用已有的序列生成新的序列)、Fetch(从磁盘上读取数据)、Store(将数据存储到磁盘)、Delete(删除)对象。此外,可以从工作文件目录中选取并双击对象,用户就可以展示和分析工作文件内的任何数据。工作文件一开始其中就包含了两个对象,一个是系数序列C(保存估计系数用),另一个残差序列RESID(实际值与拟合值之差)。小图标上标识出对象的类型,C是系数向量,曲线图是时间序列。用户选择Views对象后双击鼠标左建或直接使用EViews主窗口顶部的菜单选项,可以对工作文件和其中的对象进行一些处理。
4、保存工作成果:将工作成果保存到磁盘,点击工具条中save\输入文件名、路径\保存,或点击菜单栏中File \Save或Save as \输入文件名、路径\保存。
5、打开工作文件:我们可以打开一个已有的工作文件继续以前的工作,点击主菜单中的File \Open \ Workfile \选定文件\打开。
三、输入和编辑数据
建立或调入工作文件以后,可以输入和编辑数据。 输入数据有两种基本方法:data命令方式和鼠标图形界面方式
1、data命令方式:命令格式为:data <序列名1> <序列名2>......<序列名n>,序列名之间用空格隔开,输入全部序列后回车就进入数据编辑窗口,如图五所示。用户可以按照Excel的数据输入习惯输入数据。数据输入完毕,可以关闭数据输入窗口,点击工作文件窗口工具条的Save或点击菜单栏的File \ Save将数据存入磁盘。
2、鼠标图形界面方式——数组方式:点击Quick \ Empty Group (Edit Series), 进入数据窗口编辑窗口,点击obs 行没有数据的第一列(如图五中太阳标志处),然后输入序列名,并可以如此输入多个序列。输入数据名后,可以输入数据,方式同上。
3、鼠标图形界面方式——序列方式:点击Objects \ New object \ 选Series \ 输入序列名称\Ok,进入数据编辑窗口,点击Edit+/-打开数据编辑状态,(用户可以根据习惯点击Smpl+/-改变数据按行或列的显示形式,)然后输入数据,方式同上。
输入命令,回车 |
数据编辑窗口工具条 |
序列名称 |
输入的数据 |
(图五) |
4、编辑工作文件中已有的序列:可以按照操作Windows的习惯在工作文件主显示窗口选定一个或多个序列,点击鼠标右健打开一个或多个序列,进入数据编辑状态,可以修改数据。
四 、由组的观察查看组内序列的数据特征
(图六) |
按下数组窗口(也可以成为数组或数据编辑窗口)工具条上Views按钮,可以得到组内数据的特征,见图六。
具体介绍如下:
图七 |
图八 |
Group Members 可用于增加组中的序列;SpreadSheet以电子数据表的形式显示数据;Dated Data Table 将使时序数据以表的形式显示;Graph以各种图形的形式显示数据的;Multi Graph 以多图的形式显示组中数据;
Descriptive Stats给出组中数据的描述统计量,如均值、方差、偏度、峰度、J-B统计量(用于正态性检验)等;
Tests of equality…给出检验组中序列是否具有同方差、同均值或相同中位数的假设检验结果;N-way/One-way Tabulation…给出数组中序列观测值在某一区间的频数、频率和某一序列是否与组中其他序列独立的假设检验结果;Correlations 给出数组中序列的相关系数矩阵;Covariances给出数组中序列的斜方差矩阵;Correlogram (1)给出组内第1序列的水平序列及其差分序列的自相关函数和偏自相关函数;Cross Correlation (2)给出组内第1和第2序列的超前几期和滞后几期值之间的互相关函数; Cointegration Test执行 Johansen cointegration协整(或称为共积)检验; Granger Causality检验组内各个配对间的Granger因果关系;Lable给出数组的名称及修改时间等信息。
五、回归分析--估计消费函数
1、在经济理论指导下,利用软件包的“观察(View)”功能对数据进行“火力侦察”,观察消费性支出与可支配收入的散点图(见图七)。依据凯恩斯理论,设定理论模型:
consum= a + b (inc)
2、作普通最小二乘法估计:在主菜单选Quick\Estimate Equations,进入输入估计方程对话框, 输入待估计方程,选择估计方法—普通最小二乘法,如图八所示。点击OK进行估计,得到估计方程及其统计检验结果,如图九所示。
3、利用图九中给出的统计检验结果对模型的可靠性进行统计学检验,由统计结果可以看出该模型拟合优良,误差项不存在一阶正自相关。
4、利用图九中估计方程显示窗口中工具条View,可以显示估计方程、估计方程的统计结果、以图或表的形式显示数据的实际值、预测值和残差。
六、单方程预测
图九 |
预测是我们建立经济计量模型的目的之一,
其操作如下:进入方程估计输出窗口(可以选定一个已有的方程建打开或估计一个新方程)如图九,点击其工具栏中的Forecast打开对话框(图十),
输入序列名(Forecast name),这名称通常与方程中被解释变量的名字不同, 这样就不会混淆实际值和预测值;作为可选项,可给预测标准差随意命名[S.E(optional)],命名后,指定的序列将存储于工作文件中;用户可以根据需要选择预测区间(sample range for forecast);Dynamic 选项是利用滞后左手变量以前的预测只来计算当前样本区间的预测值,Static 选项是利用滞后左手变量的实际值来计算预测值(该选项只有在实际数值可以得到时使用),当方程中不含有滞后被解释变量或ARMA项时,这两种方法在第二步和以后各步都给出相同结果,当方程中含有滞后被解释变量或ARMA项时,这两种方法在第二步以后给出不同结果;用Output可选择用图形或数值来看预测值,或两者都用以及预测评价指标(平均绝对误差等)。将对话框的内容输入完毕,点击OK得到用户命名的预测值序列。
注意:在进行外推预测之前应给解释变量赋值。例如我们根据1980~1998年数据得到中国人均生活费支出与人均可支配收入关系的回归方程,希望预测1999、2000、2001年的人均生活费支出。为此,我们首先需要给出1999、2000、2001年人均收入可支配的数据,如果1999、2000、2001我们从历史数据中得不到1999、2000、2001年人均收入可支配的数据,就应利用其他方法估计出这些数据,把1999、2000、2001年人均收入可支配的数据(可能是估计值)输入解释变量中就可以预测出这三年的人均生活费支出。
图十 |
图十一 |
七、异方差检验
古典线性回归模型的一个重要假设是总体回归方程的随机扰动项ui同方差,即他们具有相同的方差s 2。如果随机扰动项的方差随观察值不同而异,即ui的方差为s i2, 就是异方差。检验异方差的步骤是先在同方差假定下估计回归方程,然后再对得到的的回归方程的残差进行假设检验,判断是否存在异方差。Eviews提供了怀特(White)的一般异方差检验功能。
零假设:原回归方程的误差同方差。
备择假设:原回归方程的误差异方差
我们仍利用表一数据进行分析。
操作步骤:在工作文件主显示窗口选定需要分析的回归方程 \ 打开估计方程及其统计检验结果输出窗口(见图九) \点击工具栏中的View \选Residual Tests \ White Heteroskedasticity (no cross terms) 或White Heteroskedasticity (cross terms)(图十一),可得到辅助回归方程和怀特检验统计量-即F统计量、 统计量的值及其对应的p值。由图十二中的显示结果可以看出:在1%显著水平下我们拒绝零假设,接受回归方程的误差项存在异方差的备择假设。值得重申的是:虽然图九中的信息告诉我们回归方程拟和优良,但我们还应该对其进行经济计量学检验,以确定其是否满足古典假设。
一般地,只要图十二中给出的p 值小于给定的显著水平,我们就可以在该显著水平下拒绝零假设。
图十二 |
注意:White Heteroskedasticity (no cross terms) 与White Heteroskedasticity (cross terms)选项的区别在于:在no cross terms选项下得到的辅助回归方程中不包含原回归方程左手变量的交叉乘积项作为解释变量;而cross terms选项下得到的辅助回归方程中包含原回归方程左手变量的交叉乘积项作为解释变量。在我们使用的一元回归例子中,这两个选项的作用没有区别。当我们分析多元回归模型的异方差问题时,因为所选辅助回归方程的解释变量不同,这两个选项的作用就不同了。
八、White 异方差校正功能和加权最小二乘法
1.White 异方差校正功能:我们使用表二的数据,在主菜单选Quick \Estimate Equations,进入输入估计方程对话框, 输入待估计方程 (cum in ),选择估计方法—普通最小二乘法,点击Options 按钮进入方程估计选择对话框,选择Heteroskedasticity Consistent Covariance \ White \ OK应用(见图十三)[1],回到估计方程对话框,点击OK得到校正后的回归方程(见图十四)。同学们可以比较图十四中的方程与普通最小二乘法得到的方程。
表二
中国1998年各地区城镇居民平均每人全年家庭可支配收入及交通和通讯支出
单位:人民币 元
项目 变量名
地 区 | 可 支 配 收入 | 交通和通讯支出 | 项目 变量名
地 区 | 可 支 配 收入 | 交通和通讯支出 |
|
| ||||||
in | cum | in | cum |
| ||
甘 肃 | 4009.61 | 159.60 | 新 疆 | 5000.79 | 212.30 |
|
山 西 | 4098.73 | 137.11 | 河 北 | 5084.64 | 270.09 |
|
宁 夏 | 4112.41 | 231.51 | 四 川 | 5127.08 | 212.46 |
|
吉 林 | 4206.64 | 172.65 | 山 东 | 5380.08 | 255.53 |
|
河 南 | 4219.42 | 193.65 | 广 西 | 5412.24 | 252.37 |
|
陕 西 | 4220.24 | 191.76 | 湖 南 | 5434.26 | 255.79 |
|
青 海 | 4240.13 | 197.04 | 重 庆 | 5466.57 | 337.83 |
|
江 西 | 4251.42 | 176.39 | 江 苏 | 6017.85 | 255.65 |
|
黑龙江 | 4268.50 | 185.78 | 云 南 | 6042.78 | 266.48 |
|
内蒙古 | 4353.02 | 206.91 | 福 建 | 6485.63 | 346.75 |
|
贵 州 | 4565.39 | 227.21 | 天 津 | 7110.54 | 258.56 |
|
辽 宁 | 4617.24 | 201.87 | 浙 江 | 7836.76 | 388.79 |
|
安 徽 | 4770.47 | 237.16 | 北 京 | 8471.98 | 369.54 |
|
湖 北 | 4826.36 | 214.37 | 上 海 | 8773.10 | 384.49 |
|
海 南 | 4852.87 | 265.98 | 广 东 | 8839.68 | 640.56 |
|
(数据来源:中国统计年鉴1998光盘,文件 j11c,j12c)
图十三 |
图十四 |
2、加权最小二乘法:我们使用表二的数据,在主菜单选Quick \Estimate Equations,进入输入估计方程对话框, 输入待估计方程 (cum in ),选择估计方法—普通最小二乘法,点击Options 按钮进入方程估计选择对话框,选择Weighted LS/TSLS \ 在对话框内输入用作加权的序列名称in的平方根得倒数 \ OK应用(见图十五),
回到估计方程对话框,点击OK得到加权最小二乘法回归方程(见图十六并与图十四中的方程比较)。
Eviews中进行加权最小二乘估计的过程为:选定一个与残差标准差的倒数成比例的序列作为权数,然后将权数序列除以该序列的均值进行标准化处理,将经过标准化处理的序列作为权数进行加权作最小二乘估计,这种做法不影响回归结果。但应该注意,Eviews的这种标准化处理过程对频率数据不适用。
图十五 |
图十六 |
九、一阶(高阶)序列相关校正
当线性回归模型中的随机扰动项是序列相关时,OLS估计量尽管是无偏的,但却不是有效的。当随机扰动项有一阶序列相关时,使用AR(1)可以获得有效估计量。其原理如下:
表三中的数据,设进口需求函数随机方程为
IMt= B0+ B1 GNPt+ ut (2)
IM为每年进口额, GNP每年收入的替代变量。假设误差项存在一阶自相关,则ut可以写成:
表三
我国进口支出与国内生产总值和消费者价格指数
| 国民生产总值(人民币亿元,当年价) | 进口总额(人民币亿元,当年价) | 消费价格指数(1985年=100) |
年度 | GNP | IM | CPI |
1985 | 8989.1 | 1257.8 | 100.0 |
1986 | 10201.4 | 1498.3 | 106.5 |
1987 | 11954.5 | 1614.2 | 114.3 |
1988 | 14922.3 | 2055.1 | 135.8 |
1989 | 16917.8 | 2199.9 | 160.2 |
1990 | 18598.4 | 2574.3 | 165.2 |
1991 | 21662.5 | 3398.7 | 170.8 |
1992 | 26651.9 | 4443.3 | 181.7 |
1993 | 34560.5 | 5986.2 | 208.4 |
1994 | 46670.0 | 9960.1 | 258.6 |
1995 | 57494.9 | 11048.1 | 302.9 |
1996 | 66850.5 | 11557.4 | 328.0 |
1997 | 73142.7 | 11806.5 | 337.2 |
1998 | 78017.8 | 11622.4 | 334.5 |
(数据来源::中国统计年鉴1999光盘c01、 q03和i01,)
ut = r ut-1+nt -1£r£1 (3)
其中n~N(0,s2), Cov(ni,nt) = 0, i¹j。记作ui服从AR(1)。
假定r已知,我们将方程(3)中的变量滞后一期,写为:
IMt-1= B0+ B1 GNPt-1+ ut-1 (4)
方程(4)两边同时乘以r得到:
rIMt-1=rB0+rB1GNPt-1+rut-1 (5)
将方程 (2)与方程 (4)相减并利用方程(3),得到:
IMt - rIMt-1=B0(1-r)+B1(GNPt-rGNPt-1)+nt (6)
图十七 |
图十八 |
图十八中AR(1)的系数就是r的估计值。Inverted AR Roots是残差自相关模型(3)的滞后算子多项式的根,这个根有时是虚数,但静态自回归模型的滞后算子多项式的根的模应该小于1。
如果模型(2)的误差项存在高阶自相关,形如
ut = r1 ut-1+r2 ut-2+r3 ut-3+nt -1£ri£1 i=1,2,3 (7)
我们应在图十七的估计方程对话框中输入IM C GNP AR(1) AR(2) AR(3)。如果模型(2)的误差项存在形如下式的自相关
ut = r1 ut-1+ r3 ut-3+nt -1£ri£1 i=1 ,3 (8)
我们应在图十七的估计方程对话框中输入IM C GNP AR(1) AR(3)。如果模型(2)的误差项存在形如下式的自相关
ut = r4 ut-4+nt -1£r4£1 (9)
我们应在图十七的估计方程对话框中输入IM C GNP AR(4)。这样就可以校正误差序列高阶自相关。
十、邹氏转折点检验
邹氏转折点检验的目的是检验在整个样本的各子样本中模型的系数是否相等。如果模型在不同的子样本中的系数不同,则说明该模型中存在着转折点。转折点出现的原因可能由于社会制度、经济政策的变化、社会动荡等,如固定汇率变为浮动汇率、中国的改革开放、战争等。我们可以用邹氏转折点检验来验证某点是否是转折点。这个检验使用的F是统计量和LR 统计量。
表 四
某地区1947年一季度至1957年4季度国内生产总值和投资总额数据 单位:亿美元
obs | GDP | INV | obs | GDP | INV | obs | GDP | INV |
1947-1 | 1239.5 | 43 | 1951-1 | 1504.1 | 60.4 | 1955-1 | 1742.5 | 64.7 |
1947-2 | 1247.2 | 42.3 | 1951-2 | 1548.3 | 65.4 | 1955-2 | 1758.6 | 67.9 |
1947-3 | 1255 | 43 | 1951-3 | 1585.4 | 61.7 | 1955-3 | 1778.2 | 70.8 |
1947-4 | 1269.5 | 49 | 1951-4 | 1596 | 57.3 | 1955-4 | 1793.9 | 74.2 |
1948-1 | 1284 | 49.8 | 1952-1 | 1607.7 | 58.9 | 1956-1 | 1787 | 74.2 |
1948-2 | 1295.7 | 51 | 1952-2 | 1612.1 | 51.1 | 1956-2 | 1798.5 | 73.9 |
1948-3 | 1303.8 | 51.4 | 1952-3 | 1621.9 | 52.8 | 1956-3 | 1802.2 | 75.7 |
1948-4 | 1316.4 | 49.8 | 1952-4 | 1657.8 | 55.8 | 1956-4 | 1826.6 | 76.1 |
1949-1 | 1305.3 | 43.1 | 1953-1 | 1687.3 | 56.5 | 1957-1 | 1836.4 | 77.5 |
1949-2 | 1302 | 35.6 | 1953-2 | 1695.3 | 56.2 | 1957-2 | 1834.8 | 76.8 |
1949-3 | 1312.6 | 37.8 | 1953-3 | 1687.9 | 56.1 | 1957-3 | 1851.2 | 78.5 |
1949-4 | 1301.9 | 34 | 1953-4 | 1671.2 | 51.3 | 1957-4 | 1830.5 | 68.8 |
1950-1 | 1350.9 | 43.4 | 1954-1 | 1660.8 | 51.1 |
|
|
|
1950-2 | 1393.5 | 48.6 | 1954-2 | 1658.4 | 51.6 |
|
|
|
1950-3 | 1445.2 | 53.5 | 1954-3 | 1677.7 | 54.7 |
|
|
|
1950-4 | 1484.5 | 63.9 | 1954-4 | 1698.3 | 58.8 |
|
|
|
根据表四数据建立回归方程如下:
GDP = 14.5169INV + 735.545
现在需要验证1952年4季度是不是转折点,即1952年4季度之前与之后投资对国内生产总值的贡献是否一致。操作如下:在方程估计输出窗口点击View/Stability test /Chow breakpoint test,如图十九;
进入转折点输入窗口如图二十,输入转折点日期;得到检验统计结果如图二十一。从统计结果可以看出F检验和LR 检验都拒绝零假设:1952年4季度是转折点,接受1952年4季度不是转折点的备则假设。一般地,只要图二十一中的显示的概率小于给定的显著平,如5%或1%,就可以在该显著水平拒绝原假设。
如果我们需要检验多个转折点,则可以同时输入多个转折点的时间。假如我们需要判断1952年4季度和1954年4季度是不是转折点,这时的零假设是:1952年4季度和1954年4季度都是转折点。可以验证我们拒绝零假设。如图二十二和图二十三。
图十九 |
图二十 |
图二十二 |
图二十三 |
十一、两阶段最小二乘法
在解联立方程组时,我们经常会用到两阶段最小二乘法,操作方法如下:
图二十四 |
在主菜单选Quick \Estimate Equations,选择估计方法—两阶段最小二乘法(TSLS),在估计方程对话框内, 输入待估计方程, 在工具变量窗口输入工具变量。如图二十四所示。
点击OK进行估计,就可得到估计方程及其统计检验结果。
¥29.8
¥9.9
¥59.8