Systemic Risk Modeling with Lévy Copulas

We investigate a systemic risk measure known as CoVaR that represents the value-at-risk (VaR) of a financial system conditional on an institution being under distress. For characterizing and estimating CoVaR, we use the copula approach and introduce the normal tempered stable (NTS) copula based on the Lévy process. We also propose a novel backtesting method for CoVaR by a joint distribution correction. We test the proposed NTS model on the daily S&P 500 index and Dow Jones index with in-sample and out-of-sample tests. The results show that the NTS copula outperforms traditional copulas in the accuracy of both tail dependence and marginal processes modeling.


Introduction
Value-at-risk (VaR) quantifies the level of financial risk of investments over a specific time frame, and it is the most common tool in risk management. However, VaR deals with individual institutions and therefore fails to capture risks, for example, when two financial institutions are dynamically coupled. Therefore, a systemic risk measure is required that can quantify the contribution of systemic risk of one financial institution to other institutions. To that end, the most three popular systemic risk measures are CoVaR (Adrian and Brunnermeier 2011), CoCVaR (Huang and Uryasev 2018), and the (lower) tail dependence coefficient (TDC) (Heffernan 2000). They all rely on the conditional distribution of the variable X j of one financial institution given the variable X i of another institution that is under stress. Computing CoVaR and the associated CoCVaR, CoETL, are the main measures for systemic risk as required by the FED (Tobias and Brunnermeier 2016). It is also a probabilistic version of stress testing (Basel Committee on Banking Supervision 2009). A major drawback of the stress testing when evaluating the systemic risk is that it does not determine the probabilistic structure of the stress tested scenarios. CoVaR, CoCVaR, and CoETL, provide probabilistic characterization of the systemic risk financial institutions face. One can argue that CoVaR has advantages over CoCVAR and TDC. The implementation and backtesting of CoCVaR has been criticized to be time-consuming and less straightforward (Emmer et al. 2015). Moreover, it is unnecessary to pick a more complex risk measurement than VaR because VaR (1) is sub-additive for the tails that are most commonly observed in financial applications (Danielsson et al. 2005) and (2) can be used to derive regulatory capital while operational loss data typically exhibit irregularities and complicate the mathematical modeling (Embrechts et al. 2003). Speaking of the TDC, although the TDC is a natural metric to measure tail dependence (Pourkhanali et al. 2016), CoVaR is a "prelimit" of the TDC and keeps the information of the tail dependence, even if the TDC shows 0 when the tail is thin. Furthermore, CoVaR is a measure of tail dependence regardless of the distributional tails, while the TDC might be a good measure only when the distributional tails are exactly of Pareto type (Chollete et al. 2009;Cohen et al. 2020). In summary, CoVaR is the best measure for tail dependence among the three candidates, and is mainly used in this paper to portray the systemic risk. In most of the literature, the "standard" event for CoVaR is X i = VaR α (X i ), where VaR α (X i ) is the value-at-risk of X i with a quantile α. Instead, we choose a generalized event of interest defined by X i ≤ VaR α (X i ) because the generalized definition gives more consistency with respect to concordance order between copulas (Bernardi et al. 2017;Mainik and Schaanning 2014).
Copulas represent a powerful approach to calculate CoVaR and C=°CVaR because they can efficiently describe dependence between random variables. Copulas find a wide application in quantitative finance in modeling and in minimizing tail risks and portfolio-optimizations (Braiek et al. 2020;Low et al. 2013Low et al. , 2016. A copula measures only the dependence part, leading to flexibilities of the marginal part. In the current literature, the marginal processes are usually modeled by extreme-value theory (EVT) (Di Bernardino et al. 2013;Nguyen and Samorodnitsky 2013), Markovian structures (Härdle et al. 2012;Rémillard et al. 2012), sequential Monte Carlo (Targino et al. 2015) or multivariate GARCH with normal or (skewed) Student's t innovations (Chollete et al. 2009;Garcia and Tsafack 2008;Girardi and Ergün 2013;Karimalis and Nomikos 2018;Lien et al. 2013;Reboredo and Ugolini 2015;Yang et al. 2018;Zhang and Maringer 2010). For the dependence part, a variety of copulas are developed to meet different properties (e.g., heavy tails, asymmetries, etc.). The most popular families of copulas are Archimedean copulas (Garcia and Tsafack 2008;Karimalis and Nomikos 2018;Liebscher 2008;Lien et al. 2013;Yang et al. 2018), extreme-value copulas (Jaworski 2017), elliptical copulas (Jaworski et al. 2010), and vine copulas (Braiek et al. 2020;Chollete et al. 2009;Pourkhanali et al. 2016;Zhang and Maringer 2010). A problem with these copulas is that they do not lead to semi-martingale pricing processes and thus, are useless in asset pricing (Delbaen and Schachermayer 1999;Tankov 2016). For instance, the Student's t innovation does not have an exponential moment for the price unit increment, and thus cannot be used for no-arbitrage option pricing.
Instead, we describe a systemic risk with Lévy copula family (Sato 1999) that allows for asset pricing and in particular portfolio insurance. Normal inverse Gaussian (NIG) is introduced by (Barndorff-Nielsen 1997), and NTS is introduced by (Barndorff-Nielsen and Levendorskii 2001; Barndorff-Nielsen and Shephard 2001). They both belong to Lévy family and are applied in finance (Eberlein and Madan 2009;Eberlein and Özkan 2005). To confirm that the NTS fits the marginal process sufficiently, we apply NIG as a benchmark marginal process. To verify that the NTS copula explains the dependence structure of the underlying assets or indices adequately, we compare the performance of the NTS copula with those of other copulas. Specifically, in this paper, we choose the Archimedean copulas as a benchmark copula family, for it captures a wide range of asymmetric tail dependences and allows for analytical computations of CoVaR (De Luca and Rivieccio 2012) and CoCVaR on the account of their simple generator functions. The same methodology used in our paper can readily be applied when the innovations from ARMA-GARCH are modeled with multivariate GH distribution, and in particular with an NIG distribution. The results of using (1) Multivariate GH, or (2) Multivariate NTS, or (3) other Multivariate Levy (infinitely divisible) distribution with semi-heavy distributional tails provide similar results when used in financial risk management (Schoutens 2003).
Furthermore, we also address an important question related to backtesting. Specifically, existing CoVaR backtesting methods consider the testing of CoVaR (Girardi and Ergün 2013) without testing the VaR. They simply assume that the predicted VaR is equal to the actual VaR. We propose a modified Kupiec's test that accounts for simultaneous testing of CoVaR and VaR.

Review of Literature
The recent financial crisis has stimulated financial institutions to do extensive research on systemic risk. From the definition of systemic risk given by the Group of Ten (2001), systemic financial risk is when an event triggers a loss of economic value in a substantial portion of the financial system, or even spread to the real economy. The European Central Bank defines systemic risk as a risk that "impairs the functioning of a financial system to the point where economic growth and welfare suffer materially" (Smaga 2014).
The recent literature has followed two main ways of measuring and analyzing institutional contributions to systemic risk. The first approach, called supervisory approach, relies on macroeconomic or balance sheet data to explain how the risk transmission occurred (Forbes and Rigobon 2001;Pritsker 2001) or studies specific financial institutions information not captured by the markets, such as the systemic risk measure (SRISK) (Brownlees and Engle 2017). In (Cerutti et al. 2012), a lack of useful and consistent data were pointed out and common templates for global and systematic reports for financial institutions were suggested. The second approach for measuring the systemic risk, called market-based approach, is based on high frequency time series data from publicly available markets (Segoviano Basurto and Goodhart 2009) (Giglio et al. 2016), and size_con (Giglio et al. 2016). The accuracy tests suggested that there are very marginal differences across the competing systemic risk factors. Furthermore, the first four factors of the list can be seen as 'top-bottom measures' in the sense that they aim to determine the impact of distress occurring at the level of the financial system on an individual financial institution. By contrast, the CoVaR and ∆CoVaR measures of systemic risk might be seen as 'bottom-up' measures in that they assess the impact of distress at the level of a single financial institution and the transmission of the associated risks to the entire financial system. To identify which methodology in estimating CoVaR could become a proper tool, the empirical analysis from (Bianchi and Sorrentino 2020) suggests that the GARCH model has better persistence and stability than the quantile regression method and the non-parametric (or historical) method, since it captures the time-variation of the volatility. In this paper, we adopt the CoVaR as systemic risk expressed by copulas with GARCH margins.
Methods from the existing literature on dependency of asset returns cannot be applied to option pricing. To price options by no-arbitrage arguments, a price process must be a semi-martingale (Shiryaev 1999), for example an NTS process, which can be characterized as a solution to a backward PDE. VaR and CVaR under the ARMA-GARCH with Tempered Stable Innovation have been studied in (Kim et al. 2011), and the multidimensional ARMA-GARCH-NTS was investigated in (Anand et al. 2016). Based on the previous research, we propose a CoVaR model with an NTS copula and innovations to produce semi-martingale price processes. Moreover, we propose a copula approach to compute CoCVaR.

The Methodology
In this paper, we apply a copula approach to obtain both CoVaR and CoCVaR. We note that in the existing literature, CoCVaR is only investigated based on quantile regression (Huang and Uryasev 2018). As pointed out, the copula-based method allows flexible modeling of marginal processes and the dependence structure. The steps for computing CoVaR and CoCVaR with a copula includes the following steps: 1.
estimation of the marginal process by an ARMA-GARCH-NTS model, 2.
estimation of the various copula functions by marginal innovations, and 3.
computation of the CoVaR and CoCVaR from the estimated copulas.

Copula
A copula links marginal distributions together to form a joint distribution. Often, the copulas are assumed to belong to the Archimedean family because the family allows for analytical expressions. In this paper, we apply the NTS copula since it captures dependence of extreme events (Kim and Volkmann 2013) and can be applied for option pricing.

The Archimedean Copula
A d-dimensional copula C is called Archimedean if for some generator function φ, it can be expressed as is a continuous and strictly decreasing function such that The copula density c is given by The maximum likelihood estimation method exploits Equation (2) for the estimation of the copula parameter(s). In our work, we use three different Archimedean copulas whose generator and inverse generator functions are shown in Table 1. Table 1. Archimedean copula generator functions and inverse functions.

The NTS Copula
A Lévy copula is coupled with a marginal Lévy process. The Lévy copula-based models allow for describing the sizes of simultaneous jumps of a Lévy process. An important member of the Lévy copula family is the NTS copula. The NTS distribution with the NTS copula is defined by the Gaussian subordination method (Kim and Volkmann 2013), and the one-dimensional tempered stable subordinator is T ∼ TS(α, θ) with a characteristic function where µ µ µ, β β β, γ γ γ ∈ R N , is an N-dimensional standard normal random vector distributed with a covariance matrix R = [ρ k,l ] k,l∈{1,2,...,n} , and a subordinator T ∼ TS(α, θ) that is independent of . We denote x ∼ NTS N (α, θ, β β β, γ γ γ, µ µ µ, R). If the elements of x are zero mean and with unit variance, i.e., Var(x) = (1, . . . , 1) , and γ n = 1 − β 2 n 2−α 2θ , where |β n | < 2θ 2−α , we call x a standard NTS random variable and denote it by x ∼ stdNTS N (α, θ, R). The NTS copula is constructed from the NTS distribution by where F stdNTS N is the cumulative distribution function (CDF) corresponding to the standard NTS distribution, stdNTS N , and F n is the marginal CDF of stdNTS N (α, θ, β n , 1). The F stdNTS N and F n are expressed and defined in (Kim and Volkmann 2013) in detail.
There is another family of Lévy copulas that is commonly generated from Archimedean copulas. These copulas are parametric and have generator functions. In this paper, we consider a Lévy Archimedean copula called Lévy-Clayton (L-C) copula whose generator function is defined in Table 1.

CoVaR
For a given return X i of institution i and a confidence level Compared with the traditional use of VaR, we chose to define it with a quantile because this definition allows for simplification in notation and emphasizes the risk of loss. For the return X j of institution j and a confidence level 1 − β, the CoVaR of X j |X i , formally denoted as The definition of CoVaR in fact comes from the definition of a tail dependence function (Garcia and Tsafack 2008). Let F i , F j be the marginal CDFs of X i and X j , respectively. The lower tail dependence function of X i and X j is defined as (Frahm 2004) The lower tail dependence coefficient (TDC) is simply the limit of this function when t tends to zero. More precisely, It is clear from these definitions that the CoVaR retains all the information for TDC and is better for quantifying systemic risk. Considering that the financial system is composed of n institutions and X 1 |X 2 includes indirect effects from any other institutions of the system, we can define a generalized CoVaR as (Rivera-Castro et al. 2018) 3.2.1. ∆CoVaR CoVaR itself does not give as much information as wanted about the systemic risk of a particular institution under stress. A notion that captures this better is a measure known as ∆CoVaR, which is originally defined as follows: This measure represents the difference between CoVaR j|i α,β under stress and CoVaR j|i α=0.5,β under the normal state (α = 0.5). An even better measure, which we adopt here, is the percentage change of CoVaR between stressed and normal states. It is formally defined by (Zhang 2015)

Calculation of CoVaR by a Copula
Please note that Equation (6) is equivalent to Given VaR α (X i ), i.e., P(X i ≤ VaR α (X i )) = α, we aim to find CoVaR j|i α,β that satisfies The last equation can be written in the format of a joint CDF as where F j,i (·, ·) is the joint CDF of (X j , X i ). According to Sklar's theorem (Sklar 1959), the joint CDF of two continuous variables can be represented by a copula function. Specifically, Equation (12) can be expressed by where C(·, ·) is a copula function, and F j (·) are the marginal CDFs of X i and X j , respectively. Based on Equation (14), CoVaR can be computed in four steps: 1.
estimate the marginal CDFs F i and F j from the available data, 2.
solve u from Equation (14), and 4. obtain CoVaR by inverting the marginal CDF: CoVaR j|i α,β = F −1 j (u). If the copula C(·, ·) is Archimedean copula, Equation (14) can be represented as After estimating the parameters of the generator φ, we can solve for u analytically by

CoCVaR
Recall that CVaR is the averaged VaR, defined by where f j (x) is the density function of X j . CVaR can also be backtested analogously to VaR but asymptotically (Du and Escanciano 2017 If the copula belongs to the Archimedean family, CoCVaR has an analytical expression, but still it cannot be applied for option pricing. With the NTS copula, CoCVaR is computed by the Monte Carlo approach and is time-consuming. Furthermore, for a continuous distribution CVaR 99% coincides with the expected tail loss (ETL) 99 %. On the other hand, ETL 99% is the average loss above VaR 99% whereas VaR 99.5% is the median loss above VaR 99%. As the median can be viewed as a robust version of the mean, we view VaR 99.5% as a more robust version of CVaR 99%. As discussed in the introduction section, CoVaR is sufficient and does not have any loss of information.
∆CoCVaR As with ∆CoVaR, a notion that captures systemic risk better in comparison to CoCVaR is ∆CoCVaR, which has the following definition:

The ARMA-GARCH-NTS (AGNTS) Model
For the marginal models, we consider that the market log returns (X j (t), t ≥ 0) for institution j follow an ARMA-GARCH process. Using this model, first the marginal residuals are extracted and then they are analyzed further (Anand et al. 2016). Specifically, we adopt the common assumption of ARMA(1,1)-GARCH(1,1) stochastic process: where µ j (t + 1) is the conditional mean, σ j (t + 1) is the conditional standard deviation, and j (t + 1) is the "innovation" with zero mean and unit variance. The conditional mean model is constructed by ARMA(1,1), where a, b, c ∈ R, and |a| ≤ 1. The conditional variance model is given by GARCH(1,1), where d j > 0, and k j , l j ≥ 0. We assume that the innovation ( 1 (t), 2 (t), · · · , N (t)) follows the multivariate standard NTS distribution, since the NTS distribution is more general and asymmetric compared to the Gaussian and Student's t distributions.
The NTS series are easy to simulate but this approach does not provide a closedform density function. To estimate the parameters of the ARMA-GARCH-NTS process more efficiently, an indirect inference method has to be introduced. CDF-based matching approaches, one of the indirect inference methods, have the drawback of requiring arbitrarily the frequencies of interest (Garcia et al. 2011). To avoid this problem, we adopt the quasi-likelihood-based method (Garcia et al. 2011) with Student's t innovations. We first fit the ARMA-GARCH-T process as an auxiliary model, and then fit the extracted innovations with the standard NTS. This estimation of standard NTS parameters is examined in Section 7 with the following procedure: 1.
Fit the ARMA-GARCH-T process with Student's t innovations and generate a large number of scenarios from the fitted Student's t distribution.

2.
Fit the generated scenarios with an NTS distribution and generate the same size of scenarios from the fitted NTS distribution.

3.
Use the Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) two-sample tests to confirm that these two sets of scenarios are indistinguishable.

Calculation of CoVaR and CoCVaR by AGNTS
Since the distribution of the asset returns is determined by the corresponding innovations, we first work with the CoVaR of the innovations rather than the asset returns. Considering the VaR of returns can be derived by the VaR of innovations by For β we can write, We then have Thus, once we determine the predicted CoVaR of innovations, we plug it in (20) to predict the CoVaR of log returns. In terms of linearity of integration, we obtain CoCVaR of returns also from plugging in the CoCVaR of innovations.

The Generalized Kupiec's Test
For testing P(X j ≤ CoVaR j|i α,β |X i ≤ VaR α (X i )) = β, we first consider the method from (Girardi and Ergün 2013). Let our sequence have T observations indexed from t = 1 to t = T. For institution i, once we obtain the ex-ante VaR α (X i (t)) forecast at t, we compare it with the observed return X i (t) of institution i at time t. The "VaR hit sequence" is defined by where h i (t) = 1 means an exceedance at time t. Suppose the size of the subsample {h i (t)|h i (t) = 1, t = 1, . . . , T} is N. Within this subsample, we consider the "CoVaR hit sequence" If h j|i (t) returns 1, it means both X j (t) ≤ CoVaR j|i α,t and X i (t) ≤ VaR α (X i (t)). Let the size of the subsample {h j|i (t)|h j|i (t) = 1 and h i (t) = 1, t = 1, . . . , T} be x. The unconditional Kupiec's test statistics has the form

Modification of the Generalized Kupiec's Test
The generalized Kupiec's test based on the conditional distribution assumes that the VaR estimate has no error. In other words, the uncertainly in the estimate of VaR is not accounted for in the test. It is obvious that the calculation of CoVaR is heavily related to the VaR estimate (cf. Equation (12)), and thus, the passing of the generalized Kupiec's test does not guarantee what is claimed with the test. An approach that would provide more accurate results amounts to testing the VaR and CoVaR at the same time. We include the consideration of VaR by examining the joint distribution P(X j ≤ CoVaR j|i α,β , X i ≤ VaR α (X i )) = αβ. Using the notation from Section 5.1, we can show that the modified test that achieves simultaneous testing of VaR and CoVaR has the following form: In practice, the statistic of the modified test yields a little bit different results from those of the generalized test because of the inclusion of VaR in the testing.

Goodness of Fit
For marginal processes, we use QQ-plots and density plots for graphical diagnosis, and two-sample KS and AD tests to check the quasi-likelihood-based estimators. The two-sided two-sample KS test statistic is where F 1 (x) is the proportion of sample 1 values less than or equal to x and F 2 (x) is the proportion of sample 2 values less than or equal to x. In our case, sample 1 is the generated Student's t scenarios, while sample 2 is the generated NTS scenarios. For the copula structure, we make use of several in-sample tests and error rates, including univariate KS, AD, AD squared, Mean Average Error (MAE) and Mean Squared Error (MSE). They measure the distance between the CDF of the empirical copula and the other candidate copulas. Their formal definitions are given by (Kim and Volkmann 2013) where u t = F j (X j (t)), v t = F i (X i (t)), and F j , F i are empirical marginal CDFs. C(·, ·) represents the estimated NTS or Archimedean copulas, and C D (·, ·) is the non-parametric empirical copula distribution given by the Deheuvels copula, The univariate KS statistic focuses on the superior distance while the AD statistic emphasizes the tail by a weighting denominator. The AD 2 further considers the sum of distances rather than solely the most extreme distance. MAE and MSE are both error rates, whereas the MSE stresses the large deviation by weighting a squared deviation.

Data
We tested the model on the daily S&P 500 index (SPX) and Dow Jones index (DJI) close price from 3 January 1928 to 31 December 2020. There were 24,444 samples in total. We applied the Archimedean and NTS copulas to model the systemic risk by CoVaR. The implementation is based on the R package "temStaR." 1 The scatter plot and the empirical copula plot of SPX-DJI innovation pairs estimated by kernel density are presented in Figure 1. The SPX and DJI are strongly correlated especially in the lower and upper tails. Therefore, a powerful copula must reflect the dependence structure of both tails.

Results
Backtesting is a widely used out-of-sample test, while a goodness of fit test is an in-sample test. We show in this section that NTS innovations and NTS copula perform better than traditional ones on both types of tests.
We estimate the standard NTS marginal process via the fitted unconditional Student's t distribution first. We also demonstrate Student's t, NIG, and general hyperbolic (GH) distributions as benchmarks. Figures 2, 3 and 6 show that the NIG, GH, and standard NTS fit the marginal density very closely to the empirical density and much better than the Student's t distribution.
We implement the procedures from Section 4.1 to test our applied quasi-likelihoodbased method. The QQ-plots (Figures 4 and 5) of Student's t and NTS share the same pattern, and are overlapped in density plot ( Figure 6). They both reveal that the NTS is indistinguishable from the Student's t fitting of the innovations.  We performed two-sample KS and AD tests for the Student's t and quasi-likelihoodbased NTS distributions, where the SPX and DJI data were both from 3 January 2000 to 9 December 2009 with 2500 samples for each time series. Then we repeated the tests 500 times with one-day-ahead rolling window keeping the width of 2500 samples. The p-values of the two-sample KS and AD statistics with window rolling are shown in Figures 7 and 8. The Figures imply that NTS and the Student's t distributions are quite indistinguishable. Therefore, the efficient algorithm for estimating NTS innovations is feasible and consistent.   Figure 9 illustrates the density for the empirical copula and the two-dimensional NTS copula. The x-y axis represents the CDF of marginal distributions from 0 to 1, and the z-coordinate shows the estimated copula density c(u, v). The NTS copula density captures both the upper and lower tail dependence precisely, which means NTS coupla reflects the tail behavior of empirical coupla and hence empirical joint distribution. However, the Archimedean copulas can only capture at most one tail dependence because of their properties. Table 2 supports the point that the NTS copula is superior to the Archimedean copula benchmarks, because the KS statistic of the NTS copula is the lowest. The contour plot (Figure 10) of the copula also suggests that our NTS copula fits the empirical copula much better than the Frank copula, which is the best copula among Archimedean copula benchmarks according to Table 2.

Backtesting
We applied the rolling/sliding window approach to do the backtesting with a window width equal to 500 days. After estimating all the parameters and CoVaR on day t, the window moved forward by one day, with the window width remaining unchanged. Table 3 gives the backtesting results from the modified Kupiec's test corrected by the joint distribution. There were 23,943 samples in total for backtesting, and the expected numbers of exceedances are listed in the last column as a benchmark. According to Table 2, the Frank copula is still the best copula in the Archimedean copula family in light of the number of exceedances. However, the Lévy-Clayton copula produces results far closer to the expected number than the Archimedean copula family, which means it has good performance in the tail. The NTS copula has a similar performance to the Lévy-Clayton copula but with a much more accurate fit to the global empirical copula.

CoVaR and ∆CoVaR
Figures 11 and 12 display the estimated CoVaR when α = β = 0.01 and α = β = 0.05, respectively. In particular, Figure 11 depicts the date when the exceedances of CoVaR derived from the NTS copula occurred. To make the Figure distinguishable and clearer, we only show the results of the Frank copula as a representative of the Archimedean copula family because it has the best in-sample test performance except for the NTS. The estimated CoVaRs of the Archimedean copula family is far higher than those of the NTS copula.  Figures 13 and 14 show the estimated ∆CoVaR when α = β = 0.01 and α = β = 0.05, respectively. The Archimedean copula family has a dramatic change of ∆CoVaR whereas the NTS copula is more stable. The stable values of change can be comparable with previous crisis and provide more information. We can also use the relative difference between CoCVaR and CoVaR for the early warning of crisis, where the relative difference (RD) is defined by Figures 15 and 16 show the RD derived by the NTS copula when α = β = 0.01 and α = β = 0.05, respectively. 1929-10-28 1943-04-09 1946-09-03 1948-11-03 1950-06-26 1955-09-26 1962-05-28 1986-09-11 1987-10-19 1989-10-13 1991-11-15 1993-02-16 1994-

Conclusions
We applied the NTS copula, which belongs to the Lévy copula family, for modeling CoVaR. Our results obtained from extensive backtesting and goodness of fit testing show that the NTS copula models the dependence structure more accurately than traditional copulas. We also proposed a modified Kupiec's test that accounts for the uncertainty of the estimated VaR. This translates to a test that is more stable when the data sizes are not of sufficient lengths.