试验七GARCH模型在金融数据中的应用

1、实验七（G）ARCH模型在金融数据中的应用一、实验目的理解自回归异方差（ARCH）模型的概念及建立的必要性和适用的场合。了解（G）ARCH模型的各种不同类型，如GARCH-M模型（GARCHinmean）,EGARCH模型（ExponentialGARCH）和TARCH模型（又称GJR）。掌握对（G）ARCH模型的识别、估计及如何运用Eviews软件在实证研究中实现。二、基本概念p阶自回归条件异方程ARCH（p）模型，其定义由均值方程（7.1）和条件方程方程（7.2）给出：yt=Pxt+&（7.1）222ht=var（St|Q1）=a0+a1St+a2stn+ap&（7.2）其中，表示t-1时

2、刻所有可得信息的集合，ht为条件方差。方程（7.2）表示误差项备的方差ht由两部分组成：一个常数项和前p个时刻关于变化量的信息，用前p个时刻的残差平方表示（ARCH）。广义自回归条件异方差GARCHp,q）模型可表示为：yt=Bxt+&（7.3）22ht=var（&|Q。=a0+a1stA+.+ap&.十九+一十%“（7.4）三、实验内容及要求1、实验内容：以上证指数和深证成份指数为研究对象，选取1997年1月2日2002年12月31日共6年每个交易日上证指数和深证成份指数的收盘价为样本，完成以下实验步骤：（一）沪深股市收益率的波动性研究（二）股市收益波动非对称性的研究（三）沪深股市波动溢出效

3、应的研究2、实验要求：（1）深刻理解本章的概念；（2）对实验步骤中提出的问题进行思考；（3）熟练掌握实验的操作步骤，并得到有关结果。四、实验指导（一）沪深股市收益率的波动性研究1、描述性统计（1）导入数据，建立工作组打开Eviews软件，选择File单中的NewWorkfile选项，在Workfilefrequency”框中选择undatedorirregular,在Startobservation和aEndobservation框中分另U输入1和1444,单击OK。选择File菜单中的Import-ReadText-Lotus-Excel”选项，找到要导入的名为EX6.4.xls的Excel

4、文档完成数据导入。(2)生成收益率的数据列在Eviews窗口主菜单栏下的命令窗口中键入如下命令：genrrh=log(sh/sh(-1),回车后即形成沪市收益率的数据序列rh,同样的方法可得深市收益数剧序列rz。(3)观察收益率的描述性统计量双击选取rh”数据序列，在新出现的窗口中点击ViewaDescriptiveStatistics”HistogramandStats”,则可得沪市收益率rh的描述性统计量，如图71所示：图71沪市收益率rh的描述性统计量同样的步骤可得深市收益率rz的描述性统计量。观察这些数据，我们可以发现：样本期内沪市收益率均值为0.027%,标准差为1.63%,偏度为-

5、0.146,左偏峰度为9.07,远高于正态分布的峰度值3,说明收益率rt具有尖峰和厚尾特征。JB正态性检验也证实了这点，统计量为2232,说明在极小水平下，收益率rt显著异于正态分布；深市收益率均值为-0.012%,标准差为1.80%,偏度为-0.027,左偏峰度为8.172,收益率rt同样具有尖峰、厚尾特征。深市收益率的标准差大于沪市，说明深圳股市的波动更大。2、平稳性检验再次双击选取rh序列，点击ViewUnitRootTest，出现如图72所示窗口:UnitRaotTesl.TestType&AugmentedDidciayFuierPhillips-PerronTestfenunitr

6、ootin:Leveldifference2nddifferenceIncludeintestequation:InterceptTrendandinterceptNoneLaggeddifferences:OKXcel图7-2单位根检验对该序列进行ADF单位根检验，选择?t后4阶，带截距项而无趋势项，所以采用窗口的默认选项，得到如图73所示结果：图7-3rhADF检验结果同样对rz做单位根检验后，得到如图7-4所示结果:图74rzADF检验结果在1%的显著水平下，两市的收益率rt都拒绝随机游走的假设，说明是平稳的时间序列数据。这个结果与国外学者对发达成熟市场波动性的研究一致：Pagan(19

7、96)和Bollerslev(1994)指出：金融资产的价格一般是非平稳的，经常有一个单位根(随机游走)，而收益率序列通常是平稳的。3、均值方程的确定及残差序列自相关检验通过对收益率的自相关检验，我们发现两市的收益率都与其滞后15阶存在显著的自相关，因此对两市收益率rt的均值方程都采用如下形式：rt=c+an*&(7.5)(1)对收益率做自回归在Eviws主菜单中选择QuickEstimationEquation，出现如图75所示窗口:EquatiodSpecificatianEquaboriSpecificiondependentvariablefolowedbyE$tcfregressor

8、sincludingARMAandPDLternn*.URan剪pliciequationI火eV=c1Hc2fXrhe而叫EstimationSettings:如t府日：|L5-LeastSquares(NLSandARMA)3Sample:11443一图7-5对收益率rh做自回归在Method”中选择LS(即普通最小二乘法)，然后在Estimationsettings上方空白处输入图75所示变量，单击“OK,则出现图7-6所示结果：打四1”专-CEfctiou：VITITLEDVotkfil*:OTITLED)13ilsEditQbjtetsiewtrocsquickOptions.Win

