论文（设计）中国粮食价格波动分析 基于ARCH 类模型.doc

《论文（设计）中国粮食价格波动分析 基于ARCH 类模型.doc》由会员分享，可在线阅读，更多相关《论文（设计）中国粮食价格波动分析 基于ARCH 类模型.doc（11页珍藏版）》请在三一办公上搜索。

1、中国粮食价格波动分析：基于ARCH类模型*本文得到中国社会科学院重点课题“中国粮食价格预测模型比较研究”的资助。中国社会科学院农村发展研究所各位专家、学者在双周学术交流午餐会上对本文的形成提出了建设性意见，在此表示感谢。罗万纯1 刘 锐2内容提要：了解粮食价格波动的特征对采取相应政策稳定粮食价格具有重要的现实意义。本文利用GARCH、GARCH-M、TARCH和EGARCH等ARCH类模型对粮食价格的波动、波动的非对称性进行了分析。研究表明：籼稻、粳稻、大豆价格没有显著的异方差效应；小麦和玉米价格波动有显著的集簇性；小麦市场和玉米市场没有高风险高回报的特征；小麦价格波动有非对称性，即价格上涨信

2、息引发的波动比价格下跌信息引发的波动大。本文在此基础上提出：可以利用价格波动的集簇性对未来的价格波动进行预测；要不断完善粮食市场，引导市场参与主体理性投资；要特别关注引起价格上涨的因素并采取相应措施。关键词：粮食 价格 波动 ARCH类模型一、引言近年来，中国粮食价格频繁波动。19972007年，籼稻、粳稻、小麦、大豆价格呈现相同的变化趋势，1997年3月至2003年9月价格不断下跌，但从2003年10月开始价格不断上涨。玉米价格的波动与其他品种粮食价格的波动有些差异，1997年3月至2000年4月不断下跌，但从2000年5月开始呈现在波动中不断上涨的变化趋势。粮食价格的频繁波动对生产者行为、

3、消费者行为以及宏观经济都产生了重大影响，因此，了解粮食价格波动的特征对采取相应政策稳定粮食价格具有重要的现实意义。粮食价格波动问题一直备受关注，有很多学者从不同角度进行了研究。关于粮食价格波动的特点，冯云（2008）的研究表明，粮食价格波动具有集簇性和明显的非对称性。关于粮食价格波动的影响因素，Lapp and Smith（1992）认为，粮食价格波动水平直接和间接受到宏观经济政策特别是货币政策的影响；钟甫宁（1995）强调了稳定的政策和统一的市场对避免粮食价格人为波动的重要性；柯炳生（1996）认为，农户的粮食储备及其市场反应行为是造成粮食价格波动的重要原因之一；谭江林、罗光强（2009）的

4、研究表明，通货膨胀是粮食价格波动的Granger原因。关于粮食价格波动产生的影响，石敏俊、王妍、朱杏珍（2009）的研究表明，能源价格和粮食价格上涨及其带来的饲料价格上涨和劳动力成本上升所导致的成本驱动效应对畜产品价格实际涨幅的影响在44%59%之间；对加工食品价格的影响较为显著，占加工食品价格实际涨幅的74%左右。能源价格上涨对CPI和物价总水平的拉动作用高于粮食价格波动的影响。受粮食价格上涨的影响，城镇居民净收益减少，农村居民净收益增加。何蒲明、黎东升（2009）认为，中国粮食产量与价格波动均较大，并且价格波动比产量波动更大，对国家粮食安全造成了不利影响；而粮食产量与粮食价格有密切的关系，

5、粮食价格是粮食产量变化的原因。从研究方法上来说，国外研究价格波动问题常常采用ARCH类模型。ARCH类模型能准确地模拟时间序列变量的波动，使人们能更加准确地把握波动（风险）。国内应用ARCH类模型的研究成果也比较丰富，现有的研究成果主要集中在股票市场和期货市场方面。比如，陈千里（2002）利用GARCH类模型，以上证综合指数为对象，对中国股市波动进行了实证研究，并在中国股市的背景下对集簇性和不对称性的几种经济解释进行了理论分析。华仁海、仲伟俊（2003）运用ARCH类模型对中国期货市场中期货价格、收益、交易量、波动性相互之间的关系进行了动态分析。刘宁（2004）应用ARCH类模型对上海股市的波

