计量经济学课程实验方法与案例汇编.docx

1、计量经济学课程实验方法与案例汇编计量经济学：试验方法与案例汇编（2011年03）实验一 一元线性回归模型分析某市城市居民年人均鲜蛋需求量（Y：公斤），年人均可支配收入（X：元)的例子。通过抽样调查，得到19881998年的样本观测值.其数值结果如下：解：年份1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 合计年人均收入x847.26 820。99 884。21 903.66 984.09 1035.26 1200.90 1289.77 1432.93 1538.97 1663。63 12601.67 人均需求量Y14.4 14.4

2、14。4 14。7 17.0 16。3 18。0 18。5 18。2 19.3 17。1 182。3 X2 15325548。8y2 3055。25xy 213349.96X2=15325548。8 Y2=3055。25 XY=213349.96 n=11 X=12601.67 Y=182。30，离差形式：解：（1)计算一元线性回归模型参数估计量及建立模型则直线回归方程为：(请注意：如果样本数据较大,也可以采用离差形式的公式计算出参数估计量）即:结果与用基本公式完全相等.(2)样本决定系数（或可决系数）r2=(y2-e2）/y2=1-(e2/y2)=1-(11.89/35。04）=10。339

3、3=0。6607其中，残差平方和:e2=y2- bxy=35。040。00514539.63=11。89（3）参数估计量的方差、标准差S2（a)= (e2X2)/n（n2） x2 =(11。8915325548.8/11*9*888903。8）=2.0736S2（b)=(e2)/ （n-2）x2=11.89/9888903.8=0。00000148则S(a)=1。44，S(b）=0。0012(4） 参数估计量的显著性检验T（a）=a/S(a)=10.72/1。44=7.44，T（b)=b/S（b)=0。0051/0.0012=4。25得知T（a)T/2（n2)=2。26， T（b)T/2（n2

4、）=2.26即a、b均显著不为零，说明解释变量X对Y有显著影响.（5）方程显著性检验F检验：F=（y2-e2)/（ e2/n-2)=（35.0411。89）/(11.89/112）=17.523查表得F临界值：F（1，n2）= F0。05（1，9)=5。12，故FF,通过F检验,说明Y与X线性显著。实验二 多元线性回归模型（为了便于同学们理解,现以二元线性回归模型为例）模型形式：1、参数估计量公式：2、随机误差项方差估计量3、残差平方和练习题:1、若回归模型为Y=a+b X+,并且已经根据X和Y的16组样本数据，计算出下列各值：X2=5089.84 Y2=1111。01 XY=2367.19

5、n=11 X=240。78 Y=114.22求:（1）模型参数的估计量；(2)总离差平方和、残差平方和；（3）样本决定系数； （4）参数估计量的标准差。2、已知Y和X之间存在线性因果关系,并已经获得如下数据年份XYX2Y2XY1990199119921993199419951996合计344449586776854132528343640424625111561936240133644489577672252634762578411561296160017642116934185012321666208826803192391015618请利用一元线性回归模型进行分析（参数估计、显著性检验、)

