1. 引言
经济全球化的快速发展使得风险在国家之间的传导和外溢成为一种普遍现象。在全球一体化以及金融国际化的过程中,风险溢出效应也已经超越了单一国家的金融市场,成为国际金融体系中不可忽视的因素,一个国家金融市场的波动不再仅仅受制于自身内部因素,同时也会受到其他国家金融市场的影响。
2023年3月,以硅谷银行破产为代表的美国新一轮银行破产危机给欧洲、亚洲等诸多国家银行业带来了巨大冲击,面对日渐频发的银行破产事件及银行业危机,加强对国际银行市场间风险外溢效应的研究具有重要的实践意义。从国内外研究现状来看,风险外溢的研究主要关注点还是在国内各金融子市场,而对国际间金融风险溢出效应的研究则相对较少。对国际银行市场间风险溢出效应方面的研究不足,很可能导致金融管理当局低估银行业风险,对宏观金融风险把握不准确,从而使得宏观审慎监管政策不够完善;导致银行业对国际间金融风险溢出重视程度不够,缺乏相关风险管理观念,缺少相关风险管方法和风险管理措施;导致投资者或储蓄者对银行业存在的风险认识不充分,对银行业的监督不够。鉴于此,本文利用中美银行业指数,对中美银行市场间金融风险溢出效应进行定量研究,以确定风险溢出是否存在以及其溢出强度,并分析相关风险溢出路径,从而弥补相关领域研究的不足,为中国银行业的风险管理提供依据和方法,也为监管当局的宏观审慎监管提供参考。
2. 文献综述
随着中国金融市场的深入改革和全面开放,金融市场正经历着剧烈变迁,我国金融市场与国际金融市场之间的联系日益紧密,国际金融市场间风险溢出效应的研究成为近年来金融风险研究的重点领域之一。在有关国际间风险溢出效应的实证研究中,所使用的度量方法主要包括格兰杰因果关系检验,MES,VAR向量自回归模型,CCA模型、CoVaR等,CoVaR方法是其中应用较广泛的一种。CoVaR方法由Adian和Brunnermeier (2008) [1]第一次提出,其最大的优点在于克服了以前用方差来间接测度风险的不足,国内外学者常将其与分位数回归,GARCH模型,Copula函数相结合构建不同的模型来研究风险溢出效应。
目前国内外有不少学者使用CoVaR方法研究国际间风险溢出效应。Reboredo和Ugolini (2015) [2]采用CoVaR-Copula模型分析了欧洲国家主权债务的系统性风险以及相互影响关系,Drakos和Kouretas (2015) [3]运用CoVaR模型研究了外国银行对美国系统性风险的影响程度。沈虹和邢荧(2017) [4]使用分位数回归-CoVaR方法研究了中美股票市场之间的风险溢出效应,发现中美股市之间的风险溢出效应显著且呈现不断上升的趋势。李强等(2020) [5]构建DCC-TGARCH-CoVaR模型对国际原油期市场与我国汇率、期货及股市间的风险溢出关系进行了实证,发现溢出效应显著。史嘉俊和叶李伟(2021) [6]采用广义CoVaR模型对国际主要股票市场之间的极端风险溢出效应进行检验,发现欧美成熟股票市场是风险的主要输出者,而亚洲股市是主要的风险输入者。江一帆(2021) [7]基于GARCH-Copula-CoVaR模型考察了“一带一路”沿线国家股票市场之间的风险溢出效应,结果显示各国股市之间存在正向的风险传染性。罗长青等(2021) [8]通过构建Copula-CoVaR模型,研究了中国股票市场与“一带一路”沿线国家股票股票市场之间的风险相依性,结果发现中国对各单一国家股票市场的风险溢出显著。钟雯和李保民(2021) [9]借助GARCH-Copula-CoVaR模型对中国与东盟金融市场之间风险溢出效应进行实证研究,发现相互之间存在显著的非对称的溢出效应。荆中博等(2022) [10]使用双重ΔCoVaR模型测度了跨境资本周期性波动对我国商业银行的风险溢出效应。Li Zhaohong等(2022) [11]基于GARCH-CoVaR模型分析了国际原油市场对中国汇率市场的风险溢出效应,结果显示国际原油市场与在岸汇率市场之间存在向下的风险溢出。王倩等(2022) [12]利用DCC-GARCH-CoVaR模型对中外原油期货市场之间的相依关系进行了实证,结果表面两者存在非对称的双向风险溢出。王喜平和王雪萍(2022) [13]基于时变Copula-CoVaR模型研究中欧碳市场之间的风险溢出效应,实证结果显示两者之间的风险溢出显著。张松艳和杜明娅(2022) [14]构建了DCC-GARCH-CoVaR模型来研究国际原油期货对中国农产品期货市场的风险溢出效应,发现国际原油期货市场的波动将引起我国各农产品期货市场不同程度的波动。Zhao Jing等(2023) [15]使用GARCH-EVT-Copula-CoVaR方法研究了国际石油市场对中国股票市场的风险溢出情况,实证结果显示原油市场对中国股票市场各行业均存在正向的风险溢出效应。
从以上国内外研究现状可以看出,关于国际间金融风险溢出效应的研究近年来成为了金融风险研究领域的重点之一,越来越多的学者致力于在这方面做出贡献。同时也可以看到学者们的关注点主要在两个方面,一是模型的优化,通过不断寻找更好的模型以取得更准确的实证结果,二是研究对象的细化,从研究各国整个金融市场之间的相关性,到研究各国金融子市场甚至是各国具体产业之间的风险溢出问题。尽管采用GARCH-Copula-CoVaR模型研究国际间风险溢出效应已逐渐成为一种趋势,但也不难看出使用该方法研究国际间银行业风险溢出的文献确比较少,因此本文将以中美银行业指数2018年到2022年的日收益数据为样本,构建GARCH-Copula-CoVaR模型研究美国银行业对中国银行业的风险溢出情况。
3. 模型介绍
3.1. CoVaR的定义
条件风险价值(Conditional Value at Risk,简称CoVaR)由Adrian和Brunnermeier首次提出,它将风险溢出效应纳入到了VaR框架之内,解决了传统VaR不能捕捉风险在金融机构间或金融市场之间的溢出效应这个问题。按照Adrian和Brunnermeier的定义,
表示r1关于r2的条件风险价值,即当资产r2 (一般用收益率表示,下同)处于风险水平时,资产r1所面临的风险水平,具体公式:
(3.1)
其中,
为置信水平,
表示滞后p期对应的风险价值。r2代表的是t时刻r1的总风险,它包括了无条件风险价值以及溢出风险价值,因此,风险溢出价值为:
(3.2)
因为不同金融市场之间的无条件风险价值相差可能很大,
通常不能充分反映两个金融市场之间风险溢出的强度,所以我们将其标准化:
(3.3)
3.2. GARCH模型
人们在研究经济和金融时间序列的过程中,发现不仅序列的均值存在时变的特性,方差也具有时变特性。Engel (1982) [16]提出了自回归条件异方差(ARCH)时间序列模型,规定当期的条件方差是q期滞后扰动的线性函数,ARCH模型的经济含义明确,而且可以很好的刻画市场的波动,因此,ARCH模型一经提出就受到了学术界的重视,并迅速的被应用于金融研究的各国领域。但是,ARCH模型存在滞后期数过大的问题,为解决该问题,Bollerslev [17]在ARCH基础上提出了GARCH模型,GARCH是对ARCH模型的一种发展,它有效地降低了ARCH模型中的滞后阶数,使估计变得更加准确。GARCH模型表达式如下:
(3.4)
(3.5)
(3.6)
其中,
,
,
,
,
,
,
为已知信息集
下
的条件方差,公式(3.4)代表的是条件均值方程,而公式(3.5)代表了条件方差。在GARCH模型中,一般假设
是服从正太分布的,但是金融时间序列往往并非正太分布,而是呈现尖峰厚尾的特征,因此在实证中,学者往往用t分布或GED分布代替正太分布,我们这里将采用t分布。
3.3. Copula函数
Copula函数最早是由Sklar (1959) [18]提出,它将变量的边缘分布函数连接在一起,是一种用于多元相关性分析函数,Copula函数的最大特点是可以将一个联合分布函数拆成一元的边缘分布和由Copula函数表示的相依结构两部分,从而为求多个单个变量联合分布函数提供更便捷的途径。Copula函数包含多个函数族,但在金融分析中,椭圆族和阿基米德族运用最为广泛。因此,本文将根据与实际数据的拟合情况,主要考虑从ClaytonCopula、FrankCopula、GaussianCopula以及GumbelCopula中寻找最优的拟合函数。
3.4. CoVaR计算
首先,利用GARCH-t模型分别对
和
进行拟合,获取残差序列
和
,它们均服从均值为0,方差为1,自由度为v1和v2的标准t分布。同时求出
和
的VaR (无条件风险价值):
(3.7)
其中,其中,
代表置信水平,
表示自由度为v的标准化t分布的分位数函数,
表示自由度为v的t分布的分位数函数,
为GARCH-t模型的条件标准差。
其次,根据条件分布的定义,
关于
的条件分布为:
(3.8)
由Copula函数有:
(3.9)
根据Sklar定理,可以得到联合分布函数所对应的密度函数:
(3.10)
其中,
为Copula函数的密度函数,
是边缘分布
的密度函数。因此,可将公式(3.8)简化为:
(3.11)
将
的
分位数
替代上式中的
,并令表达式恒等于
,有:
(3.12)
方程中只有一个未知数v1,求解该方程,得到v1的值。假设求解得
,则q1即为
取值
分位数(q2)的条件下,v1的
分位数,从而得到:
(3.13)
进而易求
和
。
4. 实证研究
4.1. 数据选取
本文选取2018年1月2日至2022年12月30的中证银行指数、道琼斯银行指数、纳斯达克银行指数和KBW银行指数的日收盘价作为原始数据,其中,中证银行指数数据来自于Tushare数据库,而道琼斯银行指数、纳斯达克银行指数和KBW银行指数的数据来源于英为财情网站(investing.com),对于其中的缺失数据,采取直线插值法予以填补。
4.2. 描述性统计
首先,我们通过对中证银行指数、道琼斯银行指数、纳斯达克银行指数和KBW银行指数的日收盘数据分别取对数一阶差分,从而得到四个指数的日收益率序列,为减少误差,我们对将所有结果都乘以100,即:
其中,Rt表示t日的百分比收益率,pt表示银行指数在日的收盘价格。然后,分别对四个指数的日百分比收益率序列进行描述性统计,结果见表1。
尽管四个银行指数收益率序列的偏度系数接近于0,但峰度系数均显著大于3,这说明各指数收益率序列均比正态分布更为集中,呈尖峰状态。Jarque-Bera检验结果概率值为0,即在1%显著水平下各股票指数收益率序列都显著异于正态分布,据此可以初步判断两个阶段各股票指数收益率序列均不服从正态分布。
为更好地了解四个银行指数收益率序列的特征,我们逐一绘制出他们的时间序列图,具体见图1:
Table 1. Descriptive statistics of basic information
表1. 基本信息描述性统计
序列 |
ZZYHZS |
DJ |
KBW |
NASDAQ |
Count |
1303 |
1303 |
1303 |
1303 |
Mean |
−0.011339 |
−0.002850 |
−0.004622 |
0.000991 |
Std |
1.184633 |
2.035473 |
2.074698 |
1.981586 |
Min |
−4.627665 |
−16.536890 |
−17.678360 |
−14.031583 |
Max |
8.648415 |
14.428085 |
13.816027 |
12.473007 |
Skewness |
0.502003 |
−0.265055 |
−0.371218 |
−0.199424 |
Kurtosis |
3.804713 |
11.630138 |
10.870935 |
7.579360 |
JB Statistic |
840.644978 |
7358.753886 |
6445.964656 |
3127.513369 |
JB p-value |
0.000000 |
0.000000 |
0.000000 |
0.000000 |
Figure 1. Time series changes in index returns of various banks
图1. 各银行指数收益率序列时间变化图
从上面四个图中我们可以看到,四个指数收益率序列的变动是比较平稳的,但是呈现出比较明显的波动丛集性,也就是大的波动后面往往伴随着大的波动,小的波动后面往往伴随着小的波动。此外,我们还可以发现美国三大银行指数道琼斯银行指数、KBW银行指数、NASDAQ银行指数收益率序列的时序图形状比较相似,这说明三大指数都很好地反映了美国银行业的波动,而中证银行指数收益率序列的时序图与美国三大银行指数收益率的时序图存在明显的差异,表明中国银行业收益率波动并不完全是受美国银行业风险溢出的影响,至于美国银行业对中国银行市场是否存在风险溢出以及溢出的程度多大,还需要进一步的检测。
4.3. 平稳性检验
通常情况下,在做时间序列分析之前需要进行平稳性检验,因为如果时间序列不平稳,继续进行分析就没有了意义。目前,时间序列平稳性检验的方法有多种,包括曲线图法、自相关系数检验法及单位根检验等。本文采用其中应用最广的单位根检验,而单位根检验也有多种方法,我们将采用其中的ADF检验(Augmented Dickey-Fuller test)来对四个银行指数的平稳性进行检验。检验见表2:
Table 2. ADF test results
表2. ADF检验结果
Index |
ADF Statistic |
P-value |
Critical Value (1%) |
DJ |
−10.840531 |
1.620294e−19 |
−3.435414 |
KBW |
−10.844055 |
1.588937e−19 |
−3.435414 |
NASDAQ |
−11.032208 |
5.631248e−20 |
−3.435414 |
ZZYHZS |
−24.312265 |
0.000000e+00 |
−3.435386 |
从上面的检验结果可以看出,DJ、KBW、NASDAQ、ZZYHZS的ADF检验统计量(绝对值)均明显大于各自在1%显著水平下的临界值的绝对值,因此,在1%显著水平下各指数收益率序列均是平稳的。
4.4. ARCH效应检验
ARCH (Autoregressive Conditional Heteroscedasticity)效应就是时间序列数据中存在的条件异方差性。目前用于检验ARCH效应的方法有很多,其中使用最为广泛的是由Engle提出的拉格朗日乘数检验法,即ARCH LM检验。因此,本文采用ARCH LM方法检验四个银行指数的ARCH效应,其结果见表3:
Table 3. ARCH effect test results
表3. ARCH效应检验结果
Index |
LM检验统计量 |
P-value |
DJ |
271.7175 |
0.0000 |
KBW |
248.3357 |
0.0001 |
NASDAQ |
183.5510 |
0.0000 |
ZZYHZS |
11.6132 |
0.0007 |
从检验结果中可以看到,四个银行指数收益率序列所对应的LM检验统计量的概率均小于0.01,说明在1%的显著性水平下,四个指数存在明显的ARCH效应。
4.5. Copula函数选择
Copula是一类描述随机变量之间的依赖关系的函数,其中包含多种不同的Copula函数。本文将中证银行指数分别与道琼斯银行指数、KBW银行指数、纳斯达克银行指数组合形成三个二维数据,然后分别去拟合常见的四种Copula函数,从中选取拟合最优的一个函数去描述该组数据之间的关系。首先,我们利用GARCH(1,1)-t模型对四个收益率需序列进行拟合,从拟合结果中提取出对应的独立同分布的标准残差序列,其次,再用标准残差序列去拟合GaussianCopula、ClaytonCopula、FrankCopula以及GumbelCopula,从而选出拟合最优的Copula函数。以中证银行指数和纳斯达克银行指数为例,在拟合GARCH(1,1)-t模型之后,获得两个指数的标准残差序列
和
,将
和
分别代入标准t分布函数(自由度从GARCH模型中获取),得到边缘分布序列
和
,然后将两个边缘分布序列作为观测值,运用不同的Copula函数进行估计,得到对应的AIC、BIC和HQ值,并根据AIC、BIC和HQ值最小原则选取拟合最优的Copula函数。中证银行指数和纳斯达克银行指数的拟合见表4:
Table 4. Copula fitting results
表4. Copula拟合结果
Copula |
AIC |
BIC |
HQ |
GaussianCopula |
−11.6796 |
−6.5071 |
−9.7391 |
ClaytonCopula |
−18.9882 |
−13.8157 |
−17.0477 |
FrankCopula |
−13.9961 |
−8.8237 |
−12.0556 |
GumbelCopula |
−9.9130 |
−4.7406 |
−7.9725 |
从表中可以看出,对中证银行指数和纳斯达克银行指数,ClaytonCopula函数是拟合最优的Copula函数。同理可以得到对中证银行指数和道琼斯银行指数、中证银行指数和KBW银行指数,拟合最优的Copula函数分别是ClaytonCopula函数和FrankCopula函数。
4.6. VaR和CoVaR计算
在拟合GARCH(1,1)-t模型之后,可以获取四个银行指数对应的标准残差序列的自由度,然在特定
下(5%),分别求出四个银行指数对应的残差独立同分布序列的
分为数值,再根据公式(3.7)求出中证银行指数的无条件风险价值
,将道琼斯银行指数、KBW银行指数、纳斯达克银行指数的
分为数值分别代入方程(3.12)求出中证银行指数的条件
分为数值,再由公式(3.13)求得其对应的条件风险价值
。结果见表5:
Table 5. Risk spillover strength of the three major bank indices in the United States to the CSI bank index
表5. 美国三大银行指数对中证银行指数风险溢出强度
溢出方向 |
VaR |
CoVaR |
ΔCoVaR |
%CoVaR |
DJ→ZZYHZS |
1.53σt |
1.74σt |
0.21σt |
13.73% |
KBW→ZZYHZS |
1.53σt |
1.62σt |
0.09σt |
5.89% |
NASDAQ→ZZYHZS |
1.53σt |
1.69σt |
0.16σt |
10.46% |
表5给出了在5%的置信水平下,三大银行指数代表的美国银行业对中国银行业风险溢出效应的测度结果。首先,从表中可以看出中证银行指数关于道琼斯银行业指数、KBW银行指数、纳斯达克银行指数的条件在线价值分别为1.74σt、1.62σt和1.69σt,由于CoVaR代表的是总风险水平,说明相较而言在道琼斯银行指数发生损失的情况下中国银行业发生损失的风险最大,在KBW银行指数发生损失的情况下中国银行业发生损失的风险最小;其次,从ΔCoVaR的数据可以看出,三个银行业指数对应的增量值都与CoVaR的符号相同,说明各指数对中证银行指数都产生了正向的风险溢出;最后,从%CoVaR数据看,即使是风险溢出程度最小的KBW指数也使得中证银行指数损失增长超过了5%,说明中国银行业与美国银行业关联度较高,受到美国银行业风险的影响较大。由于道琼斯银行指数代表了美国大型商业银行,KBW指数代表的是传统商业银行,而纳斯达克银行指数代表了非传统商业银行,因此据以上数据和分析可以得出,不同规模和不同类型的银行风险溢出程度不同,美国大型商业银行对我国银行业风险溢出效应最强,其次是非传统商业银行。实证结果说明,随着中国金融业深入改革和全面开放,中国银行市场与国际金融市场之间的联系日益加深,中国银行业在享受金融全球化带来的好处的同时,也面临着国际金融市场带来的风险,为保证我国金融市场的稳定,监管当局以及银行业必须对国际风险溢出加以重视。
5. 研究结论与政策建议
5.1. 研究结论
本文以美国三大银行指数和中证银行指数为研究对象,选取四个指数2018年1月1日至2022年12月30日期间的日收益数据,基于GARCH-Copula-CoVaR模型实证研究美国银行业对中国银行业的风向溢出效应,得出以下结论:
美国银行业市场对中国银行业市场存在显著的正向风险溢出效应。具体而言,当美国银行市场面临风险时,这种风险往往会传导至中国银行业市场,加剧其风险程度。进一步分析表明,美国大型商业银行对中国银行业的风险溢出效应更为显著,其影响力超过了中小型商业银行;同时,美国非传统商业银行对我国银行业市场的风险溢出效应相对较强,超越了传统商业银行的影响。这一发现揭示了国际银行业风险传递的复杂性,对于加强我国银行业风险防范与监管具有重要意义。
然而,本研究仍然存在一些不足。首先,本文虽然通过机器学习模型计算得出了美国银行业对中国银行业存在风险溢出的结论以及溢出的程度,但由于机器学习模型本身类似一个黑箱模型,并不能说明中美银行业之间风险传递的路径和机制;其次,本文仅仅研究了美国整个银行业对中国银行业的风险溢出,但是美国银行业发生风险可能并不是整个行业同时出现风险,往往是其中某一个或某几个银行先发生问题,因此就单个银行而言,多大规模的银行发生风险才会引发我国银行业受到冲击,也是有待进一步研究的另一个方面。
5.2. 政策建议
首先,要继续推进银行改革,提高风险管理能力。与西方发达国家的银行业相比,我国现代银行业起步较晚,且改革开放之前长期处于非市场化经营;改革开放之后,面对国际银行业的竞争,国内银行逐步开始改革,但目前任然存在一些问题,特别是风险管理能力较差。国内商业银行还需持续推进体制改革,完善公司治理,学习西方发达国家银行管理经验,以应对国际金融风险的冲击。
其次,要加强银行业监管力度,防止外部风险冲击。银行业对中国经济的重要性不言而喻,银行市场发生动荡,将对整个国民经济社会造成巨大破坏,在以前开放程度不高的情况下,银行市场不太容易发生风险,对风险的预测和管理也更为容易,但在金融业全面开放的今天,银行业面临全球风险事件的风险外溢,其应对风险的难度增加,在银行业增强自身风险管理能力的同时,对银行的监管也越显重要,国家应当完善银行业监管相关的法律法规,并做到与时俱进,以更好地防止银行业受到内外部风险冲击而发生金融风险。
最后,经济全球化、金融自由化仍是当经世界经济发展的大趋势,我们不能因噎废食,而是要从各个层面尽快地完善自身,增强自身竞争力和风险防范能力。从长期看,积极融入国际金融市场利大于弊,更有助于促进我国经济发展。