实验五期货最优套期保值率估计

1、实验五 期货套期保值模型一、实验项目：期货套期保值模型二、实验目的1、掌握运用时间序列模型估计中国期货交易的最优套期保值比率的方法；2、掌握评估期货套期比效果的方法；3、找到最佳的套期保值比模型。三、预备知识：（一）、关于最优套期比确定方法以空头期货保值为例1由套期保值收益方差风险达最小得到（1）用价格标准差表示风险最小套期比 单位现货相应的空头保值收益：b（k）=b(k)-b0(k)（两边求方差解出k）（2）用改变量标准差表示风险最小套期比 单位现货相应的空头保值收益: b（k）=s-kf （两边求方差解出k） 注意到（1）与（2）两种最优化方式得到有套期比k是不同的。2用收益率表示套期保值

2、比率。 空头保值收益率（V为现货市值） RH=(V-V0+D)-NF(F-F0)/V0= (V-V0+D)/V0-(NFF0/V0)(F-F0)/F0=RS-h*RF由收益率风险达最小求出套期比3 由对冲原理得到要实现期货与现货完全对冲，必须满足以下风险中性原理(现货与期货组合风险为0) Q*f +Q0*s=0 kf +s=0 k=Q/Q0=-S/F-ds/df<0(因同方向变化)上式表明，每单位现货需要k单位期货对冲其风险，负号表示交易方向要相反。 S/F 或ds/df可通过久期求出。（二）计算期货套期保值比率的相关模型虽然上述介绍的h=s/f可以求最优套期比，但是其操作性不强。首先要

3、求出三个量，然后再计算h，显然误差很大。为了减小误差，使用时间序列模型。1、简单回归模型（OLS）上述使方差风险最小求套期比的三种方法对应的三个OLS模型2 / 38OLS不足：上述三个模型假设条件是残差“独立同方差”，即在残差项具有同方差性的假设下，其回归系数即是要求的最优套期比，但是这一条件太强，在金融市场上难以满足。其中要解决最突出的两个问题（1）s与f有协整关系时，OLS所得到的结果小于最优套期比（2）三模型残差独立同方差问题。2、协整与误差修正模型（ECM）（1）期货价格序列与现货价格序列特点1）二者常常是非平稳的；2）二者具有两个经济逻辑性：二者有共同的趋势；期货到期时，二者有趋合

4、性。由此，二者存在协整关系，那么用OLS的估计量将是有偏的。Ghosh(1993)通过实证发现，当不恰当地忽略协整关系时，所计算的套期比将小于最优值。3）研究表明，使用ECM模型比OLS方法能够更有效地对冲现货头寸风险。（2）使用ECM模型计算最优套期比的两步估计法第一步：建立协整回归模型要注意在这一协整回归中保留残差：，方便第二步使用。第二步：建立误差修正模型（ECM，一般模型）：其实要建立的是ECM简单方程（*）修正误差模型（*）其中误差修正模型（*）只是模型（*）的适当变形，这两模型是等价的。与一般的修正误差模型比较要建立的修正误差模型的简单形式为ECM模型优劣：优点：考虑到s与f有协整

5、关系时，EMC模型比OLS方法能更有效地对冲现货头寸的风险。缺点：还没有解决模型残差异方差问题。3、误差修正二元GARCH模型（ECM-BGARCH模型）ECM模型虽然考虑到s与f有协整关系问题，但是还没有解决模型残差异方差问题。特别是金融资产价格、收益率等序列波动常常出现“聚集性”特征。即资产价格在大的变动后面有大的变动，在小的变动后面有小的变动。这些特征用GARCH模型来描述。注意：一元GARCH模型仅能估计单一变量的条件方差，无法估计序列之间协方差，所以要建立二元GARCH（B-GARCH）模型，常用以下两个考虑到误差修正项的模型：一是常相关系数的二元GARCH模型，二是D-BEKK模型

6、。（1）常相关系数的二元GARCH模型GARCH模型由均值方程和条件方差方程组成。1）均值方程2）GARCH方程为了简单，如果系数矩阵都用对角矩阵，则条件方差展开可得：最优套期比（2）D-BEKK模型1）均值方程2）GARCH方程如果系数矩阵都用对角矩阵，则条件方差展开可得最优套期比（三）期货套期保值比率绩效的评估以空头期货套期保值为例，一个以1单位现货多头头寸和h个单位期货空头头寸的套期保值组合，组合价值和组合利润分别为：，因此，空头套期保值组合收益率为空头套期保值（含多头）收益率的方差为若这一方差小，说明经过套期保值后，收益率稳定，保值效果好。四、实验内容（一）数据搜集和整理；（二）利用时