6、；预测当X=90时,Y的预测值.3、自己动手查找某个国家、地区或地方的收入和消费数据，模仿本章的例题,进行线性回归分析（参数估计、模型检验及解释参数经济意义)，并对结果进行讨论。4、令Y表示一名妇女生育孩子的生育率，X表示该妇女接受过教育的年数.生育率对教育年数的简单回归模型为：Y= b0+b1X问：（1）随机误差项包含哪些未引入模型中的因素，它们可能与教育水平相关吗？(2）上述简单回归模型能够准确揭示教育对生育率在其他条件不变下的影响吗？请说明。5、对于人均存款(S）与人均收入（I）之间的关系S=b0+b1I+，使用美国32年的年度数据，得到如下估计模型(括号内为标准差） S384。1050

7、.067I（151.105） （0.011） R20。538试回答:（1）b1的经济意义是什么？（2）b0和b1的符号应该怎样的？为什么？实际的符号与你的直觉是一样吗？如果有冲突的话,你可以说明原因吗？（3)你对拟合良度有什么看法?（4）对回归得到的两个参数估计量进行检验。（显著性水平为0.01）教材P98（6)6、某地区家庭人均鸡肉年消费量（Y)，家庭月平均收入X（元)，肌肉价格P1(元/kg），猪肉价格P2(元/kg)的统计资料年份YXP1P219852.783974.225.0719862.994133。815。219872。984394.035。419883.084573。955。53

8、19893。124923。735.4719903。335283。816。3719913。565603。936。9819923。646243.786。5919933.676663.846.4519943.847174.017.0019954.047683.867.3219964。038433。986。7819974。189113.977。9119984。049315.219。5419994.0710214。899。4220004.0111655.8312.3520014。2713495.7912.9920024。4114495.6711。7620034。6715756。3713.0920045.

9、0617596.1612。9820055。0119915.8912。8020065.1722586。6414。10要求：利用eviews软件,建立多元线性回归模型,并进行检验（R2，t检验和F检验）。 实验三 异方差性检验与处理一、异方差的处理方法如果预测模型通过检验证实随机项存在异方差，则可以采用两种方法对模型参数进行校正，也就是重新对其进行估计.而处理异方差,其关键问题是找出异方差的具体形式，为了简单起见，我们仅以一元线性回归模型为例，对于异方差形式提出两种假定。 一元线性回归模型为假定1：随机误差项的方差与解释变量Xi的平方成正比.即：，式中，2为常数对一元线性回归模型进行转换：用解释变

10、量Xi分别去除一元模型两端各项得：式中，=证明转换后的新模型具有等方差性.需要强调的是原模型中的常数项成了转换后模型中的系数，原模型中的系数成了转换后模型的常数项.对经过转换后的新模型应用普通最小二乘法（即OLS法)，即可得到对的回归方程：参数估计值得到，则可以建立预测模型假定2：随机误差项的方差与解释变量Xi成正比,即,式中，2为常数在此，克服异方差的方法类似于假定1，但原模型两端不是用xi去除了，而改为用去除两端各项。则经过转换后的模型：式中,同样，转换以后的新模型也克服了异方差，具有等方差性。对模型采用OLS法回归,得到因变量为，自变量为和的二元回归方程：，则原模型回归参数得到，即可写出

11、预测模型方程：二、举例：异方差判断、处理在研究某地区居民储蓄倾向时，得到了如下数据资料。判断用线性回归模型研究居民储蓄倾向时，误差项是否存在异方差，并给出处理方法.个人收入和储蓄数据n储蓄Y收入Xn储蓄Y收入X123456789101112131415162641059013112210740650343158889895077981912221702877792109954105081097911912127471349914269155221673017663185751963521163228801718192021222324252627282930311578165414001829

