时间序列分析VAR模型实验

基于VAR模型的我国房地产市场与汇率波动的因果关系————VAR模型实验2第一部分实验分析目的及方法现选取人民币对美元汇率以及商品房房价作为变量构建VAR模型。对于不满足单位根检验的序列采取对数化或差分处理，使其成为平稳序列再进行模型的拟合。对于商品房房价这一变量，由于全国各省市差异较大，故此处采用全国房地产开发业综合景气指数这一变量。此外，为了消除春节假期不固定因素带来的影响，增强数据的可比性，按照国家统计制度，从2012年起，不单独对1月份统计数据进行调查，1-2月份数据一起调查，一起发布。所以国房景气指数p这一序列缺少每年一月份的相关数据，属于非随机、不可忽略缺失，在此采用平均值填充的方法，补足数据。第二部分实验样本2.1数据来源数据来源于中经网统计数据库。具体数据见附录表。2.2所选数据变量由于我国于2005年7月实行第二次汇改，此次汇改以市场供求为基础、参考一篮子货币进行调节、有管理的浮动汇率制度取代了过去人民币汇率长达10年的紧盯美元的固定汇率体制。故本实验拟选取2005年07月到2014年10月我国以月为单位的数据。，用以上两个变量来构建VAR模型，并利用该模型进行分析预测。第四部分模型构建4.1判断序列的平稳性4.1.1汇率E序列首先绘制出E的折线图，结果如下图：3图4.1汇率E的曲线图从图中可以看出，汇率E序列较强的趋势性，由此可以初步判断该序列是非平稳的。为了减少m的变动趋势以及异方差性，先对m进行对数化处理，记为lm，其时序图如下：图4.2lm的曲线图4对数化后的趋势性减弱，但仍存在一定的趋势性，下面对lm进行一阶差分处理，去除趋势性，得到新变量dlm，观察dlm的曲线图。图4.3DLE的曲线图从图中可以看出，dle序列的趋势性基本已经消除，且新变量dle基本围绕0上下波动，因此选择形式为yt=yt-1+ut进行单位根检验：表4.1单位根输出结果NullHypothesis:DLEhasaunitrootExogenous:ConstantLagLength:2(Automatic-basedonSIC,maxlag=12)t-StatisticProb.*AugmentedDickey-Fullerteststatistic-3.0316730.0351Testcriticalvalues:1%level-3.4919285%level-2.88841110%level-2.581176*MacKinnon(1996)one-sidedp-values.AugmentedDickey-FullerTestEquationDependentVariable:D(DLE)Method:LeastSquaresDate:11/15/14Time:20:20Sample(adjusted):2005M112014M10Includedobservations:108afteradjustments5VariableCoefficientStd.Errort-StatisticProb.DLE(-1)-0.3530050.116439-3.0316730.0031D(DLE(-1))-0.5027300.115417-4.3557680.0000D(DLE(-2))-0.3115310.093265-3.3402580.0012C-0.0008880.000470-1.8875920.0619R-squared0.450240Meandependentvar1.15E-05AdjustedR-squared0.434382S.D.dependentvar0.005058S.E.ofregression0.003804Akaikeinfocriterion-8.269046Sumsquaredresid0.001505Schwarzcriterion-8.169708Loglikelihood450.5285Hannan-Quinncriter.-8.228768F-statistic28.39119Durbin-Watsonstat2.061613Prob(F-statistic)0.000000单位根统计量ADF=-3.031673小于临界值，且P为0.0351，因此该序列不是单位根过程，即该序列是平稳序列。4.1.2国房景气指数P序列首先作出P序列的时序图：图4.4P的曲线图由于每年一月份的数据缺失，故取相邻两项进行平均补全数据，得到新序列的时序图如下：6图4.5P的曲线图(补全)由上图可知，该序列P可能存在一定的趋势性和季节性，先进行单位根检验，确定改序列是否平稳。由于序列表4.2单位根输出结果NullHypothesis:PhasaunitrootExogenous:Constant,LinearTrendLagLength:3(Automatic-basedonSIC,maxlag=12)t-StatisticProb.*AugmentedDickey-Fullerteststatistic-3.9724570.0124Testcriticalvalues:1%level-4.0452365%level-3.45195910%level-3.151440*MacKinnon(1996)one-sidedp-values.由单位根检验结果可知，T值小于临界值，且P=0.0124，在5%的置信水平下，该序列不存在单位根过程。由于汇率E序列为一阶单整序列，并进行了一阶差分处理，因此样本数量减少，在下面的操作中，所有的样本序列调整为2005-08至2014-10。4.2模型参数识别7先进行VAR模型的拟合，初步选定滞后阶数为3：表4.3拟合输出结果VectorAutoregressionEstimatesDate:11/22/14Time:22:20Sample(adjusted):2005M112014M10Includedobservations:108afteradjustmentsStandarderrorsin()&t-statisticsin[]DLEPDLE(-1)0.063183-19.12274(0.09626)(14.1374)[0.65638][-1.35263]DLE(-2)0.11679815.42129(0.09604)(14.1052)[1.21615][1.09330]DLE(-3)0.24526016.39171(0.09617)(14.1243)[2.55030][1.16053]P(-1)-9.04E-051.490708(0.00066)(0.09765)[-0.13593][15.2656]P(-2)-0.000583-0.355442(0.00118)(0.17380)[-0.49226][-2.04508]P(-3)0.000346-0.160740(0.00067)(0.09872)[0.51479][-1.62821]C0.0313282.571540(0.01274)(1.87084)[2.45943][1.37454]R-squared0.2950330.979509Adj.R-squared0.2531540.978292Sumsq.resids0.00139029.99247S.E.equation0.0037100.544936F-statistic7.044848804.67678Loglikelihood454.8094-84.06138AkaikeAIC-8.2927661.686322SchwarzSC-8.1189241.860164Meandependent-0.002527100.2406S.D.dependent0.0042933.698585Determinantresidcovariance(dofadj.)4.08E-06Determinantresidcovariance3.57E-06Loglikelihood370.8871Akaikeinformationcriterion-6.609021Schwarzcriterion-6.261337再进行滞后阶数的确定：表4.4最优滞后阶数的判断VARLagOrderSelectionCriteriaEndogenousvariables:DLEPExogenousvariables:CDate:11/22/14Time:22:22Sample:2005M072014M10Includedobservations:99LagLogLLRFPEAICSCHQ0134.7784NA0.000234-2.682392-2.629965-2.6611801302.5627325.39998.57e-06-5.991165-5.833886-5.9275302329.023050.247835.45e-06-6.444909-6.182775*-6.338849*3334.37339.943949*5.30e-06*-6.472187*-6.105200-6.3237044337.45315.5997425.40e-06-6.453598-5.981758-6.2626915339.75894.0991765.60e-06-6.419372-5.842679-6.1860416345.04989.1923245.46e-06-6.445451-5.763905-6.1696967345.54840.8460765.87e-06-6.374716-5.588316-6.0565378346.73691.9687606.23e-06-6.317917-5.426663-5.9573149352.58019.4436396.01e-06-6.355154-5.359047-5.95212810353.77141.8770826.39e-06-6.298411-5.197451-5.85296111354.36490.9112796.87e-06-6.229594-5.023780-5.74172012356.46173.1346447.18e-06-6.191146-4.880479-5.660848*indicateslagorderselectedbythecriterionLR:sequentialmodifiedLRteststatistic(eachtestat5%level)FPE:FinalpredictionerrorAIC:AkaikeinformationcriterionSC:Schwarzinformationcriterion9HQ:Hannan-Quinninformationcriterion由上边可知，根据信息准则，采取少数服从多数原则，取滞后阶数为3，此外取滞后阶数为2(SC为-6.182775)或取滞后阶数为3(SC为-6.105200)时，两者SC值相差不是很大。3.3模型参数估计选取了最优滞后阶数3，进行模型的拟合。拟合结果如下：表4.5VAR(3)模型估计结果VectorAutoregressionEstimatesDate:11/22/14Time:22:23Sample(adjusted):2005M112014M10Includedobservations:108afteradjustmentsStandarderrorsin()&t-statisticsin[]DLEPDLE(-1)0.063183-19.12274(0.09626)(14.1374)[0.65638][-1.35263]DLE(-2)0.11679815.42129(0.09604)(14.1052)[1.21615][1.09330]DLE(-3)0.24526016.39171(0.09617)(14.1243)[2.55030][1.16053]P(-1)-9.04E-051.490708(0.00066)(0.09765)[-0.13593][15.26