7、间序列模型估计最优套期比；（三）评估最小方差套期比的绩效。五、实验软件环景：Eviews软件。六、实验步骤：（一）数据搜集和整理1、搜集数据期货合约在交割前两个月最活跃，使价格信息释放最为充分，更能反映期货合约的真实价值，所以中国企业多用距离交割月份较近的期货合约进行保值。表1 上海AL现货期货价格2006年4月3日至2007年4月13日数据 单位：元/吨序号期货价f现货价s序号期货价f现货价s序号期货价f现货价s1199301964081189901920016120140214202200601990082189201914016220280213803200201968083189201

8、916016320270215004200101972084191001926016420320214805206202008085191201932016520560214806205902010086189201932016620480214407205902016087185401914016720570214608202602002088185001892016820510213709203602006089184601890016920650213701020420200609018620190401702055021370112082020260911868019300171202

9、902124012205502042092190001956017220480213601321120204809318790196401732058021360142130020820941875019500174207202140015208002070095188401954017520920214301621500212409619150197201761972021440172138021120971927019920177199602124018214002134098191802002017819350209501921470212009919200200201791971020

10、950202128020580100192702032018019970210202122080216401011950020480181198302098022219702168010219670209001821981020860232270021960103205002190018319890207602423470229401042089022800184196102064025237702382010520200225001851943020640262418023500106197102166018619550201002723120229401071920020800187196

11、402008028228202240010819250207201881956019700292184021680109195302068018919500199603021720215601101929020680190195102016031208402078011119300205501911955020140322138020980112194102054019219620203603321880218001131928020500193195702036034217402140011419680206401941951020250352187021720115201602100019

12、519570201403621800217201161992021250196195302000037217902174011719480212501971959019880382114021180118200602135019819890199503920780210001192049021650199198701995040210802058012019970219002001993019900412137020980121199302180020120160202004420820207401222015021700202197601996043207002048012320010215

13、002031979019780442075020760124203102135020419740197604520560202401252099021350205197701962046203502020012620800213202061949019550472032020170127209002128020719620193604819670194401282103021320208196001948049199501968012921280213502091960019600502038020280130207402137021019580195305119650199601312082

14、021260211195701942052198001952013221320214202121951019340531988019720133212102140021319590192405419670198201342036021000214196401942055196501958013519850207002151958019380561978019620136199102034021619610194405719910197601371981020350217195601938058197101960013819780207402181955019370591971019480139

15、204602155021919550193406019510194101402037021650220195401932061197801967014120000213602211935019300621952019510142205002160022219430192706319510194401432011021000223193901924064195201944014420100211602241938019250651965016460145200502116022519240192406619560194001462012021000226192701919067195401946

16、014719810208802271941019240681941019380148197702075022819410193006919080191901491986020950229195901936070190701920015020110210502301986019440711884019020151201502115023119800196607218800190401522035021180232196101951073185901891015320380213802331977019520741842018400154203202136023419670195207518440

17、184201552037021300235761868018520156204202128023677186301860015719990213602377818650186401582030021440238791890018760159201802144023980188101886016020240214402402、建立Eviews工作文件创建工作文件并输入数据FileNewWorkfile因为数据是无观测日期的，所以选择Undated-or-irreqular栏：start：1；end：233，OK手工输入数据QuickEmpty Group在Ser01输入s列数据；在Ser02输入

18、f列数据；改变量名：点击Ser01全选第一列，在命令栏输入s；点击Ser02全选第二列，在命令栏输入f。将文件保存命名为hr（二）运用单方程时间序列模型估计最优套期比1、用OLS模型估计最优套期比（1）建立S关于F的回归方程Dependent Variable: SMethod: Least SquaresDate: 04/30/11 Time: 20:35Sample: 1 234Included observations: 234VariableCoefficientStd. Errort-StatisticProb.  C3387.983937.25643.61478

