Next Article in Journal
The Determinants of Equity Risk and Their Forecasting Implications: A Quantile Regression Perspective
Next Article in Special Issue
The Design and Risk Management of Structured Finance Vehicles
Previous Article in Journal
Down-Side Risk Metrics as Portfolio Diversification Strategies across the Global Financial Crisis
Previous Article in Special Issue
Application of Vine Copulas to Credit Portfolio Risk Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probability of Default and Default Correlations

1
Institute of Finance and Big Data, Southwest Jiaotong University, Chengdu 611756, Sichuan, China
2
Department of Finance, Oklahoma State University, Stillwater, OK 74078-4011, USA
J. Risk Financial Manag. 2016, 9(3), 7; https://doi.org/10.3390/jrfm9030007
Submission received: 30 March 2016 / Revised: 31 May 2016 / Accepted: 23 June 2016 / Published: 5 July 2016
(This article belongs to the Special Issue Credit Risk)

Abstract

:
We consider a system where the asset values of firms are correlated with the default thresholds. We first evaluate the probability of default of a single firm under the correlated assets assumptions. This extends Merton’s probability of default of a single firm under the independent asset values assumption. At any time, the distance-to-default for a single firm is derived in the system, and this distance-to-default should provide a different measure for credit rating with the correlated asset values into consideration. Then we derive a closed formula for the joint default probability and a general closed formula for the default correlation via the correlated multivariate process of the first-passage-time default correlation model. Our structural model encodes the sensitivities of default correlations with respect to the underlying correlation among firms’ asset values. We propose the disparate credit risk management from our result in contrast to the commonly used risk measurement methods considering default correlations into consideration.

1. Introduction

In practice, one of the biggest challenges of credit risk models is modeling correlations between default events, between portfolios of correlated defaultable claims and between credit quality changes. Dependence among credit default events is indispensable for determining credit risk measures used to allocate capital for solvency and to investigate systemic risk. There is a huge market for asset securitization for default swaps (for example CDO, synthetic CDO etc.) and default swaps have a basket structure with many underlying references. It is clearly important to understand the prices and credit risk of tranches of CDOs (see [1]). Risk dependence of simultaneous defaults by financial institutions is also crucial in the study of the stability of the financial system. One of the most important measures of credit risk with dependence is the probability of default under the dependence in the portfolio or the financial system. It is also important to measure the joint probability of correlated defaults and probability of portfolio default.
In this paper, we study the credit risk with dependence to develop the first step in measuring default of single firm under the asset correlation assumption. We also approach the default of multiple firms with a structural model, initiated by [2,3,4]. We also provide an improved distance-to-default with dependence, with the hope of providing a better credit rating system. Our structural model can be used to bond portfolios and other structure products.
Default correlation is essential for managing credit risk, and the structure of default correlation is crucial to price structured credit derivatives as well as in the stability of financial institutions. Lucas [5] states that general conditions can produce non-zero default correlation and specifies that general economic conditions and more specific industry and regional factors produce systematic and idiosyncratic credit risk exposure analogous to systematic and unsystematic equity price risk. Significantly biased perceptions of credit risk exposure can derive from omission (Das [6]) or underestimation (Jorion [7]) of default correlation. The most widely recognized examples are CreditMetrics; Crouhy, et al. [8] and CreditRisk + in Gundlach and Lehrbass [9] with incorrect dependence consideration due to the misunderstanding of the joint default probability. Giesecke [10] induces default correlation differently in a system of independent asset value changes through firm-specific default barriers that are themselves correlated random variables. Li [11] defines a random variable called “time-until-default”.
Unlike the n = 2 case Rebholz [12], Zhou [13], Valužis [14] and Metzler [15], which use separable variable methods to deal with initial-boundary conditions by formal series, our formulas for evaluating probabilities of (multiple) defaults are precise and computational accessible. Solutions from previous literature in closed form are given by an infinite series of modified Bessel functions. Not only with computational intensity but also with inconsistency, The default correlations and probabilities of defaults from previous literature frequently attain absolute values greater than 1 as shown in Li and Krehbiel [16]. The evaluations of probability of default and default correlation from the previous literature are not only computationally intense but also inconsistent. Zhang and Melnik [17] discuss the first passage time for multivariate jump-diffusion processes without any explicit default correlation or default probability. Li and Krehbiel [16] indicate some numerical results to the inconsistency of their default rate and default correlation. Sacerdote et al. [18] recently expressed the joint density of the first passage time in terms of the solutions of a system of Volterra-Fredholm first kind integral equation, and the solution is again given by the same method of Zhou [13] in terms of the formal power series with modified Bessel function of the first kind. Erlenmaier and Gersbach [19] regard the optimal stopping time for the default as a non-random variable to deduce the default correlation and completely ignore the main difficulty in understanding the joint distribution of first-passage-times of two correlated asset firms. Erlenmaier and Gersbach [20] provide a connection between adverse macroeconomic shocks and loan default correlations and expresses the default correlation in integral form so that the partial derivatives of the default correlation can be traced analytically. But they do not discuss the default correlation itself and the relation with the default probability. They only discuss various properties of the default correlation with respect to the distances-to-default from implicit relation.
In this paper, we specify the linkage between companies as the correlation of logarithmic changes of underlying asset values. As mentioned in Li and Krehbiel [16], there is an inconsistency between the stochastic assumptions of Merton’s firm-specific default probability model with the bivariate first passage time model of default correlation. We first evaluate the probability of single firm default with dependence by changing the driving Brownian motions among firms into a combination of driving Brownian motions through a linear algebra technique. Then we evaluate firm-specific default probabilities and default correlations under consistent stochastic assumptions in Propositions 1 and 2. To understand multiple defaults, we need to further extend the results from two joint default probability to one of n-joint default probability. By using the reflection principle in partial differential equations to solve the Kolmogorov forward equation with initial-boundary conditions, we resolve this problem to derive a closed form of the probability of multiple joint defaults in Theorem 1.
We further propose the probability of default algorithm for the dependent firms. The default algorithm in credit risk management follows from evaluating single firm default, two firms default, any k firms default simultaneously, and all firms default. For two firms case, Li and Krehbiel [16] show that the difference between KMV’s and J. P. Morgan’s approach and our probability of joint default is not trivial, among other things. Understanding default correlation and probabilities of the default state space is crucial for proper measurement of credit risk. Bond portfolio management is the quintessential example of default correlation application. Financial stability and the systemic risk are very important in the financial market. Our analysis provides a proper tool and technique to measure those risks. We propose a mixed pair default measure M D i ( t ) to incorporate the dependence by adding related default correlations. This measure should provide more information than the probability of default of a single firm from easy evaluations, and also give the credit risk estimate of the specific firm in the system.
The default probability and default correlation for more than 3 firms with constant drifts and constant diffusions can be consistently obtained in this paper. Our contribution concerns reconciliation of the assumptions of independent log asset changes with the correlated multivariate processes of the first-passage-time default correlation model. While Jarrow [21] attacks the use of Merton’s structural credit risk model as a valid source of implied default probabilities due to the violation of the model’s critical assumption that all the firm’s assets trade in frictionless markets. The mixed default measure M D i ( t ) reflects all the market information and the default probability that is directly related to the i-th firm. Even if the probability of the specified firm default is smaller than one of other firms, the mixed default measure of the specified firm would be way bigger than one of the other firms due to the larger impact of default correlation. The default correlation matrix ρ D ( t ) = ( ρ i j D ( t ) ) 1 i , j n might be useful for regulators to access the possible default contagion from one financial institution to another. The mixed pair default M D i ( t ) measures all the possible default correlation from the market to this specified firm as well as its own probability of default under the dependence by adding the row of the default correlation matrix. On the other hand, if we add the column of the default correlation matrix, the sum of default correlations on the i-th column is the total effect of dependence from the specified i-th firm to rest of firms in the market. This quantity could be used to measure the significance of the firm in the market.
The remainder of this paper is organized as follows. Section 2 presents the first-passage-time model for probability of default with dependence. Section 3 provides risk management with default correlation. Section 4 gives numerical results of the paper and compares them with previous results. Section 5 concludes our findings. We give a complete proof of main results in Appendix.

2. First-Passage-Time Model for Probability of Default and Default Correlation