12、22002017210516002250242025701720190021002300241272560426500276702830027430295602815032100325003525033500360003620038200第一步,进入软件操作系统,并创建工作文件即：FILE(文件)/NEW(新文件)/workfile(工作文件），然后根据对话框提示的数据类型和样本范围进行选择，再加以确定(OK）；第二步，输入数据.可以通过菜单或命令“data+变量名”，输入上述表中数据;第三步，作散点图。点中两组数据名称，并用graph命令或菜单操作(quick/qraph/scatter/w

13、ith regression line）根据两个变量间数据分布情况，可以认为两者之间存在明显的线性关系,但从数据的分布也可以看出异方差存在的可能性。（可通过回归残差e序列分布情况来判断模型随机项存在递增异方差）回归可以得到直线回归模型：LS C Y X下面，进一步采用GQ检验法进行确认异方差。因为X已经按照大小顺序存在，可以删除中间c=n/4=31/4=7个数据,将前后各12组数据分别构成两个子样本并进行回归,得到两个回归直线模型根据两个回归模型的残差平方和计算F统计量为：F=769899.2/144711.5=5。32给定显著性水平为0.05,查两个自由度均为12-11=10的F分布表的临界

14、值F0.05（10，10)=2.98.因为F F0.05(10，10），说明原模型存在明显的异方差。由于根据残差序列图可以看出，回归残差的绝对值有随着X线性增长的趋势，因此考虑对一元线性回归模型进行转换：用解释变量Xi分别去除一元模型两端各项得：然后再用最小二乘法，以Y/X为被解释变量,以1/X为解释变量进行参数估计,可得回归直线模型：生成两个新序列:genr y1=Y/Xgenr x1=1/X经进一步检验，模型结果得到了明显改善。第四章 异方差性：软件实现（重要归纳内容)本章的重点内容总结:异方差检验及修正一、异方差的检验方法一：图示法步骤：首先，对被解释变量和解释变量进行回归：LS Y C

15、 X 其次，在主菜单中,点击Quick/Graph 得到对话框“Series List”,在对话框中现列出横坐标X, 同时列出纵坐标Y，点击“OK”, 出现“Graph”对话框,再点击“Show Options出现对话框“Graph Options ，在左上角图形类型Graph Type 下拉列表里选择“Scatter Diagram,在右下角散点图项目下选择“Regression Line，再点击“OK”,即可得到X 、Y之间的散点图。 这是随机误差项存在异方差的初步经验证据。此外,同理也可以做出解释变量X与残差e2之间的散点图。方法二：G-Q检验法（详细内容参见教科书及本章案例）注意：该方

16、法适用于样本容量较大，且异方差为单调递增或单调递减的情况。另外，G-Q检验需要按某一被认为有可能引起异方差的解释变量观测值的大小排序，因此，可能需要对各个解释变量进行轮流试验。构造一个GQ检验的F统计量与F临界值的实现方法：操作如下: 以截面数据n=20为例(去掉中间4个观察值，即c=4）Smpl 1 8LS Y C XScalar sig1=SESmpl 13 20LS Y C XScalar sig2=SE构造G-Q统计量建立检验结果的储存向量： COEF(2）GQ第一个放置F统计量: GQ(1）=(sig22/sig12）第二个放置F临界值: GQ（2)=QFDIST（.95,6，6)以

17、上操作步骤结束后，在当前工作簿中便出现gq标志,双击此标志即可显示出F统计量和F临界值的计算结果。方法三：怀特（White）检验优点：怀特检验异方差不需要对解释变量进行排序，且对任何形式的异方差都适用。1。首先,对二元线性回归模型作普通最小二乘回归，并得到残差平方（e2）序列2。其次，作如下辅助回归方程：若 nR22（h)（注意: h为模型中卡方分布的解释变量个数,上述辅助回归模型中h=5)，则表明二元线性回归模型存在异方差性。注：在多元回归中，由于辅助回归方程中可能有太多的解释变量（x1 x2 x12 x22 x1 x2） ,从而使自由度减少, 有时可以去掉交叉项。二、异方差的修正1。加权最

18、小二乘法（Weighted Least Squares）基本思想：在采用最小二乘法时，对较小的残差平方（e2）赋予较大的权数，而较大的残差平方赋予较小的权数，对残差所提供的信息的重要程度作一番修正，以提高参数估计的精确度.注意：在同方差假定条件下,加权最小二乘法与普通最小二乘法是一致的,即模型不存在异方差情况下， 加权最小二乘法与普通最小二乘法两者是通用的.实现步骤:点击主菜单中的Quick/Estimate Equation，在Equation Specification对话框中，点击右下角的Option键，得到对话框，选中Weighted LS/Tsls前面的单选框，输入加权变量字母（残差

19、绝对值的倒数）,再点击“OK”,便回到方程对话框，再点击 “OK”即可得到加权最小二乘法的结果。点击“OK”得到回归结果课堂练习题异方差：案例分析（31个地区)（P61)省份税收Y国内生产总值X北京1435。79353。3天津438。45050.4河北618.313709.5山西430.55733。4内蒙古347。96091.1辽宁815.711023。5吉林234。75284.7黑龙江335。07065.0上海1975.512188.9江苏1894.825741。2浙江1535.418780。4安徽401。97364。2福建594。09249.1江西281.95500。3山东1308.425

20、965.9河南625.015012.5湖北434.09230.7湖南410.79200.0广东2415.531084。4广西282.75955.7海南88。01223.3重庆294。54122。5四川629。010505.3贵州211.92741.9云南378。64741.3西藏11。7342.2陕西355。55465。8甘肃142。12702.4青海43.3783.6宁夏58。8889。2新疆220。63523.2上表是2007年我国各省份全年税收Y与国内生产总值（GDP）X的数据。要求：(1)作出散点图，检验模型是否存在异方差，说明是单调递增还是单调递减异方差；（2)采用G-Q方法检验模型