9、dow151X|ViewjFr”15|jbjects|Print|Mame|Fme|Estimate|Fmcs三11St.tw|Residw|DependentVariable:RHMethod:LeastSquaresDate10727/05Time16:24SampIe(adjusted):171443includedobservations:1427afteradjustingendpointsVariableCoefficientStd.Errort-StatisticProb.c0.0002260.000431U5247E30.990RH卜0.0946970.02630113.600

10、5070.0003R-squared009015MeandependentygrD.Q0D253AdjustedR-squared0.003320S.D.dependentvar0.016346S.E.ofregression0015279Akaikeinfocriterion-5.396425Sumsquaredresid0.377654Schwarzcriterior-5.309048Loglikelihood3852.349F-statistic1296355Durbin-Watsonstat2024254ProbfF-statistic)0.000326IMFtth=c：vifrws3

11、DBnoiktV7:=antitied图7-6收益率rh回归结果(2)用Ljung-BoxQ统计量对均值方程拟和后的残差及残差平方做自相关检验：点击ViewResidualTestCorrelogram-Q-statistics”,选择10阶滞后,则可得沪市收益率rh残差项的自相关系数acf值和pacf值，如图77所示：图7-7沪市收益率rh残差项的自相关系数acf值和pacf值点击ViewResidualTestCorrelogramSquaredResiduals”,选择10阶滞后，则可得沪市收益率rh残差平方的自相关系数acf值和pacf值，如图78所示：图7-8沪市收益率rh残差平方的

12、自相关系数acf值和pacf值采用同样的方法，可得深市收益率rz的回归方程及残差、残差平方的acf值和pacf值。结果表明两市的残差不存在显著的自相关，而残差平方有显著的自相关。(3)对残差平方做线性图。对rh进行回归后在命令栏输入命令：genrres1=residA2,得到rh残差平方序列resl,用同样的方法得到rz残差平方序列res2。双击选取序列resl，在新出现的窗口中选择ViewLineGraph,得到resl的线性图如图7-9所示图7-9rh残差平方线状图同理得到rz残差平方线状图图7-10rz残差平方线状图可见&2的波动具有明显的时间可变性(timevarying)和集簇性(c

13、lustering),适合用GARCH类模型来建模。(4)对残差进行ARCH-LMTest依照步骤(1),再对rh做一次滞后15阶的回归，在出现的“Equation”窗口中点击“View”“ResidualTest”“ARCHLMTest”，选择一阶滞后，得到如图711所示结果：股OTTTLEDTorkfil*：VWTITLBJ1JFileEdit0bjectsViewProcs审aiuk0E.tionsWindowHelpViewFt白esi口卜此已七里|0七|Nam档|Ff电电工2|/5士上七口|Frueeg5t慰.士5|Re.ids|ARCHTest：F-Statislic444soi2

14、ProbabilityO.OODOOOObs*R-squared43.19040Probability0.000000TestEquation口叩end即tVariable:RESIDA2Method：LeasiSquaresDate:10/27/05Time:17:19Sanrtpiofadjusted)IB1443Includedobservations:1426afteradjustingendpointsVariableCoefficientStd.Errort-StatisticProb.C0.0002192.09E-0510,432640.0000RESlDt-U0174043D2

15、6O96666934200R-squaired0.030290M”ndependentvar0.000265AdjustedR-squared0.029609S.D.dependentvar0.000755S.E.ofregression0000743Akaikeinfocriterion11.56967Sunnsquaredresid0.000787Schwarzcriterion-11.56229Loglikelihood8251170F-Statistic44.46012Durbin-Watsonstat2.039221Prob(F-statistic)o.ooooooFathsc：DB

16、=xunw|VF=untitltd图7-11rhARCH-LMTest对rz方程回归后的残差项同样可做ARCH-LMTest,结果表明残差中ARCH应是很显著的。4、GARCHI模型建模(1) GARCH(1,1)模型估计结果点击“Quick”aEstimateEquation”,在出现的窗口中Method选项选择ARCH,可以得到如图712所示的对话框。在这个对话框中要求用户输入建立GARCH类模型相关的参数：“MeanEquationSpecification栏需要填入均值方差的形式；“ARCH-Mterm”栏需要选择ARCH-M项的形式，包括方差、标准差和不采用三种；“ARCHSpeci

