张 华1 杨 雪2 李松华2 王 菲z
（1．中国矿业大学（北京），北京市海淀区，100083：
2．华北水利水电大学，河南省郑州市，450045)
摘 要 运用GARCH类计量经济学模型，分析了煤矿安全事件舆情演化的波动性特征，基于2011年1月1日-2014年5月1 6日的煤矿事故百度搜索指数的周平均数据的实证研究发现：煤矿安全事件舆情演化的波动性呈现出聚集效应；通过对模型残差进行条件异方差即ARCH效应检验，发现煤矿事故搜索指数的波动性具有ARCH效应，因此，对煤矿安全事件舆情演化波动性的研究适用GARCH类模型，模型结果可以为煤矿安全事件舆情的预警提供先决信息。
关键词 GARCH模型 煤矿安全 舆情 波动性
中图分类号 TD-9 文献标识码 A
Zhang Hual, Yang Xue2, Li Songhua2, Wang Fei2
(l. China University of Mining&Technology, Beijing, Haidian, Beijing 100083, China;
2. North China University of Water Resources&Electric Power, Zhengzhou, Henan 450045, China)
AbstractGARCH econometric model is used to analyze volatility characteristics of public o一 pinion evolution of coal mine safety incidents. Based on an empirical study on weekly average dataof coal mine incidents from January l,2011 to May 16, 2014 in Baidu search index, it can befound that: the volatility of public opinion evolution of coal mine safety incidents shows aggrega-tion effect; the model residual error is examined by ARCH effect, and it shows that the volatilityof coal mine incidents search index has ARCH effect. Therefore, GARCH model is applicable tothe study on the volatility of public opinion evolution of coal mine safety incidents, and the modelresults can provide prerequisite information for the warning of public opinion of coal mine safetyincidents.
Key words GARCH model, coal mine safety, public opinion, volatility
煤矿安全事件舆情，即公众对煤矿事故的关注，其演化波动呈现出由初步传播到迅速扩散再到消减的生命周期性特征。同时，煤矿安全舆情的波动具有聚集性和非对称性特性，对于这些特性的把握，一般的时间序列建模的ARMA模型显得无能为力。因此，对于煤矿安全事件舆情波动性的研究，可以采用金融时间序列波动性应用较多的GARCH类模型。数据的可得性，本项研究选择了煤矿事故搜索百度指数作为煤矿安全事件舆情的替代指标，以SI，表示，样本长度为2011年1月1日-2014年5月16日，数据频率为周平均数据。舆情演化具有由初步传播（即形成）到迅速扩散（即发展）再到消退的周期性特征，煤矿安全事件的舆情也不例外、具有这一周期性特征，如图1所示。
搜索指数情况
此外，突发事件舆情演化呈现出波动聚集性特征，即一段时间内波动性较大，而在另一段时间内波动性较小。原因在于突发事件发生后，公众由于从众心理关注事件从而导致网络点击和搜索迅速增加，而关注度的上升有可能意味着新一轮的舆论波动。图1为煤矿事故搜索指数体现出煤矿安全事件舆情演化的聚集性特征，图2为煤矿事故搜索指数变化率情况，它明确显示了该聚集性特征。定义煤矿事故搜索指数变化率
1.2 变量平稳性ADF检验
为避免采用时间序列数据进行实证研究可能出现的“伪回归”现象，即回归残差非平稳、回归关系不能真实反映被解释变量和解释变量之间的均衡关系，需要对煤矿事故搜索指数序列进行平稳性检验。本部分采用应用较为广泛的ADF方法来检验序列的平稳性，检验结果如表1所示。由表1可知，煤矿事故搜索指数序列SI，的ADF检验值为-6. 236，小于l%显著性水平下的临界值，接受序列存在单位根原假设的概率p值为0. 000，根据小概率事件不易发生的原则，拒绝序列存在单位根的原假设，因此，煤矿事故搜索指数序列SI，为平稳序列。
2 煤矿事故舆情波动的ARCH效应检验
2.1 自回归阶数选择
舆情演化的生命周期呈现出自回归特性。本部分采用偏自相关系数PAC来确定自回归条件异方差ARCH模型中条件均值方程的形式，即自回归的阶数。图3的煤矿事故搜索指数SI，序列的偏自相关系数PAC的检验结果显示，条件均值方程的自回归阶数选择为1，原因是偏自相关系数PAC在滞后1阶时截尾。
煤矿事故舆情演化的一阶自回归过程可由最小二乘回归得到：
方程(l)式中圆括号中的数字为对应参数的t检验值，说明各回归系数是显著的。
2.2 ARCH效应检验
由自回归方程(l)的残差图4可以看出，残差的波动呈现聚集性，说明误差项可能存在条件异方差。
对于回归方程误差项即残差项的条件异方差特性即ARCH效应检验，通常采用Engle (1982)提出的拉格朗日乘数检验，即ARCH LM检验。
该检验通过如下辅助回归进行：
2.3 GARCH效应检验
由上文煤矿事故搜索指数SI，及其自回归残差序列呈现聚集性，以及煤矿事故搜索指数自回归残差的ARCH效应检验可知，煤矿安全事件舆情演化即煤矿事故搜索指数的残差序列存在条件异方差，则为了更好地描述煤矿安全事件舆情演化的波动性特征，本部分采用GARCH模型来建模分析。其中，GARCH(q，p)模型中的均值方程仍采取随机游走的自回归形式，条件方差方程中的ARCH项和GARCH项阶数的选择则通过不同模型的比较而得以确定。GARCH（q，p）模型的具体形式如下：
(4)式中第一个方程即为均值方程、第二个方程为条件方差方程。为确定条件方差方程中ARCH项和GARCH项的阶数，选择了GARCH(1，1)、GARCH(1，2)、GARCH(2，1)以及GARCH(2，2)4个模型进行比较，如表3所示。
最优GARCH模型的选择规则是模型的对数似然值较大而赤池和施瓦茨信息准则较小，同时参数的估计值显著，则由表3的结果可知，最优模型应为GARCH(1，1)，其估计结果如下：
(5)式GARCH(1，1)模型中的均值方程与上文(l)式的煤矿事故舆情演化的一阶自回归随机游走模型相比，具有较大的对数似然值和较小的AIC和SC值，同时GARCH(1，1)模型中条件方差方程的3个参数均为正、统计上均显著、模型估计合理，从而说明模型GARCH(1，1)较优，能够更好地模拟煤矿事故搜索指数的波动性特征，如图5所示。由图5可以看出，GARCH(1，1)模型不仅成功地模拟了煤矿事故搜索指数的每一次突变，而且对煤矿事故搜索指数突变强度的把握也较为理想，因此，GARCH(1，1)模型较为成功地刻画了煤矿安全事件舆情演化的波动性特征，可以为煤矿安全事件舆情的预警提供先决信息。
此外，GARCH(1，1)模型克服了煤矿事故搜索指数自回归残差序列所存在的条件异方差性，如表4所示：GARCH(1，1)模型残差序列异方差的ARCH I。M检验统计量接受原假设即不存在条件异方差的概率为0. 8286．置信度较高，且相伴概率也高达0. 8299，由此说明GARCH(1，1)模型克服了煤矿事故舆情演变方程残差序列所存在的ARCH效应，图6的自相关、偏自相关检验结果也支持了这一结论：自相关系数AC和偏自相关系数PAC的值都近似为O，Q统计量不显著的概率在0. 827以上。
3煤矿事故舆情波动的非对称效应检验
在资本市场上，金融资产价格的变化往往表现出这样一种现象：金融资产价格的波动对市场上的好消息和坏消息的反应不同——坏消息对金融资产价格波动性的影响要大于好消息的，即存在着非对称性效应，亦即杠杆效应。这一现象在舆情演化的波动性中也存在，如万源(2012)认为网络舆情信息演化的波动性即存在这种非对称性效应。因此，为了检验煤矿安全事件舆情演化的波动性是否存在非对称性效应，本节将应用门限ARCH模型即TARCH模型进行检验。
3.1 TARCH模型估计
一般地，TARCH模型中的条件方差方程形式如下：
表5给出了门限阈值分别取1和2的TARCH模型估计结果的相关指标，由各指标的比较可知，TARCH(1，1)模型的性能相对较好。同时，门限阈值取1和2时，对条件方差方程中的ARCH项和GARCH项阶数分别取高阶，TARCH模型估计结果的性能也不比TARCH(1，1)模型要好，因此，对于煤矿安全事件舆情演化波动性的非对称效应的检验采取了TARCH(1，1)模型。
TARCH(1，1)模型的估计结果如(7)式：
3.2 TARCH模型效果分析
由图7的TARCH模型和GARCH模型对煤矿
事故搜索指数波动性的模拟比较可知，TARCH模型比GARCH模型的表现更优，TARCH模型无论是对煤矿事故搜索指数波动性的形态还是程度的模拟都要好于GARCH模型。
同样地，TARCH(1，1)模型也克服了煤矿事故搜索指数自回归残差序列所存在的条件异方差性，如表6所示：TARCH(1，1)模型残差序列条件异方差的ARCH LM检验统计量接受原假设即不存在条件异方差的概率为0. 5827，置信度较高，且相伴概率也高达0. 5853，由此说明TARCH(1，1)模型克服了煤矿事故舆情演变方程残差序列所存在的ARCH效应。自相关、偏自相关检验结果也支持了这一结论：自相关系数AC和偏自相关系数PAC的值都近似为O，Q统计量不显著的概率基本都在0. 857以上。
4 结语
本研究运用GARCH类计量经济学模型分析了煤矿安全事件舆情演化的波动性特征，基于2011年1月1日-2014年5月16日的煤矿事故百度搜索指数的周平均数据的实证研究发现：
(l)基于对煤矿事故搜索指数及其变化率的统计分析发现，煤矿安全事件舆情演化的波动性呈现出聚集效应。
(2)关于煤矿安全事件舆情演化波动性的研究适用GARCH类模型。舆情演化具有自回归特性．基于煤矿事故搜索指数建立了煤矿安全事件舆情演化的一阶自回归随机游走模型，并对模型残差进行了条件异方差即ARCH效应检验，发现煤矿事故搜索指数的波动性具有ARCH效应，因此，对煤矿安全事件舆情演化波动性的研究适用GARCH类模型。
(3) TARCH模型对煤矿安全事件舆情演化波动性的把握优于GARCH模型。实证结果表明．GARCH(1，1)模型和TARCH (l.1)模型都成功地刻画了煤矿安全事件舆情演化的波动性特征，并均克服了煤矿事故搜索指数自回归方程残差序列存在的条件异方差性。相比较而高．TARCH(1，1)模型考虑了好消息和坏消息对煤矿安全事件舆情演化波动性影响的不同，发现好消息能够降低煤矿安全事件舆情演化的波动性；同时TARCH(1，1)模型对煤矿安全事件舆情演化波动性的形态和程度的模拟都优于GARCH(1．I)模型。
