序列相關(guān)性的檢驗(yàn)與修正

1、序列相關(guān)性的檢驗(yàn)與修正案例：書本P115進(jìn)口與國(guó)內(nèi)生產(chǎn)總值的關(guān)系。一檢驗(yàn)準(zhǔn)備工作：建立工作文件，導(dǎo)入數(shù)據(jù)。采用OLS方法建立進(jìn)口方程。在命令框輸入：equationEqOl.lsmcgdpDependentVariable:MMethod:LeastSquaresDate:11/17/10Time:23:00Sample:19782001Includedobservations:24CoefficienlStd.Errort-StatiSticProb.c152.905716.07S4S3.31S37600031GDP0J2Q3-9+Q.00101420.116S1O.OQQOR.-squar

2、ed0.943440臨朗dependentvarS26.9&42AdjustedR-squared0.946096S.D.dependsntvar667.4365S.E.ofregression154.9600Akaikeinfocriterion13.00J37Sumsquaredresid52S2774Schwarzcriterion13.10204Loglikelihood-154.0464Hannan-Quinncriter.13.02991F-statistic404.6860Durbin-Watsonstat0.627922Prob(F-statistic0.000000建立殘差序

3、列在命令框輸入：seriese=resid建立殘差序列的滯后一期序列在命令框輸入：serieseag1=resid(-1)方法1：利用兩個(gè)殘差序列畫圖、觀察。方法2：查看回歸方程的DW值=0.628,存在序列相關(guān)。方法3：LM檢驗(yàn)在命令框輸入：equationEq02.lsecgdpe(-1)e(-2)在命令框輸入：Scalarlm1=obs(e(-2)*eq02.r2可得LMl=15.006在命令框輸入：scalarchi1=qchisq(0.95,2)可得chi1=5.99可以判定模型存在2階序列相關(guān)。簡(jiǎn)便方法：在方程eq01窗口中點(diǎn)擊View/ResidualTest/SeriesCor

4、relationLMTest,并選擇滯后期為2,則會(huì)得到如下圖所示的信息。注:LM計(jì)算結(jié)果與上面有差異，因?yàn)檫@里的輔助回歸所采用的resid(-1)、resid(-2)的缺失值用0補(bǔ)齊。Breusch-GodfteySerialCorrelationLMTestF-statislic19.52905Prob.FC2.20)0.0000ObsR-squared15.37241Prob.Chi-Square(20.0004TestEqua:ion:DependentVariable:RE3IDMetnoa：LeastsquaresDate.11/1B/10Time:00.04Sample:197S

5、2001Includedobsrvatiang：24Presamplemissingvalulaggedreidualset1ozero.CoefficientStd.Errort-StalisticProbC6.59475028561860.2303940.3197GDP-0.00D3440.300633-0.5040090.6190RESID(-1)1.095930.1755246.2304590.0000RESID(-2)-07857760.212016-3.692270O.OOUR-squared0.0G1350Meandependentvar-1.52E-13AdjustedR-sq

6、uared0.61D553S.D.dlependlentvar151.5539S.E.ofregression94.57026AkaikeinfDcriterion12.00774Sumsquaredresid179900.9Schwaizcritarion12.29409Loglikslihood-1J1.D529H甘nnan-Quinncriter.12.13983F-statis1ic12.01937urbin-Watsanstat1.S73142Prob(F-statistic)0.000061檢驗(yàn)是否存在更高階的序列相關(guān)。繼續(xù)在命令框輸入:equationEq03.lsecgdpe(

7、-1)e(-2)e(-3)在命令框輸入：Scalarlm2=obs(e(-3)*eq03.r2可得LM2=14.58在命令框輸入：scalarchi2=qchisq(0.95,3)可得chi2=7.185仍然存在序列相關(guān)性，但由于e(-3)的參數(shù)不顯著，可認(rèn)為不存在3階序列相關(guān)。在方程eq01窗口中點(diǎn)擊View/ResidualTest/SeriesCorrelationLMTest，并選擇滯后期為3,則會(huì)得到如下圖所示的信息。Breusch-GodfreySerialCorrelationLMTest:F-statistic1237576Prob.F(3,19)0.0001Obs*R-squ

8、ared15.87561Prob.Chi-Square(30.0012TestEquation:DependentVariable:RESIDMethod:LeastSquaresDate:11/13/10Time:00:11Sample:19782001Includedobservations:24Presamplemissingvaluelaggedresidualssettozero.CoefficientStd.Errort-StatisticProb.C6.69163829.319490.22S2320.3219GDP-0.0003490.000703-0.4967170.6251R

9、ESID(-1)1.1078380.24S95545411510.0002RESID(-2)-0.8192950444735-1.3422100.0811RESID(-3)0.0322970.3733510.0065070.9320R-squared0.661404Meandepe-1.52E-13AdjustedR-squared0.590217S.D.dependentvar151.5539S.E.ofregression97.01614Akaikeinfocriterion12.17068Sumisquaredresid178830.5Schwarzcriterion12.4-1611L

10、oglikelihood-141.0482Hannan-Quinncriter.12.23500F-statistic9.281822urbin-Watsonstat1.380605Prob(F-statistic0.000247顯然，LM檢驗(yàn)的結(jié)果拒絕原假設(shè)(無序列相關(guān))，表明存在序列相關(guān)性。二序列相關(guān)性的修正與補(bǔ)救廣義差分法就是廣義最小二乘法(GLS)，但損失了部分樣本觀測(cè)值，損失的數(shù)量依賴于序列相關(guān)性的階數(shù)(如一階序列相關(guān)，至少損失1個(gè)樣本值)。在實(shí)際操作中，往往基于廣義差分法完成估計(jì)。根據(jù)隨機(jī)擾動(dòng)項(xiàng)相關(guān)系數(shù)估計(jì)方法的不同,可以分為C-0迭代法和Durbin兩步法。(1)Durbin兩步

11、法第一步，估計(jì)隨機(jī)擾動(dòng)項(xiàng)的相關(guān)系數(shù)根據(jù)前面檢驗(yàn)可知存在二階序列相關(guān)，因此設(shè)定方程為M*+pM+pM+*GDP+*GDP+*GDP,t01t12t21t2t13t2t在命令框輸入lsmcm(-1)m(-2)gdpgdp(-1)gdp(-2)即得到隨機(jī)擾動(dòng)項(xiàng)相關(guān)系數(shù)的估計(jì)值，結(jié)果記為eq04DependentVariable:MMethod:LeastSquaresate:11/22/10Time:17:03Sample(adjusled):19802001Includedobservations:22alteradjustmentsCoefficientStdErrort-StaiisticPr

12、ob.C78.0857544.4625617562130.09S20.932490.14125066424660.0000-0.46S6520.265736-1.76326S0.0969GDP0.0547770.0093115.SQ2762ci.oooaGDP(-1)-0.0963950.013567-5.1913770.0001GDP(-2)0.0535830.0101065.3020790.0001R-squared0.991335Meandependentvar890.0591AdjustedR-squared0.98B627S.D.dependentvar661.6499S.E.ofr

13、egression70.56198Akaikeinfocriterion11.57736Sumsquaredreaid79663.88Schwarzcriterion11.37542Loglivelihood-121.3565Hannan-Quinncriter.11.64796F-statistic366.0372Djrbin-Watsonstat2.306915Prob(F-statistic)0.000000第二步，進(jìn)行差分變換，然后對(duì)應(yīng)書本（4.2.20）公式進(jìn)行回歸。M*二M(0.938M0.469M)ttt1t2GDP*=GDP(0.983GDP0.469GDP)ttt1t2在命令

14、框中輸入Seriesm_star=m-(0.938*m(-1)-0.469*m(-2)Seriesgdp_star=gdp-(0.938*gdp(-1)-0.469*gdp(-2)將M*對(duì)GDP*回歸。在命令框輸入ttLsm_starcgdp_star結(jié)果記為eq05第三步，還原（參考書本4.2.20式）*1垐-p1278.161-0.9380.469，147.19這一結(jié)果與eq01有差別。(2)C-O迭代法原理：參見教材113頁(yè)。操作：非常簡(jiǎn)單。根據(jù)前面檢驗(yàn)可知存在二階序列相關(guān)，因此設(shè)定方程為M，GDPpAR(1)pAR(2)t01t12t在命令框中輸入Equationeq06.lsmcgd

15、par(1)ar(2)DependentVariable:LIMethod:LeastSquaresate:11/22/10Time:17:29Sample沮djusted:19302001Includedobservations:22afteradjustmentsConvergenceachievedafter5iterationsCoefficientStd.Errort-StatisticProb.c169.321044.390073.8143090.0013GDP0.0197920.00107318.45-2500.0000AR(1)1.1031770.1-31246&1142270

16、.0000AR(2)-0.3011940.221892-3.6107360.0020R-squared0.932325Meandependsntvmr890.0591AdjustedR-squared0.979379S.D.dependentvar661.6499S.E.ofregression95.01304Akaikzeinfocriterion12.10387Sumsqua.redresid162494.6Schwarzcriterion12.30724Loglikelihood-129.1976Hannan-Quinncriter.12.15560F-statistic333.4596

17、urbin-Watsonstat1.353364Prob(F-statistic0.000000InvertedARRoots.55+.70i.55.7Oi從DW值來看，已經(jīng)不存在序列相關(guān)性，并且每項(xiàng)AR的回歸系數(shù)都具有統(tǒng)計(jì)顯著性。假設(shè)采用1階廣義差分，在命令框中輸入Equationeq07.1smcgdpar(1)DependentVariable:MMethod:LeastSquaresDate:11/22/10Time:T7:37Sample(adjusted):19792Q01Includedobservations:23afteradjustmentsConvergenceachievedafte1iterationsCoefficientStd.Errort-StatisticProb.C163.&43斗142.31671.1491560.2640GDP0.0214200.0025080.0000ARC1)072-80930.1B12864.0162B00.0007R-squared0.970150Meandependentvar858.1739AdjustedR.-squared0.967165S.D.dependentvar664.2776S.E.ofregression120.699Aaikeinfocriterion12.&4Q12Sums