6、动进行了分析。唐衍伟、陈刚、张晨宏（2004）运用ARCH类模型对中国铜、大豆、小麦期货市场的波动和有效性进行了研究。总的来说，关于粮食价格波动和ARCH类模型应用的研究成果非常丰富，为本文的研究提供了参考和借鉴。但是，现有研究还有进一步拓展的空间：一是目前关于粮食价格波动的描述性分析较多，但计量分析较少；二是利用ARCH类模型研究粮食价格波动的成果很少，冯云（2008）的研究没有分具体粮食品种，从实际情况看，不同粮食品种价格的波动存在差异，因而有必要细分品种来研究。为丰富和深化现有研究，本文将利用ARCH类模型对粮食价格的波动、波动的非对称性进行分析，主要验证以下几个问题：粮食价格波动是否具

7、有集簇性？粮食市场是否有高风险高回报的特征？粮食价格波动是否具有非对称性？在本文中，粮食品种指籼稻、粳稻、小麦、玉米、大豆；集簇性指大的价格变化往往跟随着大的价格变化，小的价格变化往往跟随着小的价格变化；非对称性指价格下跌信息引发的波动和价格上涨信息引发的波动不一样大。二、研究方法（一）波动分析本文首先对价格序列和价格收益率序列进行描述性分析和单位根平稳性检验，然后设定均值方程并进行ARCH-LM检验，最后建立GARCH模型和GARCH-M模型。GARCH模型主要检验波动是否具有集簇性；GARCH-M模型主要检验粮食市场是否有高风险高回报的特征。1.(G)ARCH模型。自回归条件异方差（aut

8、oregressive conditional heteroskedasticity，简写为“ARCH”）模型由Engle（1982）提出，由两个方程组成： （1） （2）（1）式称为均值方程，其中，为被解释变量，在本文中表示粮食价格收益率，为解释变量，在本文中只包括的滞后项；（2）式称为方差方程，其中，表示在时刻的条件方差，在方差方程中它被定义为残差滞后项的加权平方和。（2）式中，1，以确保条件方差。为ARCH项，如果ARCH项高度显著，说明粮食价格收益率具有显著的波动集簇性。Bollerslev（1986）提出了（2）式的拓展形式，即广义自回归条件异方差（generalized autor

9、egressive conditional heteroskedasticity，简写为“GARCH”）模型。在ARCH模型的方差方程（2）中加入条件方差自身的滞后项就得到GARCH模型： （3）（3）式中，为ARCH项，为GARCH项，和分別为它们的滞后阶数，如果ARCH项和GARCH项都高度显著，说明粮食价格收益率具有显著的波动集簇性。为保证条件方差非负，一般要求系数和，但这个系数的非负性要求只是保证模型有意义的充分条件而非必要条件。GARCH模型将波动来源划分为两部分：变量过去的波动和外部冲击，而和则分别反映了它们对本期波动的作用强度。模型系数之和的大小反映了波动的持续性，当它小于1时，

10、说明冲击的影响会逐渐消失，当它大于1时，说明冲击的影响不但不会消失，反而会扩散。相对于ARCH，GARCH模型的优点在于：可以用较为简单的GARCH模型来代表一个高阶ARCH模型，从而使模型的识别和估计都变得比较容易。2.(G)ARCH-M模型。(G)ARCH-M（GARCH-in-mean）模型由Engle，Lilien and Robins（1987）提出。在(G)ARCH模型的均值方程（1）中加入就转化为(G)ARCH-M模型： （4）（4）式中，是条件标准差的一个倍数，若为正数，就意味着市场参与主体因风险增加而要求更高的收益，该参数用来验证粮食市场是否有高风险高回报的特征。 （二）波动

