全球股市风险传染:来自频域多层连通网络的证据
T
T
^("T") { }^{\text {T}}
Zisheng Ouyang
a
a
^("a ") { }^{\text {a }} , Xuewei Zhou
a,*,
a,*,
^("a,*, ") { }^{\text {a,*, }} , Yongzeng Lai
b
b
^("b ") { }^{\text {b }}
a
a
^("a ") { }^{\text {a }} 商学院,湖南师范大学,长沙,410081,中国
b
b
^(b) { }^{\mathrm{b}} 数学系,威尔弗里德·劳里埃大学,滑铁卢,安大略省,N2L 3C5,加拿大
文章信息
JEL 分类:
代码
G14
G15
G32 关键词: 金融风险 频域 多层连接性网络 全球股市
摘要
多层连接网络是揭示金融风险传染机制的有前景的工具。本文构建了频域中的多层连接网络,以考察 2006 年 1 月 4 日至 2022 年 6 月 30 日期间全球股市之间的风险传染。我们分析了单层连接网络的全球效率、平均连接强度和网络密度,并研究了层间连接网络的平均重叠度、网络相关系数和网络参与系数。我们观察到:(i) 全球股市之间的风险传染在短期、中期和长期表现出不同的行为,(ii) 在金融压力期间,中期和长期的连接性显著增加,而短期连接性显著下降,(iii) 亚洲股市是主要的风险接收者,但它们接收的风险在频率上是异质的。我们的工作为研究全球风险传染提供了新的视角,并为投资者和监管者提供了有价值的知识。
1. 引言
随着全球经济一体化的加深,各国金融市场之间出现了密集的连接网络。信息和风险的传播在根本上依赖于这些连接(Broner & Ventura, 2016),例如 2008 年的全球金融危机、2011 年的欧洲债务危机以及 2020 年的 COVID-19 大流行。事实上,自 2008 年全球金融危机以来,各国监管机构加强了对系统性风险的防范和识别。但金融系统中实体(例如金融市场和金融机构)之间的联系通常是复杂和多维的(Wang, Wan et al., 2023),这增加了风险识别的难度。
系统风险的一个常用衡量指标是条件风险价值(CoVaR),它基于“自下而上”的视角,反映了金融机构对市场的风险贡献(Adrian & Brunnermeier, 2016; Borri & Giorgio, 2022; Morelli & Vioto, 2020)。另一个常用的衡量指标是边际预期损失(MES),它关注金融机构在市场危机时对市场的风险贡献。它具有可加性优势(Acharya et al., 2017; Duan et al., 2021; Ouyang et al., 2021)。金融机构的资本短缺程度对系统风险有重要影响。Brownlees 和 Engle(2017)提出了一种基于财务报表数据的 SRISK 风险指数,但金融的频率
声明数据较低。它往往无法及时反映系统性风险的演变。尽管 CoVaR、MES 和 SRISK 等方法被广泛使用,但它们忽视了金融系统中的连通性。因此,它们无法找到金融风险的传播路径,这影响了系统性风险测量的准确性和真实性。
网络分析有潜力帮助监测金融市场的连通性并评估系统脆弱性(Bianchi et al., 2019)。因此,一些文献应用复杂网络理论研究系统性风险。Billio et al.(2012)通过格兰杰因果网络测量了 100 家金融机构的系统性风险和连通性。然而,格兰杰因果网络无法识别每条边的强度。Diebold 和 Yilmaz(2014)基于方差分解构建了溢出指数模型以避免这个问题。这种方法在系统性风险中有许多应用(Ando et al., 2022; Benkraiem et al., 2022; Demirer et al., 2018; Dong et al., 2022; Grillini et al., 2022; Guo & Hou, 2022; Liu & Huang, 2022; Wang et al., 2018)。监管机构在极端事件期间比在正常情况下更关注风险传染。因为极端风险更具破坏性和传染性,它们也严重威胁金融系统的稳定性。Hautsch et al.(2015)使用风险价值(VaR)构建了一个极端风险网络,以检查金融机构之间的系统性风险传染。Hardle et al. (2016) 基于 CoVaR 测量了金融机构之间的尾部连通性。总的来说,关于互联金融系统网络的研究相对丰富。然而,有趣的是,上述研究主要集中在时间域的网络上,忽略了频域中的连通性信息。
事实上,金融市场之间的连通性和风险传染随着频率而变化(Huang et al., 2023; Jian et al., 2023)。此外,从投资组合管理的角度来看,短期投资者(投机者)更关注高频率下股票市场的连通性(即短期连通性),而长期投资者(基金经理和机构投资者)则关注低频率下的关系(即长期连通性)。因此,金融连通性是市场参与者在不同投资期限下的投资行为的结果,这导致了股票市场之间的频率连通性(Bhanja et al., 2023; Luo et al., 2021; Sun et al., 2020)。此外,频率连通性对政策制定者也具有指导意义。因为政策制定者致力于长期市场发展,他们可以根据长期连通性制定政策(Jian et al., 2023; Sun et al., 2020)。 理解金融市场在不同频率下的连通机制对国际投资者的投资决策和政策制定者的风险防范具有重要意义。
为了分析金融数据的频率动态,现有文献开发了傅里叶变换(FT)、小波分析(WA)和经验模态分解(EMD)等方法。傅里叶变换可以将时域信号转换为频域信号(Luo et al., 2021; Olbrys & Mursztyn, 2019)。然而,它仅在时域和频域之间建立一对一关系,分离时间和频率。小波分析可以将时间序列数据分解为不同的频率,广泛应用于金融风险管理(Bales, 2022; Choi, 2022; Qiao et al., 2020; Ren et al., 2022)。小波分析严格依赖于数据的稳定性(Luo et al., 2021),但金融市场数据通常是非平稳的,这限制了小波分析的应用范围。为了克服上述方法的局限性,一些文献应用经验模态分解进行时频分析(Jiang et al., 2022; Jin et al., 2021; Nava et al., 2018; Zhu et al., 2020)。
尽管频率分析已被应用于研究金融系统之间的连通性和风险传染(Liu & Jiang, 2020; Luo et al., 2021; Mensi et al., 2019; Qiao et al., 2020; Sun et al., 2020),但这些研究仅在系统层面上考察金融实体之间的频率机制。因为它们主要关注不同频率下整体金融连通性的变化,忽视了金融实体之间网络拓扑的频率动态。网络分析是一种有前景的工具,可以帮助识别风险传染渠道并发现系统重要节点,这对微观审慎监管具有重要意义。然而,我们发现很少有研究将频率分析和网络分析结合起来,以调查金融连通性和风险传染,关于全球股市之间风险传染的研究更是寥寥无几。
由于上述缺点,我们主要做以下创新工作。首先,基于动态网络(Barunik & Ellington, 2020),我们在频域中构建多层连接网络。该方法考察金融系统在不同频率下的风险传染,如短期、中期或长期。我们将频域中的多层连接网络应用于全球股市,旨在研究全球风险传染的频率机制。其次,我们构建全球效率、平均连接强度和网络密度,以考察单层连接网络的风险传播效率、风险连接强度和节点密度。我们还应用平均重叠度、网络相关系数和网络参与系数来考察多层连接网络之间的结构相似性。第三,我们对全球股市的风险传染进行静态和动态分析,以研究风险传染的频率动态。 最后,我们探讨了全球金融危机和 COVID-19 大流行期间的多层连接网络,以发现在这两场全球危机下网络结构是否以及如何发生变化。
我们的工作与 Ren 等(2022)和 Jian 等(2023)密切相关,他们研究了全球股市之间连通性的频率异质性。但我们的工作与这些最相关的论文有所不同:(1)最相关的论文基于单层网络视角研究全球股市之间的连通性,而没有考虑短期、中期和长期连通网络的层间差异。例如,风险传播效率、风险传染结构和风险传染强度在不同频率网络上是否存在差异?此外,全球股市在不同频率网络上是否发挥相同的作用?这些问题需要通过多层网络来回答。(2)大多数相关论文的动态分析基于滚动窗口,这对异常值非常敏感(Lee et al., 2023)。此外,大多数相关论文没有探讨全球股市之间的连通性在危机前后如何变化。澄清全球股市在危机前后的连通性对国际风险防范具有重要的指导意义。 (3) 最相关的文献发现全球股市之间的连通性在频率上是异质的,但他们没有研究网络结构如何随频率变化。
我们的工作有三个贡献。 (1) 我们通过 Barunik 和 Ellington(2020)提出的动态网络构建频域中的多层连通性网络(例如,短期、中期和长期),旨在研究全球股市在频域中的风险传染网络拓扑。此外,我们通过平均重叠度和网络相关系数来测量短期、中期和长期连通性网络的相似性。 (2) 此外,我们应用 TVP-VAR 模型在频域中构建多层连接网络。与滚动窗口方法相比,TVP-VAR 模型具有不受异常值影响、无需设置滚动窗口以及不损失观测值的优点(Liu & Huang, 2022)。因此,我们获得了一个动态的每日网络,以研究全球股市之间的风险传染,有效捕捉金融连接性和风险传染的动态。 (3) 我们通过所提出的方法得出了有趣的结论。首先,我们发现短期、中期和长期的连通性网络在平静时期高度相关,但在危机期间则相反。其次,值得注意的是,在 2008 年的全球金融危机、2011 年的主权债务危机和 2020 年的 COVID-19 大流行期间,我们观察到全球股市之间的长期连通性急剧增加。第三,我们注意到亚洲国家是主要的风险接收者,但他们所接收的风险往往是异质的。
本文的其余部分组织如下。第二节介绍所使用的方法。第三节将我们的方法应用于全球股市。最后,在第四节中我们总结本文。
2. 方法论
2.1. 测量动态连通性矩阵
现有的连通性网络主要通过溢出指数模型(Demirer et al., 2018; Diebold & Yilmaz, 2014; Wang et al., 2021)构建,但该模型存在一些局限性:(1)溢出指数模型通过滚动窗口获得动态连通性网络,这会丢失部分样本信息,并且对窗口宽度的选择更加敏感(Dai & Zhu, 2022; Huang et al., 2023; Zhou et al., 2023);(2)溢出指数模型基于时间域框架,忽略了连通性网络在频率域中的重要风险信息(Hung & Vo, 2021; Jian et al., 2023; Qiao et al., 2022)。因此,它无法揭示连通性网络的频率演变。值得注意的是,尽管 BK 溢出模型(Barunik & Krehlik, 2018)也可以检查频率异质性,但其动态是通过滚动窗口实现的。鉴于上述缺点,Barunik 和 Ellington(2020)将 TVP-VAR 模型与频率域思想相结合,提出了频率域 TVP-VAR 模型。该模型揭示了瞬态和持久的连通性,并阐明了金融系统中的频率依赖性。 与滚动窗口方法相比,TVP-VAR 模型具有不受异常值影响、无需设置滚动窗口以及不损失观测值的优点(Liu & Huang, 2022)。几篇论文通过频域 TVP-VAR 模型研究了金融连通性中的频率异质性(Gong et al., 2022;Ouyang & Zhou, 2023;Wang, Wei et al., 2023)。
参考 Barunik 和 Ellington (2020),我们考虑一个时间序列
X
t
,
T
=
(
X
t
,
T
1
,
…
,
X
t
,
T
N
)
T
X
t
,
T
=
X
t
,
T
1
,
…
,
X
t
,
T
N
T
X_(t,T)=(X_(t,T)^(1),dots,X_(t,T)^(N))^(T) \boldsymbol{X}_{t, T}=\left(\boldsymbol{X}_{t, T}^{1}, \ldots, \boldsymbol{X}_{t, T}^{N}\right)^{T} ,其中
t
t
t t 指的是离散时间索引,
T
T
T T 是一个额外的索引,表示时间序列
X
t
,
T
X
t
,
T
X_(t,T) \boldsymbol{X}_{t, T} 的局部近似的锐度。我们假设
X
t
,
T
X
t
,
T
X_(t,T) \boldsymbol{X}_{t, T} 遵循滞后阶数为
p
p
p p 的 TVP-VAR。
X
t
,
T
=
Φ
1
(
t
/
T
)
X
t
−
1
,
T
+
⋯
+
Φ
p
(
t
/
T
)
X
t
−
p
,
T
+
ϵ
t
,
T
X
t
,
T
=
Φ
1
(
t
/
T
)
X
t
−
1
,
T
+
⋯
+
Φ
p
(
t
/
T
)
X
t
−
p
,
T
+
ϵ
t
,
T
X_(t,T)=Phi_(1)(t//T)X_(t-1,T)+cdots+Phi_(p)(t//T)X_(t-p,T)+epsilon_(t,T) \boldsymbol{X}_{t, T}=\boldsymbol{\Phi}_{1}(t / T) \boldsymbol{X}_{t-1, T}+\cdots+\boldsymbol{\Phi}_{p}(t / T) \boldsymbol{X}_{t-p, T}+\boldsymbol{\epsilon}_{t, T}
其中
ϵ
t
,
T
=
σ
−
1
/
2
(
t
/
T
)
η
t
,
T
ϵ
t
,
T
=
σ
−
1
/
2
(
t
/
T
)
η
t
,
T
epsilon_(t,T)=sigma^(-1//2)(t//T)eta_(t,T) \boldsymbol{\epsilon}_{t, T}=\sigma^{-1 / 2}(t / T) \boldsymbol{\eta}_{t, T} 与
η
∼
N
I
D
(
0
,
I
M
)
.
I
M
η
∼
N
I
D
0
,
I
M
.
I
M
eta∼NID(0,I_(M)).I_(M) \boldsymbol{\eta} \sim N I D\left(0, \boldsymbol{I}_{M}\right) . \boldsymbol{I}_{M} 是单位矩阵。
ϕ
i
(
t
/
T
)
(
i
=
1
,
2
,
…
,
p
)
ϕ
i
(
t
/
T
)
(
i
=
1
,
2
,
…
,
p
)
phi_(i)(t//T)(i=1,2,dots,p) \boldsymbol{\phi}_{i}(t / T)(i=1,2, \ldots, p) 是时间变化的自回归系数。在固定时间点
u
=
t
/
T
u
=
t
/
T
u=t//T u=t / T 的邻域内,过程
X
t
,
T
X
t
,
T
X_(t,T) \boldsymbol{X}_{t, T} 由平稳过程
X
~
t
(
u
)
X
~
t
(
u
)
tilde(X)_(t)(u) \tilde{\boldsymbol{X}}_{t}(u) 表示。
X
~
t
(
u
)
=
Φ
1
(
u
)
X
~
t
−
1
(
u
)
+
⋯
+
Φ
p
(
u
)
X
~
t
−
p
(
u
)
+
ϵ
t
X
~
t
(
u
)
=
Φ
1
(
u
)
X
~
t
−
1
(
u
)
+
⋯
+
Φ
p
(
u
)
X
~
t
−
p
(
u
)
+
ϵ
t
tilde(X)_(t)(u)=Phi_(1)(u) tilde(X)_(t-1)(u)+cdots+Phi_(p)(u) tilde(X)_(t-p)(u)+epsilon_(t) \tilde{\boldsymbol{X}}_{t}(u)=\boldsymbol{\Phi}_{1}(u) \tilde{\boldsymbol{X}}_{t-1}(u)+\cdots+\boldsymbol{\Phi}_{p}(u) \tilde{\boldsymbol{X}}_{t-p}(u)+\boldsymbol{\epsilon}_{t}
它可以转化为
V
M
A
(
∞
)
V
M
A
(
∞
)
VMA(oo) V M A(\infty) 表示
X
~
t
(
u
)
=
∑
h
=
−
∞
∞
Ψ
(
u
,
h
)
ϵ
t
−
h
X
~
t
(
u
)
=
∑
h
=
−
∞
∞
Ψ
(
u
,
h
)
ϵ
t
−
h
tilde(X)_(t)(u)=sum_(h=-oo)^(oo)Psi(u,h)epsilon_(t-h) \tilde{\boldsymbol{X}}_{t}(u)=\sum_{h=-\infty}^{\infty} \boldsymbol{\Psi}(u, h) \epsilon_{t-h}
根据 Barunik 和 Krehlik(2018)的研究,我们知道
Ψ
(
u
)
e
−
i
ω
=
∑
h
e
−
i
ω
h
Ψ
(
u
,
h
)
Ψ
(
u
)
e
−
i
ω
=
∑
h
e
−
i
ω
h
Ψ
(
u
,
h
)
Psi(u)e^(-i omega)=sum_(h)e^(-i omega h)Psi(u,h) \boldsymbol{\Psi}(u) e^{-i \omega}=\sum_{h} e^{-i \omega h} \boldsymbol{\Psi}(u, h) ,其中
h
h
h h 是预测范围,
i
=
−
1
i
=
−
1
i=sqrt(-1) i=\sqrt{-1} 。在频率
ω
ω
omega \omega 下,
X
~
t
(
u
)
X
~
t
(
u
)
tilde(X)_(t)(u) \tilde{X}_{t}(u) 的谱密度可以定义为
V
M
A
(
∞
)
V
M
A
(
∞
)
VMA(oo) V M A(\infty) 滤波序列的傅里叶变换,如下所示。
S
x
(
u
,
ω
)
=
∑
h
=
−
∞
∞
E
[
X
~
t
+
h
(
u
)
X
~
t
T
(
u
)
]
e
−
i
ω
h
=
{
Ψ
(
u
)
e
−
i
ω
}
Σ
(
u
)
{
Ψ
(
u
)
e
+
i
ω
}
T
S
x
(
u
,
ω
)
=
∑
h
=
−
∞
∞
E
X
~
t
+
h
(
u
)
X
~
t
T
(
u
)
e
−
i
ω
h
=
Ψ
(
u
)
e
−
i
ω
Σ
(
u
)
Ψ
(
u
)
e
+
i
ω
T
S_(x)(u,omega)=sum_(h=-oo)^(oo)E[ tilde(X)_(t+h)(u) tilde(X)_(t)^(T)(u)]e^(-i omega h)={Psi(u)e^(-i omega)}Sigma(u){Psi(u)e^(+i omega)}^(T) \boldsymbol{S}_{x}(u, \omega)=\sum_{h=-\infty}^{\infty} E\left[\tilde{\boldsymbol{X}}_{t+h}(u) \tilde{\boldsymbol{X}}_{t}^{T}(u)\right] e^{-i \omega h}=\left\{\boldsymbol{\Psi}(u) e^{-i \omega}\right\} \Sigma(u)\left\{\boldsymbol{\Psi}(u) e^{+i \omega}\right\}^{T}
接下来,计算第
k
k
k k 个变量对第
j
j
j j 个变量在频带
d
d
d d 上的风险传染,其中
d
=
(
a
,
b
)
:
a
,
b
∈
(
−
π
,
π
)
d
=
(
a
,
b
)
:
a
,
b
∈
(
−
π
,
π
)
d=(a,b):a,b in(-pi,pi) d=(a, b): a, b \in(-\pi, \pi) ,
a
<
b
a
<
b
a < b a<b 。
θ
(
u
,
d
)
j
,
k
=
σ
k
k
−
1
∫
a
b
(
[
Ψ
(
u
)
e
−
i
ω
Σ
(
u
)
]
j
,
k
)
2
d
ω
∫
−
π
π
[
{
Ψ
(
u
)
e
−
i
ω
}
Σ
(
u
)
{
Ψ
(
u
)
e
+
i
ω
}
T
]
j
,
j
d
ω
θ
(
u
,
d
)
j
,
k
=
σ
k
k
−
1
∫
a
b
Ψ
(
u
)
e
−
i
ω
Σ
(
u
)
j
,
k
2
d
ω
∫
−
π
π
Ψ
(
u
)
e
−
i
ω
Σ
(
u
)
Ψ
(
u
)
e
+
i
ω
T
j
,
j
d
ω
theta(u,d)_(j,k)=(sigma_(kk)^(-1)int_(a)^(b)([Psi(u)e^(-i omega)Sigma(u)]_(j,k))^(2)d omega)/(int_(-pi)^(pi)[{Psi(u)e^(-i omega)}Sigma(u){Psi(u)e^(+i omega)}^(T)]_(j,j)d omega) \boldsymbol{\theta}(u, d)_{j, k}=\frac{\sigma_{k k}^{-1} \int_{a}^{b}\left(\left[\boldsymbol{\Psi}(u) e^{-i \omega} \boldsymbol{\Sigma}(u)\right]_{j, k}\right)^{2} d \omega}{\int_{-\pi}^{\pi}\left[\left\{\boldsymbol{\Psi}(u) e^{-i \omega}\right\} \boldsymbol{\Sigma}(u)\left\{\boldsymbol{\Psi}(u) e^{+i \omega}\right\}^{\mathrm{T}}\right]_{j, j} d \omega}
表 1 动态连通性矩阵。
X
t
/
T
1
X
t
/
T
1
X_(t//T)^(1) X_{t / T}^{1}
X
t
/
T
2
X
t
/
T
2
X_(t//T)^(2) X_{t / T}^{2}
⋯
⋯
cdots \cdots
X
t
/
T
N
X
t
/
T
N
X_(t//T)^(N) X_{t / T}^{N}
X
t
/
T
1
X
t
/
T
1
X_(t//T)^(1) X_{t / T}^{1}
θ
~
(
u
,
d
i
)
1
,
1
θ
~
u
,
d
i
1
,
1
tilde(theta)(u,d_(i))_(1,1) \tilde{\theta}\left(u, d_{i}\right)_{1,1}
θ
~
(
u
,
d
i
)
1
,
2
θ
~
u
,
d
i
1
,
2
tilde(theta)(u,d_(i))_(1,2) \tilde{\theta}\left(u, d_{i}\right)_{1,2}
⋯
⋯
cdots \cdots
θ
~
(
u
,
d
i
)
1
,
N
θ
~
u
,
d
i
1
,
N
tilde(theta)(u,d_(i))_(1,N) \tilde{\theta}\left(u, d_{i}\right)_{1, N}
X
t
/
T
2
X
t
/
T
2
X_(t//T)^(2) X_{t / T}^{2}
θ
~
(
u
,
d
i
)
2
,
1
θ
~
u
,
d
i
2
,
1
tilde(theta)(u,d_(i))_(2,1) \tilde{\theta}\left(u, d_{i}\right)_{2,1}
θ
~
(
u
,
d
i
)
2
,
2
θ
~
u
,
d
i
2
,
2
tilde(theta)(u,d_(i))_(2,2) \tilde{\theta}\left(u, d_{i}\right)_{2,2}
⋯
⋯
cdots \cdots
θ
~
(
u
,
d
i
)
2
,
N
θ
~
u
,
d
i
2
,
N
tilde(theta)(u,d_(i))_(2,N) \tilde{\theta}\left(u, d_{i}\right)_{2, N}
⋮
⋮
vdots \vdots
⋮
⋮
vdots \vdots
⋮
⋮
vdots \vdots
⋱
⋱
ddots \ddots
⋮
⋮
vdots \vdots
X
t
/
T
N
X
t
/
T
N
X_(t//T)^(N) X_{t / T}^{N}
θ
~
(
u
,
d
i
)
N
,
1
θ
~
u
,
d
i
N
,
1
tilde(theta)(u,d_(i))_(N,1) \tilde{\theta}\left(u, d_{i}\right)_{N, 1}
θ
~
(
u
,
d
i
)
N
,
2
θ
~
u
,
d
i
N
,
2
tilde(theta)(u,d_(i))_(N,2) \tilde{\theta}\left(u, d_{i}\right)_{N, 2}
⋯
⋯
cdots \cdots
θ
~
(
u
,
d
i
)
N
,
N
θ
~
u
,
d
i
N
,
N
tilde(theta)(u,d_(i))_(N,N) \tilde{\theta}\left(u, d_{i}\right)_{N, N}
X_(t//T)^(1) X_(t//T)^(2) cdots X_(t//T)^(N)
X_(t//T)^(1) tilde(theta)(u,d_(i))_(1,1) tilde(theta)(u,d_(i))_(1,2) cdots tilde(theta)(u,d_(i))_(1,N)
X_(t//T)^(2) tilde(theta)(u,d_(i))_(2,1) tilde(theta)(u,d_(i))_(2,2) cdots tilde(theta)(u,d_(i))_(2,N)
vdots vdots vdots ddots vdots
X_(t//T)^(N) tilde(theta)(u,d_(i))_(N,1) tilde(theta)(u,d_(i))_(N,2) cdots tilde(theta)(u,d_(i))_(N,N) | | $X_{t / T}^{1}$ | $X_{t / T}^{2}$ | $\cdots$ | $X_{t / T}^{N}$ |
| :--- | :--- | :--- | :--- | :--- |
| $X_{t / T}^{1}$ | $\tilde{\theta}\left(u, d_{i}\right)_{1,1}$ | $\tilde{\theta}\left(u, d_{i}\right)_{1,2}$ | $\cdots$ | $\tilde{\theta}\left(u, d_{i}\right)_{1, N}$ |
| $X_{t / T}^{2}$ | $\tilde{\theta}\left(u, d_{i}\right)_{2,1}$ | $\tilde{\theta}\left(u, d_{i}\right)_{2,2}$ | $\cdots$ | $\tilde{\theta}\left(u, d_{i}\right)_{2, N}$ |
| $\vdots$ | $\vdots$ | $\vdots$ | $\ddots$ | $\vdots$ |
| $X_{t / T}^{N}$ | $\tilde{\theta}\left(u, d_{i}\right)_{N, 1}$ | $\tilde{\theta}\left(u, d_{i}\right)_{N, 2}$ | $\cdots$ | $\tilde{\theta}\left(u, d_{i}\right)_{N, N}$ |
需要注意的是,
θ
(
u
,
d
)
j
,
k
θ
(
u
,
d
)
j
,
k
theta(u,d)_(j,k) \theta(u, d)_{j, k} 是传统方差分解到时变频带的自然分解。假设
d
i
(
iin
{
1
,
2
,
…
,
s
}
)
d
i
(
iin
{
1
,
2
,
…
,
s
}
)
d_(i)(iin{1,2,dots,s}) d_{i}(\operatorname{iin}\{1,2, \ldots, s\}) 是一个特定的频带,将不同频带中的风险传染相加,以获得时域中的风险传染。
θ
(
u
)
j
,
k
=
∑
i
=
1
s
θ
(
u
,
d
i
)
j
,
k
θ
(
u
)
j
,
k
=
∑
i
=
1
s
θ
u
,
d
i
j
,
k
theta(u)_(j,k)=sum_(i=1)^(s)theta(u,d_(i))_(j,k) \boldsymbol{\theta}(u)_{j, k}=\sum_{i=1}^{s} \boldsymbol{\theta}\left(u, d_{i}\right)_{j, k}
为了方便解释,我们对风险传染进行了标准化
θ
~
(
u
,
d
i
)
j
,
k
=
θ
(
u
,
d
i
)
j
,
k
∑
k
=
1
N
θ
(
u
)
j
,
k
×
100
θ
~
u
,
d
i
j
,
k
=
θ
u
,
d
i
j
,
k
∑
k
=
1
N
θ
(
u
)
j
,
k
×
100
tilde(theta)(u,d_(i))_(j,k)=(theta(u,d_(i))_(j,k))/(sum_(k=1)^(N)theta(u)_(j,k))xx100 \tilde{\boldsymbol{\theta}}\left(u, d_{i}\right)_{j, k}=\frac{\boldsymbol{\theta}\left(u, d_{i}\right)_{j, k}}{\sum_{k=1}^{N} \boldsymbol{\theta}(u)_{j, k}} \times 100
此外,我们获得了在时间
u
=
t
/
T
u
=
t
/
T
u=t//T u=t / T 和频段
d
i
d
i
d_(i) d_{i} 下
N
N
N N 变量的动态连通性矩阵
C
(
u
,
d
i
)
C
u
,
d
i
C(u,d_(i)) C\left(u, d_{i}\right) (见表 1)。
2.2. 频域中的多层网络
基于动态连通性矩阵
C
(
u
,
d
i
)
C
u
,
d
i
C(u,d_(i)) \boldsymbol{C}\left(u, d_{i}\right) ,我们可以获得频域中的多层连通性网络
Ω
=
(
G
(
u
,
d
1
)
,
G
(
u
,
d
2
)
,
…
,
G
(
u
,
d
i
)
)
Ω
=
G
u
,
d
1
,
G
u
,
d
2
,
…
,
G
u
,
d
i
Omega=(G(u,d_(1)),G(u,d_(2)),dots,G(u,d_(i))) \Omega=\left(G\left(u, d_{1}\right), G\left(u, d_{2}\right), \ldots, G\left(u, d_{i}\right)\right) ,其中
G
(
u
,
d
i
)
=
G
(
V
,
C
(
u
,
d
i
)
)
,
V
=
{
1
,
2
,
…
,
N
}
G
u
,
d
i
=
G
V
,
C
u
,
d
i
,
V
=
{
1
,
2
,
…
,
N
}
G(u,d_(i))=G(V,C(u,d_(i))),V={1,2,dots,N} G\left(u, d_{i}\right)=G\left(V, C\left(u, d_{i}\right)\right), V=\{1,2, \ldots, N\} 是节点集合,
C
(
u
,
d
i
)
C
u
,
d
i
C(u,d_(i)) C\left(u, d_{i}\right) 是边集合或频带
d
i
d
i
d_(i) d_{i} 的动态连通性矩阵。因此,
C
(
u
,
d
i
)
C
u
,
d
i
C(u,d_(i)) \boldsymbol{C}\left(u, d_{i}\right) 中的元素
θ
~
(
u
,
d
i
)
j
,
k
θ
~
u
,
d
i
j
,
k
tilde(theta)(u,d_(i))_(j,k) \tilde{\boldsymbol{\theta}}\left(u, d_{i}\right)_{j, k} 是频带
d
i
d
i
d_(i) d_{i} 的连通性网络中的一个有向加权边,表示从第
k
k
k k 个变量到第
j
j
j j 个变量的相应连通性强度。由于方差分解获得了一个完全连通性网络,这将掩盖一些重要信息。因此,为了去除每一层中的噪声信息,我们使用阈值方法来过滤完全连通性网络,并保留最重要的连通性信息。参考 Ren 等 (2022),我们取矩阵
C
(
u
,
d
i
)
C
u
,
d
i
C(u,d_(i)) C\left(u, d_{i}\right) 的平均值作为阈值,并去除低于阈值的边,以获得过滤后的连通性网络。值得注意的是,多层网络可以大致分为两种类型:多层连通性网络和互联多层网络 (Wang, Wan 等, 2023)。 多层连接网络在每一层具有相同的节点集,而在互联的多层网络中,不同层的节点集可以不同。因此,本文中的多层网络指的是多层连接网络,即每一层的节点集是相同的。最后,本文中频域的多层连接网络被定义为每一层的节点集相同,但每一层的节点在不同频率下连接。
2.3. 网络度量
2.3.1. 全球效率
根据王等人(2021)的研究,我们计算了频域中多层连接网络的全球效率(GE),它反映了金融风险在连接网络中的传播效率。层
d
i
d
i
d_(i) d_{i} 的 GE 定义为
GE
d
i
=
1
N
(
N
−
1
)
∑
j
,
k
=
1
,
j
≠
k
N
1
l
j
,
k
d
i
GE
d
i
=
1
N
(
N
−
1
)
∑
j
,
k
=
1
,
j
≠
k
N
1
l
j
,
k
d
i
GE^(d_(i))=(1)/(N(N-1))sum_(j,k=1,j!=k)^(N)(1)/(l_(j,k)^(d_(i))) \mathrm{GE}^{d_{i}}=\frac{1}{N(N-1)} \sum_{j, k=1, j \neq k}^{N} \frac{1}{l_{j, k}^{d_{i}}}
其中
l
i
j
l
i
j
l_(ij) l_{i j} 是从节点
k
k
k k 到
j
j
j j 在带
d
i
d
i
d_(i) d_{i} 上的最短路径长度。更大的 GE 意味着更有效的风险传播。
2.3.2. 平均连通性强度
我们引入了频域中多层连通网络的平均连通强度(ACS)的计算,它反映了不同连通网络的连通强度。层
d
i
d
i
d_(i) d_{i} 的 ACS 定义为
ACS
d
i
=
1
N
∑
j
,
k
=
1
,
j
≠
k
N
θ
~
(
u
,
d
i
)
j
,
k
ACS
d
i
=
1
N
∑
j
,
k
=
1
,
j
≠
k
N
θ
~
u
,
d
i
j
,
k
ACS^(d_(i))=(1)/(N)sum_(j,k=1,j!=k)^(N) tilde(theta)(u,d_(i))_(j,k) \operatorname{ACS}^{d_{i}}=\frac{1}{N} \sum_{j, k=1, j \neq k}^{N} \tilde{\boldsymbol{\theta}}\left(u, d_{i}\right)_{j, k}
其中
θ
~
(
u
,
d
i
)
j
,
k
θ
~
u
,
d
i
j
,
k
tilde(theta)(u,d_(i))_(j,k) \tilde{\boldsymbol{\theta}}\left(u, d_{i}\right)_{j, k} 表示节点
k
k
k k 对
j
j
j j 在频带
d
i
d
i
d_(i) d_{i} 上的冲击强度。
2.3.3. 网络密度
为了表征网络中节点之间相互连接边的密度,我们引入了网络密度(ND)。层
d
i
d
i
d_(i) d_{i} 的 ND 定义为
ND
d
i
=
2
L
d
i
N
(
N
−
1
)
ND
d
i
=
2
L
d
i
N
(
N
−
1
)
ND^(d_(i))=(2L^(d_(i)))/(N(N-1)) \mathrm{ND}^{d_{i}}=\frac{2 L^{d_{i}}}{N(N-1)}
其中
L
d
i
L
d
i
L^(d_(i)) L^{d_{i}} 是层
d
i
d
i
d_(i) d_{i} 的实际连接边数。更高的 ND 表明网络更复杂。
2.3.4. 平均重叠度
为了阐明频域中多层连通网络的边结构相似性,我们构建了平均重叠度(AOD)指标。层
d
i
d
i
d_(i) d_{i} 的 AOD 定义为
AOD
=
1
K
∑
j
,
k
=
1
,
j
≠
k
N
∑
i
=
1
s
a
j
,
k
d
i
AOD
=
1
K
∑
j
,
k
=
1
,
j
≠
k
N
∑
i
=
1
s
a
j
,
k
d
i
AOD=(1)/(K)sum_(j,k=1,j!=k)^(N)sum_(i=1)^(s)a_(j,k)^(d_(i)) \mathrm{AOD}=\frac{1}{K} \sum_{j, k=1, j \neq k}^{N} \sum_{i=1}^{s} a_{j, k}^{d_{i}}
其中
a
j
,
k
d
i
=
sign
(
θ
~
(
u
,
d
i
)
j
,
k
)
,
K
=
∑
j
,
k
=
1
,
j
≠
k
N
[
1
−
∏
i
=
1
s
(
1
−
a
j
,
k
d
i
)
]
a
j
,
k
d
i
=
sign
θ
~
u
,
d
i
j
,
k
,
K
=
∑
j
,
k
=
1
,
j
≠
k
N
1
−
∏
i
=
1
s
1
−
a
j
,
k
d
i
a_(j,k)^(d_(i))=sign(( tilde(theta))(u,d_(i))_(j,k)),K=sum_(j,k=1,j!=k)^(N)[1-prod_(i=1)^(s)(1-a_(j,k)^(d_(i)))] a_{j, k}^{d_{i}}=\operatorname{sign}\left(\tilde{\boldsymbol{\theta}}\left(u, d_{i}\right)_{j, k}\right), K=\sum_{j, k=1, j \neq k}^{N}\left[1-\prod_{i=1}^{s}\left(1-a_{j, k}^{d_{i}}\right)\right] 。如果
θ
~
(
u
,
d
i
)
j
,
k
≠
0
,
a
j
,
k
d
i
=
1
θ
~
u
,
d
i
j
,
k
≠
0
,
a
j
,
k
d
i
=
1
tilde(theta)(u,d_(i))_(j,k)!=0,a_(j,k)^(d_(i))=1 \tilde{\boldsymbol{\theta}}\left(u, d_{i}\right)_{j, k} \neq 0, a_{j, k}^{d_{i}}=1 。当
s
s
s s 层的边结构相同时,AOD 的值为
s
s
s s ,而如果每个边仅存在于一个层中,则其值为 1。
2.3.5. 网络相关系数
AOD 仅测量频域中多层连接网络的边缘结构相似性,无法检查边缘权重(连接强度)的相似性。因此,参考皮尔逊相关系数的做法,我们计算网络相关系数(NCC)。层
d
i
d
i
d_(i) d_{i} 的 NCC 定义为
NCC
=
∑
t
=
1
T
(
x
t
d
i
−
x
¯
d
i
)
(
x
t
d
j
−
x
¯
d
j
)
∑
t
=
1
T
(
x
t
d
i
−
x
¯
d
i
)
2
∑
t
=
1
T
(
x
t
d
j
−
x
¯
d
j
)
2
NCC
=
∑
t
=
1
T
x
t
d
i
−
x
¯
d
i
x
t
d
j
−
x
¯
d
j
∑
t
=
1
T
x
t
d
i
−
x
¯
d
i
2
∑
t
=
1
T
x
t
d
j
−
x
¯
d
j
2
NCC=(sum_(t=1)^(T)(x_(t)^(d_(i))- bar(x)^(d_(i)))(x_(t)^(d_(j))- bar(x)^(d_(j))))/(sqrt(sum_(t=1)^(T)(x_(t)^(d_(i))- bar(x)^(d_(i)))^(2))sqrt(sum_(t=1)^(T)(x_(t)^(d_(j))- bar(x)^(d_(j)))^(2))) \mathrm{NCC}=\frac{\sum_{t=1}^{T}\left(x_{t}^{d_{i}}-\bar{x}^{d_{i}}\right)\left(x_{t}^{d_{j}}-\bar{x}^{d_{j}}\right)}{\sqrt{\sum_{t=1}^{T}\left(x_{t}^{d_{i}}-\bar{x}^{d_{i}}\right)^{2}} \sqrt{\sum_{t=1}^{T}\left(x_{t}^{d_{j}}-\bar{x}^{d_{j}}\right)^{2}}}
其中
x
t
d
i
x
t
d
i
x_(t)^(d_(i)) x_{t}^{d_{i}} 是通过将过滤后的网络转换为列向量获得的,
x
¯
d
i
x
¯
d
i
bar(x)^(d_(i)) \bar{x}^{d_{i}} 是
x
t
d
i
x
t
d
i
x_(t)^(d_(i)) x_{t}^{d_{i}} 的均值。
2.3.6. 网络参与系数
单个节点在多层网络中可能扮演不同的角色,例如,节点
k
k
k k 在层
d
i
d
i
d_(i) d_{i} 中具有更高的风险,但在层
d
j
d
j
d_(j) d_{j} 中不一定具有相同的风险水平。因此,为了衡量节点风险水平的离散程度,定义了以下网络参与系数 (NPC):
NPC
k
out
=
s
s
−
1
[
1
−
∑
i
=
1
s
(
out
k
d
i
sout
k
)
2
]
NPC
k
in
=
s
s
−
1
[
1
−
∑
i
=
1
s
(
i
n
k
d
i
sin
k
)
2
]
NPC
k
out
=
s
s
−
1
1
−
∑
i
=
1
s
out
k
d
i
sout
k
2
NPC
k
in
=
s
s
−
1
1
−
∑
i
=
1
s
i
n
k
d
i
sin
k
2
{:[NPC_(k)^("out ")=(s)/(s-1)[1-sum_(i=1)^(s)((" out "_(k)^(d_(i)))/(" sout "_(k)))^(2)]],[NPC_(k)^("in ")=(s)/(s-1)[1-sum_(i=1)^(s)((in_(k)^(d_(i)))/(sin_(k)))^(2)]]:} \begin{aligned}
& \mathrm{NPC}_{k}^{\text {out }}=\frac{s}{s-1}\left[1-\sum_{i=1}^{s}\left(\frac{\text { out }_{k}^{d_{i}}}{\text { sout }_{k}}\right)^{2}\right] \\
& \mathrm{NPC}_{k}^{\text {in }}=\frac{s}{s-1}\left[1-\sum_{i=1}^{s}\left(\frac{i n_{k}^{d_{i}}}{\sin _{k}}\right)^{2}\right]
\end{aligned}