21、是否存在异方差?若存在异方差，采用加权最小二乘法(WLS)进行模型估计（即修正异方差） （3）采用怀特（White）检验法检验模型是否存在异方差性（4）对修正异方差后所建立的方程进行检验（经济学检验、统计学检验：拟合良度、变量显著性检验）（5）解释方程斜率的经济含义。实验四 自相关性检验与处理一、广义差分法：基本步骤：1、建立模型：LS Y C X2、检验自相关：德宾沃森检验（DW检验)3、若存在自相关，则采取广义差分法进行修正：操作步骤：(1)进行广义差分变换（见P134之（6。7）式），生成两个新序列,软件操作：Genr Y2=Y-*Y（-1）,ok(按回车键）Genr X2=X-*X(-

22、1)，ok（按回车键）(2）以上述两个新生成的序列为变量,建立新的回归模型（见P134之公式（6.8）求得b0,b1参数估计值,其中， ,（3）对新模型再次进行自相关检验，方法同第2步二、Durbin(杜宾）两步法第一,估计自相关系数，见P135其中，前面的回归系数恰好为第二，再用的估计值进行广义差分（同广义差分法）自相关举例：下表资料为某国消费（C）和可支配收入（I）(单位：百万美元）年份消费（C）可支配收入（I）198719881989199019911992199319941995199619971998113781301215263168731776418857200742143922

23、833242052530727020116171329715790180171931420198215122312424724261752721928915（1)试求其残差,并利用图示法检验随机项是否存在自相关；（2)若存在自相关,应用迭代法对原模型自相关进行处理。解:第一步，求出模型估计值的残差e，并检验随机项自相关（利用图示法)数据录入后，在命令行输入LS Y C X，运行显示结果。 （4。003） (77。657) R2=0.9983 d=0.933第二步,生成Y的估计值Y和残差值 e两个新数列:genr y1=984.6433+0。8892Xgenr e=Y-Y1第三步，绘制出残差数列

24、趋势图,可以看出随机项呈现正自相关。以上所做的三个步骤是为了检验模型随机项自相关问题,接下来是对自相关进行克服或处理。采用迭代法对原模型进行处理首先输入命令：LS Y C X AR(1）可以得到(即AR（1)对应的参数值）或者采用残差值e 与其滞后值e（-1）进行回归，也可以得到然后，采取广义差分法生成两个新数列，即，genr Y2=Y-*Y(1）genr X2=X*X（-1)LS y2 c x2得到新模型：Y2=A+ b1X2其中，b0=A/1-， b1=b1*则克服自相关以后得到模型：Y=b0+b1X练习与思考题1、若一个含有三个解释变量的多元线性回归模型,用样本容量为100的数据回归得到

25、DW值为1.05，能得出什么结论（即是否可以表明存在一阶自相关）？若DW值等于2。35， 结论又是怎样?2、上述资料为日本1971-1990年20年间税收Y与国民生产总值X数据。P141（6）：年份税收（Y)国民生产总值（X)197127.9181。9197231.6198.3197336.6207。7197436.0207。3197532.1215.6197634。6224。3197736.4235.0197842.0247.1197945.1260.6198048.5268.8198150.9277.4198253。2287.2198355。9295。8198458.9309.119856

26、2。2324.0198666。2333。3198773。6349。8198880。4370.6198984.9387。5199090。0407。2（1）依据上述资料，对回归模型Y=+X进行OLS估计，并读出t值、R2以及DW.（2）检验所估计的模型是否存在自相关性。 实验五 多重共线性检验与处理一、多重共线性举例:某地区19912000年服装消费量、可支配收入、流动资产、服装价格、相关商品价格等资料，试分析变量之间的多重共线性。年份服装消费(Y）可支配收入(X1)流动资产（X2)服装价格(X3)相关商品价格（X4)19918.482.917。1929419929.688。021。3939619

27、9310。499.925。19697199411。4105.329。09497199512.2117.734.0100100199614.2131。040。0101101199715.8148.244。0105104199817.9161。849。0112109199919。3174.251。0112111200020。8184.753.0112111利用Eviews 软件来完成整个分析过程.第一步，进入软件操作系统，并创建工作文件即：FILE（文件）/NEW(新文件)/workfile（工作文件），然后根据对话框提示的数据类型和样本范围进行选择，再加以确定（OK）；第二步,输入数据。可以通过菜单或命令“data+变量名,输入上述表中数据；第三步，作散点图。点中两组数据名称，并用graph命令或菜单操作（quick/qraph/scatter/with regression line)第四步，在建立了工作文