Abstract
One of the key questions in credit dependence modelling is the specfication of the copula function linking the marginals of default variables. Copulae functions are important because they allow to decouple statistical inference into two parts: inference of the marginals and inference of the dependence. This is particularly important in the area of credit risk where information on dependence is scant. Whereas the techniques to estimate the parameters of the copula function seem to be fairly well established, the choice of the copula function is still an open problem. We find out by simulation that the tcopula naturally arises from a structural model of credit risk, proposed by Cossin and Schellhorn (2007). If revenues are linked by a Gaussian copula, we demonstrate that the tcopula provides a better fit to simulations than does a Gaussian copula. This is done under various specfications of the marginals and various configurations of the network. Beyond its quantitative importance, this result is qualitatively intriguing. Student's tcopulae induce fatter (joint) tails than Gaussian copulae ceteris paribus. On the other hand observed credit spreads have generally fatter joint tails than the ones implied by the Gaussian distribution. We thus provide a new statistical explanation why (i) credit spreads have fat joint tails, and (ii) financial crises are amplified by network effects.
1. Introduction
One of the key questions in credit dependence modelling is currently the specification of the copula function linking the marginals of default variables. Several books have been written on copulae as well as their application to finance, for example, Cherubini et al. [1], Embrechts et al. [2], Joe [3], and Nelsen [4]. We refer to these books for an exposition of copula theory. The main application of copulae seems to be the following. In several domains, and in particular in credit risk, there is often a lack of information to model the dependence between random variables. However, there is sufficient information to model the marginals of each random variable. It is in this particular context thatcopulae functions are important for dependence modelling. The estimation of the joint distribution can thus be decoupled into estimation of the marginals, which is more robust, and estimation of the copula function. This estimation method, coupled with maximum likelihood (MLE) estimation, is often called inference for the margins, or IFM (see, e.g., Cherubini et al. [1, Section ]).
Whereas the techniques to estimate the parameters of the copula function seem to be fairly well established, the choice of the copula function is still an open problem. Among several copulae, researchers in finance seem to have employed mostly the Gaussian, , and Archimedean (family of) copulae. One of the advantages of the Gaussian copula is ease of simulation, while the Archimedean copula offers some advantages for inference. There has been widespread use of Gaussian copulas by rating agencies lately.
The initial goal of this research was to help practitioners in the selection of a particular copula function. More specifically, our main goal was to compare tcopulae with the Gaussian copula as a way to model counterparty risk. We found out by simulation that the tcopula naturally arises from a structural model of credit risk, proposed by Cossin and Schellhorn [5], and henceforth abbreviated CS model. The CS model links operating revenue to the credit spreads of firms in a network economy. If revenues are linked by a Gaussian copula, we demonstrate that the tcopula provides a better fit to simulations than does a Gaussian copula. This is done under various specifications of the marginals and various configurations of the network. Beyond its quantitative importance, this result is qualitatively intriguing. It has been recognized by various researchers (see, e.g., Bluhm [6]), that tcopulae induce fatter (joint)^{1} tails than Gaussian copulae ceteris paribus. On the other hand, the finance literature has abundantly documented in the last twenty years that credit spreads have generally fatter joint tails than implied by the Gaussian distribution. We thus provide a new statistical explanation why credit spreads have fat joint tails: even if the driver of credit spreads (namely revenue) has normal tails, because of counterparty relationships, credit spreads have fatter joint tails. Thus, financial crises are amplified by network effects.
This work also confirms the plausibility of the CS model. Indeed, under the hypothesis that (i) the (CS) model is valid and (ii) revenue has normal tails, we conclude that credit spreads have fat (joint) tails. Since fat joint tails are empirically observed, this seems to indicate that the (CS) model does a good job at explaining contagion. In this article, we do not dwell on this positive feature of our model. Indeed most credit risk models attempt to explain nonnormal risk of contagion, and a model which does not exhibit this feature would not be very interesting.
The curious reader may wonder what is the point of doing a statistical inference on a deterministic model. It is tempting to allude to Einstein's criticism of quantum mechanics “God does not play with dice”. While there may be a purely deterministic description of microphysical phenomena beyond quantum mechanics, the latter theory proved to be fruitful for many years. We could argue that the same applies in credit risk. While a better description of credit risk if probably is given by models like the (CS) model, the lack of identifiability of model parameters (due to lack of disclosure of exposure parameters ) forces practitioners to resort to statistical models. The point is that the statistical model must be indistinguishable from benchmark deterministic models.
Our article is composed of three parts. We first expose the models of network economies that we test. We then describe the statistical methodology. The statistical methodology is split into two parts: parameter estimation and hypothesis testing. In the hypothesis testing section, we compare the two nullhypotheses: Gaussian copula, and tcopula. Finally, we present the results.
2. Models of Network Economies
We first describe briefly the (CS) model in general and then introduce the types of network economies that we will study. We then discuss the distributional assumptions.
2.1. The CS Model
In the CS model, each firm is (potentially) at the same time a borrower from another (or several) firm(s) in the network and a lender to another (or several) firm(s). Our main assumption is not only that the quantity of debt is fixed, but also the network of lending and borrowing is fixed. In other terms, firms have preferred lenders, and the amount borrowed from them does not change with time. In addition to their borrowing and lending function, each firm produces goods and distributes dividends to its shareholders, in a manner similar to Leland's [7] model. We show an example of such a network economy in Figure 1.
We call the long run operational revenue rate of firm . In the CS model, this rate is not directly observable. At each time the market calculates an estimator of . For the sake of brevity in this document, we will call simply the production revenue. The total revenue of firm consists of its operational revenue plus debt payments from its borrowers. The main assumption of our model is that network relationships are fixed: debt payments from firm to firm are, at all times, proportional to the total revenue of firm . In other terms, there are constants so that
The payout ratio expresses what percentage of total expenses of firm is distributed to bondholders and equityholders. It is modelled as a geometric Brownian motion with relative drift and volatility . A payout ratio superior to one corresponds to a recapitalization of the firm. We assume that the payout ratio is independent of the level of expenses. Payout ratios across firms can be correlated. Upon default, firm incurs a loss rate of . Theorem 2.1 was proved in Cossin and Schellhorn [5]:
Theorem 2.1. In steady state, the values of equity and debt are, for finite , and not a bankruptcy time:
The parameter scales the production revenue intensity, as is common in diffusion approximations. We refer to Cossin and Schellhorn [5] for a full definition of . Roughly speaking, a large corresponds to a large cash account, which acts as a (partial) buffer against cash flow risk. In this paper, we make the assumption that is very large, for simplicity. The functions and are fully analytical and equal (mutatis mutandis) to the formulae for the price of debt and equity in Leland's [7] model.
We want to show that the copula function of the counterparty risk premium (which we define in (2.8)) is equal to the copula function of total revenue. This occurs in the case of a stationary debt structure a la Leland and Toft [8], where the distribution of principal is uniform between . In such an environment, a constant bankruptcy level is optimal. It can also be seen that liquidity risk is also decoupled from counterparty risk. More formally, let be the total value of debt. There is now a new function so that Differentiating (2.3) with respect to , we see that the value of debt with maturity is equal to On the other hand, we can rewrite the lefthandside as where is the total principal and the credit spread of firm for maturity . For simplicity we write Thus, equating (2.4) with (2.5), and taking we have Our main goal in this article is to assess the impact of counterparty risk. We are thus more interested in the credit risk coming from than the credit risk coming from in the last expression. We thus define the counterparty risk premium and a bankruptcy risk premium :
The counterparty risk premium stems from the network term and the bankruptcy risk premium comes from the idiosyncratic variable , which represents the payout ratio. While the bankruptcy risk premium is different for debt and for equity, the counterparty risk premium is the same for both securities. Assembling (2.8), (2.9), and (2.7), we see that we can decompose the credit spread into counterparty risk premium and bankruptcy risk premium:
It is well known that the copula function is unchanged under monotone transformations, such as the transformation (2.8) between total revenue and counterparty risk premium . We thus showed the following fact.
Fact 2.1. The copula function of the counterparty risk premium is the same as the copula function of total revenue.
Finally, we perform only a static analysis; that is, we analyze only the counterparty risk premium viewed as a random variable. What remains to be specified are then the network configurations and the distributional assumptions.
2.2. Networks
We analyze 3 different network configurations, that is possible relationships (networks) between the production revenue and the total revenue :(i)a “triangular”, or “full” network,(ii)a “star” network,(iii)a “series” network.
2.2.1. Triangular or “Full” Network
In this type of network, firm makes loans to firms . Firm loans to firms , and so on. The network equations are then (see also Figure 2)
We analyze the specific case for . Also, all the firms have the same size, that is,
2.2.2. Series Network
In a series network, each firm borrows from the “next” firm . It is illustrated in Figure 3: where .
2.2.3. Star Network
A star network is an idealization of an economy with a small group of large “middle firms” who both lend to a set of “start firms” and borrow from a set of “end firms.” Start and end firms all have the same revenue size, that is, . All large firms also have the same revenue size, .
We choose a particular network with 7 start firms, 2 middle firms, and 3 end firms. A typical example is the US car industry, with Ford and GM as “middle firms”, dealers as end firms, and suppliers as start firms. Figure 4 illustrates this type of network. Since only the production revenue of a (say) start firm times lambda matters in determining the impact of a start firm on a middle firm, it is needless to vary lambda in this network to account for different relative impact of start firms on middle firms. A different relative revenue suffices.
2.3. Distributional Assumptions
Our main assumption is that the dependence between production revenues is described by the Gaussian copula. Since copulae are invariant under monotone transformations of variables, the copula of the logarithm of revenue is also Gaussian under this hypothesis. We chose the following variancecovariance matrix for :
with . This choice was fairly arbitrary, but we found out in other tests (not reported here) that other values of give us similar results. We test different types of marginal distributions for production revenue : (i)exponential with various rates, namely, 1, 5, and 20(ii)uniform with various support,(iii)gamma with various shapes and scales,(iv)quadratic where ,(v)cubic where .
Exponential, gamma, and uniform random variables are standard choices for modelling positive random variables (see, e.g., Lando [9]).
3. Statistical Methodology
As stated earlier, the main goal in this paper was to compare tcopula with a Gaussian copula as a way to model counterparty risk. Since the Gaussian copula is a special case of a tcopula, namely a tcopula with an infinite number of degrees of freedom, we try to fit a tcopula to our simulated data. The calibrated number of degrees of freedom will be a good indicator whether a nonGaussian tcopula is a better choice than the Gaussian copula. We then expose our methodology for hypothesis testing.
3.1. Parameter Estimation
Assuming that the dependence of the firms' counterparty risk premia is the Student's tcopula, we conduct our estimation analyses based on the IFM method discussed in Chapter 5.3 of the book Copula Methods in Finance by Cherubini et al. [1]. Cherubini et al's method is composed of two following steps: (1)infer the parameter(s) of the marginals,(2)infer the parameter(s) of the Student's tcopula.
Our analyses focus on the second step: estimating the parameter(s) of the copula . We bypass the first step because we take the empirical distribution of the marginals as given. The Student's tcopula has two parameters, namely, the correlation matrix and the degrees of freedom . For simplicity, we use the method of moments to first infer the correlation matrix, and use maximum likelihood to estimate the degrees of freedom : where
Let be the total number of firms in the network and the total number of revenue scenarios. The following algorithm constructs the network economy and estimates the parameter vector .(1)Generate a table of total revenues based on the relationship specified in each network discussed in the previous section: (2)Calculate, for each scenario , (3)Transform the variable into the variable for each scenario : (4)Calculate the covariance matrix of : (5)Maximize over the degrees of freedom the log likelihood of the Student's tcopula density: where
3.2. Hypothesis Testing
In the section, we describe our methodology to determine whether or not we made the right assumption, that the Student's copula fits our counterparty risk premium (simulated) data better than the Gaussian copula. The goodnessoffit test we use to compare the empirical distribution with the hypothesized cumulative distribution is called the KolmogorovSmirnov test. We perform two KS tests: () empirical versus normal and () empirical versus Student's .
3.2.1. First KS Test: Empirical versus Normal
The first KS test rests on the fact that if the copula of the counterparty risk premium is normal then the distribution of the variable (defined below) should be chisquare with degrees of freedom equal to the number of the firms in the network . For a given network, we conduct the first KS test as follows.
The hypotheses are We follow the methodology discussed in Malevergne and Sornette [10].
() Find the empirical distribution of and calculate as in the parameter estimation step: where
() Calculate the Gaussian variables via the following transformation:
() Determine the covariance matrix of the Gaussian variables
() Calculate the variable
() Find the empirical distribution of , namely,
() Calculate the Kolmogorov distance:
() Verify the result by using the kstest2 (Kolmogorov test to compare the distribution of two samples) command in Matlab.
3.2.2. Second KS Test: Empirical versus
As previously mentioned, the second KS test will allow us to compare the empirical distribution of the copula to the theoretical Student's distribution. We will then compare the results of the second KS test to those of the first KS test. We proceed with the second KS test as follows:
Unlike in the Gaussian copula case, there is no analytical expression for the density of the theoretical KolmogorovSmirnov statistic in the tcopula case. Thus we generate by simulation (over a very large number of scenarios ) the distribution . The overall algorithm is as follows.
() Generate the “theoretical” data by Conditional Monte Carlo (CMC) simulation following the method adapted from that of Aas et al. [11] and find the empirical distribution of the theoretical data. Please refer to the appendix for the Conditional Monte Carlo (CMC) simulation algorithm.
() Calculate :
() Calculate the variables
() Determine the covariance matrix of the variables .
() Calculate the variable :
() Find the theoretical distribution of , namely, :
() Calculate the Kolmogorov distance
() Verify the result by using the kstest2 command in Matlab.
4. Results
4.1. Parameter Estimation
The principal result that we obtain is that, for counterparty risk premia, any Student's tcopula results in a better fit than the Gaussian copula.
We show the loglikelihood defined on the righthandside of (3.1) as a function of the number of degrees of freedom of the Student's distribution. In our simulation, we use various numbers of scenarios: 50, 100, 150, 200, 300, 500, 700 to 900 scenarios. For the triangular and series networks, the number of firms varies from 2, 3, 5, 10, 20, 50, 100 to 200. In all these cases, the loglikelihood function is uniformly decreasing. In other words, the degrees of freedom that maximize the loglikelihood are finite numbers in all the cases.
We show hereafter a subset of our simulation results. In Figures 5 to 16 are shown the results for the triangular, or “full” network. In Figure 17, are shown the results for the series network, while in Figure 18 we show results for the star network. The benchmark marginal distribution is exponential. For the triangular network, we tried also all the marginal distributions specified earlier in order to verify that the marginal distributions of the firms' production revenues have no effect on the dependence structure of the their counterparty risk premia.
It is worth mentioning again that the loglikelihood function is uniformly decreasing across all samples. As mentioned above, the Gaussian distribution is a Student's distribution with an infinite number of degrees of freedom. Thus, from our parameter estimation step, we conclude that any Student's distribution results in a better fit than the Gaussian distribution. This fact will be confirmed when we perform hypothesis testing for goodnessoffit check.
4.1.1. “Full” or Triangular Networks
4.1.2. Other Networks
4.2. Hypothesis Testing
4.2.1. First KS Test: Empirical versus Normal
We analyze the Kolmogorov statistic which represents the maximum distance between the empirical distribution and the distribution.
As in the parameter estimation step, we varied our revenue scenarios from 50, 100, 150, 200, 300, 500, 700 to 900 scenarios and the number of firms from 2, 3, 5, 10, 20, 50, 100 to 200 in each of our networks. As the number of firms in the triangular network increases, the Kolmogorov statistic increases, thus allowing us to reject the null hypothesis that the copula of the counterparty risk premium is Gaussian if the number of firms is sufficiently large. We also use Matlab's kstest2 function to verify the reliability of our Kolmogorov statistic . Note that our Kolmogorov statistic takes into account only the maximum difference value while the Matlab's kstest2 function takes into account the difference values at all data points between the empirical distribution and the theoretical distribution (, in this case). The Matlab's kstest2 results are consistent with our conclusion. Specifically, when the number of firms in the network exceeds 10, we reject the null hypothesis in favor of the alternative hypothesis at a 5% significance level. In plain words, when the number of firms exceeds 10 in the network, the copula that captures the dependence of the counterparty risk premium is not Gaussian.
Similar results hold for the series and “star” networks. In the series network, when the number of firms in the network exceeds 30, we reject the null hypothesis that the copula is Gaussian. We also reject the null hypothesis in the “star” network of 12 firms. In sum, when the number of the firms increases, we are more likely to reject the null hypothesis. This leads us to conclude that as the number of firms increases, the copula that captures the dependence of the counterparty risk premium is more likely to be nonGaussian.
4.2.2. Second KS Test: Empirical versus
Note that the copula is likely to be nonGaussian, is it a Student's tcopula? Which copula, Gaussian or Student's , is more likely? The second KS test will allow us to compare the empirical distribution of the copula to the theoretical Student's distribution.
As before, we did the test on the triangular or “full” network, the series network, and the “star” network. For each type of network, we varied the number of revenue scenarios from 50, 100, 150, 200, 300, 500, 700 to 900 scenarios and the number of firms from 2, 3, 5, 10, 20, 50, 100 to 200 firms. The results we obtained indicate that we cannot reject the hypothesis that the copula is Student's in all cases at the 5% significance level. In other words, we do not have sufficient evidence to reject that the copula capturing the dependence of the counterparty risk premium is Student's .
4.2.3. PValue
For networks of 2–5 firms, the values indicate that we cannot reject the null hypothesis in neither the first nor the second KolmogorovSmirnov tests. However, the values resulting from the first KS test (which tests whether the empirical distribution comes from the normal distribution) are lower than the values resulting from the second KS test (which tests whether the empirical distribution comes from the Student's distribution) in all cases. For networks of more than 10 firms, the values indicate that we can reject the null hypothesis in the first KS test but not the second KS test. In other words, in networks of more than 10 firms, we have sufficient evidence to conclude that the empirical distribution does not come from the normal distribution but we do not have sufficient evidence to conclude that the empirical distribution does not come from the Student's distribution.
4.2.4. Test Power
It is well known that the onesided KolmogorovSmirnov test is superior in many aspects to the traditional goodnessoffit test: see, for instance, Massey [12]. The latter paper provides a figure, from which the power of that test is given. For twosided KolmogorovSmirnov tests, much less is known (see however Milbrodt and Strasser [13]). Given the complexity of calculating the power of that test, we doubt that many practioners engage in it, especially since more powerful tests probably do not exist.
5. Conclusion
In our study, we use the copula method to model the dependence between the counterparty risk premia of various firms. Specifically, we study the impact of the dependence among production revenues on the dependence among counterparty risk premia in several different network economies.
Taking as given that the dependence between the production revenues was a Gaussian copula, we generate the production revenues by simulation. Then we calculate the total revenues from the simulated production revenues. Two main tests—the parameter estimation and hypothesis testing—are carried out on each network setup. In the first test, the parameter of the copula being estimated is the number of degrees of freedom. The estimation results obtained from each network indicate that the Student's tcopula is likely to be the copula capturing the relationship between the firms' counterparty risk premia, reaffirming our assumption. The results from our second test—the KolmogorovSmirnov hypothesis testing—confirms our conclusion.
Appendices
Conditional Monte Carlo (CMC) simulation follows the method adapted from that of Aas et al. [11]
A. Generation of the Theoretical Data
To generate the “theoretical” sample data from the Student's copula, we use conditional Monte Carlo (CMC) simulation. Provided that multivariate data can be modelled using a set of paircopulae which act on two variables at a time, we generate samples from the Student's copula by way of conditional Monte Carlo simulation.
Aas et al. [11] provide a way to simplify the copula function in order to derive an efficient simulation algorithm as follows. We define
We can then calculate Now, set Then Since the final expression for is The inverse function is given by Using the simplified paircopulae functions and shown previously, the general algorithm for sampling from the “theoretical” Student's copula becomes.(1)Sample where independent standard uniform.(2)Given that and , we can calculate as follows: where
B. Monte Carlo Simulation Algorithm for the “Full” Network
As before, let be the total number of firms in the network and let be the total number of revenue scenarios for each firm in the network. and are the functions derived previously. The detailed algorithm of our simulation is as in Algorithm 1.

Acknowledgment
The authors thank Fitch Ratings for their financial and intellectual support. All errors are ours.
Endnotes
 By this, we mean that the dependence between tail events is stronger for the tcopula than for the Gaussian copula.