11、非对称性分析本文用门槛ARCH（threshold ARCH，简写为“TARCH”）模型和指数GARCH（exponential generalized autoregressive conditional heteroskedasticity，简写为“EGARCH”）模型对波动的非对称性进行分析，两个模型的估计结果可以相互验证。1.TARCH模型。TARCH模型由Rabemananjara and Zokoian（1993）提出，其条件方差方程为： （5）（5）式中，是虚拟变量，当时，否则，。此模型中，价格上涨信息（）对条件方差的影响为，而价格下跌信息（）的影响为。如果，表明波动具有非对称性

12、。当时，表明价格下跌信息引发的波动比价格上涨信息引发的波动大；当时，表明价格上涨信息引发的波动比价格下跌信息引发的波动大。2.EGARCH模型。EGARCH模型由Nelson（1991）提出，其条件方差方程为： （6）（6）式中，价格上涨信息（）对的影响为，价格下跌信息（）的影响为。如果，表明波动具有非对称性。当时，表明价格上涨信息引发的波动比价格下跌信息引发的波动大；当时，表明价格下跌信息引发的波动比价格上涨信息引发的波动大。三、数据描述性统计分析本文使用的是1997年3月至2007年12月的集贸市场月度价格（元/公斤）数据，数据来源于中国农产品价格调查年鉴（20042008年，历年）国家统

13、计局农村社会经济调查司：中国农产品价格调查年鉴（20042008年，历年），中国统计出版社。价格收益率以相邻月份粮食价格的对数一阶差分表示，即，其中，和分别表示第t月和第t-1月的价格。价格收益率的描述性统计量见表1。籼稻、粳稻、小麦、玉米、大豆价格收益率的峰度都高于正态分布的峰度值3，说明价格收益率具有尖峰和厚尾特征；JB正态性检验也证实了价格收益率显著异于正态分布。本文用Eviews3.0来进行计量分析。表1 价格收益率基本统计量籼稻粳稻小麦玉米大豆平均值0.00150.00110.00010.00320.0022标准差0.03190.03410.02810.03100.0251偏度1.9

14、5161.66860.2095-0.12131.0825峰度14.353114.22846.28904.90207.6186JB正态性检验774.6908737.526959.086919.7606139.8495观察数129129129129129从价格收益率变化图（图1图5）可见，籼稻、粳稻、小麦、玉米、大豆价格存在波动的集聚现象和异方差效应，为此，要对各个序列进行ARCH-LM检验，来判断它们是否存在异方差效应。图1 籼稻价格收益率图2 粳稻价格收益率图3 小麦价格收益率图4 玉米价格收益率图5 大豆价格收益率四、模型估计结果 单位根平稳性检验结果表明：价格序列都非平稳，而价格收益率序列

15、都平稳。ARCH-LM检验结果表明：籼稻、粳稻、大豆价格异方差效应不显著。对小麦而言，选择滞后阶数为5阶时，检验概率p值小于5%，说明残差序列存在异方差效应。至少存在5阶的ARCH效应。这就意味着必须估计很多个参数，而这很难精确做到，可以用一个低阶的GACH模型代替，以减少待估参数个数。对玉米而言，选择滞后阶数为1阶时，检验概率p值小于5%，说明残差序列存在异方差效应。根据ARCH-LM检验结果，本文对存在显著异方差效应的小麦和玉米市场建立ARCH类模型。（一）小麦小麦市场的ARCH类模型估计结果如下（详见表2）：1.GARCH模型估计结果。小麦价格收益率条件方差方程中，和都在10%的水平下显

16、著，表明价格收益率序列具有显著的波动集簇性。和之和为0.68，小于1，因此，过去的波动对未来的影响逐渐消失。2. GARCH-M模型估计结果。小麦价格收益率均值方程中估计值为0.01，但不显著，这反映了小麦市场没有高风险高回报的特征。表2 小麦市场ARCH类模型估计结果GARCHGARCH-MTARCHEGARCH0.01（0.11）0.19（1.54）0.21（2.01）0.18（1.70）0.18（1.72）0.00（3.92）0.00（3.65）0.00（3.98）-2.4（-3.86）0.19*（1.87）0.20（1.81）0.35（2.31）-0.34*（-2.09）0.49*（4

