TeX 62.4%
Python 37.5%
1% ============================================================================2\section{Methodology}3\label{sec:methodology}4% ============================================================================56\subsection{Realized Volatility}78We construct daily realized volatility from intraday 5-minute log returns. For asset $i$ on day $t$, let $P_{i,t,j}$ denote the closing price of the $j$-th 5-minute bar, where $j = 1, \ldots, M_t$ and $M_t$ is the number of intraday observations. The 5-minute log return is:9\begin{equation}10 r_{i,t,j} = \ln\left(\frac{P_{i,t,j}}{P_{i,t,j-1}}\right).11 \label{eq:logreturn}12\end{equation}1314The daily realized variance is computed as the sum of squared intraday returns \citep{Andersen2003,BarndorffNielsen2002}:15\begin{equation}16 \RV_{i,t} = \sum_{j=1}^{M_t} r_{i,t,j}^2.17 \label{eq:rv}18\end{equation}1920The annualized daily realized volatility is then:21\begin{equation}22 \sigma_{i,t}^{\RV} = \sqrt{252 \cdot \RV_{i,t}}.23 \label{eq:rv_annualized}24\end{equation}2526We impose a minimum threshold of $M_t \geq 20$ intraday observations per day to ensure reliable realized volatility estimation. Days with fewer observations are excluded. Under standard assumptions, as the sampling frequency increases ($M_t \to \infty$), $\RV_{i,t}$ converges to the integrated variance:27\begin{equation}28 \RV_{i,t} \xrightarrow{p} \int_0^1 \sigma_{i,t}^2(s)\, ds.29 \label{eq:rv_convergence}30\end{equation}3132The 5-minute sampling frequency represents a widely accepted compromise between capturing sufficient intraday variation and avoiding microstructure noise contamination \citep{Liu2015}.3334\subsection{Vector Autoregression}3536Let $\mathbf{y}_t = (\RV_{1,t}, \RV_{2,t}, \ldots, \RV_{N,t})'$ denote the $N \times 1$ vector of daily realized volatilities for $N = 10$ assets. We model the joint dynamics using a VAR($p$) process:37\begin{equation}38 \mathbf{y}_t = \boldsymbol{\mu} + \sum_{\ell=1}^{p} \boldsymbol{\Phi}_\ell \mathbf{y}_{t-\ell} + \boldsymbol{\varepsilon}_t,39 \label{eq:var}40\end{equation}41where $\boldsymbol{\mu}$ is an $N \times 1$ vector of intercepts, $\boldsymbol{\Phi}_\ell$ are $N \times N$ coefficient matrices, and $\boldsymbol{\varepsilon}_t \sim (0, \boldsymbol{\Sigma})$ is the vector of reduced-form innovations with covariance matrix $\boldsymbol{\Sigma}$. We set $p = 4$ lags throughout. The VAR($p$) admits the moving average representation:42\begin{equation}43 \mathbf{y}_t = \boldsymbol{\mu}^* + \sum_{s=0}^{\infty} \boldsymbol{\Psi}_s \boldsymbol{\varepsilon}_{t-s},44 \label{eq:vma}45\end{equation}46where the $N \times N$ coefficient matrices $\boldsymbol{\Psi}_s$ are computed recursively:47\begin{equation}48 \boldsymbol{\Psi}_s = \sum_{\ell=1}^{p} \boldsymbol{\Phi}_\ell \boldsymbol{\Psi}_{s-\ell}, \quad s = 1, 2, \ldots49 \label{eq:psi_recursion}50\end{equation}51with $\boldsymbol{\Psi}_0 = \mathbf{I}_N$ and $\boldsymbol{\Phi}_\ell = \mathbf{0}$ for $\ell > p$.5253\subsection{Generalized Forecast Error Variance Decomposition}5455Following \citet{DieboldYilmaz2012}, we employ the generalized variance decomposition of \citet{Koop1996} and \citet{Pesaran1998}, which is invariant to the ordering of variables. The $H$-step-ahead generalized forecast error variance decomposition is defined as:56\begin{equation}57 \theta_{ij}^g(H) = \frac{\sigma_{jj}^{-1} \sum_{h=0}^{H-1} (\mathbf{e}_i' \boldsymbol{\Psi}_h \boldsymbol{\Sigma} \mathbf{e}_j)^2}{\sum_{h=0}^{H-1} \mathbf{e}_i' \boldsymbol{\Psi}_h \boldsymbol{\Sigma} \boldsymbol{\Psi}_h' \mathbf{e}_i},58 \label{eq:gfevd}59\end{equation}60where $\sigma_{jj}$ is the $j$-th diagonal element of $\boldsymbol{\Sigma}$, $\mathbf{e}_i$ is the $i$-th column of $\mathbf{I}_N$, and $H = 10$ days. Since the rows of $\theta_{ij}^g(H)$ do not necessarily sum to unity under the generalized framework, we normalize:61\begin{equation}62 \tilde{\theta}_{ij}^g(H) = \frac{\theta_{ij}^g(H)}{\sum_{j=1}^{N} \theta_{ij}^g(H)},63 \label{eq:gfevd_normalized}64\end{equation}65such that $\sum_{j=1}^{N} \tilde{\theta}_{ij}^g(H) = 1$ and $\sum_{i,j=1}^{N} \tilde{\theta}_{ij}^g(H) = N$.6667\subsection{Spillover Measures}6869From the normalized GFEVD matrix $\tilde{\boldsymbol{\Theta}}^g(H) = [\tilde{\theta}_{ij}^g(H)]$, we define four spillover measures.7071\paragraph{Total Spillover Index.}72\begin{equation}73 \TSI(H) = \frac{\sum_{\substack{i,j=1 \\ i \neq j}}^{N} \tilde{\theta}_{ij}^g(H)}{N} \times 100.74 \label{eq:tsi}75\end{equation}7677\paragraph{Directional Spillovers TO Others.} The volatility spillover transmitted by market $j$ to all other markets:78\begin{equation}79 S_{\cdot j}^g(H) = \sum_{\substack{i=1 \\ i \neq j}}^{N} \tilde{\theta}_{ij}^g(H) \times 100.80 \label{eq:to_spillover}81\end{equation}8283\paragraph{Directional Spillovers FROM Others.} The volatility spillover received by market $i$ from all other markets:84\begin{equation}85 S_{i \cdot}^g(H) = \sum_{\substack{j=1 \\ j \neq i}}^{N} \tilde{\theta}_{ij}^g(H) \times 100.86 \label{eq:from_spillover}87\end{equation}8889\paragraph{Net Directional Spillovers.} The net volatility spillover for market $i$:90\begin{equation}91 S_i^{\text{net}}(H) = S_{\cdot i}^g(H) - S_{i \cdot}^g(H).92 \label{eq:net_spillover}93\end{equation}94A positive (negative) value indicates that market $i$ is a net transmitter (receiver) of volatility.9596\subsection{Rolling Window Analysis}9798To capture time-varying dynamics, we estimate the VAR model and compute the GFEVD using a rolling window of $w = 200$ trading days, advancing by one day:99\begin{equation}100 \TSI_t(H) = \TSI(H; \mathbf{y}_{t-w+1}, \ldots, \mathbf{y}_t).101 \label{eq:rolling_tsi}102\end{equation}103104\subsection{Granger Causality Tests}105106For realized volatility series $\RV_{i,t}$ and $\RV_{j,t}$, we test whether lagged values of $\RV_{j,t}$ have predictive power for $\RV_{i,t}$ beyond own lags \citep{Granger1969}. The null hypothesis is tested using an $F$-statistic:107\begin{equation}108 F = \frac{(\text{RSS}_R - \text{RSS}_U) / p}{\text{RSS}_U / (T - 2p - 1)},109 \label{eq:granger_f}110\end{equation}111where $\text{RSS}_R$ and $\text{RSS}_U$ are the residual sums of squares from the restricted and unrestricted models.112113\subsection{Network Representation}114115We represent the volatility connectedness structure as a weighted directed graph $\mathcal{G} = (\mathcal{V}, \mathcal{E}, \mathcal{W})$, where the vertex set $\mathcal{V}$ corresponds to the $N$ assets, the edge set $\mathcal{E}$ contains directed edges $(j, i)$ where the GFEVD contribution $\tilde{\theta}_{ij}^g(H)$ exceeds the median of all off-diagonal elements, and edge weight equals the magnitude of the pairwise spillover. Node size is proportional to $S_{\cdot j}^g(H)$, and nodes are colored by asset class.116