19、80.0004F0.8491880.04669418.186210.0000R-squared0.587730    Mean dependent var20414.70Adjusted R-squared0.585953    S.D. dependent var1035.698S.E. of regression666.4354    Akaike info criterion15.85027Sum squared resid1.03E+08  

20、60; Schwarz criterion15.87981Log likelihood-1852.482    F-statistic330.7382Durbin-Watson stat0.296083    Prob(F-statistic)0.000000图1 S关于F回归方程 （1）t=（3.614788）（18.18621）p=（0.0004） （0.0000）ft系数的p值接近0，回归系数是显著的。回归结果得到每单位现货用0.849188单位期货进行空头保值，即最优套期比是0.849188。结论

21、1：由现货价S关于期货价F回归模型得到的套期比是：0.84918。评价：1）虽然模型（1）系数显著，但模型精度R2=0.587730离1较远，精度不太高。而且不能排除模型（1）是伪回归。2）这一结论只能保证在保值策略实施前（建模的样本内），模型（1）在一定程度上是有效的，不能保证在策略实施期（样本外）模型同样有效，所以使用这一结论进行套期保值需要注意到这些情况。（2）建立st关于ft的回归方程在工作文件窗口的命令区，生成差分序列，以及st，ft序列：GENR ds=s-s(-1)GENR df=f-f(-1)建立st与ft的OLS简单回归模型最小二乘估计的命令OLS：ds c dfDepend

22、ent Variable: DSMethod: Least SquaresDate: 04/30/11 Time: 20:40Sample (adjusted): 2 234Included observations: 233 after adjustmentsVariableCoefficientStd. Errort-StatisticProb.  C0.10853723.082110.0047020.9963DF0.5588040.0755777.3938740.0000R-squared0.191373    Mean dep

23、endent var-0.515021Adjusted R-squared0.187872    S.D. dependent var390.9655S.E. of regression352.3308    Akaike info criterion14.57556Sum squared resid28675647    Schwarz criterion14.60519Log likelihood-1696.053    F-sta

24、tistic54.66937Durbin-Watson stat2.638061    Prob(F-statistic)0.000000图2 S关于F的回归方程（含常数项）常数项概率很大，接受常数为0的假设，重新定义方程：OLS：ds dfDependent Variable: DSMethod: Least SquaresDate: 04/30/11 Time: 20:44Sample (adjusted): 2 234Included observations: 233 after adjustmentsVariableCoefficientStd

25、. Errort-StatisticProb.  DF0.5588030.0754137.4098920.0000R-squared0.191373    Mean dependent var-0.515021Adjusted R-squared0.191373    S.D. dependent var390.9655S.E. of regression351.5707    Akaike info criterion14.56698Sum square

26、d resid28675650    Schwarz criterion14.58179Log likelihood-1696.053    Durbin-Watson stat2.638060图3 S关于F的回归方程（不含常数项） （2）t=（7.409892）p=（0.0000） ft系数的p值接近0，回归系数是显著的，但每单位现货用0.558803单位期货进行空头保值，即最优套期比是0.558803。可见，分别用套期比公式得到有结果k是不同的：，结论2：由现货价差分S关于期货价差分F回归模型得到的套期比是：0

27、.558803。评价：1）虽然这一模型系数显著，但模型精度R2=0.191373，精度非常低。而且也不能排除模型（2）是伪回归。2）结论2只能保证在保值策略实施前（建模的样本内），S与F在一定程度上满足模型（2），不能保证在策略实施期（样本外）模型（2）同样有效。3）结论2与结论1相比，结论1是保证在保值策略实施前（建模的样本内），S与F在一定程度上满足模型（1）；结论2是保证在保值策略实施前（建模的样本内），S与F在一定程度上满足模型（2）。4）差分模型一般用于分析短期波动情况，所以模型（2）在不顾伪回归下，也只用于动态套期保值。2、用ECM模型估计最优套期比（1）对f和s分别进行平衡性检验

28、在F页面上，选ViewCorrelogramLevel，滞后期空格处填写24（用234除以10近似）Date: 04/30/11 Time: 20:55Sample: 1 234Included observations: 234AutocorrelationPartial CorrelationAC  PAC Q-Stat Prob       .|*|       .|*|10.9460.946211.990.0

29、00       .|*|       .|. |20.8960.018403.210.000       .|*|       .|. |30.849-0.007575.380.000       .|* |   &#