17、.53）0.50（4.27）0.46（4.59）0.23（1.46）0.14*（1.71）注：括号中的数值为z统计值；*表示在10%水平下显著，*表示在5%水平下显著，*表示在1%水平下显著。3.TARCH和EGARCH模型估计结果。在TARCH模型中，的估计值小于零，且在5%的水平下显著。在EGARCH模型中，的估计值大于零，且在10%的水平下显著。这说明，小麦市场中价格上涨信息引发的波动比价格下跌信息引发的波动要大，小麦价格波动具有显著的非对称性。（二）玉米玉米市场的ARCH类模型估计结果如下（详见表3）：1.ARCH模型估计结果。玉米价格收益率条件方差方程中在1%水平下高度显著，表明玉米

18、价格波动具有显著的集簇性。2. ARCH-M模型估计结果。玉米价格均值方程中估计值为0.04，但不显著，这反映了玉米市场不存在高风险高回报的特征。3.TARCH和EGARCH模型估计结果。在TARCH模型中，的估计值大于零，但不显著；在EGARCH模型中，的估计值小于零，但不显著。这说明，玉米价格波动没有显著的非对称性。表3 玉米市场ARCH类模型估计结果ARCHARCH-MTARCHEGARCH0.04（0.60）0.43（4.08）0.40（3.52）0.39（3.51）0.39（4.08）0.00（8.38）0.00（8.47）0.00（8.7）-7.48（-52.02）0.43*（2.

19、89）0.38（2.69）0.28（1.67）0.15（0.61）0.57（3.04）-0.12（-1.32）注：括号中的数值为z统计值；*表示在10%水平下显著，*表示在5%水平下显著，*表示在1%水平下显著。 总的来说，ARCH类模型的估计结果表明：小麦、玉米价格波动具有显著的集簇性；小麦、玉米市场没有高风险高回报的特征；小麦价格波动具有显著的非对称性，即价格上涨信息引发的波动大于价格下跌信息引发的波动，玉米价格波动没有显著的非对称性。关于产生波动集簇性的原因，目前有两种解释：第一种是将波动与宏观经济形势联系在一起，认为信息产生过程的序列相关性导致了波动的集簇性；第二种是认为市场上存在看法

20、不同的风险厌恶交易者，Harris and Raviv（1993）和Shalen（l993）都提出了看法不同会产生波动集簇性的理论模型。从中国粮食市场的状况来看，笔者更倾向于第二种解释，因为中国粮食市场与宏观经济走向的联系比较弱，随着粮食市场化政策的实施，市场参与主体不断增多，人们对信息的理解方式有很大的不同，这导致价格波动集簇性的产生。对波动非对称性的解释，French，Schwert and Stambaugh（1987）和Campell and Hentschel（1992）提出了波动反馈效应：在价格上涨信息的冲击下，增加的波动带来降低当期价格的效应与信息引起的价格上涨相抵消，从而使波动

21、减小；在价格下跌信息的冲击下，增加的波动带来降低当期价格的效应与信息引起的价格下降相叠加，使波动增大。但是，粮食市场比较特殊，和股市存在较大差异，人们更关心价格上涨信息，存在“跌价无人问，涨价多头管”的现象，这导致价格上涨信息引发的波动比价格下跌信息引发的波动更大。五、总结和讨论本文通过对籼稻、粳稻、小麦、玉米、大豆价格波动进行分析，得出以下结论：籼稻、粳稻、大豆价格没有显著的异方差效应。小麦和玉米价格波动具有显著的集簇性。小麦市场和玉米市场没有高风险高回报的特征。小麦价格波动具有非对称性，即价格上涨信息引发的波动比价格下跌信息引发的波动大。与现实相联系，模型估计结果表明：小麦、玉米市场大的价

