Modeling the Dependency Structure of Integrated Intensity Processes

This paper studies an important issue of dependence structure. To model this structure, the intensities within the Cox processes are driven by dependent shot noise processes, where jumps occur simultaneously and their sizes are correlated. The joint survival probability of the integrated intensities is explicitly obtained from the copula with exponential marginal distributions. Subsequently, this result can provide a very useful guide for credit risk management.


Introduction
Modeling and evaluation of credit risk have been popular subjects for the last decade. They have become crucial but complicated mathematical issues. Such problems are often subtle because of the complex structures of credit risk and the correlation structures of the overall economy. It is well known that there are two types of mathematical models of credit risk: the structural approach and the intensity-based approach (see Ma and Kim [1] for a general reference). Our current paper concerns the latter.
Since Jarrow and Turnbull [2], one credit risk modeling development has been based on the default intensity of a Poisson process. The Poisson process counts the number of events and the time that these events occur in a given time interval. Compared to the Poisson process, a key feature of a Cox process (see Cox [3] and Lando [4]), also known as a doubly stochastic Poisson process, is that its arrival rate is also stochastic. It is well known that the Cox process can be applied in diverse areas such as insurance, finance, neuroscience, astrophysics etc. The properties of the Cox process are described in detail in Grandell [5]. One of the intensity processes within the Cox process which can be used to measure the impact of rare events is the shot noise process (refer to Klüppelberg and Mikosch [6]).
Joint events impact financial markets significantly because of a disproportionately large amount of joint defaults by different counterparties. One of the best ways to understand the dependence structure of joint defaults is by using copulas. A copula is a function that links univariate distributions on [0, 1] n to a joint multivariate distribution on [0, 1]. The existence of such a function is guaranteed by the well-known Sklar theorem (See Nelsen [7] for more details). Choe et al. [8] firstly developed a theoretical framework addressing the joint distribution and provided a general equation for time-dependent copulas related to stochastic processes.
A copula approach has been developed to deal with default correlations. Li [9] introduced the one-factor Gaussian copula and showed how to build a multivariate distribution of survival times for the valuation of credit derivatives. To relax the assumption of the Gaussian distribution in the one-factor Gaussian copula, student t-copula [10,11] is commonly used. Also, as shown by Ma and Kim [12], the copula approach enables modeling of default correlation under an intensity-based framework, where the default intensities are driven by dependent jumps. Consequently, given by an intensity-based credit risk model, we focus on the dependence structure of intensities via copula in this study.
In our previous study, we derived a generalized theoretical result for the joint Laplace transform under the Farlie-Gumbel-Morgenstern (FGM) copula, and then applied this to credit default swap (CDS) rates. To extend the dependency structure of the default intensities, we replace the FGM copula by the extended FGM copula (refer to [13]). Thus, we obtain explicitly the joint Laplace transform of integrated intensities, which promptly provides us with the joint survival (or default) probability, by using the copula with exponential marginal distributions. This is then applied to the pricing of CDS rates. Our numerical study shows that the CDS rate has a diverse dependency structure from a financial point of view, depending on the chosen parameters.
The remainder of this paper is structured as follows. In Section 2, we describe the shot noise intensity processes and the extended FGM copula with exponential marginal distributions, respectively. We calculate the joint survival probability of the integrated intensities in Section 3. The relevant jointly asymptotic distributions are obtained in Section 4. In Section 5, we numerically investigate the functional behavior of the CDS rates in terms of the relevant parameters. We present our concluding remarks in Section 6.