17、fication”栏需要选择ARCH和GARCH项的阶数，以及估计方法包括GARCH、TARCH和EGARCH等等；“VarianceRegressors”栏需要填如结构方差的形式，由于Eviews默认条件方差方程中包含常数项，因此在此栏中不必要填入“C”。我们现在要用GARCH(1,1)模型建模，以沪市为例，只需要在“MeanEquationSpecification”栏输入均值方差“RHCRH(-15)”，其他选择默认即可，得到如图713和图714所示的结果。图712EquationSpecification窗口置JEViet,-EqasGQiL：QSTITLEDWTITLED-JOIX|

18、口EiltEditQbjctsfie*ftuickOttioftsIindawHelpiew工口史写IClbjwcit|FrintFarecaEt|ftats|ResidsDependentVariable;RHMethod,ML-ARCHDate10/27/05Time17:29Sample(adjusted)：171443Includedobservations1427afteradjuslingendpointsConvergenceachievedafter25iterationsCoefficientSidErrorz-StatisticProbc-2.21E-OB0.000203-

19、B.0077950.9938RH(-16)0.0693020.020200293574500033VarianceEquationC872E-061.45E-066.0035430.0000ARCH(1)0.176514001524211,59078口.000。GARCH(1)08073910.01370356920G000000Rsquared0.007543Meandependentvar0.000253AdjustedR-squared0.004751S.D.dependentvar0.016348S.E.ofregression0.016309Akaikeinfocriterion-6

20、.645039Sumsquaredresid0.379216Schwarzcriterion5.626597Loglikelihood4032.735F-statistic27D1026Durbin-Watsonstat2.023588Prob(F-statistic)0.029225MFath=Xtvisws3IiB=rmm胖=uMirUd会EYiis-Equation.UITITUDTorkfilt.ITHTITLEDI寸FileEditDbjfistsViwPrcsgiukOs.tioiliXWindowHalp.|国|巽View|FrocslObjFrintiHamelFrre工uEs

21、timate|Forets5tl5lti|HrsiAeDependentVariable:RZMethod:ML.ARCHDate:10/27AJ5Time:17:28Sample(adjusted):171443includedobservations.1427afteradjustingendpoinlsConvergenceachievedafter29iterationsCodhciernSidErrorz-StatisticProbc-UOOO57SO.ODD327-1.77(3677a.O7G6RZ(-15)0.0575420.0219142.S2576BO.OOSSVarianc

22、eEquationc6.32E-061.29E-064.906562o.oooaARCH(1)1221S30.009763124G695,0000GARCHnBE3BMO.OD923193.56S23a.ooooRsquared0,005931Meandependentwar-0.000173AdjustedR-squaredU.M3135SD.depsndeiitvar0.017933S.E.ofregression0.0179KAkwikeinfocrrterioni-5.47E727Sumsquaredresid0456386Schwarzcriterion6.456285Loglike

23、lihood3912.646F-statiistic2.120990Durbin-Watsonstat1.91436Prob(F-statistic)0.075965二|IEPtth=c：i：views3DB=nDnWF=untitl：UHTITUDOFile旦查LObjectsViewProcsauickOfitifnsWindowHel*5DependantVairiabh:RZMethod:ML-ARCHDale：10/2705Time：17:36.Sample(adjusted):171443Includedobservatioins:1427afteradjustingendpoin

24、tsConvergenceachievedafter29iterationsCoefficintStd,Errorz-StaiiisticPrab.c-,07840.000355-2.20681302730.056272002224J2.52973D0114VananceEquationc5.74E-06123F-064B71746ooonARCH(1)U.Q96SD20.0D97519.927374O.OOOU(RESIDJJiARCHfl)0.M7B69,0143873.3272760.0009GARCH(1)0.967877n.00075499,140710000R-squared.OS

25、2D5Meandependentvar-0.000173AdjustedR-squared0.001705S.D.dependentvar0.017933S.Eofregression.17S10AKaikeinfocriterion-5.47F775Sumsquaredresidl,456219Schwarzent&rion-5.456645Loglikelihood3915.1D6F-statistiic1.406995Path=c：Sifivier3DB=monoWF-untitlsd图7-18深市收益率TARCH(1,1)模型估计结果在TARCH中，qdt二项的系数估计值都大于0,而且

26、都是显著的。这说明沪深股市中坏消息引起的波动比同等大小的好消息引起的波动要大，沪深股市都存在杠杆效应。2、EARCH模型估计结果rz在图7-12的“ARCHSpecification”下拉列表中选择“EGARCH则可得到rh、的EGARCH莫型估计结果，分别如下图7-19和图7-20所示。lEVievs-TEquation.UTTITLEDTorkfil*：VHTLTLZDI11FileEdlitOLj*ct-EVioas童uiukOjtionsWindowHlpVi”|FrtiCEOhectwPrintHstim昼tc|Fcirecest|StatwKisidm|DependentVaria

27、ble:RHMethod:ML，ARCHDale:10/27JO5Time.17:43Sample(adjusted).171443Includedobservations.1427afteradjuMingendpoint;Convergenceachievedafter57iterationsCoefficientStd.Erraiz-StaiisticPn)b.cRH(15)-0.000276054325o.oonsQj0.020970*,9080022.5906690.36390.0096VarianceEquationCIRES1/5QRGARCH1RES/SQRGARCH(1)EG