22、格波动后面往往跟随着大的价格波动，小的价格波动后面经常跟随着小的价格波动，这说明，小麦、玉米价格波动在一定程度上是可以预测的。小的价格波动影响较小，但大的价格波动会对粮食产业发展、居民生活、宏观经济产生较大的影响，为稳定粮食市场，大的价格波动不可不防。粮食市场没有体现出高风险高回报、低风险低回报的特征，这说明，中国粮食市场的大部分交易者在做决策时非理性因素大于理性因素，中国粮食市场有待进一步发展和完善。对小麦市场来说，价格上涨信息引发的波动比价格下跌信息引发的波动大。为稳定小麦市场，要特别关注引起小麦价格上涨的因素并采取相应措施。本文只是利用ARCH类模型对不同品种粮食价格波动问题进行了初步分

23、析，未来的研究可以从以下几个方面进行深化：一是波动的溢出效应；二是ARCH类模型的修正；三是ARCH类模型在粮食价格预测中的应用。参考文献1.冯云：中国粮食价格波动的实证分析，价格月刊2008年第2期。2.钟甫宁：稳定的政策和统一的市场对我国粮食政策的影响，中国农村经济1995年第7期。3.柯炳生：中国农户粮食储备及其对市场的影响，中国农村观察1996年第6期。4.谭江林、罗光强：粮食价格波动与通货膨胀关系的实证研究，价格月刊2009年第3期。5.石敏俊、王妍、朱杏珍：能源价格波动与粮食价格波动对城乡经济关系的影响基于城乡投入产出模型，中国农村经济2009年第5期。6.何蒲明、黎东升：基于粮食

24、安全的粮食产量和价格波动实证研究，农业技术经济2009年第2期。7.陈千里：中国股市波动集簇性和不对称性研究，湖北大学学报（自然科学版）2002年（第24卷）第3期。8.华仁海、仲伟俊：我国期货市场期货价格收益、交易量、波动性关系的动态分析，统计研究2003年第7期。9.唐衍伟、陈刚、张晨宏：我国期货市场的波动性与有效性基于三大交易市场的实证分析，财贸研究2004年第5期。10.刘宁：对上海股票市场波动性的ARCH研究，兰州大学学报（自然科学版）2004年（第40卷）第6期。 11.Bollerslev, T.: Generalized Autoregressive Conditional H

25、eteroskedasticity, Journal of Econometrics, 31: 307-328,1986.12.Campell, J. Y. & Hentschel, L.: No News is Good News: an Asymmetric Model of Changing Volatility in Stock Returns, unpublished manuscript, Princeton University, 1990.13.Engle, R. F.: Autoregressive Conditional Heteroskedasticity with Es

26、timates of the Variance of U.K. Inflation, Econometrica, 50: 987-1008, 1982.14.Engle, R. F.; Lilien, D. M. & Robins, R. P.: Estimating Time Varying Risk Premia in the Term Structure: the ARCH-M Model, Econometrica , 55: 391-407, 1987.15.French, K. R.; Schwert, G. W. & Stambaugh, R. E.: Expected Stoc

27、k Returns and Volatility, Journal of Financial Economics, 19: 3-29, 1987.16.Harris, M.; Raviv, A.:Differences of Opinion Make a Horse Race, The Review of Financial Studies, 6(3): 473-506,1993.17.Lapp, J. S. & Smith, V. H.: Aggregate Sources of Relative Price Variability among Agricultural Commoditie

28、s, American Journal of Agricultural Economics, 74: 1-9, 1992.18.Nelson, D. B.: Conditional Heteroskedasticity in Asset Returns: a New Approach, Econometrica, 59: 347-70,1991.19.Rabemananjara, R. & Zakoian, J. M.: Threshold ARCH Models and Asymmetries in Volatility, Journal of Applied Econometrics, 8

29、: 31-49,1993.20.Shalen ,C. T.: Volume, Volatility and the Dispersion of Beliefs, The Review of Financial Studies, 6(2): 405-434,1993.（作者单位：1中国社会科学院农村发展研究所；2农业部农村经济研究中心）（责任编辑：小 林）Editors note: Judson Jones is a meteorologist, journalist and photographer. He has freelanced with CNN for four years, cov

30、ering severe weather from tornadoes to typhoons. Follow him on Twitter: jnjonesjr (CNN) - I will always wonder what it was like to huddle around a shortwave radio and through the crackling static from space hear the faint beeps of the worlds first satellite - Sputnik. I also missed watching Neil Arm