30、160;   .|. |40.798-0.050728.360.000       .|* |       .|. |50.7560.051866.320.000       .|* |       .|* |60.7260.092994.010.000   

31、    .|* |       .|* |70.7060.0901115.20.000       .|* |       *|. |80.676-0.0921226.90.000       .|* |      &#

32、160;.|. |90.6510.0201331.00.000       .|* |       .|. |100.625-0.0171427.30.000       .|* |       .|. |110.594-0.0211514.80.000     

33、60; .|* |       *|. |120.558-0.0781592.20.000       .|* |       *|. |130.518-0.0581659.20.000       .|* |       .|. |140

34、.480-0.0171717.20.000       .|* |       *|. |150.438-0.0611765.50.000       .|* |       .|. |160.4030.0221806.60.000       .|

35、* |       .|. |170.3760.0421842.50.000       .|* |       *|. |180.343-0.0871872.60.000       .|* |       .|. |190.3150.021189

36、8.10.000       .|* |       *|. |200.283-0.0631918.80.000       .|* |       .|. |210.250-0.0311935.00.000       .|* | 

37、60;     *|. |220.210-0.0601946.50.000       .|* |       .|* |230.1820.0761955.10.000       .|* |       .|. |240.152-0.0391961.20.000图4

38、F序列相关分析图从图4的F序列自相关系数（AC）没有很快趋近0，说明序列F是非平稳的。因为期货价格（资产价格序列）往往有一定的趋势和截距，所以对ADF单位根检验时，选择同时具有趋势项和常数项的模型。滞后项p要精确确定就是AIC准则，粗略确定由系统默认。由上面分析，选择模型进行单位检验（Unit Root Test）。假设；备择假设。在工作文件窗口，选定变量F，双击它，在F页面上，点击ViewUnit Root TestADF，表示已经进入扩展的DF检验。在Test for unit root in中，选择Level(对水平变量进行单位根检验，检验系数对应的项Ft-1)Tren and inte

39、ret(含漂移项和时间趋势项)，其它选系统默认。Null Hypothesis: F has a unit rootExogenous: Constant, Linear TrendLag Length: 0 (Automatic based on SIC, MAXLAG=14)t-Statistic  Prob.*Augmented Dickey-Fuller test statistic-2.759187 0.2141Test critical values:1% level-3.9979305% level-3.42922910% level-3.13809

40、2*MacKinnon (1996) one-sided p-values.Augmented Dickey-Fuller Test EquationDependent Variable: D(F)Method: Least SquaresDate: 04/30/11 Time: 20:58Sample (adjusted): 2 234Included observations: 233 after adjustmentsVariableCoefficientStd. Errort-StatisticProb.  F(-1)-0.0627440.022740-2.7591

41、870.0063C1300.498471.05942.7607940.0062TREND(1)-0.3714370.316004-1.1754170.2410R-squared0.032185    Mean dependent var-1.115880Adjusted R-squared0.023769    S.D. dependent var306.0687S.E. of regression302.4093    Akaike info criterion14.274

42、23Sum squared resid21033818    Schwarz criterion14.31867Log likelihood-1659.948    F-statistic3.824393Durbin-Watson stat2.024316    Prob(F-statistic)0.023233图5 F序列单位根检验期货价格F序列的ADF检验统计量观察值为t=-2.759187，比概率1%、5%和10%对应的三个临界值都大。对应的概率0.2141也比1%、5

43、%和10%都大。所以这次ADF检验接受F非平稳的原假设，即认为F是非平稳的。对F序列一次差分进行ADF检验，与上不同的是在Test for unit root in中，选择1st difference(对F序列的一次差分进行单位根检验)，其它都相同。Null Hypothesis: D(F) has a unit rootExogenous: Constant, Linear TrendLag Length: 0 (Automatic based on SIC, MAXLAG=14)t-Statistic  Prob.*Augmented Dickey-Fuller tes

44、t statistic-15.81343 0.0000Test critical values:1% level-3.9981045% level-3.42931310% level-3.138142*MacKinnon (1996) one-sided p-values.Augmented Dickey-Fuller Test EquationDependent Variable: D(F,2)Method: Least SquaresDate: 04/30/11 Time: 21:06Sample (adjusted): 3 234Included observations: 2

