I’m trying to estimate CoVaR using bivariate DCC GARCH in R. The concept of CoVaR is the dependence adjusted of VaR, which was first introduced by Adrian and Brunnermeier (2011). However, this original definition of CoVaR presented some limitations, which called for a modified version proposed by Girardi and Tolga Ergun (2013). For given confidence levels α and β,
Girardi et al. suggests a three-step procedure is used to estimate CoVaR.
Step 1: VaR of each institution i is obtained by estimating a univariate GARCH (1,1) models for each time period. (This part I already know how to do using the package
rugarch— make forecast return and volatility and find the quantile.)
Step 2: For the return of institution i and j, set up a bivariate GARCH model with DCC specification to estimate the pdf of (Xi,Xj). (I’m also fine with this part — I know how to solve for the parameters using the
Step 3: If we start from the definition of CoVaRi|jβ,t above, by definition Pr(Xit≤VaRiα,t)=α so
Once VaRiα,t and pdf of (Xi,Xj) have been estimated in previous two steps, CoVaR can be obtained by numerically solving the equation,
where ft(Xi,Xj) is the bivariate density of (Xi,Xj).
This is the part that got me stumped. How can I obtain the bivariate density from the DCC GARCH parameters? And after I have the density figured out, how can I solve for the integral bound to get CoVaR?
Any input is appreciated, thanks!
1. Tobias Adrian and Markus K. Brunnermeier. CoVaR. Technical report, National Bureau of Economic Research, 2011.
2. Giulio Girardi and A. Tolga Ergun. Systemic risk measurement: Multivariate GARCH estimation of CoVaR. Journal of Banking & Finance, 37(8):3169–3180, 2013.