28、ARCH-0.512770 279B14-0D51846 9639030057778002221G0.011154O005960-0.07406712,59365-4.64311516171940ooooOOODOoooonO0000R-squaredlAdjustedR-squaredS.E.ofregressionSumsquaredresidlLoglikelihood.0(E3900.Q0Q8930.01632403786554050.076MeandependentvarS.D.dependentvarAkaikeimfbcriterionSchwarzcriterionF-stat

29、istic.002530.01634B5.6679415.B45ei11.B27619二1Path二c:tevide:；nolleWF-untitZ-edEViews-EquatiQii：UBTITLEBTorkfile：VBTITLED-1x!EiLeRdi1QbjadroctQuickOjtiOXikgild。*Ulp_|ffjx|Yi*|Pro=|口b.ct.|FrintKtmt|FrtztEi+ft+lFwyt|StHiR”idi|DependentVariable:RZ一Method:ML-ARCHDate:10/27/05Time:17:41SamplB(adjusted)-171

30、443includedobseivations】427afteradjustingendpomtsConwergencanotachievedafter100iterationsCoefficientStd.Errorz-StatisticProb.C-0.000775000340-22776340227RZ(-15)0D49J52002237322103260.0271VarianceEquationC-0.3713630.045840-3.1012130.0000PEc|/S?RC-APLHfl23664501587114.92301,OQOORES/SQR(GARCH(1)-003206

31、90DO0Q9S3.60277S00003EGARCH0.97B4B60.005116190.38190.0000R-squared0.004B09Meandependentvar-0.000173AdjustedR-squared0.001307S.D.dependentvar0.017933S.Eofregression017922Akaikeinfocriterion与487199Sumsquaredresid0456401Schwancnierion-5,65060Loglikelihood3921116Fatalistic1.373207Path=c：eviews3DB=noneWF

32、=untitled图7-20深市收益率EGARCH(1,1膜型估计结果在EGARCH中，告三项的系数估计值都小于零。在估计结果中沪市为-0.051846,深市为-0.032059,而且都是显著的，这也说明了沪深股市中都存在杠杆效应。(三)沪深股市波动溢出效应的研究当某个资本市场出现大幅波动的时候，就会引起投资者在另外的资本市场的投资行为的改变，将这种波动传递到其他的资本市场。这就是所谓的“溢出效应”。例如9.11恐怖袭击后，美国股市的大震荡引起欧洲及亚洲股市中投资者的恐慌，从而引发了当地资本市场的大动荡。接下来我们将检验深沪两市之间的波动是否存在“溢出效应”。1、检验两市波动的因果性(1)提取

33、条件方差重复前面GARCH-MI型建模的步骤，选择主菜单栏“Procs”下的“MakeGARCVarianceSeries，得到rh回归方程残差项的条件方差数据序列GARCHQ1同样的步骤rz回归方程残差项的条件方差数据序列GARCH02(2)检验两市波动的因果性在“WorMle”中同时选中“GARCH01和“GARCH02,右击，选择“Open”一“AsGroup”，在弹出的窗口中点击ViewGrangerCausality”,并选择滞后阶数5,得到如图721所示结果。图7-21Granger因果检验可见，我们不能拒绝原假设：上海的波动不能因果深圳的波动。但是可以拒绝原假设：深圳的波动不能因

34、果上海的波动。这初步证明沪深股市的波动之间存在溢出效应，且是不对称，单向的，表明是由于深圳市场的波动导致了上海市场的波动，而不是相反。2、修正GARCH-罐!型在沪市GARCH-M模型的条件方差方程中加入深市波动的滞后项，应该会改善估计结果。在EquationSpecification”窗口中，按图7-22示输入如下变量，即在模型的条件方差方程中加入了深市波动的滞后项。图7-22修正GARCH-罐!型点击“OK,则得到加入滞后项GARCH0后沪市GARCH-MI型重新估计的结果，如图723所不。-.51lX9EVievs-LEquation：UITITLEDVorkfile：CMATTER6.4.FileEdit0bjctsViewFrees