45、32 after adjustmentsVariableCoefficientStd. Errort-StatisticProb.  D(F(-1)-1.0438060.066008-15.813430.0000C3.41852440.788420.0838110.9333TREND(1)-0.0436580.301593-0.1447570.8850R-squared0.521985    Mean dependent var-0.991379Adjusted R-squared0.517810   &

46、#160;S.D. dependent var443.0252S.E. of regression307.6361    Akaike info criterion14.30856Sum squared resid21672554    Schwarz criterion14.35313Log likelihood-1656.793    F-statistic125.0324Durbin-Watson stat2.001017    

47、Prob(F-statistic)0.000000图6 F序列一次差分单位根检验从图6看到，期货价格F一次差分序列的ADF检验统计量观察值为t= -15.81343，比概率1%、5%和10%对应的三个临界值都小。对应的概率0.0000也比1%、5%和10%都小。所以这次ADF检验拒绝F一次差分序列非平稳的原假设。即认为F一次差分序列是平稳的。所以FI（0），因此FI（1）。同理检验得到SI（0），因此SI（1）。（2）进行S和F的协整检验由于S和F都是一阶单整的，满足协整检验的前提。由前面已用OLS方法建立了S关于F的回归方程：t=（3.614788）（18.18621）p=（0.0004）

48、（0.0000）根据协整检验要求，还要检验残差是否平稳，先生成残差变量：GENR e=resid在工作文件窗口，选定变量e，双击e，在e页面上，点击ViewUnit Root TestADF，选择Level(检验的模型为被检验变量是e的差分，检验系数对应的项是et-1)None(不含常数项，不含时间项)。取系统默认p=4（因为没有做p=？AIC最小值），所以在lag lengthuser specifi中取4（表示差分滞后项数取4，即p=4）。得到结果：Null Hypothesis: D(E) has a unit rootExogenous: NoneLag Length: 2 (Auto

49、matic based on SIC, MAXLAG=4)t-Statistic  Prob.*Augmented Dickey-Fuller test statistic-14.14951 0.0000Test critical values:1% level-2.5751445% level-1.94222410% level-1.615772*MacKinnon (1996) one-sided p-values.Augmented Dickey-Fuller Test EquationDependent Variable: D(E,2)Method: Le

50、ast SquaresDate: 04/30/11 Time: 21:12Sample (adjusted): 7 234Included observations: 228 after adjustmentsVariableCoefficientStd. Errort-StatisticProb.  D(E(-1)-2.3848930.168550-14.149510.0000D(E(-1),2)0.6600240.1243015.3098930.0000D(E(-2),2)0.1869350.0651312.8701140.0045R-squared0.784291&#

51、160;   Mean dependent var1.614791Adjusted R-squared0.782374    S.D. dependent var748.7820S.E. of regression349.3100    Akaike info criterion14.56287Sum squared resid27453931    Schwarz criterion14.60799Log likelihood-1657.167

52、    Durbin-Watson stat2.049519图7 S关于F协整回归残差的单位根检验从图7看到，S关于F协整回归残差的ADF检验统计量观察值为t= -14.14951，比概率1%，5%、10%对应的两个临界值都小，对应的概率 0.0000也比1%，5%、10%小。ADF检验得到拒绝残差序列非平稳的原假设。即认残差序列是平稳的，即残差eI（0）。所以在1%的概率水平下，S与F序列存在协整关系，其协整方程为：t=（3.614788）（18.18621）p=（0.0004） （0.0000）（3）建立误差修正模型（ECM）由上面得知，S

53、与F序列存在协整关系。建立误差修正模型可分析向长期均衡状态调整的非均衡动态调整过程。原来协整模型形式如下（*）变成为误差修正模型（*）其中与一般的修正误差模型比较要建立的修正误差模型的简单形式为前面已生成了ECM所需要的序列变量GENR E=RESIDGENR DS=S-S(-1)GENR DF=F-F(-1)最小二乘估计命令建立修正误差模型OLS：DS C DF E(-1)得到回归结果为Dependent Variable: DSMethod: Least SquaresDate: 04/30/11 Time: 21:16Sample (adjusted): 4 234Included observations: 231 after adjustmentsVariableCoefficientStd. Errort-StatisticProb.  C0.05392122.384660.0024090.9981DF0.5715310.07