Problem formulation
Definition 1 Let (O, F, P) be a probability space, filtered by a nondecreasing right-continuous family {F t , t 2 [0, T]} of sub-σ-fields of F. Let N t be a Poisson process adapted to F. And put λ t be a non-negative stochastic process independent of N t , adapted to F and assumed to be rightcontinuous such that (i.e. R t 0 l s ds < 1 almost surely). If for all 0 t 1 t 2 and u 2 R then N t is call a Cox process with intensity λ t , where F l t ¼ sfl s ; s tg. Recall that Eq (1) satisfies Consider the process L t ¼ R t 0 l s ds, then from Eq (2) we can solve that Eq (3) implies that the problem of finding the probability generating function of N t is equal to the problem of finding the moment generating function of Λ t . A CDS is a financial instrument for swapping the credit risk. The protection buyer of a CDS pays a premium for effectively insuring against a debt default. He receives a lump sum payment if the debt instrument is defaulted. The protection seller of a CDS receives regularly payments from the protection buyer. If the debt instrument defaults they have to pay the agreed amount to the protection buyer. Then, the credit risk of the reference entity is transferred from the protection buyer to the protection seller. Especially, in the worst case about the protection buyer, the reference entity and the protection seller are influenced by common external events. This paper therefore study the issue of the default correlation between protection seller and reference entity.
The quantitative modeling for credit risk is difficult because default events are rare and involve significant losses. In financial market, the default events can be affected by the external events, such as the subprime crisis of 2007 and firms' bankruptcy. To capture the effects of external events, we consider the default intensities within the Cox processes. Intensity processes, which are affected by the external events, can give rise to default events. Generally, these external events result in a simultaneous positive jump in intensity processes. As time goes by, these intensity processes decrease to prevent default events after the arrival of an external event; the decreases continue until another external event occurs. To model this, we consider the default intensity processes as shot noise processes of the following form: is a sequence of dependent and identically distributed random variables with distributions F (s) (x (s) ) and F (re) (x (re) ) and fX ðbÞ j g is a sequence of independent and identically distributed random variables with distribution F (b) (x (b) ). The point processes M t and M 0 t represent the total number of events up to time t with frequency ρ and ρ 0 , respectively. We assume that M t and M 0 t are independent of fX ðsÞ j ; X ðreÞ j g and fX ðbÞ j g, respectively. Note that shot noise intensity processes have the common point process M t to model simultaneous positive jumps between the protection seller and reference entity.
To model the problem of the correlation between jump sizes for the protection seller and reference entity, we use the extended FGM copula, developed by Huang and Kotz [14], in the following form: where α > 0, β > 0, and γ > 0. Then, the extended FGM copula with exponential marginal distributions, used to model simultaneous positive jumps between the reference entity and the protection seller, is given by The joint survival probability In the following, we derive the joint survival probability of the shot noise intensity processes with a extended FGM copula. The piecewise deterministic Markov process (PDMP) theory developed by Davis [15] is a powerful mathematical tool in this context. Using the theory, we derive the joint Laplace transform of the distribution of the vector ðL ðsÞ t ; L ðreÞ t Þ that provides us with the joint survival (or default) probability, where ðL ðsÞ t ; L ðreÞ t Þ, where where F(x (s) ) and F(x (re) ) denote the distributions of the jump sizes, respectively, dC(F(x (s) ), F (x (re) )) implies the joint probability density function of jump sizes using extended FGM copula. This is well-defined process if f is differentiable with respect to each independent variable and the following holds: We consider shot noise intensity processes where jumps occur simultaneously and their sizes are correlated.
Lemma 1 For arbitrary constants γ (s) , γ (re) , κ (s) and κ (re) such that γ (s) ! 0, γ (re) ! 0, κ (s) ! 0 and κ (re) ! 0, the process respectively. Then Lemma 1, namely, the process given by (Eq 6) is a martingale that yields the equality (Eq 9). We can get the joint Laplace transforms of the distributions ðL ðsÞ t ; L ðreÞ t Þ and ðl ðsÞ t ; l ðreÞ t Þ at time t, respectively, from Lemma 2. We merely need the joint distribution ðL ðsÞ t ; L ðreÞ t Þ. However, the joint distribution ðl ðsÞ t ; l ðreÞ t Þ at time t would be used when calculating the joint Laplace transform of the distribution ðL ðsÞ t ; L ðreÞ t Þ. If we set ψ (s) = ψ (re) = 0 in Lemma 2, then we get promptly Proof. Let t 2 = t 1 and t 1 = t 0 in (Eq 11), and then put t 0 ! −1. If we apply (Eq 4) and (Eq 7) to (Eq 11), then we have by direct computation. Now let us take the expectation of Eq (12). Then, the conditional expectation on the left-hand side (LHS) of the equation is changed into the expectation; this is the same as the LHS of (Eq 13). For the right-hand side (RHS), if we choose ψ (s) and ψ (re) as c ðsÞ ¼ x ðsÞ d ðsÞ 1 À e Àd ðsÞ ðt 2 Àt 1 Þ and c ðreÞ ¼ x ðreÞ d ðreÞ 1 À e Àd ðreÞ ðt 2 Àt 1 Þ ; and substitute the resultant ψ (s) and ψ (re) into the expectation on the RHS of Eq (14), the we get the RHS of (Eq 13). Note that (Eq 13) yields the joint survival (or default) probability by setting ξ (s) = ξ (re) = 1.

