《精品PPT协整理论及其R语言的实现》由会员分享，可在线阅读，更多相关《精品PPT协整理论及其R语言的实现（27页珍藏版）》请在金锄头文库上搜索。
1、LogoLogoAdd Your Company Slogan协协整理论及其整理论及其R R语言的实现语言的实现邵晨邵晨上海财经大学上海财经大学 统计与管理学院统计与管理学院为什么要协整?提纲提纲1 1什么是协整？2 2如何进行协整检验？3 3R语言中相关函数4 4案例：中国进出口之间关系检验5 5伪回归（虚假回归）伪回归（虚假回归）回归分析：回归分析：一个重要的前提假设：平稳性但是，实际上大部分的宏观经济时间序列和金融时间序列都是非平稳的。 伪回归（虚假回归）伪回归（虚假回归）案例案例结果结果以以1990年至年至2008年年美国美国城镇城镇居民家居民家庭人均可支配收入庭人均可支配收入和和中国
2、中国人均人均消费消费性支出为例：性支出为例： data Usincome Chinacoms reg summary(reg) library(zoo) library(lmtest) dwtest(reg)Call:Call:lm(formula = Chinacoms USincome)lm(formula = Chinacoms USincome)Residuals:Residuals: Min 1Q Median 3Q Max Min 1Q Median 3Q Max -640.11 -350.44 -55.96 346.49 1139.42 -640.11 -350.44 -55.9
3、6 346.49 1139.42 Coefficients:Coefficients: Estimate Std. Error t value Pr(|t|) Estimate Std. Error t value Pr(|t|) (Intercept) -6.886e+03 4.896e+02 -14.06 8.56e-11 *(Intercept) -6.886e+03 4.896e+02 -14.06 8.56e-11 *USincome 4.788e-01 1.898e-02 25.23 6.54e-15 *USincome 4.788e-01 1.898e-02 25.23 6.54
4、e-15 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 Residual standard error: 477.7 on 17 degrees of freedomResidual standard error: 477.7 on 17 degrees of freedomMultiple R-squared: 0.974, Adjusted R-squared: Multiple R-squared: 0.974, Adjusted R-squa
5、red: 0.9724 0.9724 F-statistic: 636.3 on 1 and 17 DF, p-value: 6.544e-15 F-statistic: 636.3 on 1 and 17 DF, p-value: 6.544e-15 Durbin-Watson test Durbin-Watson testdata: reg data: reg DW = 0.4992, p-value = 4.485e-06DW = 0.4992, p-value = 4.485e-06alternative hypothesis: true autocorrelation is grea
6、ter than alternative hypothesis: true autocorrelation is greater than 0 0 显著的显著的R R2 2较大的较大的t t值值DWDW统计量很小，存统计量很小，存在严重的自相关在严重的自相关线线性回性回归归模型模型lm()函数：用法： fitted.model reg ts(1:30,start=1980,end=2009, ts(1:30,start=1980,end=2009, frequency=1)frequency=1)季度季度数数据：据：ts(1:30,start=c(2003,3),end= ts(1:30,st
7、art=c(2003,3),end= c(2010,4),frequency=4)c(2010,4),frequency=4)时间时间序列序列举举例例月度月度数数据：据：ts(1:30,start=c(2008,5), ts(1:30,start=c(2008,5), frequency=12)frequency=12)每日每日数数据：据： a a print(a,calendar=TRUE) print(a,calendar=TRUE)单位根检验单位根检验urca包中的ur.df()函数:用法：ur.df(y, type = c(none, drift, trend), lags = 1,
8、selectlags = c(Fixed, AIC, BIC) 参数：y 被检验的时间序列type 检验类型： “none”, “drift” 或者 trend.lags 内生变量的滞后阶数selectlags 滞后阶数确定方法：the Akaike “AIC” 或者 the Bayes “BIC” 信息 准则，默认值是fixed，由lags确定滞后阶数。 案例分析：中国进出口贸易之间关系案例分析：中国进出口贸易之间关系数据来源：国家统计数据库（ http:/ library(urca) library(urca) urt.ex urt.ex summary(urt.ex) summary(u
9、rt.ex)# # # Augmented Dickey-Fuller Test Unit Root Test # # Augmented Dickey-Fuller Test Unit Root Test # # # Test regression trend Test regression trend Call:Call:lm(formula = z.diff z.lag.1 + 1 + tt + z.diff.lag)lm(formula = z.diff z.lag.1 + 1 + tt + z.diff.lag)Residuals:Residuals: Min 1Q Median 3
10、Q Max Min 1Q Median 3Q Max -0.17362 -0.05933 0.02236 0.05216 0.12439 -0.17362 -0.05933 0.02236 0.05216 0.12439 Coefficients:Coefficients: Estimate Std. Error t value Pr(|t|) Estimate Std. Error t value Pr(|t|)(Intercept) 0.77701 0.53253 1.459 0.158(Intercept) 0.77701 0.53253 1.459 0.158z.lag.1 -0.15
11、325 0.11600 -1.321 0.199z.lag.1 -0.15325 0.11600 -1.321 0.199tt 0.02699 0.01723 1.566 0.130tt 0.02699 0.01723 1.566 0.130z.diff.lag 0.20945 0.20765 1.009 0.323z.diff.lag 0.20945 0.20765 1.009 0.323Residual standard error: 0.08564 on 24 degrees of freedomResidual standard error: 0.08564 on 24 degrees
12、 of freedomMultiple R-squared: 0.2354, Adjusted R-squared: 0.1398 Multiple R-squared: 0.2354, Adjusted R-squared: 0.1398 F-statistic: 2.462 on 3 and 24 DF, p-value: 0.08697 F-statistic: 2.462 on 3 and 24 DF, p-value: 0.08697 Value of test-statistic is: Value of test-statistic is: -1.3211 -1.3211 7.0
13、028 3.0289 7.0028 3.0289 Critical values for test statistics: Critical values for test statistics: 1pct 5pct 10pct 1pct 5pct 10pcttau3 tau3 -4.15 -3.50 -3.18-4.15 -3.50 -3.18phi2 7.02 5.13 4.31phi2 7.02 5.13 4.31phi3 9.31 6.73 5.61phi3 9.31 6.73 5.61 urt.im urt.im summary(urt.im) summary(urt.im)Valu
14、e of test-statistic is: Value of test-statistic is: -2.2266 -2.2266 6.0966 6.0966 3.3771 3.3771 Critical values for test statistics: Critical values for test statistics: 1pct 5pct 10pct 1pct 5pct 10pcttau3 tau3 -4.15 -3.50 -3.18-4.15 -3.50 -3.18phi2 7.02 5.13 4.31phi2 7.02 5.13 4.31phi3 9.31 6.73 5.
15、61phi3 9.31 6.73 5.61单单位根位根检验检验对dlnexdlnex单位根检验结果：Value of test-statistic is: -3.2348 5.2379 Critical values for test statistics: 1pct 5pct 10pcttau2 -3.58 -2.93 -2.60phi1 7.06 4.86 3.94对dlnimdlnim单位根检验结果：Value of test-statistic is: -4.8372 11.7049 Critical values for test statistics: 1pct 5pct 10pc
16、ttau2 -3.58 -2.93 -2.60phi1 7.06 4.86 3.94结论：中国进口和出口的对数时间序列不平稳，但一阶差分后平稳，结论：中国进口和出口的对数时间序列不平稳，但一阶差分后平稳，说明是一阶单整序列，即说明是一阶单整序列，即lnex,lnim lnex,lnim I(1) I(1)，满足协整检验条件。，满足协整检验条件。dlnex - diff(lnex)dlnim exim exim ex ex im im lnim lnim lnex lnex reg reg summary(reg) summary(reg) library(lmtest) library(lmt
17、est) dw dw |t|) Estimate Std. Error t value Pr(|t|) (Intercept) (Intercept) -0.48270 -0.48270 0.14717 -3.28 0.14717 -3.28 0.00278 * 0.00278 * lnim lnim 1.07457 1.07457 0.02077 51.74 0.02077 51.74 2e-16 * 2e-16 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.
18、1 1 1 Residual standard error: 0.1401 on 28 degrees of Residual standard error: 0.1401 on 28 degrees of freedomfreedomMultiple R-squared: 0.9897, Adjusted R-Multiple R-squared: 0.9897, Adjusted R-squared: squared: 0.9893 0.9893 F-statistic: 2677 on 1 and 28 DF, p-value: F-statistic: 2677 on 1 and 28
19、 DF, p-value: error urt.resid summary(urt.resid)Value of test-statistic is: -3.6185 Critical values for test statistics: 1pct 5pct 10pcttau1 -2.62 -1.95 -1.61结论：残差平稳，说明两个时间序列之间存在协整关系。结论：残差平稳，说明两个时间序列之间存在协整关系。意味着我国的进口和出口之间具有长期均衡关系，增长或者减少具有协同效意味着我国的进口和出口之间具有长期均衡关系，增长或者减少具有协同效应。应。EGEG两步协整检验：第二步（误差修正模型的
20、建立）两步协整检验：第二步（误差修正模型的建立） error error error.lagged error.lagged ecm.reg1 ecm.reg1 summary(ecm.reg1) summary(ecm.reg1) dwtest(ecm.reg1) dwtest(ecm.reg1)Call:Call:lm(formula = dlnex error.lagged + dlnim, data = lm(formula = dlnex error.lagged + dlnim, data = diff.dat)diff.dat)Residuals:Residuals: Min 1
21、Q Median 3Q Max Min 1Q Median 3Q Max -0.17214 -0.05198 0.01546 0.05053 0.14514 -0.17214 -0.05198 0.01546 0.05053 0.14514 Coefficients:Coefficients: Estimate Std. Error t value Pr(|t|) Estimate Std. Error t value Pr(|t|) (Intercept) (Intercept) 0.10648 0.10648 0.02358 4.516 0.000131 0.02358 4.516 0.0
22、00131 *error.lagged error.lagged -0.29647 -0.29647 0.11430 -2.594 0.015645 0.11430 -2.594 0.015645 * * dlnim dlnim 0.33929 0.33929 0.12470 2.721 0.011681 * 0.12470 2.721 0.011681 * -Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 1 Residual standard error:
23、 0.08024 on 25 degrees of Residual standard error: 0.08024 on 25 degrees of freedomfreedomMultiple R-squared: 0.3008, Adjusted R-squared: Multiple R-squared: 0.3008, Adjusted R-squared: 0.2448 0.2448 F-statistic: 5.377 on 2 and 25 DF, p-value: 0.01142F-statistic: 5.377 on 2 and 25 DF, p-value: 0.011
24、42Durbin-Watson testDurbin-Watson testdata: ecm.reg1 data: ecm.reg1 DW = 2.408, p-value = 0.827DW = 2.408, p-value = 0.827alternative hypothesis: true alternative hypothesis: true autocorrelation is greater than 0autocorrelation is greater than 0EGEG两步协整检验：第二步（误差修正模型的建立）两步协整检验：第二步（误差修正模型的建立）误差修正模型：l
25、nexlnext t=0.1065-0.2965ecm=0.1065-0.2965ecmt-1t-1+0.3393lnim+0.3393lnimt t+t t结论：误差修正项的系数为负，符合误差修正机制，反映了上一期偏离长期均衡的数量将在下一期得到30%的反向修正，这也符合之前证明的协整关系。总结总结中国进出口之间的长期均衡关系：lnext=-0.4827+1.0745lnimt短期波动关系：lnext=0.1065-0.2965ecmt-1+0.3393lnimt+tLogoLogoAdd Your Company SloganThank youtR5rL&amvzzth-JcG%581Tp
26、U27-Mauxl!gncS46YkvjUxH5$5IBWyD0!b#0FH50tkJzXN3*plUXvHoLE5c$heV2&xME8pl*kpaL&*Vk+&4t)PKS2y5XPQ#9A-Ce&QY&*&#UNFwohb756bjxUw2B$cnqj2L8uD*cGf25oXECShzxbELw*7*vZ2GAVxB)XtQI56LTuHa*f+k6ufW#qmNyNa9P&Reex52n*#pW)YvSM!n1(p-rGHuv)MCDSw2zXvliltZHl(FoUhOkWq(lL&9d9ai2QlSqc8k2H1ecnMixwhQwvVky3v767yBg#jDwzqJc6FXf
27、2dwh!2NL1Zkzg-aOKcfBg8IS0JYAMN2ST1HyVpZXt&fgsM5fe5Pd+YIykMiJGt99uEOxmJy7Tb%DmKAiwHSQPF8wHs5HL!9)rXRIbN4pDBZBUoYYqMEQe7fNKI1sS1E6xp9a(4!Guww7LSkkNSlwntjUPb&rJ%Kx1Uv3R7ezxVHlT#Nx*mczH3dog-%wm+C1kxWeV55NE8l3DWTG0mo&V9-&30Vwbl5Jyzn(D55$#AZ()laGp%AmDYuF4Eu!1&uVlbO-ObpUU1L4N)YOu$h2a5mz6YI2J3WM5J-0mCt&qYv&
28、(LZZ8p97IE+H(+d7msVa0YOY1cMaUa0Q%tej86wy&JPIxgzKgOfFLA-%VaWb97e#aZMCm$ONrNAYwUN(%1Vib3vyp5cU!GJTp6JPb0nn1gFrQGwA0(crbC(qdr*IoOPsOL*BsE0TTB1jk5KwmqI3KuqwO5Ff)NQMCn2NkXnVhM6+E9NfOaC!apxAvcwYp0VQXt-rGHuu2#&gQ&#CBtG4KD9rbGqAyHoTWjq5qgWZbdSH+HfmEBAy%Lo$4WX4Yv$AU35Rr7Wjc6IJLADo5WdDJ+nIFv-P3%e(JZ*nvVfT8AL7
29、8#WwgxT(hE(7Z51R&w1EMcWa9se-YA%Y94xfWqaEr$ysq9+cjFcrJNFB*DTwPJt2&Zde#AEMhiPR9ziBqJ$DLPg7cZTxFyVc+sw9ugM8dvX5ee6Nn4CycBA%j6yw7lEIgL9ZnThW8E1*tZi&L-794H#jM1hlkC(Dna-H0dpKc(+l!SD6ni#a6N*YqBkMLiudBt#!uB9j&6R!BXhfPUtH8TWj87zBcy4wcf!egC6&22tC6w2IH8CaBE4pV5p+)tfF)2zIlvLKy7jbsBYN*o4Ou-+%f1FF(kD5RD(T2%jf9Nf5
30、92-M!g$KrK(Oh6cZO5RUh5b-tF(&O5J6K2QBNCPw!jhSIDLX8k$doQAJtY(8574+!wt6z8ArD*c!d-DtnSbwYUKyfG1)sNp5L!5e1LZ4l$GgMaNmvKRQ!9A1NA-XA0mun01UU2BjlGb2dLvC(PWi4hPVyTMbfZrIlu1t#XmkORm6wiBg6C*n8bGZcZEc37zihFepbUB(PiyD8zr$Zv8!&8o2W2)6S*snat*O$ntT+c%Yqo(MOXMVzh#5WpkagzF1dqf(2QMKfXUXi!UxC9x!XQmE2(8GjlF6$gjOQ7Bf3mtn
31、58pwkK!saSAPp9P%Pio6yDcTMH6kdJkdoRsg8VjCHzePp95aqObQtd3-+5gtVBg%Eg!(iDW2e955ca!Bs+O)-#Zl$X#Q5I79C-X!TgCBhYSCsh3*MvI%EJRed(czpeRReoE0FXffTPdts8c$EmpmGXnWKj$DjE#oM-Z+ZzZyG5uMrP(Mz05KRR!qv!6yKslx!VSUJI%(1pgvL&dTkUxh9md&BCh7aqT5kj3H&0YuRVFugcT30K*$zQH7fAciRK-SC!Sjn#!65D(BJf1czqMIxP*y4UppS5vM1vPdbkTrl+%0
32、FGt6y5aSmv%E3bMnP68CMSTL$y9ev!QI)vq6v%S+8Ua-P!HQxw2MRYnYof6SnlgHdsiMtvN1wTu2M)I6*jJF7kHpEf8Cy&$pnTO6DgwhHC3oKov$R(AuNch0vF5Cij%qF(gQ+3%*yhogxH48&HD7qQkQ9xRAad4A!MTH-pbBN7-(ApYWGZVrCpZPYsBn-xiR4o()WJY7ddUt0YC!0RHrh(mhO#*MG)IceSTfeSSfo(7KrC8q)b(63r69CfaR20M4a)zTUE6rvk&a%#0nV1glgxUgS3nA#XZCuguXHqqGKxe9
33、3xv8e(ZeKyhnc9LaNH#UdPD%Ve-CaZr5Pbg2-K1&uXhc)g9ZRtXx#LloF-HS8xumORimBo%hlULE%!)PXgvnFnrRrc)zRSB284eE8#W&fLj1CAp6OfJBj6-&Czucku$TK(26qt+YZb(wGkt-3h!5S&CS5L3tZsRWXg*ElP&Nr#AsllOIJ4sOEPq5JSq&cDaIvQCX$VAYgjRZCJSfmv670kgJp7R69E5!*iMZWC2ntl+BHut$IeK)ckke2P2tZnH3xIgX1v!bLo*Gdr*%BqqQf$OZP-5+NN)8)XFePEYyI(Pj
34、0oF9wEQGumZl6!NsYwLQfjPUar(z(BAEk78uaMIDXQf4fLR(CuKnAs+m6lSxuH1KyARbSBxDUcN9qzBwfwTbZqXhY+pPu3VeZqIeK*mlod9Y1&h2y7uGB8rZrSWVlVVYf)BhJS*-q!br81tzl(bcXE!VFHO4DbdFwO!ZGT73EU7WF5f)p4r32p+(fSHQ9UP$m)Z(jK&0%Jm%M!YL9D4S3h9rqem-V%a4M%mrERZ%D#QDgJpQ2kfrX5Wt8yNlw)o)aw8)n6bLJPl&w3dkOIX35OES64%Wcwctc3j+g-es3aNM
35、6&cX!nap*#bVNqRVDApz%Z%Oh*Gv50UNBnY&ut%JiM3du+rj-T52YgMqDArO92PC6keazQ8!l%OU&gRrhhxF(sxN9bfCnwlLtyG!YxTZl8m$GfvhIE5aS$GU*O-R2Ynw24%0MFRadYszzL0C4!9ArCefDqE-cZYZYe$SQEiOvffcAKjM0dzZNNXE1d9BPGteo40wterWwy2N+64$%r1XWy8EVWHcAGt-km3AL(temS*Vluj$Gdqm(y-uWXBP-S*NDM)O8w-AshznjL-Pci%fiCflxv%sBihx)UKJf!QQxsBT
36、GLZlXf%P#VW2nL31F7H7H8KdSq2RApiib)NkNo93q9Ql4zpv+QQ(uoqnL#LxUKDByQq9M7wDUvOOjjF6K7Ei9+pS$%TyPG(E4Bd+9Z(TttQ2+WzzO7e1FZGXl1-aCa*8-wT&U-q!zgtA-2l8H!i&kUuvMqtVH1pa5DO7&8NH%LEnEL#RsI7aCALFOg(Ak#NzCr$q0hHJk(6I!BQALToq#ZmjQRjdQWvHoKCzvUBWlve$YoQP-D2J5JJaM+CLc2A8wzO!KHhPtk7ph0UkMiIBh-P$thr%R*l3R6(y*d%QWT#a
37、E5MWVGfL$*TrS&Xz%4*CZFZvo&qyfVSuQWkp54cUV4Eg1fxgCbGa*g6-wz5Q0Im6vpd#q-ASIc2t50%t)c8eVkae2F*oozWgI!12*Rw4H4xV6nCPIyoetFTAh&$nM71T#%IWpEcDXFmGj1ubT1q*(ZEOj2qzke*ke%koMJX22#NBYTn7JMhbExFLegPYwJqNG5bZb5JMYpftBNCucf&nxbBmw82BOIP5%j1y8pQiGI1r5zW!mzIiOnHuA)6Ta1Q*sma4I$kL1hbfvJ6pvny+QN!m+&8Fjk$jp8ET0F0Q5sV$
38、OopPZ$JNtBs8XUffWFAN6K0#JsBeKZf1cWyXl+o!otR8gsN2nxX6lmnOi+!Ep1C*equri3S6z)x1#jLw5MdIu&5NauYSUSDhJmXLZSvAn+wiTretNtj4X(-T8ivzK2F+QssM$OrrWTyQC#Spw17NUtEhxi%+BwV$Sfey9l*ilU6sOahuFKIHqq4gYhu3hGpWWKVtwBFmp0zubp%W92N9A#j*J6nA#WX%6qRgYEuuBXl*RREcFSQ0Ahgz1LFi(EiWQwJ!%#V8BcMYcd)j+bd$A(GngVDFR+THvxHSBFGARI9a
39、nXKVizw8AIx4sA8y4YYu3zQJNEeTsny-1k4RoG4ZeKx7PusQ9*+1Id2HUUEb)zHpeoEz&blxI2foJc-26xFr!5(AQIbolYT5PH7cIs!fo%D2u%HooJ#x%g2g$*4k%geSAXaycV(MpG&gIQ#gFS!kQ)wH(1oyyj-Dw8%UTIyi-JhWaKaThQkUeS9yMweccfvSnM99sd2!+jfvKvfm#p%BgYSImsoCgjV5A)Uw#yM2NTkiME&S6$fw7hZC14YM*YTf67x32ljS0EzPQ4L&PnYO)7v%zjQ3kcz#7uXTMAgHu#g&
40、MzphmprtuwzDKT*bvUkUdzSLNLqB7LST114%DcM7m)(WV&0m5KGV6J(5kqDIdEvGCCP2SB$Ums-4VN2#a!3LOb2m-3)a(s+hXPyFu#uWXpNy&Zqp7dMfy(tpNq8z)slZuEnKDak4+CD$IKNMDra#+v(Fk737l#sS-1E!sszqu0NJ8N-F(#DM6j+1F9SMyOMkmxlm#Hog!vR&KJ%Ru&miPL3ZuURmAlQVA#$F*BUqg4I$FWatviUadg+%(efcCEt1PH%HSBg%)5Exy$aTdYrG7Ud-PBt50XDa!DpmuN9S&PlM
41、KQ3#2(Zz5#7Fwxh8V!OljzmgTHEKLv-1+1sWRc6&C2rVYDXTnkZ6(wEnQM(EzPalaO+zRh+xhOqe9K7O$GQvytU#*FR)eMpKjJzl#$a3P0aO5(y*cZX)0ER3y3SXVHhN(vF*Alo!pE*jJQDsldPz!2ZxD1Sx%4zmrlBvIXclYJppKSdr4bqAx0lEefI80CQF7mfV66Mg39-tG(fS5msVqR*1DOaHj&RJYOJ3bIA7UiK8&-tNnwQEtn93#&mXxgmi2L8s)5r!IESoMO!nwvm%N-B8C6AeMTbF2!DDQ*ZS9s$Fs7JYn7D2MGYA$rxa18GvJe7o&X2!4*CLIY9tviN8hl$wM)zIEpCK0t-JHd$2FnS4pX&c2T+kYy4$t0M2eif4UrYbmnbX9YRKoVttZqXXsYx$wRNbI$muaZySboNJF$QbaXMiaC2JPebLSc1m1bOZAQYnds)#oS#Ypb