31、strong step foot on the moon and the first space shuttle take off for the stars. Those events were way before my time.As a kid, I was fascinated with what goes on in the sky, and when NASA pulled the plug on the shuttle program I was heartbroken. Yet the privatized space race has renewed my childhoo

32、d dreams to reach for the stars.As a meteorologist, Ive still seen many important weather and space events, but right now, if you were sitting next to me, youd hear my foot tapping rapidly under my desk. Im anxious for the next one: a space capsule hanging from a crane in the New Mexico desert.Its l

33、ike the set for a George Lucas movie floating to the edge of space.You and I will have the chance to watch a man take a leap into an unimaginable free fall from the edge of space - live.The (lack of) air up there Watch man jump from 96,000 feet Tuesday, I sat at work glued to the live stream of the

34、Red Bull Stratos Mission. I watched the balloons positioned at different altitudes in the sky to test the winds, knowing that if they would just line up in a vertical straight line we would be go for launch.I feel this mission was created for me because I am also a journalist and a photographer, but

35、 above all I live for taking a leap of faith - the feeling of pushing the envelope into uncharted territory.The guy who is going to do this, Felix Baumgartner, must have that same feeling, at a level I will never reach. However, it did not stop me from feeling his pain when a gust of swirling wind k

36、icked up and twisted the partially filled balloon that would take him to the upper end of our atmosphere. As soon as the 40-acre balloon, with skin no thicker than a dry cleaning bag, scraped the ground I knew it was over.How claustrophobia almost grounded supersonic skydiverWith each twist, you cou

37、ld see the wrinkles of disappointment on the face of the current record holder and capcom (capsule communications), Col. Joe Kittinger. He hung his head low in mission control as he told Baumgartner the disappointing news: Mission aborted.The supersonic descent could happen as early as Sunday.The we

38、ather plays an important role in this mission. Starting at the ground, conditions have to be very calm - winds less than 2 mph, with no precipitation or humidity and limited cloud cover. The balloon, with capsule attached, will move through the lower level of the atmosphere (the troposphere) where o

39、ur day-to-day weather lives. It will climb higher than the tip of Mount Everest (5.5 miles/8.85 kilometers), drifting even higher than the cruising altitude of commercial airliners (5.6 miles/9.17 kilometers) and into the stratosphere. As he crosses the boundary layer (called the tropopause), he can

40、 expect a lot of turbulence.The balloon will slowly drift to the edge of space at 120,000 feet (22.7 miles/36.53 kilometers). Here, Fearless Felix will unclip. He will roll back the door.Then, I would assume, he will slowly step out onto something resembling an Olympic diving platform.Below, the Ear

41、th becomes the concrete bottom of a swimming pool that he wants to land on, but not too hard. Still, hell be traveling fast, so despite the distance, it will not be like diving into the deep end of a pool. It will be like he is diving into the shallow end.Skydiver preps for the big jumpWhen he jumps

42、, he is expected to reach the speed of sound - 690 mph (1,110 kph) - in less than 40 seconds. Like hitting the top of the water, he will begin to slow as he approaches the more dense air closer to Earth. But this will not be enough to stop him completely.If he goes too fast or spins out of control,

43、he has a stabilization parachute that can be deployed to slow him down. His team hopes its not needed. Instead, he plans to deploy his 270-square-foot (25-square-meter) main chute at an altitude of around 5,000 feet (1,524 meters).In order to deploy this chute successfully, he will have to slow to 1

44、72 mph (277 kph). He will have a reserve parachute that will open automatically if he loses consciousness at mach speeds.Even if everything goes as planned, it wont. Baumgartner still will free fall at a speed that would cause you and me to pass out, and no parachute is guaranteed to work higher tha

45、n 25,000 feet (7,620 meters).It might not be the moon, but Kittinger free fell from 102,800 feet in 1960 - at the dawn of an infamous space race that captured the hearts of many. Baumgartner will attempt to break that record, a feat that boggles the mind. This is one of those monumental moments I will always remember, because there is no way Id miss this.