本文在多种金融资产的收益率服从混合不确定多元分布的假设下,建立鲁棒投资组合模型。首先,使用GARCH-EVT模型描述单一金融资产收益率分布的厚尾、异方差特征,然后运用Copula理论描述多个收益率之间的相依结构,建立了描述多种金融资产收益率分布的GARCH-EVT-Copula模型。最后,运用不同的Copula函数构建了收益率的混合不确定多元分布集合,并在WCVaR-Copula鲁棒投资组合模型基础上,建立了GARCH-EVT-Copula-WCVaR鲁棒投资组合模型。通过与经典的均值方差模型、以正态分布刻画边缘分布的Normal-Copula-WCVaR鲁棒模型的实证比较,可以发现,在股价出现极端情况的股灾期间和收益率大幅波动期间,GARCH-EVT-Copula-WCVaR鲁棒模型所建立的投资组合的回报都要高于均值方差模型和Normal-Copula-WCVaR鲁棒模型。 In this paper, a robust portfolio model is established under the assumption that the underlying distribution is mixed and uncertain. Firstly, we use the GARCH-EVT model to describe the fat tail and heteroscedasticity characteristics of single financial asset return. Then we use Copula to describe the dependence structure between the yields, and establish the GARCH-EVT-Copula model. Finally, a mixture of uncertain multivariate distributions of returns is constructed using different Copula functions, and a robust portfolio model of GARCH-EVT-Copula-WCVaR is estab-lished on the basis of WCVaR-Copula robust model. Compared to the Normal-Copula-WCVaR robust model and the classical mean variance model in the experimental study, the return of GARCH-EVT-Copula-WCVaR robust model is higher than that of Normal-Copula-WCVaR robust model and mean variance model during the stock market crisis and the period when the volatility of return fluctuates.
陆家涛
北京理工大学数学与统计学院,北京
收稿日期:2018年2月5日;录用日期:2018年2月20日;发布日期:2018年2月28日
本文在多种金融资产的收益率服从混合不确定多元分布的假设下,建立鲁棒投资组合模型。首先,使用GARCH-EVT模型描述单一金融资产收益率分布的厚尾、异方差特征,然后运用Copula理论描述多个收益率之间的相依结构,建立了描述多种金融资产收益率分布的GARCH-EVT-Copula模型。最后,运用不同的Copula函数构建了收益率的混合不确定多元分布集合,并在WCVaR-Copula鲁棒投资组合模型基础上,建立了GARCH-EVT-Copula-WCVaR鲁棒投资组合模型。通过与经典的均值方差模型、以正态分布刻画边缘分布的Normal-Copula-WCVaR鲁棒模型的实证比较,可以发现,在股价出现极端情况的股灾期间和收益率大幅波动期间,GARCH-EVT-Copula-WCVaR鲁棒模型所建立的投资组合的回报都要高于均值方差模型和Normal-Copula-WCVaR鲁棒模型。
关键词 :CVaR,GARCH,EVT,Copula理论,鲁棒投资组合
Copyright © 2018 by author and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
1952年Markowitz提出的均值方差模型奠定了现代投资组合理论的基础,在资产配置领域取得了广泛的应用。但是,该模型存在对参数估计非常敏感的不足,Black与Litterman (1992) [
在鲁棒投资组合领域,往往通过构建资产收益率的不确定分布来建立鲁棒投资组合模型。Zhu和Fukushima (2009) [
本文在多种金融资产收益率服从混合不确定多元分布的假设下建立鲁棒投资组合模型。首先,运用GARCH-EVT模型描述单一金融资产收益率分布的厚尾、异方差特征。极值理论主要研究金融资产收益的尾部分布,可以很好的描述收益分布的非正态厚尾特征。同时,GARCH模型在刻画金融资产收益率的异方差方面得到了广泛应用。利用极值理论处理极端情形的优势,并结合GARCH模型描述单个金融资产收益率的异方差特征,可以实现对单一金融资产收益率分布的有效描述。Bali和Neftci (2003) [
其次,运用Copula理论描述多种金融资产收益率之间的非对称性和尾部非线性相关性。基于Copula理论构建的多元分布函数,可以有效描述随机变量之间的相关程度,并且能够反映它们之间的相依结构。Jondeau和Rockinger [
将GARCH模型、EVT以及Copula理论结合起来,建立Copula-GARCH-EVT模型可以有效捕捉金融资产收益率的厚尾、异方差以及收益率之间的非对称、非线性相关的特点。李秀敏、史道济(2007) [
综上,本文将运用GARCH-EVT模型刻画单一金融资产的收益率分布,捕捉金融资产收益分布的厚尾、异方差特征;同时通过Copula理论建立金融资产收益率的多元分布,捕捉金融资产收益率之间的非对称、非线性特征。最后,通过使用不同的Copula函数构建收益率的混合不确定多元分布集合,并在WCVaR理论框架下使得最坏情况下的CVaR达到最小来建立鲁棒投资组合模型。
GARCH模型由条件均值方程和条件方差方程组成。GARCH(p,q)模型的一般结构如下:
{ r t = μ t + ε t ε t = σ t ⋅ z t σ t 2 = α 0 + ∑ i = 1 p α i σ t − i 2 + ∑ j = 1 q β j ε t − j 2 α 0 > 0 , α i ≥ 0 , β j ≥ 0 , ∑ α i + β j < 1 (1)
其中, μ t 为t时刻的均值方程,一般采用ARMA(m,n)模型, z t 是独立同分布变量,表示新息过程。各种GARCH模型的区别通常在于条件方差方程采取的形式不同或者是对 z t 的分布假设不同,通常假设 z t 服从正态分布或学生t分布。为了刻画金融资产收益的厚尾特征,本文采用极值分布来描述 z t 。
假设随机变量序列 x 1 , x 2 , ⋯ , x n 独立同分布,来自于同一个未知分布F,用 u 表示某一充分大的阈值。假定大于 u 的样本个数为 N u , X 1 , X 2 , ⋯ , X n 表示大于阈值 u 的样本观测值, Y 1 , Y 2 , ⋯ , Y n 表示相应的超额差值, Y i = X i − u ( i = 1 , 2 , ⋯ , N u ) ,定义 F u 为大于阈值 u 的随机变量的差额损失分布,则
F u ( y ) = P r ( X − u < y | X > u ) = F ( u + y ) − F ( u ) 1 − F ( u ) = F ( x ) − F ( u ) 1 − F ( u )
由Pickands定理,对于充分大的阈值 u ,令 y = x − u ,则超额损失分布函数 F u ( y ) 收敛于如下的广义Pareto分布:
G ξ , β ( y ) = { 1 − ( 1 + ξ β y ) − 1 ξ , ξ ≠ 0 1 − exp ( − y β ) , ξ = 0
其中, ξ 是形状参数, β 是尺度参数。
本文使用对数收益率来表示金融资产的收益, R i , t = ln P i , t − ln P i , t − 1 ,其中 P i , t 和 R i , t 分别表示第i种资产在t时刻的价格和收益率。
使用GARCH-EVT模型刻画第i种资产的收益率分布,首先,根据对数收益率序列 ( R i , t − m + 1 , ⋯ , R i , t − 1 , R i , t ) 建立ARMA-GARCH模型,得到条件均值序列的估计值 ( μ ^ i , t − m + 1 , ⋯ , μ ^ i , t − 1 , μ ^ i , t ) 和条件标准差序列的估计值 ( σ ^ i , t − m + 1 , ⋯ , σ ^ i , t − 1 , σ ^ i , t ) 。由GARCH模型(1)可以得到新息序列:
( z i , t − m + 1 , ⋯ , z i , t − 1 , z i , t ) = ( R i , t − m + 1 − μ ^ i , t − m + 1 σ ^ i , t − m + 1 , ⋯ , R i , t − 1 − μ ^ i , t − 1 σ ^ i , t − 1 , R i , t − μ ^ i , t σ ^ i , t )
其次,采用两端服从GDP分布,中间服从标准正态分布的形式来拟合新息序列。
F ( z ) = { T u L T { 1 + ξ L u L − z β L } − 1 ξ L , z < u L Φ ( z ) , u L ≤ z ≤ u R 1 − T u R T { 1 + ξ R z − u R β R } − 1 ξ R , z > u R (2)
其中, ξ L 是左尾形状参数, β L 是左尾尺度参数, u L 是左尾阈值, T u L 为低于左尾阈值的样本数, ξ R 是右尾形状参数, β R 是右尾尺度参数, u R 是右尾阈值, T u R 为高于右尾阈值的样本数, Φ ( z ) 表示标准正态分布的分布函数。
Copula函数可以分为椭圆族Copula和阿基米德族Copula。椭圆族Copula包括正态Copula函数以及t-Copula函数。椭圆Copula适合描述具有对称相依结构的数据,但是在金融领域,许多金融数据表现出非对称、非线性相关的特点。Genest和Mackay (1986) [
利用GARCH-EVT模型(1)、(2)对投资组合中各金融资产的对数收益率分别构建边缘分布 F i ( z i ) ,然后根据金融资产的相依结构,选择合适的Copula函数 c ( ⋅ ) ,构建多种资产收益率的联合分布 F ( z 1 , ⋯ , z n ) = c ( F 1 ( z 1 ) , ⋯ , F n ( z n ) ) 。
本文通过GARCH-EVT-Copula模型刻画多种金融资产收益率的多元分布,然后运用不同的Copula函数构建资产收益分布的混合不确定集合,在WCVAR鲁棒投资组合模型基础上建立稳健的投资组合模型。
令 x ∈ R n 表示随机收益率向量, w ∈ R n 表示决策向量, f ( w , x ) 表示相应的损失函数。假设 x 服从密度函数为 p ( . ) 的连续分布。
给定决策向量 w ∈ W ,损失 f ( w , x ) 不超过阈值 α 的概率为
Ψ ( w , α ) ≜ ∫ f ( w , x ) ≤ α p ( x ) d x
对于固定的决策向量 w ∈ W 及置信水平 β ,用 V a R β ( w ) 及 C V a R β ( w ) 表示相应的风险度量VaR与CVaR,则
V a R β ( w ) ≜ min { α ∈ R : Ψ ( w , α ) ≥ β }
C V a R β ( w ) ≜ 1 1 − β ∫ f ( w , x ) ≥ V a R β ( w ) f ( w , x ) p ( x ) d x
Rockafellar和Uryasev (2003) [
F β ( w , α ) ≜ α + 1 1 − β ∫ x ∈ R n [ f ( w , x ) − α ] + p ( x ) d x (3)
C V a R β ( w ) = min a ∈ R F β ( w , α )
Kakouris和Rustem (2014) [
结论1:设随机向量 u ∈ [ 0 , 1 ] n , F ( x ) = ( F 1 ( x 1 ) , ⋯ , F n ( x n ) ) 表示边缘累计分布的集合,令 u = F ( x ) ,损失函数为 f ˜ ( w , u ) ,假设 u 服从Copula密度函数为 c ( ⋅ ) 的连续分布。则对于置信度为 β 的VaR及CVaR定义为
V a R β ( w ) ≜ min { α ∈ R : c ( u | f ˜ ( w , u ) ≤ α ) ≥ β }
C V a R β ( w ) ≜ 1 1 − β ∫ f ˜ ( w , u ) ≥ V a R β ( w ) f ˜ ( w , u ) c ( u ) d u
所以,根据(3)式,求 C V a R β ( w ) 可简化为求解如下最小化问题:
Y β ( w , α ) = α + 1 1 − β ∫ u ∈ I n [ f ˜ ( w , u ) − α ] + c ( u ) d u
C V a R β ( w ) = min a ∈ R Y β ( w , α )
本文使用GARCH-EVT模型来刻画收益率的边缘分布, u = F ( x ) = ( F 1 ( x 1 ) , ⋯ , F n ( x n ) ) ,其中 F i ( x i ) 表示由GARCH-EVT模型(2)描述的第i种资产的收益率分布;然后假设Copula函数 c ( u ) 属于一个Copula函数集合 C = ( c 1 ( u ) , ⋯ , c l ( u ) ) ,通过不同的Copula函数 c i ( u ) 建立相应的GARCH-EVT-Copula模型得到刻画收益率的混合不确定多元分布集合 H ( x ) 。通过使此不确定集合下的风险度量CVaR在最坏情况(worst-case CVaR)下达到最优,建立鲁棒投资组合模型。
描述资产收益率向量 x 的混合不确定多元分布集合 H ( x ) :
H ( x ) ≜ { ∑ i = 1 l λ i c i ( F ( x ) ) : λ = ( λ 1 , ⋯ , λ l ) ∈ Λ , c i ∈ C , i = 1 , ⋯ , l }
其中 Λ = { λ = ( λ 1 , ⋯ , λ l ) | ∑ i = 1 l λ i = 1 , λ i ≥ 0 } 。
所以,在混合不确定多元分布集合 H ( x ) 下的风险度量CVaR有如下形式:
C V a R β ( w ) = min a ∈ R Y β ( w , α )
Y β ( w , α ) = α + 1 1 − β ∫ u ∈ I n [ f ˜ ( w , u ) − α ] + ∑ i = 1 l λ i c i ( u ) d u = ∑ i = 1 l λ i Y β i ( w , α )
当收益率的分布是不确定的,可以通过定义不确定分布下最坏的CVaR (worst-case CVaR)来建立鲁棒投资组合模型。根据Zhu和Fukushima (2009) [
定义:对于Copula函数集合C、置信度 β 及固定的决策向量 w ∈ W ,定义WCVaR:
W C V a R β ( w ) ≜ sup c ( . ) ∈ C C V a R β ( w )
结论2:对于任意的 w 及不确定分布集合 H ( x ) ,
W C V a R β ( w ) = min a ∈ R max λ ∈ Λ Y β ( w , α )
所以,在混合不确定集合 H ( x ) 下最小化 W C V a R β ( w ) 的鲁棒投资组合模型为:
min w ∈ W W C V a R β ( w ) = min w ∈ W min a ∈ R max λ ∈ Λ Y β ( w , α )
该模型可以简化成如下形式:
min ( w , α , θ ) ∈ W × R × R { θ : ∑ i = 1 l λ i Y β i ( w , α ) ≤ θ , λ ∈ Λ }
Zhu和Fukushima (2009) [
min ( w , α , θ ) ∈ W × R × R { θ : Y β i ( w , α ) ≤ θ , i = 1 , ⋯ , l } (4)
综上,在Zhu和Fukushima (2009) [
min w ∈ R n , v ∈ R l , α ∈ R , θ ∈ R θ
s . t . w ∈ W , v ∈ R m , α ∈ R , θ ∈ R
α + 1 S i ( 1 − β ) ( 1 i ) T v i ≤ θ , i = 1 , ⋯ , l
v k i ≥ f ˜ ( w , u k i ) − α , k = 1 , ⋯ , S i , i = 1 , ⋯ , l
v k i ≥ 0 , k = 1 , ⋯ , S i , i = 1 , ⋯ , l
其中 u k i 表示从第i个Copula函数 c i 建立的多元分布中抽取的第k个样本, v = ( v 1 , ⋯ , v l ) ∈ R m , m = ∑ i = 1 l S i , 1 i = ( 1 , ⋯ , 1 ) ∈ R S i , W 表示决策向量的限制集合, W = { e T w = 1 , w i ≥ 0 , E ( u T w ) ≥ γ } ,其中
本文从Wind数据库中提取2004/1/1至2016/12/30共3158个交易日的上证180指数的能源行业指数和工业行业指数数据,构建二元投资组合,图1为两种指数在样本期的走势图。使用2004/1/2至2014/12/31的对数收益率数据作为训练样本,使用2015/1/5至2016/12/30的对数收益率数据作为验证样本。训练样本数据包括2008年上证综指从5272点跌至1850点的股灾时期,验证样本数据中包含2015年下半年上证综指从5178点跌至3539点的股灾时期。从图2中可以发现,在2008年和2015年日对数收益率序列
图1. 上证180的能源和工业行业指数数据
图2. 上证180的能源和工业行业指数日对数收益率
波动明显增大,展现出了异方差的特点。
表1列出了能源和工业行业指数日对数收益率的描述统计量以及正态分布的Jarque-Bera检验结果。从表中的数据可知,对于全样本、训练样本以及验证样本,两种行业指数的日对数收益率都是尖峰的,而且Jarque-Bera检验拒绝了序列服从正态分布的假设。
为了检验本文建立的以GARCH-EVT为边际分布的GARCH-EVT-Copula MCVaR鲁棒模型(G-E-Copula-WCVaR)的效果,我们使用经典的均值方差模型(MV-Model)和以正态分布为边际分布的Normal-Copula WCVa鲁棒模型(N-Copula-WCVaR)作为对照,比较三种模型的累计对数收益率。在两种鲁棒投资组合模型中,本文使用t-Copula、Gumbel Copula以及Clayton Copula来构建资产收益率的混合不确定二元分布集合,并进行2000次模拟实验。关于权重向量的限制集合W,本文不允许卖空。
静态投资组合策略以训练样本(2004/1/2~2014/12/31)建立投资组合模型,得到相应的投资组合权重,
行业 | 均值 | 标准差 | 峰度 | Jarque-Bera | |
---|---|---|---|---|---|
全样本期 2004/1/2~2016/12/30 | 能源 | 0.0001 | 0.0213 | 5.756 | 0.0000 |
工业 | 0.0002 | 0.0199 | 6.461 | 0.0000 | |
训练样本期 2004/1/2~2014/12/31 | 能源 | 0.0002 | 0.0209 | 5.366 | 0.0000 |
工业 | 0.0003 | 0.0186 | 6.250 | 0.0000 | |
验证样本期 2015/1/5~2016/12/30 | 能源 | −0.0006 | 0.0230 | 6.955 | 0.0000 |
工业 | −0.0002 | 0.0258 | 5.483 | 0.0000 |
表1. 日对数收益率序列的描述统计量及JB检验
并在整个验证期(2015/1/5~2016/12/30)固定该投资组合权重不变。
上证180能源行业指数收益率分布的GARCH-EVT模型:
{ r 1 , t = − 0.0313 r 1 , t − 1 − 0.0311 r 1 , t − 2 − 0.0048 r 1 , t − 3 + 0.0336 r 1 , t − 4 + 0.0714 ϵ 1 , t − 1 + ε 1 , t ε 1 , t = σ 1 , t ⋅ z 1 , t , z 1 , t ~ F 1 ( z ) σ 1 , t 2 = 0.0575 σ 1 , t − 1 2 + 0.9338 ε 1 , t − 1 2
F 1 ( z ) = { 0.0959 ⋅ ( 1 − 0.1629 ⋅ − 1.2076 − z 0.7225 ) 1 0.1629 , z < − 1.2076 Φ ( z ) , − 1.2076 ≤ z ≤ 1.2172 0.9000 ⋅ ( 1 + 0.0325 ⋅ z − 1.2172 0.6157 ) − 1 0.0325 , z > 1.2172
上证180工业行业指数收益率分布的GARCH-EVT模型:
{ r 2 , t = − 0.6649 r 2 , t − 1 + 0.0018 r 2 , t − 2 + 0.7050 ϵ 2 , t − 1 + ε 2 , t ε 2 , t = σ 2 , t ⋅ z 2 , t , z 2 , t ~ F 2 ( z ) σ 2 , t 2 = 0.0558 σ 2 , t − 1 2 + 0.9352 ε 2 , t − 1 2
F 2 ( z ) = { 0.0997 ⋅ ( 1 + 0.0825 ⋅ − 1.2190 − z 0.5814 ) − 1 0.0825 , z < − 1.2190 Φ ( z ) , − 1.2190 ≤ z ≤ 1.2074 0.9003 ⋅ ( 1 + 0.0184 ⋅ z − 1.2074 0.5169 ) − 1 0.0184 , z > 1.2074
从表2中三种模型结果的比较来看,三种模型在收益及夏普比率上没有明显差异,而且在测试期间发生了股灾的极端情况,三种模型的日均收益率和累计收益率均为负。考虑到静态投资组合策略没有及时对组合权重进行调整,不能充分体现出模型间的差异,下面,本文将构建动态的投资组合策略对三种不同模型进行比较分析。
动态投资组合策略在验证期间的每一交易日结束后,将该日数据加入到训练样本中,滚动地建立投资组合模型,共得到488个交易日的每日投资组合权重。
图3给出了三种不同模型所构建的动态投资组合在验证期2015~2016年的累计回报率,GARCH-EVT-Copula-WCVaR鲁棒模型的累计收益率要明显高于Normal-Copula-WCVaR模型和MV-Model;由于鲁棒模型的保守性,在2015上半年的牛市中,两种鲁棒模型的累计收益率均未超过经典的均值方差模型的累计收益率;由于鲁棒模型稳健性,在包含2015年下半年股灾的整个测试期间,两种鲁棒模型的回报都高于经典的均值方差模型,且GARCH-EVT-Copula-WCVaR模型的累计对数收益率最大。
表3给出了三种模型所建立的动态投资组合的收益情况,从平均收益、总收益率、95%的VaR及夏普比率比较,GARCH-EVT-Copula-WCVaR鲁棒模型都是最优的。
图4给出了股灾期间GARCH-EVT-Copula-WCVaR鲁棒模型和Normal-Copula-WCVaR鲁棒模型的累计收益率。从图4中可以看出,2015年下半年股灾期间,GARCH-EVT-Copula-WCVaR鲁棒模型所建立的投资组合收益率要略高于用正态分布刻画边缘分布的Normal-Copula-WCVaR鲁棒模型;GARCH-EVT-Copula-WCVaR鲁棒模型的平均对数收益率及累计对数收益率分别为−0.0027、−0.3934,Normal-Copula-WCVaR鲁棒模型的平均对数收益率及累计对数收益率分别为−0.0029、−0.4295。这体现出GARCH-EVT模型能够捕捉到收益率分布厚尾特征的优势。
图5分别展示了2016年1月至2016年5月间,两种指数的对数收益率,可以看出1月至2月两种指数对数收益率的波动明显高于3月至5月。能源指数的对数收益率在两个阶段的标准差分别为0.0318
权重 w | 日均收益率 | 累计收益率 | 夏普比率 | |
---|---|---|---|---|
MV model | ( 0.2404 , 0.7596 ) T | −2.98 × 10−4 | −0.1456 | −0.0122 |
Normal-Copula-WCVaR | ( 0.4271 , 0.5729 ) T | −3.72 × 10−4 | −0.1815 | −0.0157 |
GARCH-EVT-Copula-WCVaR | ( 0.3688 , 0.6312 ) T | −3.49 × 10−4 | −0.1703 | −0.0146 |
表2. 三种模型构建的静态投资组合策略
平均对数收益率 | 累计收益 | VaR0.95 | 标准差 | 夏普比率 | |
---|---|---|---|---|---|
MV model | −0.30 × 10−3 | −0.142 | −0.045 | 0.0241 | −0.0120 |
Normal-Copula-WCVaR | −0.11 × 10−3 | −0.055 | −0.046 | 0.0246 | −0.0046 |
GARCH-EVT-Copula-WCVaR | 0.06 × 10−3 | 0.030 | −0.041 | 0.0237 | 0.0026 |
表3. 三种模型构建的投资组合收益比较
图3. 三种模型的累计对数收益率
图4. 股灾期间两种鲁棒模型的累计对数收益率
图5. 2016/1~2016/5上证180能源指数及工业指数对数收益率
图6. 2016/1~2016/5两种鲁棒模型的累计对数收益率
和0.0135,工业指数的对数收益率在两个阶段的标准差分别为0.0350和0.0143。两种指数的对数收益率在两个阶段都展示出了异方差的特征。
从图6中可以看出在2016/1~2016/5两种指数对数收益率的波动明显发生变化期间,GARCH-EVT-Copula-WCVaR模型的累计对数收益率要明显高于仅用正态分布刻画边缘分布的Normal-Copula -WCVaR模型。这体现出GARCH-EVT模型能够捕捉到收益率异方差特征的优势。
本文先使用GARCH-EVT模型来刻画单一金融资产收益率的分布,然后使用Copula理论建立多种资产收益率的多元分布,并通过运用三种不同的Copula函数构建了金融资产收益率的混合不确定多元分布集合。最后,在该混合不确定多元分布集合上定义WCVaR风险度量,建立了鲁棒投资组合模型。通过与均值方差模型和仅用正态分布刻画收益率边缘分布的Copula-WCVaR鲁棒模型的实证对比,在资产价格出现极端情况的股灾期间及资产收益率波动发生变化期间,本文建立的GARCH-EVT-Copula-WCVaR鲁棒模型都要优于Normal-Copula-WCVaR鲁棒模型和均值方差模型。
陆家涛. 基于GARCH-EVT-Copula的WCVaR鲁棒投资组合模型 Robust Portfolio Selection with GARCH-EVT-Copula[J]. 统计学与应用, 2018, 07(01): 54-64. http://dx.doi.org/10.12677/SA.2018.71008