2. 华侨大学 现代应用统计与大数据研究中心,福建 厦门 361021
2. Research Center of Applied Statistics and Big Data, Huaqiao University, Xiamen 361021, China
股指期货是基于股票指数的金融衍生品,不仅能够对现货资产进行对冲来降低系统风险,还能作为一种投机套利工具来丰富资产组合,获得良好收益。股指期货与现货市场存在一定的关联,正常市场条件下双方表现出同涨同跌的情况[1]。现货组合与股指期货表现出的趋同程度越大,套期保值的效果越明显。标的指数是一篮子股票,在实际操作中,是无法使用所有现货股票来构建组合的。此外,随机或者人工指定少数股票来构建投资组合,不仅无法充分把握整体信息,又会产生一定的模拟误差,增加投资风险。因此,用较少的股票来构建现货投资组合,降低同步买卖的成本,并反映该投资组合与标的指数之间的相关性,是本文的研究关键。
目前,针对股指期货的套期保值研究有不少的研究和应用成果[2-6],其中通过指数复制和组合追踪来研究股指期货的套期保值也日益得到国内外研究者的关注,采用新的量化研究方法成为研究指数复制问题的一个方向。苏治等[7]建立混合整数线性规划模型并引入内核搜索分析框架,通过实证分析发现,增强型内核搜索法在成分股很大时才能够得到高质量的解,考虑现实交易成本特征的模型具有更好的稳健性,在指数追踪时投资组合的动态调整具有一定的必要性。倪禾[8]提出一种基于启发式遗传算法的寻优方案,通过最大效用函数来寻找一个最为经济的组合,该组合拥有最少的资产数、较少的权重调整次数和尽可能的接近或超越标的指数收益的优势,具有较强的实用性。胡春萍等[9]构建时间加权SVN的指数优化复制模型,不仅能够考虑时间因素对历史追踪误差的影响,还能够提高追踪组合的追踪效果。刘睿智等[10]使用自适应套索方法构造稀疏投资组合,使得该组合对指数的追踪效果良好。Chen等[11]设计了一种稳健的跟踪市场指数投资组合模型,该模型旨在利用0-1规划找到组合和标的指数之间的最大相似性,且在实验中,该模型得到的投资组合具有更低的跟踪误差和投资风险。Guastaroba等[12]为了平衡最小跟踪误差和超额收益的关系,设计出混合整数线性规划方法,利用内核搜索方法寻求最优解。该方法可以引入真实交易的特征,获得的最优组合也具有较好的表现,相比于传统的启发式方法具有更优的性能。Filippi等[13]针对指数追踪问题中获得超额收益和最小跟踪误差双的目标问题,提出了双目标整数线性规划模型。该模型旨在获得的组合在指数追踪上具有最小的追踪误差,并且能够获得超额收益。
时间序列聚类在金融领域的应用越来越受到研究者的关注,对探究金融市场的发展规律、把握市场信息起到重要的作用。Dose等[14]利用基于随机优化技术的时间序列聚类来增强指数跟踪,通过层次聚类得到不同的簇,每个簇中选择股票的子集并确定相应的权重后作为投资组合。该方法表明,通过聚类方法要比随机选择股票组合具有更少的噪声和更好的稳健性。Nanda等[15]利用K-means聚类将股票划分成不同的簇,并从不同簇中选择一些股票作为组合。由于组合来自不同的簇,分散投资使得风险最小化。柴尚蕾等[1]利用基于独立主成分分析和模糊C均值聚类的两阶段优化方式构建现货组合,相比随机方法得到更低的跟踪误差,提高了组合对大盘动荡的抵抗力。Lemieux等[16]分别探讨了通过3种传统聚类方法获得的投资组合上在实际交易中的应用,分析不同聚类技术是如何影响分析师对不同风险组合的看法。
时间序列聚类在金融领域得到了充分的利用,但传统时间序列聚类方法一般先指定初始簇中心,也不能充分反映其空间组织联系。社区发现是复杂网络研究中的重要研究工具,根据某种规则将网络中的节点划分为若干个社区,每个社区内部的节点连接较为紧密,而社区之间连接稀疏。因此,社区发现与聚类是不谋而合的。目前结合社区发现探索时间序列聚类的有关成果仍旧不足[17],标签传播[18]作为简便的划分算法,可以在相连的节点中传播有用的信息。将标签传播应用到时间序列聚类分析,是一种新的研究视角。此外,动态时间弯曲是一种度量准确性高的度量方法,是各个领域理论和应用研究的热点。综上,本文提出一种基于标签传播的时间序列聚类方法来对股指期货套期保值策略进行研究。利用动态时间弯曲构建股票池网络,将每只股票看作网络中独立的节点;通过标签传播方法将网络中的节点划分成不同的社区,实现时间序列聚类;构建最小追踪误差模型,优化每只股票在组合中的权重。
1 相关理论基础标签传播是社区发现领域的重要方法之一,其简单高效的思想使其得到了广泛的关注。动态时间弯曲是泛化能力强的度量方法,是各个领域的研究热点。为了充分理解标签传播方法计算原理,对标签传播和动态时间弯曲的基本原理进行阐述,明确两者在本研究中的重要作用。并给出了优化权重过程所使用的最小追踪误差模型。
1.1 标签传播标签传播(label propagation algorithm, LPA)是根据网络的局部信息结构,利用节点的连接关系自动传播信息,最终得到社区的划分结果。LPA不需要事先指定社区规模、个数、优化目标函数,算法简单又容易实现,同时具有接近线性的时间复杂度。因此,获得了社区发现领域的广泛关注。设无向无权网络
算法 标签传播算法(label propagation algorithm, LPA)
输入 网络
输出
1)
2)
3) 打乱节点顺序,获取打乱顺序的节点集
4) 对于各个节点
5) 若
动态时间弯曲(dynamic time warping, DTW)为一种鲁棒性强的度量方法,最早用于语音识别[19]。与欧氏距离相比,DTW能够弯曲时间轴达到不等长时间序列的度量,充分反映时间序列的形态[20]。特别地,在时间数据挖掘[21]领域中,动态时间弯曲具有较为广泛的应用。利用欧氏距离构建两条时间序列
${\rm{DTW}}(S,Q) = \mathop {\rm min}\limits_P (\sum\limits_{h = 1}^H {{p_h}} )$ | (1) |
式中:
${{\textbf{γ}}} (i,j) = D(i,j) + {\rm min}\left\{ {\begin{array}{*{20}{c}} \!\!\!\!\!\!\!\!\!\! {{\textbf{γ}} (i,j - 1)} \\ {{\textbf{γ}} (i - 1,j - 1)} \\ \!\!\!\!\!\!\!\!\!\! {{\textbf{γ}} (i - 1,j)} \end{array}} \right.$ | (2) |
式(2)表示当前的累积代价为当前距离加上相邻3个最小的累积代价。最后
为得到最优现货组合,从标的指数的
${r_t} = \ln (\sum\limits_{i = 1}^n {{p_{i,t}} \times {w_i}} ) - \ln (\sum\limits_{i = 1}^n {{p_{i,t - 1}} \times {w_i}} )$ | (3) |
标的指数的收益率为
${R_{_t}} = \ln ({P_t}) - \ln ({P_{t - 1}}) = \ln (\frac{{{P_t}}}{{{P_{t - 1}}}})$ | (4) |
式中:
${\rm TE} = \sqrt {\frac{1}{n}\sum\limits_{t = 1}^{{T}} {{{({r_t} - {R_t})}^2}} } $ | (5) |
优化现货投资组合需要结合以下约束:一是资本的预算;二是投资比例的最低要求和最大限制,表示对风险的控制,也可根据投资者的风险偏好来制定。综合目标函数和约束条件,METM如下:
$\begin{aligned} {\rm{min}}{\rm TE}{\rm{ = }}\qquad\quad\quad\quad\quad\quad\qquad\\ \sqrt {\displaystyle\frac{1}{n}\displaystyle\sum\limits_{t = 2}^T {\left[ {\ln (\sum\limits_{i = 1}^n {{p_{i,t}} \times {w_i}} ) - \ln (\sum\limits_{i = 1}^n {{p_{i,t - 1}} \times {w_i}} ) - \ln (\frac{{{P_t}}}{{{P_{t - 1}}}})} \right]} }\\ {\rm{s.t.}}\displaystyle\sum\limits_{i = 1}^n {{w_i} = 1,\alpha \leqslant {w_i} \leqslant \beta ,i = 1,2,···,n} \quad\quad\quad \end{aligned}$ | (6) |
聚类是数据挖掘重要任务之一,根据一定规则将数据划分为若干个簇,簇内的对象保持着高度相似性,簇之间的对象尽可能不同。在金融领域中,聚类对于板块分析、投资组合分析有着重要的意义。LPA的便捷高效,应用到时间序列聚类则是一种新兴尝试。为了充分反映时间序列的网络空间结构,并且能够根据时间序列之间的关系相互影响自动实现聚类,提出一种基于标签传播的时间序列聚类方法,并将其应用到股指期货套期保值优化策略中。
算法 标签传播时间序列聚类(time series clustering based on label propagation, TCLP)
输入 时间序列数据集;
输出 时间序列聚类结果。
1) 对时间序列进行标准化处理;
2) 利用DTW度量每条时间序列之间的距离;
3) 将时间序列视为节点,指定距离阈值
4) 利用LPA对时间序列网络进行划分,实现聚类。
为消除量纲对度量的影响,对时间序列进行标准化处理。例如,由于公司资本大小不同,经营业绩不同等原因,使得不同股票价格千差万别,若没有对股票价格进行标准化处理,那么时间序列相似性度量结果则不准确,使得两条形态很相似的时间序列之间的距离很大。步骤2) 是时间序列相似性度量,DTW可以实现数据点“一对多”匹配,从而实现两条不等长时间序列的度量。由于股票数据存在可能停牌、数据缺失、数据错误等原因,利用DTW可以有效地解决数据清洗之后带来的序列不等长问题。使用LPA之前,需要创建好时间序列网络,而构建时间序列网络的方法通常有两种,即
为验证方法的有效性,通过对Keogh教授提供CBF数据集[22]进行检验。CBF,即Cylinder-Bell-Funnel,它是一种人工数据集,每个类的数据都是服从标准正态分布的噪声加上一个不同于每个类的偏移量。CBF共有3个类别,每个类别时间序列形态各不相同,可以直观地体现聚类效果。度量每条时间序列之间的DTW距离,以时间序列
Download:
|
|
图1(a)左侧给出CBF三种形态相异的时间序列,右侧是对应形态的时间序列的ID。图1(b)展示新方法的聚类效果,用不同颜色代表不同簇,发现新方法能够成功将数据集划分为3个簇。具体分析图1(b),尽管节点3与节点1相连且真实情况也同属一类,但是节点3与更多的绿色节点相连,导致在标签传播过程中节点3接受了绿色节点传过来的标签,被划分到了绿色的簇中。同理可分析节点13和节点24。度量方法和网络构建方式影响着时间序列空间网络布局和最后的聚类结果。然而,通过构建时间序列网络并利用标签传播方法实现聚类,为聚类分析提供了一种新的研究模式。
3.2 实证分析为了检验方法的真实有效性,使用真实股票数据来进行实证分析。采用金融行业的股票数据,采用2014年1月2日至2014年12月31日沪深300股成分股的日收盘价作为实验数据,数据从锐思数据库下载获得。为保证后续挖掘的质量,对数据进行清洗,即剔除7天以上未开盘股票,默认值为该股票的平均收盘价。数据清洗后剩余265只股票数据。
对每只股票数据做标准化处理:
${Y_i} = \frac{{{X_i} - \mu }}{\delta }$ | (7) |
式中:
利用基于TCLP时间序列聚类方法,对实验股票数据集进行聚类。TCLP不必先指定簇数,能依据度量公式的特点自行划分数据集。本实验在构建股票网络结构之前,先利用DTW度量每条时间序列的距离,接着由前50相似序列距离之和的均值确定距离阈值,最后利用标签传播方法对股票网络结构进行划分,划分结果如图2。
Download:
|
|
TCLP-DTW将股票数据集划分为了4个形态不同的股票簇,并且簇内的股票时间序列表现出类似形态,而不同股票簇之间的形态有着明显的区别,说明聚类效果较好。
从每个股票簇中选定追踪组合成分股,并通过METM计算每只追踪成分股的权重,得出追踪误差。首先,从每个股票簇中确定成分股。选择能呈现该股票簇走势的几只股票,在建立组合时表征了整体情况,分散风险。
各个股票簇的股票依据在沪深300指数中的比重从大到小排列,采用前两只比重大的股票,获得8只股票,选取结果如图3所示。
Download:
|
|
可以发现,各个簇选择股票的走势是相似的,而不同股票簇的成分股走势也不同,充分体现
其次,利用METM计算每只成分股的权重。观察4种投资约束条件下TCLP-DTW得到的现货组合追踪误差情况。
如表1所示,得到了TCLP-DTW构建优化现货组合的TE。尽管不同约束条件下能够获得最小的TE,但是部分股票的权重占比很大,不利于风险分散。此外,尽管前面两种约束条件得到的追踪误差并非最低的,然而贵州茅台等股票价格较高,因此在组合中所占权重相对少一些。TCLP-DTW构建优化现货组合是比较倾向于低投资比例的,并且股价较高的成分股占比较少,有利于投资者控制成本风险。
由于K-means聚类方法简单、高效,成为金融领域中应用得较多的聚类方法。为对比K-means和TCLP-DTW创建投资组合的效果,利用K-means将股票池划分成4个簇,同时对比同一条件下双方的
Download:
|
|
从图4中发现,K-means方法在投资比例不设限的条件下同样得到最低的TE。由于投资比例不设限,就能在更大的范围中搜索最优解。对比TLPC-DTW、K-means和随机这3种选股方式所构建优化组合的TE。随机选股为对于各个约束比例模拟追踪10次,每回随机抽取8只股票创建组合。利用METM得到股票的权重,以10次TE的均值当做该种投资约束下的TE。尽管K-means在初始化时随机选取簇中心,但整个聚类过程是一个迭代优化的过程,使得每个簇中对象尽可能相似,簇间的对象尽可能相异,这与随机选取方式有着本质的区别。为了比较聚类方式构建最优现货组合和随机构建最优现货组合,利用
${\rm T{E}_{others}} = \frac{{{E_{{\rm TLPC - NSM}}} - {E_{{\rm others}}}}}{{{E_{{\rm others}}}}}$ | (8) |
作为相对追踪误差率。
从表2给出的TE以及TER发现,TCLP-DTW构建的现货组合在3种投资比例约束条件下现货组合的TE较于K-means、随机方法的TER均为负数,表明误差均得到了改进。K-means只对等长的时间序列聚类,若股票长度不同,则先预处理达到等长效果。K-means的簇中心是簇内的对象的平均值,不是实际中的股票。TCLP-DTW是根据股票的相互影响进行划分,并能从中得到股票连接关系。随机方式的随机性容易存在收益波动大的情况。通过对比发现,TCLP-DTW构建的组合得到更小和稳定的TE,为研究套期保值提供新的视角。
4 结束语本文提出一种标签传播时间序列聚类方法,使用动态时间弯曲能够较好地度量时间序列之间的距离,结合距离阈值构建反映时间序列之间关联关系的网络,再利用标签传播来实现新方法。该方法能够较好地用于股票聚类,用来选择代表股票以确定投资组合,并通过最小误差追踪模型优化组合中的股票权重。实证分析得出,本文方法对追踪误差有一定优化,为进一步了解市场规律、提高投资效率提供一定的技术支撑。
[1] |
柴尚蕾, 郭崇慧, 徐旭. 股指期货套利中的最优现货组合构建策略研究[J]. 运筹与管理, 2012, 21(2): 154-161. CHAI Shanglei, GUO Chonghui, XU Xu. Optimal spot portfolio construction strategy in stock index futures arbitrage[J]. Operations research and management science, 2012, 21(2): 154-161. DOI:10.3969/j.issn.1007-3221.2012.02.023 (0) |
[2] |
郑尊信. 股指期货持有成本模型的修正与比较[J]. 哈尔滨工业大学学报, 2009, 41(2): 248-250. ZHENG Zunxin. Modification and comparison of cost-of-carry model of index futures[J]. Journal of Harbin Institute of Technology, 2009, 41(2): 248-250. DOI:10.3321/j.issn:0367-6234.2009.02.063 (0) |
[3] |
韩立岩, 任光宇. 基于已实现二阶矩预测的期货套期保值策略及对股指期货的应用[J]. 系统工程理论与实践, 2012, 32(12): 2629-2636. HAN Liyan, REN Guangyu. Hedging strategy with futures based on prediction of realized second moment: An application to stock index futures[J]. Systems engineering—theory and practice, 2012, 32(12): 2629-2636. DOI:10.3969/j.issn.1000-6788.2012.12.005 (0) |
[4] | HOU Yang, LI S. Hedging performance of Chinese stock index futures: an empirical analysis using wavelet analysis and flexible bivariate GARCH approaches[J]. Pacific-basin finance journal, 2013, 24: 109-131. DOI:10.1016/j.pacfin.2013.04.001 (0) |
[5] | SU E D. Stock index hedging using a trend and volatility regime-switching model involving hedging cost[J]. International review of economics and finance, 2017, 47: 233-254. DOI:10.1016/j.iref.2016.10.016 (0) |
[6] | TRABELSI N, NAIFAR N. Are Islamic stock indexes exposed to systemic risk? Multivariate GARCH estimation of CoVaR[J]. Research in international business and finance, 2017, 42: 727-744. DOI:10.1016/j.ribaf.2017.07.013 (0) |
[7] |
苏治, 蔡腾腾, 马泽伟. 一种改进的不完全指数复制方法[J]. 数量经济技术经济研究, 2013, 30(6): 149-160. SU Zhi, CAI Tengteng, MA Zewei. An improved solution for incomplete index tracking problem[J]. The journal of quantitative & technical economics, 2013, 30(6): 149-160. (0) |
[8] |
倪禾. 基于启发式遗传算法的指数追踪组合构建策略[J]. 系统工程理论与实践, 2013, 33(10): 2645-2653. NI He. Heuristic genetic algorithm for optimizing an index tracking portfolio[J]. Systems engineering—theory and practice, 2013, 33(10): 2645-2653. DOI:10.12011/1000-6788(2013)10-2645 (0) |
[9] |
胡春萍, 薛宏刚, 徐凤敏. 基于时间加权SVM的指数优化复制模型与实证分析[J]. 系统工程理论与实践, 2014, 34(9): 2193-2201. HU Chunping, XUE Honggang, XU Fengmin. An stock index replicating model based on time weighted SVM and it’s empirical analysis[J]. System engineering—theory and practice, 2014, 34(9): 2193-2201. (0) |
[10] |
刘睿智, 周勇. 指数跟踪投资组合与多信息下指数可预测性——基于Adaptive LASSO和ARIMA-ANN方法[J]. 系统工程, 2015, 33(4): 1-7. LIU Ruizhi, ZHOU Yong. The portfolio of index tracing and index predictability under multi-information—Based on adaptive LASSO and ARIMA-ANN method[J]. Systems engineering, 2015, 33(4): 1-7. DOI:10.3969/j.issn.1009-6744.2015.04.001 (0) |
[11] | CHEN Chen, KWON R H. Robust portfolio selection for index tracking[J]. Computers and operations research, 2012, 39(4): 829-837. DOI:10.1016/j.cor.2010.08.019 (0) |
[12] | GUASTAROBA G, SPERANZA M G. Kernel search: An application to the index tracking problem[J]. European journal of operational research, 2012, 217(1): 54-68. DOI:10.1016/j.ejor.2011.09.004 (0) |
[13] | FILIPPI C, GUASTAROBA G, SPERANZA M G. A heuristic framework for the bi-objective enhanced index tracking problem[J]. Omega, 2016, 65: 122-137. DOI:10.1016/j.omega.2016.01.004 (0) |
[14] | DOSE C, CINCOTTI S. Clustering of financial time series with application to index and enhanced index tracking portfolio[J]. Physica A: statistical mechanics and its applications, 2005, 355(1): 145-151. DOI:10.1016/j.physa.2005.02.078 (0) |
[15] | NANDA S R, MAHANTY B, TIWARI M K. Clustering Indian stock market data for portfolio management[J]. Expert systems with applications, 2010, 37(12): 8793-8798. DOI:10.1016/j.eswa.2010.06.026 (0) |
[16] | LEMIEUX V, RAHMDEL P S, WALKER R, et al. Clustering techniques and their effect on portfolio formation and risk analysis[C]//Proceedings of the International Workshop on Data Science for Macro-Modeling. Snowbird, UT, USA, 2014: 1–6. (0) |
[17] | FERREIRA L N, ZHAO Liang. Time series clustering via community detection in networks[J]. Information sciences, 2016, 326: 227-242. DOI:10.1016/j.ins.2015.07.046 (0) |
[18] | RAGHAVAN U N, ALBERT R, KUMARA S. Near linear time algorithm to detect community structures in large-scale networks[J]. Physical review E: covering statistical, nonlinear, biological, and soft matter physics, 2007, 76(3): 036106. (0) |
[19] | SAKOE H, CHIBA S. Dynamic programming algorithm optimization for spoken word recognition[J]. IEEE transactions on acoustics, speech, and signal processing, 1978, 26(1): 43-49. DOI:10.1109/TASSP.1978.1163055 (0) |
[20] |
李海林, 梁叶. 基于数值符号和形态特征的时间序列相似性度量方法[J]. 控制与决策, 2017, 32(3): 451-458. LI Hailin, LIANG Ye. Similarity measure based on numerical symbolic and shape feature for time series[J]. Control and decision, 2017, 32(3): 451-458. (0) |
[21] | LI H. Accurate and efficient classification based on common principal components analysis for multivariate time series[J]. Neurocomputing, 2016, 171: 744-753. DOI:10.1016/j.neucom.2015.07.010 (0) |
[22] | CHEN Yanping, KEOGH E, HU Bing, et al. The UCR time series classification archive[EB/OL]. 2015. (2015–07-01)[2015–12-01]. http://www.cs.ucr.edu/~eamonn/time_series_data/. (0) |