In the basic model for default correlations based on the first-passage-time model of Black and Cox [4], a firm defaults when its asset value breaches the default barrier. We assume that there are n credit entities with a collection of n firms (the i-th credit entity is referred to as the i-th firm). The i-th firm defaults whenever this firm’s assets falls below the i-th firm specific threshold.
Assumption 1. 
The dynamics of total asset values of both firms is given by
d ln V = μ d t + Ω d W ( t )
where ln V T = ( ln V 1 ( t ) , , ln V n ( t ) ) , μ T = ( μ 1 , , μ n ) are constant drift terms for each firm and W T ( t ) is the independent standard n-dimensional Brownian motion. The log-asset covariance matrix is given by Ω · Ω = ρ i j σ i σ j with a constant correlation ρ i j = C o r r ( d ln V i , d ln V j ) for 1 i , j n .
Assumption 2. 
For each firm i, there exists a time-dependent value C i ( t ) such that the firm continues to operate and meet its contractual obligations if V i ( t ) > C i ( t ) . The firm defaults on all of its obligations immediately if V i ( t ) C i ( t ) , and some form of corporate restructuring takes place, for 1 i n .
For simplicity, we assume the drift vector and variance matrix of the total asset values’ logarithmic change are constant. This is the same setup as Zhou [13] for n = 2 . For simplicity, we focus on the the default threshold level C i ( t ) = e λ i t K i for constant λ i and K i for i = 1 , 2 , , n as proposed by Black and Cox [4]. Black and Cox [4] interpret C i ( t ) as the smallest possible value required by the safety covenant of a debt contract. One can extend to a random barrier and random interest rates as studied in Nielsen et al. [22], Leland [23] and Leland and Toft [24].
If the default status are described as Bernoulli random variables,
D i ( t ) = 1 if the firm i defaults by  t 0 otherwise . 1 i n
then we have, for 1 i n ,
D i ( t ) = { τ i t } , τ i = min { t 0 : V i ( t ) e λ i t K i }
as the first-passage-time that the i-th firm reaches its default threshold.
Define the default correlation C o r r ( D i ( t ) , D j ( t ) ) = ρ i j D ( t ) between the firms i and j over the time interval [ 0 , t ] as
C o r r ( D i ( t ) , D j ( t ) ) = C o v ( D i ( t ) , D j ( t ) ) V a r ( D i ( t ) ) · V a r ( D j ( t ) ) = E [ D i ( t ) · D j ( t ) ] - E [ D i ( t ) ] · E [ D j ( t ) ] V a r ( D i ( t ) ) · V a r ( D j ( t ) )
Default correlation is central to determining the joint default probability that is the probability of multiple defaults. If P ( D i ( t ) = 1 ) = P i , then
P ( D i ( t ) = 1 and D j ( t ) = 1 ) = P i P j + ρ i j D ( t ) ( 1 - P i ) ( 1 - P j ) P i P j P ( D i ( t ) = 1 or D j ( t ) = 1 ) = 1 - ( 1 - P i ) ( 1 - P j ) - ρ i j D ( t ) P i P j ( 1 - P i ) ( 1 - P j )
We need to first derive the probability of i-th firm default under Assumptions 1 and 2.
Proposition 1. 
Under Assumptions 1 and 2, the probability of i-th firm default is given by
P ( D i ( t ) = 1 ) = P i = N ( - ln ( V i , 0 / K i ) Σ i t - μ i - λ i Σ i t )
+ ( V i , 0 K i ) 2 ( λ i - μ i ) σ i 2 N ( - ln ( V i , 0 / K i ) Σ i t + μ i - λ i Σ i t )
and the distance-to-default at the time t for the i-th firm is given by
d d i ( t ) = ln ( V i , 0 / K i ) + ( μ i - λ i - 1 2 Σ i 2 ) ( T - t ) Σ i T - t
where Σ i 2 is the i-th eigenvalue of Ω · Ω T and N ( · ) is the cumulative probability distribution function for a standard normal random variable.
The proof of Proposition 1 is given in Appendix.
Then, the i-th firm-specific probability P ( D i ( t ) = 1 ) = P ( τ i t ) of default given by Proposition 1 is derived under our Assumptions 1 and 2 with underlying asset correlations. The probability of i-th firm default is given by
p i ( t ) = P ( D i ( t ) = 1 ) = N ( - ln ( V i , 0 / K i ) σ i t - μ i - λ i σ i t )
+ ( V i , 0 K i ) 2 ( λ i - μ i ) σ i 2 N ( - ln ( V i , 0 / K i ) σ i t + μ i - λ i σ i t )
It is derived without the correlation Assumption 1 for each individual firm. We have P i ( t ) = p i ( t ) if and only if ρ i j = 0 for all i j . The probability p i ( t ) of the i-th firm default is obtained assuming independence of the stochastic process in Merton (1974). Hence the firm-specific default probabilities (3) are not consistent with the default correlation of Assumptions 1 and 2, as shown in Li and Krehbiel [16]. The distance-to-default d d i ( t ) is also different from the distance-to-default with independent asset processes:
d i ( t ) = ln ( V i , 0 / K i ) + ( μ i - λ i - 1 2 σ i 2 ) ( T - t ) σ i T - t
Proposition 2. 
The probability that either i-th firm or j-th firm defaults and the probability that both firms default under Assumptions 1 and 2 given { u i , λ i , σ i , V i , 0 , K i , ρ } are given by
P ( D i ( t ) = 1 o r D j ( t ) = 1 ) = 1 - F ( b i , b j , t ) = 1 - N 2 , ρ i j ( b i σ i t , b j σ j t )
+ N 2 , - ρ i j ( - b i σ i t , b j σ j t ) - N 2 , ρ i j ( - b i σ i t , - b j σ j t ) + N 2 , - ρ i j ( b i σ i t , - b j σ j t )
and the default correlation of the i-th firm and the j-th firm is
ρ i j D ( t ) = P i + P j - P i P j - P ( D i ( t ) = 1 o r D j ( t ) = 1 ) P i ( 1 - P i ) P j ( 1 - P j )
where b i = - ln [ K i V i ( 0 ) ] + ( μ i - λ i ) t , P i is given in Proposition 1 for 1 i , j n , F ( b i , b j , t ) = P ( τ i j > t ) with τ i j = τ i τ j = min ( τ i , τ j ) and the probability of the standard bivariate normal distribution on y 1 x 1 and y 2 x 2 is denoted by
N 2 , ρ ( x 1 , x 2 ) = - x 2 - x 1 1 2 π 1 - ρ 2 e - y 1 2 - 2 ρ y 1 y 2 + y 2 2 2 ( 1 - ρ 2 ) d y 1 d y 2
and n 2 , ρ ( y 1 , y 2 ) = 1 2 π 1 - ρ 2 e - y 1 2 - 2 ρ y 1 y 2 + y 2 2 2 ( 1 - ρ 2 ) is the standard bivariate normal distribution density function denoted by N 2 ( 0 , 1 , ρ ) for random variables Y i ( i = 1 , 2 ) with mean zero and variance 1 and correlation ρ = C o v ( Y 1 , Y 2 ) .
By Proposition 2, the probability that both i-th firm and j-th firm default under Assumptions 1 and 2 is given by
P ( D i ( t ) = 1 and D j ( t ) = 1 ) = P i + P j - 1 + N 2 , ρ i j ( b i σ i t , b j σ j t )
- N 2 , - ρ i j ( - b i σ i t , b j σ j t ) + N 2 , ρ i j ( - b i σ i t , - b j σ j t ) - N 2 , - ρ i j ( b i σ i t , - b j σ j t )
which is definitely different from
P ( V i ( t ) e λ i t K i and V j ( t ) e λ j t K j ) = N 2 , ρ i j ( b i σ i t , b j σ j t )
where (6) is used to evaluate the default correlation in Crouhy et al. [8] (also see Bielecki and Rutkowski [25], p. 117). The essential difference of these evaluations is discussed in [16], Section 2.1.
The proof of Proposition 2 follows from the same method in [16]. Formula (6) evaluates the probability of joint default on the horizon date, while (5) evaluates the joint default probability prior to and inclusive of the horizon date. Zhou’s analogous default probability is achieved by solving the Kolmogorov forward equation with initial-boundary conditions, with separable variable method and formal series of eigenfunctions of the boundary value problem. Zhou’s closed-form solution is given by an infinite series involving the double integral of Bessel functions. The series solution is formal and not necessarily convergent because the domain of the boundary value problem is unbounded. Numerical evidence for those inconsistencies is illustrated in [16].
Theorem 1. 
The evaluation of the probability that either firm defaults under Assumptions 1 and 2 given { μ i , λ i , σ i , V i , 0 , K i , ρ i j } 1 i , j n is given by
P ( τ t ) = P ( D 1 ( t ) = 1 o r D 2 ( t ) = 1 , o r D n ( t ) = 1 ) = 1 - F ( b 1 , b 2 , , b n , t )
F ( b 1 , b 2 , , b n , t ) = i = 1 2 n ( - 1 ) s i g n Q i N Ω Q i Ω Q i T ( s i g n ( z 1 ) b 1 , s i g n ( z 2 ) b 2 , , s i g n ( z n ) b n )
where τ = min { τ 1 , , τ n } , ( s i g n ( z 1 ) , , s i g n ( z n ) Q i with s i g n ( z i ) = ± 1 , s i g n Q i is the sum of negative 1 in ( s i g n ( z 1 ) , , s i g n ( z n ) , Ω Q i = d i a g ( s i g n ( z 1 ) , , s i g n ( z n ) · Ω , and N A ( x 1 , , x n ) = - x n - x 1 1 ( 2 π t ) n / 2 det A exp ( - 1 2 t z T A z ) d z 1 d z n .
We solve the Kolmogorov forward equation with initial-boundary conditions using the reflection principle in partial differential equations. This is an extension of n = 2 case in [16]. We first solve the Kolmogorov forward equation with initial condition from the constant coefficients in the partial differential operator. The solution gives the n-dimensional Gaussian distribution function that satisfies all three boundary conditions automatically, except the boundary on the finite regions x i = b i for 1 i n . This is the main difficulty that led Rebholz [12], Zhou [13] and Valužis [14] to use the separable variable methods to obtain a formal series solution, even in n = 2 case. In contrast, we use the reflection principle first to shift the boundary region as the boundary of the first quadrant, and then extend the solution to the whole plane R n by odd extension. After identifying the involution through the odd extension and using the Dirac function properties, the probability of at-least one default among all n firms is obtained in Theorem 1. For the completeness of the result, we present the detailed proof in Appendix. The probability of at-least one default in Theorem 1 is given by 2 n terms with standard multivariate normal distribution function.
For all the firms under Assumptions 1 and 2, we can use Propositions 1 and 2 and Theorem 1 to determine any possible combination of default events. We outline the basic scheme and hope this will be useful in determining systemic risk.
  • Probability of single firm default We first use Proposition 1 to evaluate the probability P ( D i ( t ) = 1 ) = P i ( t ) of i-th firm default at the time t. The probability P i of the i-th firm default incorporates with other firms through Assumptions 1 and 2, and this is essentially different from the probability p i of the i-th firms default in (3) where the i-th firm is independent of the rest firms.
  • Probability of two firms default For any pair of two firms, we first evaluate the probability P ( D i ( t ) = 1 or D j ( t ) = 1 ) that either i-th firm or j-th firm defaults by Proposition 2. There are n 2 many pairs of two firms from those n firms. Similarly for k firms, there are n k many k firms from those n firms.
    Hence, the probability P ( D i ( t ) = 1 and D j ( t ) = 1 ) = P i j ( t ) that both i-th firm and j-th firm default at the time t is given by
    P ( D i ( t ) = 1 and D j ( t ) = 1 ) = P i + P j - P ( D i ( t ) = 1 or D j ( t ) = 1 )
  • Probability of three firms default For any triple firms ( i , j , k ) , we apply the sub-matrix taking the i , j , k rows and columns and the probability P ( D i ( t ) = 1 or D j ( t ) = 1 or D k ( t ) = 1 ) that either i-th firm or j-th firm or k-th firm defaults by Theorem 1 for n = 3 . Hence, the probability P ( D i ( t ) = 1 and D j ( t ) = 1 and D k ( t ) = 1 ) = P i j k ( t ) that both i-th firm and j-th firm and k-th firm default at the time t is given by
    P i j k ( t ) = - P i ( t ) - P j ( t ) - P k ( t ) + ( P i j ( t ) + P i k ( t ) + P j k ( t ) ) + P ( D i ( t ) = 1 or D j ( t ) = 1 or D k ( t ) = 1 )
  • Probability of k firms default For k 2 , we apply Theorem 1 for all k × k sub-matrix from Ω Ω T to evaluate the probability P ( j = 1 k D i j ( t ) = 1 ) , and then evaluate the probability P i 1 i 2 i k ( t ) that k firms (from i 1 -th to i k -th firms) both default at the time t from the following identity and previous steps,
    P ( j = 1 k D i j ( t ) = 1 ) = j = 1 k P i j ( t ) - j 1 , j 2 P i j 1 i j 2 ( t ) + + ( - 1 ) k - 1 P i 1 i 2 i k ( t )
  • Probability of all n firms default Repeat the previous step for 2 k n - 1 . Evaluate the probability P ( τ t ) that either firm defaults among n firms by Theorem 1. Therefore, the probability P 12 n ( t ) that all n firms default at the time t is derived from the following identity.
    P ( τ t ) = i = 1 n P i ( t ) - i , j P i j ( t ) + + ( - 1 ) n - 2 P i 1 i 2 i n - 1 ( t ) + ( - 1 ) n - 1 P 12 n ( t )

3. Risk Management with the Default Correlation

A portfolio manager would not only be concerned with the default of a single firm but also be worried with the probability of multiple defaults in the portfolio. Similarly, a regulator would have to estimate the probability of multiple defaults of a number of financial institutions (banks) if there is no bailout to let a default of a single financial institution (the too big to fail puzzle). Modeling credit risk from a portfolio perspective is indispensable if derivatives on credit portfolios are priced and risk measures of credit portfolios are to be computed. New regulatory rules make it mandatory for financial institutions to extend their risk measures from a single-contract setting to a portfolio perspective.
An effective credit risk measurement in a portfolio consists of three important components: (1) the probability of default for each individual firm over various investment horizons; (2) the joint probability of default between any two firms in the portfolio during various investment horizons; and (3) the magnitude of financial loss in the case of default. As noted in the literature, the precise quantity of the default correlations is the most crucial and the most challenging part of credit risk analysis in a portfolio.
The structural credit risk models assume that all of the firms assets trade and their values are observable. The major drawback of the structure credit risk models is that both firms’ assets are not tradable and their volatilities are unobservable. See Jarrow [26] for more debates between structural models and reduced form models. In practice, one can calibrate the observable stock prices and stock price’s volatility to estimate the firms’ asset values and volatilities, and further minimize the error between the market observed stock price and volatility and the model’s asset value and volatility. In industry practice, the structural models are still used incorrectly as a theoretical model for estimating the firm’s default probability. Reduced form model is adapted to estimate default probabilities from historical default data. Jarrow [21] shows that the implied default probability form a structural credit risk model and the default correlations obtained from credit risk copula models lead to mis-specified estimates.

3.1. CreditRisk + Incorporating Default Correlations

Credit Suisse Financial Products (CSFP) introduced the CreditRisk + in October 1997. CreditRisk + applies techniques from actuarial mathematics in an enhanced way to compute the probabilities for portfolio loss levels. Among other fundamental ideas in CreditRisk + , the linear relationship between the systemic risk factors and the probabilities of default is assumed; the defaults of obligators are assumed independent; correlations among obligators are assumed implicitly due to common risk factors which drive the probability of defaults.
The original CreditRisk + does not allow correlations in modelling default events. Akkaya et al. [27] extend the CreditRisk + that is able to model default correlations among segments while preserving the analytical solution for the loss distribution. The portfolio loss X is given by A S 1 A μ A , where 1 A is the default indicator for obligor A with P [ 1 A = 1 ] = p A and μ A is the exposure net of recovery for a single-segment portfolio S. For the N different and correlated segments S 1 , , S N , the total loss is given by
X = k = 1 N X k = k = 1 N A S k 1 A μ A
The expected loss and variance of X are given by
E L ( X ) = k = 1 N A S k p A μ A , V a r ( X ) = k = 1 N V a r ( X k ) = k = 1 N σ S k 2 E L 2 ( X k ) + k l C o r r ( S k , S l ) σ S k σ S l E L ( X k ) E L ( X l ) + k = 1 N A S k p A μ A 2
The variance V a r ( X ) contains a term k l C o r r ( S k , S l ) σ S k σ S l E L ( X k ) E L ( X l ) which is the default correlation measured in the extended CreditRisk + in [27]. The analytic solution for the default correlation depends on the assumption on the portfolio loss distributions. The multivariate structure with nontrivial correlation leads a pragmatic simplification to reduce into a single-segment model.
Due to the non-uniqueness in the parameter estimation, one has to assess the impact of the choice of parameters on the overall loss distribution. This reduced-form approach to default correlation does not reflect any underlying structure of the assets. The conditionally independent defaults (CID) approach creates default dependence among firms through the dependence of the firm’s intensity processes on a common set of state variables. Contagion models extend the CID model to account for the empirical information of default clustering and to include the existence of credit risk contagion mechanisms between firms.

3.2. CreditMetrics Incorporating Default Correlations

From the CreditMetrics TM Technical Document, we examined several rating changes and defaults in order to establish that such correlations indeed exist. The direct method to estimate joint rating change likelihoods is to test credit rating time series across many firms which are synchronized in time as a statistical approach. The advantage of this method is that it is independent of the underlying process and the joint distribution shape and the underlying correlation. The disadvantage is that it treats all firms with a same credit rating identical. To improve this method, the CreditMetrics TM estimates credit quality correlation through bond spreads. Whenever bond price histories are available, it is reasonable to estimate a certain type of credit correlation by extracting credit spreads from the bond prices and then estimating the correlation in the movements of credit spreads. This adapts a reduced-form approach. The biggest drawback of this approach is that the bond spread data is notoriously scarce, especially for low credit quality issues.
For the structural model, CreditMetrics uses the correlation between equity returns as a proxy for the correlation of asset returns. The drawback of this replacement is overlooking discrepency between equity and asset correlation. It is better than using a fixed correlation as well as much more readily availability than credit spreads or actual joint rating changes. Due to the scarcity of data for many obligors and the impossibility of storing a huge correlation matrix, one has to rely on correlations within a set of indices and mapping scheme to build the obligor-by-obligor correlations from the index correlations.
Copula methods have been extensively used in both industry and academia to assess the joint default probability of groups of obligors. However, the choice of the copula and its calibration is still an issue to debate. The copula reduced form model separates the modelling and estimation of individual default probabilities from the modelling and calibration of the estimation of credit risk dependence (the copula). The inconsistency among the reduced form methods is also unavoidable. The copula model for pricing collateralized debt obligations (CDO’s) and implied default correlations is misused. Jarrow [21] points out that the default correlation obtained from credit risk copula models for computing VaR measures leads to misspecified estimates.

3.3. Proposed Mixed Pair Defaults Incorporating Default Correlations

Default dependence, mostly represented by default correlation, is one of most important features in credit derivatives pricing, hedging and risk management. Despite the inconsistency of the individual default probability and the default correlation, the default correlation measured in Theorem 1 gives the proper indication of the default dependence under Assumptions 1 and 2. We define a mixed pair for the defaults incorporating default correlations.
Under Assumptions 1 and 2, we define a mixed pair of defaults for each firm by ( P i ( t ) , ρ i j D ( t ) ) for 1 i , j n . The mixed pair of defaults indicates the measure of the default probability and the default correlation under the realistic correlated asset assumption. Both measures should play an important role in credit risk management. We can simply use the sum of the default probability and the default correlation for the default measure under Assumptions 1 and 2:
M D i ( t ) = P i ( t ) + j i ρ i j D ( t ) , 1 i n
where M D i ( t ) is the mixed default measure for the i-th firm and the default correlation ρ i j D ( t ) between the i-th firm and the j-th firm is given in Proposition 2. A single default correlation between two firms can be extended to a default correlation matrix. Theoretically the mixed default measure can be larger even the individual firm specific default probability is quite lower due to the larger default correlation from the underlying asset correlation between this firm and the other highly “toxic” firm. The mixed default measure M D i ( t ) reflects all the market information and the default probability that is directly related to the i-th firm. Even if P i ( t ) < P k ( t ) , the mixed default measure M D i ( t ) > M D k ( t ) due to the impact of default correlations.
The default correlation matrix ρ D ( t ) = ( ρ i j D ( t ) ) 1 i , j n might be useful for regulators to access the possible default mitigation from one financial institution to another. The mixed pair default M D i ( t ) measures all the possible default correlation from the market to this specified firm as well as its own probability of default under the dependence by adding the row of the default correlation matrix. On the other hand, if we add the column of the default correlation matrix, the quantity j i ρ j i D ( t ) is the total effect of dependence from the specified i-th firm to the rest of firms in the market.
The individual asset risk is usually characterized by its return variance in the modern investment theory. The covariance matrix for a portfolio can be a measure of a risk for the portfolio. Similarly, the default correlation matrix (or default covariance) is the key to managing a bond portfolio or other portfolios. One of main risk management tasks is to keep the default risk and the default correlation under control. The ideal situation would spread possible defaults out rather than cluster, if the default cannot be avoided in the risk management. The mixed default measure for each firm or a single defaultable bond would provide an important measure in risk management for fixed income portfolios, bank management, insurance industry and capital management, etc. More empirical analysis on the mixed default measure will be given in a future study.

4. Numerical Analysis

In this section, we show the numerical analysis of our main results. For computation simplification, we assume λ i = μ i ( i = 1 , 2 , , n ) throughout this section.
Credit quality parameters, V K , for the numerical analysis of Propositions 1 and 2 and Theorem 1 are selected from implied standardized distances to default estimated from the Standard and Poor’s 2010 annual cumulative default study, Credit Week: 20 April 2011. Implied standardized distance to default, Z i based on p i ( t ) in (3), can be estimated from the study’s cumulative default rates P ( Z , t ) to the historical default rate A ˜ i ( t ) , under the maintained hypothesis, λ i = μ i by minimizing the sum of squared errors. We have Z i = l n ( V i , 0 / K i ) σ i and
Z = a r g m i n Z t ( P ( Z , t ) t - A ˜ ( t ) t ) 2
Empirical analysis of default correlations can be done by using firms’ data ( V i , 0 , K i , ρ , μ i , λ i ) to compute the default correlation by Proposition 2. We use the Standard and Poor’s 2010 annual cumulative default study to first estimate the standard distance to default from default data.
Table 1 is borrowed from Li and Krehbiel [16], and presents the standardized distance to default for S & P credit ratings implied by the cumulative default probabilities reported in the 2010 cumulative default study.
Assuming the annualized standard deviation of logarithmic change of asset market value is σ 1 = 0.4 , the S & P BBB rated bonds with a ratio of initial asset value to default threshold V 1 , 0 K 1 = 11 , and σ 2 = 0.3 for bonds rated BB by S & P with V 2 , 0 K 2 = 5 , and σ 3 = 0.2 for bonds rated AA V 3 , 0 K 3 = 33 . The constant correlation matrix of the logarithms of asset values is given by parameters ( ρ 12 , ρ 13 , ρ 23 ) :
ρ = σ 1 2 ρ 12 σ 1 σ 2 ρ 13 σ 1 σ 3 ρ 12 σ 1 σ 2 σ 2 2 ρ 23 σ 2 σ 3 ρ 13 σ 1 σ 3 ρ 23 σ 2 σ 3 σ 3 2 = 0.16 0.12 ρ 12 0.08 ρ 13 0.12 ρ 12 0.09 0.06 ρ 23 0.08 ρ 13 0.06 ρ 23 0.04
Table 2 illustrates the probability of default for each asset under the correlation changes. We evaluate the probability of default for each individual firm under Assumptions 1 and 2 in Proposition 1, and compare the probability of default from Merton (1974) under the isolated stochastic process for t = 1 . For simplicity, we have μ i = λ i for i = 1, 2, 3. The ratio V 1 , 0 K 1 = 11 for the BBB rated bonds with σ 1 = 0.4 , V 2 , 0 K 2 = 5 for the BB rated bonds with σ 2 = 0.3 and V 3 , 0 K 3 = 33 for the AA rated bonds with σ 3 = 0.2 .
Under the asset correlation changes, we see how the eigenvalues of the covariance matrix varying. These simple numerical results illustrate the important role of the underlying asset correlations, and the probability of default P 1 behaves more complex than the probability of p 1 . Note that p i in Merton [3] does not vary at all for i = 1, 2, 3. But the probability of default P i is very sensitive to the correlations ρ 12 , ρ 13 , ρ 23 .
We compute the distance-to-default for each individual firm under Assumptions 1 and 2, and compare the distance-to-default from [3] under the isolated stochastic process for t = 1 . For simplicity, we have μ i = λ i for i = 1, 2, 3. The ratio V 1 , 0 K 1 = 11 for the BBB rated bonds with σ 1 = 0.4 , V 2 , 0 K 2 = 5 for the BB rated bonds with σ 2 = 0.3 and V 3 , 0 K 3 = 33 for the AA rated bonds with σ 3 = 0.2 .
Table 3 illustrates the distance-to-default in Proposition 1 under Assumptions 1 and 2. The distance-to-default is the key to the rating of the underlying firm or bonds. We see the changes of d d i through the variations of ρ 12 , ρ 13 and ρ 23 . The distance-to-default d d 1 increases dramatically when the correlation ρ 23 of other two firms increases from 0.5 to 0.25 , and decreases rapidly when ρ 23 increases from 0.25 to 0.5 . This is the essence of correlation. As the correlation of other firms varies, the distance-to-default as well as the credit rating of the first firm may change dramatically. We see the big difference between d d 3 and d 3 . Without considering the underlying asset correlations, the third firm or bonds has the distance-to-default d 3 = 17.383 . Its actual distance-to-default d d 3 with the underlying asset correlation into consideration is much smaller than d 3 . Hence, the third firm (AA rated bonds) is overrated if we have the correlation input.
We evaluate the time series of the probability of default for each individual firm under Assumptions 1 and 2 in Proposition 1 as well as the time series of the distance-to-default, and compare the probability of default and the distance-to-default from [3] under the isolated stochastic process for t = 0.5 , 1 , 1.5 , 2 , 2.5 , 3 and ρ 12 = 0.5 , ρ 13 = 0.25 and ρ 23 = 0.75 . For simplicity, we have μ i = λ i for i = 1 , 2 , 3 . The ratio V 1 , 0 K 1 = 11 for the BBB rated bonds with σ 1 = 0.4 , V 2 , 0 K 2 = 5 for the BB rated bonds with σ 2 = 0.3 and V 3 , 0 K 3 = 33 for the AA rated bonds with σ 3 = 0.2 .
Table 4 illustrates that three firms or bonds positively correlated have their probabilities of default and distances-to-default changing as t varies. The probability of default P 3 for the third AA rated bonds increases clearly, and P 1 stays stagnantly. The distances-to-default for both three assets decrease as t increases, as well as the difference d d i ( t ) - d i ( t ) gets smaller as t varies.
Compute F ( b 1 , b 2 , b 3 , t ) by Theorem 1 for n = 3 . We have ( s i g n ( z 1 ) , s i g n ( z 2 ) , s i g n ( z 3 ) ) = ( ± 1 , ± 1 , ± 1 ) Q i , and
Ω Q i Ω Q i T = σ 1 2 s i g n ( z 1 ) s i g n ( z 2 ) ρ 12 σ 1 σ 2 s i g n ( z 1 ) s i g n ( z 3 ) ρ 13 σ 1 σ 3 s i g n ( z 1 ) s i g n ( z 2 ) ρ 12 σ 1 σ 2 σ 2 2 s i g n ( z 2 ) s i g n ( z 3 ) ρ 23 σ 2 σ 3 s i g n ( z 1 ) s i g n ( z 3 ) ρ 13 σ 1 σ 3 s i g n ( z 2 ) s i g n ( z 3 ) ρ 23 σ 2 σ 3 σ 3 2
For instance with the quadrant ( 1 , - 1 , - 1 ) Q i , we have s i g n Q i = 2 and
Ω Q i Ω Q i T = σ 1 2 - ρ 12 σ 1 σ 2 - ρ 13 σ 1 σ 3 - ρ 12 σ 1 σ 2 σ 2 2 ρ 23 σ 2 σ 3 - ρ 13 σ 1 σ 3 ρ 23 σ 2 σ 3 σ 3 2
Hence ( - 1 ) s i g n Q i N Ω Q i Ω Q i T ( s i g n ( z 1 ) b 1 , s i g n ( z 2 ) b 2 , s i g n ( z 3 ) b 3 )
= ( - 1 ) 2 N Ω Q i Ω Q i T ( b 1 , - b 2 , - b 3 ) = N n , A i ( b 1 σ 1 t , - b 2 σ 2 t , - b 3 σ 3 t )
where A i is the matrix Ω Q i Ω Q i T with σ 1 = σ 2 = σ 3 = 1 and N n , A i is the standard multivariate normal distribution with A i . Hence, we have F ( b 1 , b 2 , b 3 , t ) given by the following.
N Ω Q 1 Ω Q 1 T ( b 1 , b 2 , b 3 ) - N Ω Q 2 Ω Q 2 T ( - b 1 , b 2 , b 3 ) + N Ω Q 3 Ω Q 3 T ( - b 1 , - b 2 , b 3 ) - N Ω Q 4 Ω Q 4 T ( b 1 , - b 2 , b 3 )
+ N Ω Q 5 Ω Q 5 T ( - b 1 , b 2 , - b 3 ) - N Ω Q 6 Ω Q 6 T ( - b 1 , - b 2 , - b 3 ) + N Ω Q 7 Ω Q 7 T ( b 1 , - b 2 , - b 3 ) - N Ω Q 8 Ω Q 8 T ( b 1 , b 2 , - b 3 )
= N n , A 1 ( b 1 σ 1 t , b 2 σ 2 t , b 3 σ 3 t ) - N n , A 2 ( - b 1 σ 1 t , - b 2 σ 2 t , b 3 σ 3 t ) + N n , A 3 ( - b 1 σ 1 t , - b 2 σ 2 t , b 3 σ 3 t )
- N n , A 4 ( b 1 σ 1 t , - b 2 σ 2 t , b 3 σ 3 t ) + N n , A 5 ( - b 1 σ 1 t , b 2 σ 2 t , - b 3 σ 3 t ) - N n , A 6 ( - b 1 σ 1 t , - b 2 σ 2 t , - b 3 σ 3 t )
+ N n , A 7 ( b 1 σ 1 t , - b 2 σ 2 t , - b 3 σ 3 t ) - N n , A 8 ( b 1 σ 1 t , b 2 σ 2 t , - b 3 σ 3 t )
where we give a complete formula for n = 3 that the probability of either three firms default in the above.
The following is the algorithm for the numerical evaluation in practice for n = 3 case.
  • Probability of Default with underlying correlations With the explicit numbers, we first compute P 1 ( t ) , P 2 ( t ) and P 3 ( t ) by Proposition 1 as in Table 2, and evaluate the distance-to-default d d 1 ( t ) , d d 2 ( t ) and d d 3 ( t ) for each three firms as in Table 3. (One can compare these values with p i ( t ) and d i ( t ) from previous evaluation to adjust their credit ratings with the consideration of the underlying correlations under Assumptions 1 and 2.)
  • Probabilities of either two firms default at the time t From the three firms with BBB, BB and AA rated by the S & P, we choose any pair (there are three choices) and evaluate the chosen two firms either one will be default at the time t by Proposition 2. Hence we obtain
    P ( D 1 ( t ) = 1 or D 2 ( t ) = 1 ) , P ( D 1 ( t ) = 1 or D 3 ( t ) = 1 ) , P ( D 2 ( t ) = 1 or D 3 ( t ) = 1 )
  • Probability of two firms both default at the time t By the previous two steps, we can evaluate both i-th firm and j-th firm default at the time t,
    P 12 ( t ) = P 1 ( t ) + P 2 ( t ) - P ( D 1 ( t ) = 1 or D 2 ( t ) = 1 )
    P 13 ( t ) = P 1 ( t ) + P 3 ( t ) - P ( D 1 ( t ) = 1 or D 3 ( t ) = 1 )
    P 23 ( t ) = P 2 ( t ) + P 3 ( t ) - P ( D 2 ( t ) = 1 or D 3 ( t ) = 1 )
  • Probability of both three firms default at the time t This is to analyze the systemic risk by the probability of both firms in the considered system to default at the time t. We first evaluate either firm from the three firms default by Theorem 1,
    P ( D 1 ( t ) = 1 or D 2 ( t ) = 1 or D 3 ( t ) = 1 )
    and hence the probability P 123 ( t ) is given by
    P 123 ( t ) = P ( D 1 ( t ) = 1 or D 2 ( t ) = 1 or D 3 ( t ) = 1 ) - P 1 ( t ) - P 2 ( t ) - P 3 ( t ) + P 12 ( t ) + P 13 ( t ) + P 23 ( t )
  • Default Correlations and Mixed defaults We can now get the default correlation matrix ( ρ i j D ( t ) ) by Proposition 2. Now we can formulate the mixed pair matrix in the following,
    P 1 ( t ) ρ 12 D ( t ) ρ 13 D ( t ) ρ 12 D ( t ) P 2 ( t ) ρ 23 D ( t ) ρ 13 D ( t ) ρ 23 D ( t ) P 3 ( t )
    By adding the rows of the mixed pair matrix, we obtain the mixed default measures
    M D 1 ( t ) = P 1 ( t ) + ρ 12 D ( t ) + ρ 13 D ( t )
    and M D 2 ( t ) , M D 3 ( t ) to access the other firms effected on the firm. The quantity ρ 12 D ( t ) + ρ 13 D ( t ) is the total effect of dependence from the specified first firm to the market of rest firms. The quantity ρ 12 D ( t ) + ρ 23 D ( t ) is the total effect of dependence from the specified second firm to the market of rest firms. The quantity ρ 13 D ( t ) + ρ 23 D ( t ) is the total effect of dependence from the specified third firm to the market of rest firms.

5. Conclusions

We give the answer to this issue by providing explicit formulas of the probability of multiple defaults and default correlations in this paper under the dependence. Merton [3] derives the probability of a single firm default under the isolated asset firm value hypothesis. By using a linear algebra technique, we change the driving Brownian motions into uncorrelated ones and apply the formula of Merton [3] to derive the probability of single firm default under the correlated asset values. In order to solve the probabilities of multiple defaults, we use the reflection principle in PDE to solve the Kolmogorov forward equation with initial-boundary conditions to evaluate the probability that any collection of some firms from the choice of all firms default in our main result Theorem 1. Then we propose an algorithm to determine probabilities of any multiple choice of firms defaults.
We further propose the mixed pair defaults to understand all the possible default correlations from the market to effect the specific firm. The sum of the i-th column of the default correlation matrix provide the quantity of the i-th firm to effect the market from the default correlation perspective. The mixed default pair would provide an important measure in credit risk management for fixed income portfolio, financial stability and insurance industry.
We also point out that the difference between the probability of single firm default under dependence in our result and the one used in structural models presently, as well as the difference between the joint default probability in this paper and the one used by KMV’s and J. P. Morgan’s. By fully considering the underlying asset correlations, the probability that multiple firms default can be measured precisely with parameters with respect to the time horizon. Our results can be applied in broad pricing and risk management related to credit risk.
The probability of multiple defaults and the default correlations are important factors in evaluating the credit risk of portfolios or correlated firms in credit markets. The financial industry and regulators recognize the importance of controlling default risks and default correlations for the stability of financial markets. A closed-form formula for evaluating default correlations in credit risk analysis is absolutely necessary for both practical and theoretical purposes. The importance and significance of the probability of multiple defaults and default correlations has been addressed extensively.
no

Acknowledgments

I would like to thank Tim Krehbiel for many helpful discussions related to this topic, and also thanks for Oliver Li and Yuecheng Jia for excellent research assistance.

Conflicts of Interest

The author declares no conflict of interests.

Appendix A

Proof of Proposition 1. 
We use the symmetric matrix property in Ω · Ω T to incorporate with Assumption 1 by standard Linear Algebra diagonalization. Under Assumptions 1 and 2, we first derive the default probability for each individual firm then, we use the first-passage-time to get the default probability and new distance-to-default corresponding to eigenvalues of the covariance matrix. Note that Ω · Ω T is a real symmetric matrix if and only if there is a spectral decomposition
Ω · Ω T = i = 1 n Σ i 2 q i q i T
where Ω · Ω T q i = Σ i 2 q i due to the covariance matrix’s positivity and { q 1 , , q n } is an orthnormal basis such that the orthnormal basis can form an orthogonal matrix Q such that Q T Ω · Ω T Q = Σ 1 2 0 0 0 0 0 0 Σ n 2 . Hence, the new covariance matrix can be treated as uncorrelated diffusions as in [3]. The eigenvalues of the covariance matrix Ω · Ω T = ( ρ i j σ i σ j ) 1 i , j n are given by
Σ 1 2 , Σ 2 2 , , Σ n 2
where Σ i = σ i if and only if ρ i j = 0 for all i j . If n = 2 , we have
Σ 1 2 = 1 2 σ 1 2 + σ 2 2 + ( σ 1 2 - σ 2 2 ) 2 + 4 ρ σ 1 σ 2 , Σ 2 2 = 1 2 σ 1 2 + σ 2 2 - ( σ 1 2 - σ 2 2 ) 2 + 4 ρ σ 1 σ 2
Then, the i-th firm specific probability of default P ( D i ( t ) = 1 ) = P ( τ i t ) is given by
P i ( t ) = P ( D i ( t ) = 1 ) = N ( - ln ( V i , 0 / K i ) Σ i t - μ i - λ i Σ i t )
+ ( V i , 0 K i ) 2 ( λ i - μ i ) σ i 2 N ( - ln ( V i , 0 / K i ) Σ i t + μ i - λ i Σ i t )
The underlying asset values of these firms are now uncorrelated through the new driving Brownian motions with respect to the orthnormal basis. Hence, one can adapt Merton’s structural model for the isolated firm asset value to find the probability of default in (A2) for the i-th firm. ☐
Proof of Theorem 1. 
The evaluation of the probability that either firm defaults with default correlation is reduced to the following probability.
P ( D 1 ( t ) = 1 or D 2 ( t ) = 1 , or D n ( t ) = 1 )
= P ( τ 1 t or τ 2 t , or τ n t ) = P ( τ t )
where τ = min { τ 1 , τ 2 , , τ n } and τ i = min { t 0 : V i ( t ) e λ i t K i } for i = 1 , 2 , , n .
If λ i = μ i , then we can set X i ( t ) = - ln [ e - λ i t V i ( t ) / V i ( 0 ) ] and b i = - ln [ K i / V i ( 0 ) ] . The default condition V i ( t ) e λ i t K i is equivalent to X i ( t ) b i .
If λ i μ i , then set X i ( t ) = - ( λ i - μ i ) t - ln [ e - λ i t V i ( t ) / V i ( 0 ) ] and b i ( t ) = - ln [ e λ i t K i e μ i t V i ( 0 ) ] , the default is equivalent to X i ( t ) b i ( t ) .
For both cases, we have d X ( t ) = - Ω d W ( t ) with X ( t ) T = ( X 1 ( t ) , , X n ( t ) ) and n-dimensional independent Brownian motion W ( t ) T = ( W 1 ( t ) , W n ( t ) ) .
Let f ( x 1 , x 2 , , x n , t ) be the transition probability density of a particle in the region { ( x 1 , x 2 , , x n ) R n : x i < b i 1 i n before the time t. Therefore, we have
F ( b 1 , b 2 , , b n , t ) = - b 1 - b 2 - b n f ( x 1 , x 2 , , x n , t ) d x 1 d x 2 d x n = P ( τ > t ) = 1 - P ( τ t ) = 1 - P ( D 1 ( t ) = 1 or D 2 ( t ) = 1 , or D n ( t ) = 1 )
Following the standard Kolmogorov forward equation for the transition probability density f ( x 1 , x 2 , , x n , t ) , we have the classical problem to determine the density from the PDE, i.e., the transition probability density is the solution of the Kolmogorov forward equation
1 2 ( i = 1 n σ i 2 2 f x i 2 + 2 1 i j n ρ i j σ i σ j 2 f x i x j ) = f t , x i < b i , 1 i n
subject to the following boundary conditions:
f ( - , x 2 , , x n , t ) = f ( x 1 , , - , , x n , t ) = f ( x 1 , , x n - 1 , - , t ) = 0 f ( x 1 , x 2 , , x n , 0 ) = δ ( x 1 ) δ ( x 2 ) δ ( x n ) F ( b 1 , b 2 , , b n , t ) 1 , t > 0 f ( b 1 , x 2 , , x n , t ) = = f ( x 1 , , x n - 1 , b n , t ) = 0
where δ ( x ) is the Dirac’s delta function with - δ ( x ) d x = 1 . The equation has a solution
f ( x 1 , x 2 , , x n , t ) = 1 ( 2 π t ) n / 2 det ( Ω Ω T ) 1 / 2 exp ( - 1 2 t x T ( Ω Ω T ) - 1 x )
as a time-varying Gaussian distribution with the symmetric matrix Ω Ω T , where x T = ( x 1 , , x n ) . It is straightforward to verify that the multivariate normal distribution function f ( x 1 , x 2 , , x n , t ) satisfies the Kolmogorov forward equation and the first three boundary conditions, except f ( b 1 , x 2 , , x n , t ) = = f ( x 1 , , x n - 1 , b n , t ) = 0 .
For completeness, we present the complete proof of the simple reflection principle of the partial differential equation for this problem. This settles the correlated default probability from the first-passage-time approach.
Let P u = ( t - 1 2 ( i = 1 n σ i 2 2 x i 2 + 2 1 i j n ρ i j σ i σ j 2 x i x j ) ) u = 0 be the partial differential equation with u ( x , 0 ) = 0 , u ( - , x 2 , t ) = = f ( x 1 , , - , t ) = 0 , u ( b 1 , x 2 , , x n , t ) = u ( x 1 , , x n - 1 , b n , t ) = 0 . Let X i = - x i + b i ( 1 i n ) be a coordinate shift, and g ( X 1 , X 2 , , X n , t ) = u ( x 1 , x 2 , , x n , t ) be the function in terms of the new variables X i for 1 i n . Thus, we have
P g ( X 1 , X 2 , , X n , t ) = 0 , 0 X i < , 1 i n
g ( X 1 , X 2 , , X n , 0 ) = δ ( X 1 ) δ ( X 2 ) δ ( X n )
g ( + , X 2 , t ) = = g ( X 1 , , + , t ) = 0
g ( X 1 , , 0 | X i , , X n , t ) = 0 , 0 X i < , 1 i n
Now, we use the reflection principle to solve the above boundary valued PDE. Define h ( X 1 , , X n , t ) be the odd extension of g ( X 1 , , X n , t ) as the following:
h ( y 1 , y 2 , , y n , t ) = ( - 1 ) i = 1 n s i g n - y i g ( s i g n ( y 1 ) y 1 , , s i g n ( y n ) y n , t )
where s i g n - ( y i ) = - 1 if s i g n ( y i ) = - 1 and s i g n - ( y i ) = 0 if s i g n ( y i ) = 1 , and the odd extension of g runs over the 2 n quadrants in the Euclidean plane R n with respect to the quadrant where ( s i g n ( y 1 ) , , s i g n ( y n ) ) belongs to. For example when n = 2 , the odd extension of g in first quadrant Q 1 can be expressed by
h ( y 1 , y 2 , t ) = g ( y 1 , y 2 , t ) if   ( y 1 , y 2 ) Q 1 - g ( - y 1 , y 2 , t ) if   ( y 1 , y 2 ) Q 2 g ( - y 1 , - y 2 , t ) if   ( y 1 , y 2 ) Q 3 - g ( y 1 , - y 2 , t ) if   ( y 1 , y 2 ) Q 4
where Q i is the i-th quadrant in the Euclidean plane R 2 . Similarly, we define an odd extension for the initial value for δ ( y 1 , y 2 , , y n ) accordingly:
δ ( y 1 , y 2 , , y n ) = ( - 1 ) i = 1 n s i g n - y i δ ( s i g n ( y 1 ) y 1 ) δ ( s i g n ( y n ) y n )
For example when n = 2 ,
δ ( y 1 , y 2 ) = δ ( y 1 ) δ ( y 2 ) if   ( y 1 , y 2 ) Q 1 - δ ( - y 1 ) δ ( y 2 ) if   ( y 1 , y 2 ) Q 2 δ ( - y 1 ) δ ( - y 2 ) if   ( y 1 , y 2 ) Q 3 - δ ( y 1 ) δ ( - y 2 ) if   ( y 1 , y 2 ) Q 4
Then, we solve the initial value problem on the whole plane R n :
P h ( y 1 , y 2 , , y n , t ) = 0 , t > 0 , ( y 1 , , y n ) R n ; h ( y 1 , , y n , 0 ) = δ ( y 1 , y 2 , y 2 , , y n )
The problem can be solved through the standard heat kernel on R 2 ,
h ( y 1 , , y n , t ) = R n f ( y 1 - z 1 , y 2 - z 2 , , y n - z n , t ) δ ( z 1 , z 2 , , z n ) d z 1 d z 2 d z n
where f ( z 1 , z 2 , , z n , t ) = 1 ( 2 π t ) n / 2 det ( Ω Ω T ) 1 / 2 exp ( - 1 2 t z T ( Ω Ω T ) - 1 z ) . Due to the odd extension, h ( y 1 , , y n , t ) = 0 whenever y i = 0 for some i. Therefore g ( y 1 , y 2 , , y n , t ) is simply the restriction of h ( y 1 , y 2 , , y n , t ) on the first quadrant Q 1 . For ( y 1 , y 2 , , y n ) Q 1 ,
g ( y 1 , y 2 , , y n , t ) = h ( y 1 , y 2 , , y n , t ) = R n f ( y 1 - z 1 , y 2 - z 2 , , y n - z n , t ) δ ( z 1 , z 2 , , z n ) d z 1 d z 2 d z n
By using the definition of δ ( z 1 , z 2 , , z n ) on different regions to decompose the integral into 2 n terms corresponding to each quandrant, we obtain g ( y 1 , y 2 , , y n , t )
= R n f ( y 1 - z 1 , y 2 - z 2 , , y n - z n , t ) δ ( z 1 , z 2 , , z n ) d z 1 d z 2 d z n = f ( y 1 - z 1 , y 2 - z 2 , , y n - z n , t ) ( - 1 ) i = 1 n s i g n - z i δ ( s i g n ( z 1 ) z 1 ) δ ( s i g n ( z n ) z n ) d z 1 d z 2 d z n = Q 1 ( - 1 ) i = 1 n s i g n - z i f ( y 1 - s i g n ( z 1 ) z 1 , y 2 - s i g n ( z 2 ) z 2 , , y n - s i g n ( z n ) z n , t ) δ ( z 1 ) δ ( z 2 ) δ ( z n ) d z 1 d z 2 d z n
where the second identity follows from the odd extension of δ ( z 1 , z 2 , , z n ) from the first quadrant to all the plane R n , the third identity from the change of variables to switch all different quadrants into the first quadrant Q 1 . Therefore the solution for the boundary conditions is given by u ( x 1 , x 2 , , x n , t )
= g ( X 1 , X 2 , , X n , t ) = g ( - x 1 + b 1 , - x 2 + b 2 , , - x n + b n , t ) = Q 1 ( - 1 ) i = 1 n s i g n - z i f ( - x 1 + b 1 - s i g n ( z 1 ) z 1 , , - x n + b n - s i g n ( z n ) z n , t ) δ ( b 1 - z 1 ) δ ( b n - z n ) d z 1 d z n = i = 1 2 n ( - 1 ) s i g n Q i - b 1 - b n f ( - x i + s i g n ( z i ) y i + b i ( 1 - s i g n ( z i ) ) , t ) δ ( y 1 ) δ ( y n ) d y 1 d y n = i = 1 2 n ( - 1 ) s i g n Q i - b 1 - b n f Ω Q i ( s i g n ( z i ) x i + b i ( 1 - s i g n ( z i ) ) - y i , t ) δ ( y 1 ) δ ( y n ) d y 1 d y n = i = 1 2 n ( - 1 ) s i g n Q i f Ω Q i ( s i g n ( z 1 ) x 1 + b 1 ( 1 - s i g n ( z 1 ) ) , , s i g n ( z n ) x n + b n ( 1 - s i g n ( z n ) ) , t )
where the first two equalities follow from the coordinate shift, the third from the previous identity with the substitution y i = - x i + b i ( 1 i n ), sign Q i is the number of negative 1 for the representative ( ± 1 , ± 1 , , ± 1 ) in Q i , the fourth identity from the variable change y i = b i - z i ( 1 i n ) and ( s i g n ( z 1 ) , , s i g n ( z n ) is the representative of Q i , the fifth identity from the even function property of f, where f Ω Q i is the function f ( s i g n ( z 1 ) z 1 , , s i g n ( z n ) z n , t ) which is same as replacing Ω by
diag ( s i g n ( z 1 ) , , s i g n ( z n ) ) Ω = Ω Q i
and the last from the Dirac delta function property.
Thus we have F ( b 1 , b 2 , , b n , t )
= - b n - b 1 u ( x 1 , x 2 , , x n , t ) d x 1 d x n = i = 1 2 n ( - 1 ) s i g n Q i - b n - b 1 f Ω Q i ( s i g n ( z 1 ) x 1 + b 1 ( 1 - s i g n ( z 1 ) ) , , s i g n ( z n ) x n + b n ( 1 - s i g n ( z n ) ) ) d x 1 d x n = i = 1 2 n ( - 1 ) s i g n Q i N Ω Q i Ω Q i T ( s i g n ( z 1 ) b 1 , s i g n ( z 2 ) b 2 , , s i g n ( z n ) b n )
where the last equality follows from various changing variables and N A ( x 1 , , x n ) = - x n - x 1 1 ( 2 π t ) n / 2 det A exp ( - 1 2 t z T A z ) d z 1 d z n .
For example n = 2 , the representative ( s i g n ( z 1 ) , s i g n ( z 2 ) ) = ( ± 1 , ± 1 ) for the quadrants with s i g n ( z 1 ) = ± 1 and s i g n ( z 2 ) = ± 1 has
( 1 , 1 ) Q 1 , ( - 1 , 1 ) Q 2 , ( - 1 , - 1 ) Q 3 , ( 1 , - 1 ) Q 4
s i g n Q 1 = 0 , s i g n Q 2 = 1 , s i g n Q 3 = 2 , s i g n Q 4 = 1
Ω Q 1 = Ω , Ω Q 2 = d i a g ( - 1 , 1 ) Ω , Ω Q 3 = d i a g ( - 1 , - 1 ) Ω , Ω Q 4 = d i a g ( 1 , - 1 ) Ω
Hence F ( b 1 , b 2 , t ) = i = 1 2 2 ( - 1 ) s i g n Q i N Ω Q i ( s i g n ( z 1 ) b 1 , s i g n ( z 2 ) b 2 ) . Thus, we have
F ( b 1 , b 2 , t ) = N 2 , ρ ( b 1 σ 1 t , b 2 σ 2 t ) - N 2 , - ρ ( - b 1 σ 1 t , b 2 σ 2 t ) + N 2 , ρ ( - b 1 σ 1 t , - b 2 σ 2 t ) - N 2 , - ρ ( b 1 σ 1 t , - b 2 σ 2 t )
where N 2 , ρ ( x 1 , x 2 ) = - x 2 - x 1 1 2 π 1 - ρ 2 e - y 1 2 - 2 ρ y 1 y 2 + y 2 2 2 ( 1 - ρ 2 ) d y 1 d y 2 is the probability of the standard bivariate normal distribution on y 1 x 1 and y 2 x 2 , and 1 2 π 1 - ρ 2 e - y 1 2 - 2 ρ y 1 y 2 + y 2 2 2 ( 1 - ρ 2 ) is the standard bivariate normal distribution density function N 2 ( 0 , 1 , ρ ) for random variable Y i with mean zero and variance 1 and correlation ρ = C o v ( Y 1 , Y 2 ) . ☐

References

  1. D. Duffie, and N. Gârleanu. “Risk and valuation of collateralized debt obligation.” Financ. Anal. J. 57 (2001): 41–59. [Google Scholar] [CrossRef]
  2. F. Black, and M. Scholes. “The pricing of options and corporate liabilities.” J. Polit. Econ. 81 (1973): 637–659. [Google Scholar]
  3. R.C. Merton. “On the pricing of corporate debt: The risk structure of interest rates.” J. Financ. 29 (1974): 449–470. [Google Scholar] [CrossRef]
  4. F. Black, and J.C. Cox. “Valuing Corporate Securities: Some effects of bond indenture provisions.” J. Financ. 31 (1976): 351–367. [Google Scholar] [CrossRef]
  5. D.J. Lucas. “Default correlation and credit analysis.” J. Fixed Income 4 (1995): 76–87. [Google Scholar] [CrossRef]
  6. S.R. Das, G. Fong, and G. Geng. “Impact of correlated default risk on credit portfolios.” J. Fixed Income 11 (2001): 9–19. [Google Scholar] [CrossRef]
  7. P. Jorion. “Risk management lessons from the credit crisis.” Eur. Financ. Manag. 15 (2009): 923–933. [Google Scholar] [CrossRef]
  8. M. Crouhy, D. Galai, and R. Mark. “A comparative analysis of current redit risk models.” J. Bank. Financ. 24 (2000): 59–117. [Google Scholar] [CrossRef]
  9. M. Gundlach, and F. Lehrbass, eds. CreditRisk+ in the Banking Industry. Berlin/Heidelberg, Germany: Springer, 2004.
  10. K. Giesecke. “Correlated default with incomplete information.” J. Bank. Financ. 28 (2004): 1521–1545. [Google Scholar] [CrossRef]
  11. D.X. Li. “On default correlation: A copula function approach.” J. Fixed Income 9 (2000): 43–54. [Google Scholar] [CrossRef]
  12. J.A. Rebholz. “Planar Diffusions with Applications to Mathematical Finance.” Ph.D. Thesis, University of California, Berkeley, CA, USA, 1994. [Google Scholar]
  13. C. Zhou. “An analysis of default correlations and multiple defaults.” Rev. Financ. Stud. 14 (2001): 555–576. [Google Scholar] [CrossRef]
  14. M. Valužis. “On the probabilities of correlated defaults: A first passage time approach.” Nonlinear Anal. Model. Control 13 (2008): 117–133. [Google Scholar]
  15. A. Metzler. “On the first passage problem for correlated Brownian motion.” Stat. Probab. Lett. 80 (2010): 277–284. [Google Scholar] [CrossRef]
  16. W. Li, and T. Krehbiel. “An improved approach of default probabilities and default correlations with consistency.” Int. J. Theor. Appl. Financ. 19 (2016): 1650036. [Google Scholar]
  17. D. Zhang, and R.V. Melnik. “First passage time for multivariate jump-diffusion processes in finance and other areas of applications.” Appl. Stoch. Models Bus. Ind. 25 (2009): 565–582. [Google Scholar] [CrossRef]
  18. L. Sacerdote, M. Tamborrino, and C. Zucca. “First passage times of two-dimensional correlated diffusion processes: Analytical and numerical methods.” , 2012. [Google Scholar]
  19. U. Erlenmaier, and H. Gersbach. Default Probabilities and Default Correlations. Frankfurt am Main, Germany: Deutsche Bank Research, 2001. [Google Scholar]
  20. U. Erlenmaier, and H. Gersbach. “Default correlations in the merton model.” Rev. Financ. 18 (2014): 1775–1809. [Google Scholar] [CrossRef]
  21. R. Jarrow. “Risk management models: Construction, testing, usage.” J. Deriv. 18 (2011): 89–98. [Google Scholar] [CrossRef]
  22. T.N. Nielsen, J. Saá-Requejo, and P. Santa-Clara. “Default risk and interest rate risk, the term structure of default spreads.” Unpublished work. 1993. [Google Scholar]
  23. H. Leland. “Corporate debt value, bond covnants, and optimal capital structure.” J. Financ. 49 (1994): 1213–1252. [Google Scholar] [CrossRef]
  24. H. Leland, and K. Toft. “Optimal capital structure, endogenous bankruptcy, and the term structure of credit spreads.” J. Financ. 51 (1996): 987–1019. [Google Scholar] [CrossRef]
  25. T.R. Bielecki, and M. Rutkowski. Dredit Risk: Modelling, Valuation and Hedging. Berlin/Heidelberg, Germany: New York, NY, USA: Springer-Verlag, 2002. [Google Scholar]
  26. R. Jarrow. “Credit market equilibrium theory and evidence: Revisiting the structural vs. reduced form credit risk model debate.” Financ. Res. Lett. 8 (2011): 2–7. [Google Scholar] [CrossRef]
  27. N. Akkaya, A. Kurth, and A. Wagner. “Incorporating Default correlations and severity variations.” In CreditRisk+ in the Banking Industry. Edited by M. Gundlach and F. Lehrbass. Heidelberg/Berlin, Germany: Springer, 2004, pp. 129–152. [Google Scholar]
Table 1. Implied values Z from 2010 Cumulative Default Study.
Table 1. Implied values Z from 2010 Cumulative Default Study.
AAAAAABBBBBBCINVSPECALL
Z8.9168.7697.5695.9803.9022.5151.1387.0242.7824.457
V / K 353321115321736
Table 2. Eigenvalues and Probability of Default with underlying asset correlations.
Table 2. Eigenvalues and Probability of Default with underlying asset correlations.
ρ 12 ρ 13 ρ 23 Σ 1 Σ 2 Σ 3 P 1 p 1 P 2 p 2 P 3 p 3
−0.75−0.75−0.750.06250.06256.258.67E-222.04E-090.000%5.73E-0516.19%2.87E-05
−0.75−0.75−0.50.06581.91873.57818.93E-212.04E-0924.527%5.73E-056.45%2.87E-05
−0.75−0.75−0.250.00081.08564.28860.00E+002.04E-0912.242%5.73E-059.13%2.87E-05
−0.75−0.7500.003714.24630.00E+002.04E-0910.752%5.73E-058.97%2.87E-05
−0.75−0.750.250.00320.56254.80930.00E+002.04E-093.188%5.73E-0511.08%2.87E-05
−0.75−0.750.50.16030.52.33972.11E-092.04E-092.284%5.73E-052.23%2.87E-05
−0.75−0.750.750.06250.06256.258.67E-222.04E-090.000%5.73E-0516.19%2.87E-05
00.250.750.043913.20612.50E-302.04E-0910.752%5.73E-055.08%2.87E-05
00.50.750.00972213.6153221.22E-1302.04E-0910.752%5.73E-056.59%2.87E-05
00.750.750.00368414.2464840.00E+002.04E-0910.752%5.73E-058.97%2.87E-05
0−0.250.750.043913.20612.50E-302.04E-0910.752%5.73E-055.08%2.87E-05
0−0.50.750.00972213.6153221.22E-1302.04E-0910.752%5.73E-056.59%2.87E-05
0−0.750.750.00368414.2464840.00E+002.04E-0910.752%5.73E-058.97%2.87E-05
0.50.250.750.03980.5962934.3388892.81E-332.04E-093.714%5.73E-059.32%2.87E-05
0.50.50.750.06250.3301654.7323658.67E-222.04E-090.509%5.73E-0510.80%2.87E-05
0.50.750.750.0256960.255.4741961.36E-502.04E-090.129%5.73E-0513.51%2.87E-05
0.50.250.50.1655680.56253.3970183.79E-092.04E-093.188%5.73E-055.78%2.87E-05
0.50.50.50.250.2541.62E-062.04E-090.129%5.73E-058.04%2.87E-05
0.50.750.50.06250.3301654.7323658.67E-222.04E-090.509%5.73E-0510.80%2.87E-05
Table 3. Comparing distances-to-default.
Table 3. Comparing distances-to-default.
ρ 12 ρ 13 ρ 23 d d 1 d 1 d d 2 d 2 d d 3 d 3
−0.75−0.75−0.759.4675.7956.3135.2150.14917.383
−0.75−0.75−0.59.2205.7950.4695.2150.90317.383
−0.75−0.75−0.2584.7645.7951.0245.2150.65317.383
−0.75−0.75039.3915.7951.1095.2150.66617.383
−0.75−0.750.2542.3615.7951.7715.2150.49817.383
−0.75−0.750.55.7895.7951.9235.2151.52117.383
−0.75−0.750.759.4675.7956.3135.2150.14917.383
00.250.7511.3405.7951.1095.2151.05717.383
00.50.7524.2705.7951.1095.2150.88817.383
00.750.7539.4745.7951.1095.2150.66617.383
0−0.250.7511.3405.7951.1095.2151.05717.383
0−0.50.7524.2705.7951.1095.2150.88817.383
0−0.750.7539.4745.7951.1095.2150.66617.383
0.50.250.7511.9205.7951.6985.2150.63717.383
0.50.50.759.4675.7952.5145.2150.52017.383
0.50.750.7514.8795.7952.9695.2150.32517.383
0.50.250.55.6905.7951.7715.2150.97617.383
0.50.50.54.5465.7952.9695.2150.74817.383
0.50.750.59.4675.7952.5145.2150.52017.383
Table 4. Time-Series Default Probability and Distance-to-Default.
Table 4. Time-Series Default Probability and Distance-to-Default.
t P 1 p 1 P 2 p 2 P 3 p 3 d d 1 d 1 d d 2 d 2 d d 3 d 3
0.501.15E-170.0001356.34E-090.2544332.1E-3516.9288.33610.3525.5491.63712.221
101.02E-090.0069532.87E-050.4203281.15E-1812.0205.7956.5733.8240.6378.541
1.504.92E-070.027540.0005090.5105534.76E-139.8144.6504.7573.0400.0956.892
201.12E-050.0563220.002220.5687973.18E-108.4993.9563.5912.562−0.2865.898
2.507.49E-050.0878140.0054680.6102851.62E-087.6023.4752.7392.229−0.5855.212
34.2E-2650.0002690.119160.0100890.6417462.25E-076.9393.1152.0691.977−0.8354.700

Share and Cite

MDPI and ACS Style

Li, W. Probability of Default and Default Correlations. J. Risk Financial Manag. 2016, 9, 7. https://doi.org/10.3390/jrfm9030007

AMA Style

Li W. Probability of Default and Default Correlations. Journal of Risk and Financial Management. 2016; 9(3):7. https://doi.org/10.3390/jrfm9030007

Chicago/Turabian Style

Li, Weiping. 2016. "Probability of Default and Default Correlations" Journal of Risk and Financial Management 9, no. 3: 7. https://doi.org/10.3390/jrfm9030007

APA Style

Li, W. (2016). Probability of Default and Default Correlations. Journal of Risk and Financial Management, 9(3), 7. https://doi.org/10.3390/jrfm9030007

Article Metrics

Back to TopTop