Pricing CDS rates
To determine the CDS rates, we take a generalized Cox-Ingersoll-Ross model [16] as the interest rate process r t for a zero-coupon default-free bond. We also use a theoretical result Eq (13) for the pricing of a zero-coupon defaultable bond with zero recovery and a deterministic payoff value of 1, paid at t k when a reference entity defaults in (t k−1 , t k ] and the protection seller survives up to t k for 0 k K. Then, the CDS rates, denoted by S, are given by where π, N, and k n denote the recovery rate, the number of payment dates for the protection buyer, and the dates of these payments (i.e., t k n , n = 1, 2, Á Á Á, N), respectively. Refer to Ma and Kim [12].
We calculate the price of the CDS rates numerically and investigate the functional behavior of the CDS rates with respect to the relevant parameters. The parameter values used to calculate the CDS rates are Table 1 shows that the CDS rate increases with the value of the copula parameter p. It also is approximately converging to 0.42 as the values of p are increased. In the financial market, stabler protection sellers must receive higher CDS rates, whereas protection buyers pay lower CDS rates with respect to a reference entity with a relatively high credit rating. The increasing behavior shown in Table 1 is due to the greater impact of the protection seller, but this is a minor issue because the CDS rate is monotonically increasing only slightly as p increases. The calculations of CDS rates caused by changes in the values of δ (s) and α for the protection seller and δ (re) and β for the reference entity are shown in Tables 2 and 3, respectively, for p = 0.1. Table 2 shows that the CDS rates increase when the values of δ (s) and α increase, respectively. However, Table 3 shows that the values of the CDS rate decrease when the values of δ (re) and β increase, respectively. The parameters δ (s) and δ (re) are the exponential decay rates (that is, how quickly the entity can prevent a default after the arrival of a primary event) for the protection seller and reference entity, respectively. In addition, the parameters α and β describe the magnitude of jump sizes for the protection seller and reference entity, respectively. The monotonic behavior with respect to the protection seller and reference entity shown in Tables 2 and  3, respectively, is an obvious effect from a financial point of view. Tables 2 and 3 confirm how fast the CDS rate varies. From the above investigation, we note that the CDS rate has a diverse dependency structure regarding the relevant parameters. This may be monotonic behavior from a financial point of view.

Final remarks
Based on the shot noise intensity processes, we studied the important issue of default correlation in credit risk modeling. To model this, we used a multivariate shot noise intensity process, where jumps occur simultaneously and their sizes are correlated. The joint Laplace transforms of the integrated intensities, which promptly provides us with the joint survival (or default) probability, were obtained from the extended FGM copula with exponential marginal distributions. Also, we applied this result to CDS rates and then our numerical results suggested that CDS rates have a diverse dependency structure with respect to the relevant parameters from a financial point of view. Our result provides a guide for credit risk management. More comparative studies on Gaussian and Student-t copulas and applying our theoretical result to market data for calibration (refer to [13] and [17]) would be helpful to complete the present work. These remain topics for future research.