# A Theoretical Argument Why the t-Copula Explains Credit Risk Contagion Better than the Gaussian Copula

Authors: Didier Cossin; Henry Schellhorn; Nan Song; Satjaporn Tungsong
Journal: Advances in Decision Sciences
 (2010)

---

## 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 t-copula 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 t-copula 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 t-copulae induce fatter (joint) tails than Gaussian copulaeceteris 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.


---

## Body

## 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 that-copulae 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 5.3]).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,t, 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 wide-spread 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 t-copulae with the Gaussian copula as a way to model counterparty risk. We found out by simulation that the t-copula 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 t-copula 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 t-copulae 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 null-hypotheses: Gaussian copula, and t-copula. 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.Figure 1
An example of a network economy with 3 firms.We callν¯i the long run operational revenue rate of firm i. In the CS model, this rate is not directly observable. At each time t the market calculates an estimator ν̂i(t) of ν¯i. For the sake of brevity in this document, we will call ν̂i simply the production revenue. The total revenueα̂i of firm i 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 k to firm i are, at all times, proportional to the total revenue of firm k. In other terms, there are constants λki so that(2.1)α̂i(t)=ν̂i(t)+∑kα̂k(t)λki.Thepayout ratioδi expresses what percentage of total expenses of firm i is distributed to bondholders and equityholders. It is modelled as a geometric Brownian motion with relative drift μi and volatility σi. 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 i incurs a loss rate of wi. Theorem 2.1 was proved in Cossin and Schellhorn [5]:Theorem 2.1.
In steady state, the values of equityS and debt B are, for finite N, and t not a bankruptcy time:
(2.2)Si(t,ω)=α̂i(t,ω)s(δi(t,ω);δi(0),μi,σi,wi,r)+O(1N),Bi(t,ω)=α̂i(t,ω)b(δi(t,ω);δi(0),μi,σi,wi,r)+O(1N).The parameterN scales the production revenue intensity, as is common in diffusion approximations. We refer to Cossin and Schellhorn [5] for a full definition of N. Roughly speaking, a large N corresponds to a large cash account, which acts as a (partial) buffer against cash flow risk. In this paper, we make the assumption that N is very large, for simplicity. The functions s and b 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 [t,t+T]. 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 Bi(t,T,ω) be the total value of debt. There is now a new function bLT so that
(2.3)Bi(t,T,ω)=α̂i(t,ω)biLT(δi(t,ω),T;δi(0),μi,σi,wi,r)+O(1N).
Differentiating (2.3) with respect to T, we see that the value of debt with maturity [T,T+dT] is equal to
(2.4)∂∂TBi(t,T,ω)dT=α̂i(t,ω)∂∂TbiLT(δi(t,ω),T;δi(0),μi,σi,wi,r)dT+O(1N).
On the other hand, we can rewrite the left-hand-side as
(2.5)∂∂TBi(t,T,ω)dT≡LiTexp(-(r+csi(t,T))T)dT,
where Li is the total principal and csi(T) the credit spread of firm i for maturity T. For simplicity we write(2.6)φi(t,T,ω)=∂biLT∂T|δi(t,ω),T.
Thus, equating (2.4) with (2.5), and taking N→∞, we have(2.7)r+csi(t,T,ω)=1T[ln(Liα̂i(t,ω)T)-lnφi(t,T,ω)].
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 α̂i than the credit risk coming from φi in the last expression. We thus define the counterparty risk premiumcrpi and a bankruptcy risk premiumbrpi:(2.8)crpi=1Tln(Liα̂i(t,ω)T),(2.9)brpi=-1Tlnφi(t,T,ω)-r.The counterparty risk premium stems from the network termα̂i and the bankruptcy risk premium comes from the idiosyncratic variable δi, 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:(2.10)crpi+brpi=csi.It is well known that the copula function is unchanged under monotone transformations, such as the transformation (2.8) between total revenue α̂i and counterparty risk premium crpi. 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ν̂i and the total revenue α̂i:(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, firmm makes loans to firms 1,2,…,m-1. Firm m-1 loans to firms 1,2,…,m-2, and so on. The network  equations are then (see also Figure 2) (2.11)α1(ω)=ν1(ω),α2(ω)=λ12*ν1(ω)+ν2(ω),α3(ω)=λ13*ν1(ω)+λ23*ν2(ω)+ν3(ω),⋮αm(ω)=λ1m*ν1(ω)+λ2m*ν2(ω)+⋯+λm-1,m*νm-1(ω)+νm(ω).Figure 2
A Triangular or “Full” network.We analyze the specific caseλi,j=0.5 for i<j. Also, all the firms have the same size, that is, E[νi]=1.


#### 2.2.2. Series Network

In a series network, each firmm-1 borrows from the “next” firm m. It is illustrated in Figure 3: (2.12)α1(ω)=ν1(ω),α2(ω)=0.5α1(ω)+ν2(ω),α3(ω)=0.5α2(ω)+ν3(ω),α4(ω)=0.5α3(ω)+ν4(ω),⋮αm(ω)=0.5αm-1,m(ω)+νm(ω),
where λm-1,m=0.5.Figure 3
A Series network.


#### 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,E[νi]=1. All large firms also have the same revenue size, E[νi]=5.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. Figure4 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.Figure 4
The “Star” network.


### 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 variance-covariance matrix forν:(2.13)Var[ννt]=[1ρ⋯ρρ1··ρρ·ρ1]withρ=0.5. 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 F for production revenue ν̂i: (i)
exponential with various rates, namely, 1, 5, and 20(ii)
uniform with various support,(iii)
gamma with various shapes and scales,(iv)
quadraticF(x)=x2 where x∈[0,1],(v)
cubicF(x)=x3 where x∈[0,1].Exponential, gamma, and uniform random variables are standard choices for modelling positive random variables (see, e.g., Lando [9]).


## 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.Figure 1
An example of a network economy with 3 firms.We callν¯i the long run operational revenue rate of firm i. In the CS model, this rate is not directly observable. At each time t the market calculates an estimator ν̂i(t) of ν¯i. For the sake of brevity in this document, we will call ν̂i simply the production revenue. The total revenueα̂i of firm i 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 k to firm i are, at all times, proportional to the total revenue of firm k. In other terms, there are constants λki so that(2.1)α̂i(t)=ν̂i(t)+∑kα̂k(t)λki.Thepayout ratioδi expresses what percentage of total expenses of firm i is distributed to bondholders and equityholders. It is modelled as a geometric Brownian motion with relative drift μi and volatility σi. 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 i incurs a loss rate of wi. Theorem 2.1 was proved in Cossin and Schellhorn [5]:Theorem 2.1.
In steady state, the values of equityS and debt B are, for finite N, and t not a bankruptcy time:
(2.2)Si(t,ω)=α̂i(t,ω)s(δi(t,ω);δi(0),μi,σi,wi,r)+O(1N),Bi(t,ω)=α̂i(t,ω)b(δi(t,ω);δi(0),μi,σi,wi,r)+O(1N).The parameterN scales the production revenue intensity, as is common in diffusion approximations. We refer to Cossin and Schellhorn [5] for a full definition of N. Roughly speaking, a large N corresponds to a large cash account, which acts as a (partial) buffer against cash flow risk. In this paper, we make the assumption that N is very large, for simplicity. The functions s and b 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 [t,t+T]. 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 Bi(t,T,ω) be the total value of debt. There is now a new function bLT so that
(2.3)Bi(t,T,ω)=α̂i(t,ω)biLT(δi(t,ω),T;δi(0),μi,σi,wi,r)+O(1N).
Differentiating (2.3) with respect to T, we see that the value of debt with maturity [T,T+dT] is equal to
(2.4)∂∂TBi(t,T,ω)dT=α̂i(t,ω)∂∂TbiLT(δi(t,ω),T;δi(0),μi,σi,wi,r)dT+O(1N).
On the other hand, we can rewrite the left-hand-side as
(2.5)∂∂TBi(t,T,ω)dT≡LiTexp(-(r+csi(t,T))T)dT,
where Li is the total principal and csi(T) the credit spread of firm i for maturity T. For simplicity we write(2.6)φi(t,T,ω)=∂biLT∂T|δi(t,ω),T.
Thus, equating (2.4) with (2.5), and taking N→∞, we have(2.7)r+csi(t,T,ω)=1T[ln(Liα̂i(t,ω)T)-lnφi(t,T,ω)].
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 α̂i than the credit risk coming from φi in the last expression. We thus define the counterparty risk premiumcrpi and a bankruptcy risk premiumbrpi:(2.8)crpi=1Tln(Liα̂i(t,ω)T),(2.9)brpi=-1Tlnφi(t,T,ω)-r.The counterparty risk premium stems from the network termα̂i and the bankruptcy risk premium comes from the idiosyncratic variable δi, 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:(2.10)crpi+brpi=csi.It is well known that the copula function is unchanged under monotone transformations, such as the transformation (2.8) between total revenue α̂i and counterparty risk premium crpi. 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ν̂i and the total revenue α̂i:(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, firmm makes loans to firms 1,2,…,m-1. Firm m-1 loans to firms 1,2,…,m-2, and so on. The network  equations are then (see also Figure 2) (2.11)α1(ω)=ν1(ω),α2(ω)=λ12*ν1(ω)+ν2(ω),α3(ω)=λ13*ν1(ω)+λ23*ν2(ω)+ν3(ω),⋮αm(ω)=λ1m*ν1(ω)+λ2m*ν2(ω)+⋯+λm-1,m*νm-1(ω)+νm(ω).Figure 2
A Triangular or “Full” network.We analyze the specific caseλi,j=0.5 for i<j. Also, all the firms have the same size, that is, E[νi]=1.


### 2.2.2. Series Network

In a series network, each firmm-1 borrows from the “next” firm m. It is illustrated in Figure 3: (2.12)α1(ω)=ν1(ω),α2(ω)=0.5α1(ω)+ν2(ω),α3(ω)=0.5α2(ω)+ν3(ω),α4(ω)=0.5α3(ω)+ν4(ω),⋮αm(ω)=0.5αm-1,m(ω)+νm(ω),
where λm-1,m=0.5.Figure 3
A Series network.


### 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,E[νi]=1. All large firms also have the same revenue size, E[νi]=5.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. Figure4 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.Figure 4
The “Star” network.


## 2.2.1. Triangular or “Full” Network

In this type of network, firmm makes loans to firms 1,2,…,m-1. Firm m-1 loans to firms 1,2,…,m-2, and so on. The network  equations are then (see also Figure 2) (2.11)α1(ω)=ν1(ω),α2(ω)=λ12*ν1(ω)+ν2(ω),α3(ω)=λ13*ν1(ω)+λ23*ν2(ω)+ν3(ω),⋮αm(ω)=λ1m*ν1(ω)+λ2m*ν2(ω)+⋯+λm-1,m*νm-1(ω)+νm(ω).Figure 2
A Triangular or “Full” network.We analyze the specific caseλi,j=0.5 for i<j. Also, all the firms have the same size, that is, E[νi]=1.


## 2.2.2. Series Network

In a series network, each firmm-1 borrows from the “next” firm m. It is illustrated in Figure 3: (2.12)α1(ω)=ν1(ω),α2(ω)=0.5α1(ω)+ν2(ω),α3(ω)=0.5α2(ω)+ν3(ω),α4(ω)=0.5α3(ω)+ν4(ω),⋮αm(ω)=0.5αm-1,m(ω)+νm(ω),
where λm-1,m=0.5.Figure 3
A Series network.


## 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,E[νi]=1. All large firms also have the same revenue size, E[νi]=5.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. Figure4 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.Figure 4
The “Star” network.


## 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 variance-covariance matrix forν:(2.13)Var[ννt]=[1ρ⋯ρρ1··ρρ·ρ1]withρ=0.5. 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 F for production revenue ν̂i: (i)
exponential with various rates, namely, 1, 5, and 20(ii)
uniform with various support,(iii)
gamma with various shapes and scales,(iv)
quadraticF(x)=x2 where x∈[0,1],(v)
cubicF(x)=x3 where x∈[0,1].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 t-copula with a Gaussian copula as a way to model counterparty risk. Since the Gaussian copula is a special case of a t-copula, namely a t-copula with an infinite number of degrees of freedom, we try to fit a t-copula to our simulated data. The calibrated number of degrees of freedom will be a good indicator whether a nonGaussian t-copula 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 t-copula, 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)θ1 of the marginals,(2)
infer the parameter(s)θ2 of the Student's t-copula.Our analyses focus on the second step: estimating the parameter(s) of the copulaθ2. We bypass the first step because we take the empirical distribution of the marginals as given. The Student's t-copula has two parameters, namely, the correlation matrix R 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 ν:(3.1)ν̂=argmaxν∑ω=1ΩlncR,ν(F1(αω,1),F2(αω,2),…,Fn(αω,m)),
where
(3.2)cR,ν(u1,u2,…,um)=|R|-1/2Γ((ν+2)/2)Γ(ν/2)(Γ(ν/2)Γ((ν+1)/2))2(1+(1/ν)ςtR-1ς)-(ν+2)/2∏j=1m(1+ςj2/ν)-(ν+1)/2,ςj=tν-1(uj).Letm 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:
(3.3)α11,α21,…,αΩ1,α12,α22,…,αΩ2,⋮α1m,α2m,…,αΩm.(2)
Calculate, for each scenarioω,
(3.4)Ûω1=F̂1(αω,1)=1Ω∑i,k=1Ω1{αi,1<αk,1},Ûω2=F̂2(αω,2)=1Ω∑i,k=1Ω1{αi,2<αk,2},⋮Ûωm=F̂2(αω,m)=1Ω∑i,k=1Ω1{αi,m<αk,m}.(3)
Transform the variableÛ into the variable ζ for each scenario ω:
(3.5)ζω1=tν-1(Ûω1),ζω2=tν-1(Ûω2),⋮ζωm=tν-1(Ûωm).(4)
Calculate the covariance matrix ofζ:
(3.6)Remp=Cov(ζω1,ζω2,…,ζωm).(5)
Maximize over the degrees of freedomν the log likelihood of the Student's t-copula density:
(3.7)maxν∑ω=1Ωlnc(F1(αω,1),F2(αω,2),…,F2(αω,n))=maxln(|Remp|-1/2Γ((ν+2)/2)Γ(ν/2)(Γ(ν/2)Γ((ν+1)/2))2(1+(1/ν)ςitRemp-1ςi)-(ν+2)/2∏j=1n(1+ςi,j2/ν)-(ν+1)/2),
where
(3.8)ςω,1=tν-1(Ûω1),ςω,2=tν-1(Ûω2),⋮ςω,n=tν-1(Ûωn).


### 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 goodness-of-fit test we use to compare the empirical distribution with the hypothesized cumulative distribution is called the Kolmogorov-Smirnov test. We perform two K-S tests: (1) empirical versus normal and (2) empirical versus Student's t.


#### 3.2.1. First K-S Test: Empirical versus Normal

The first K-S test rests on the fact that if the copula of the counterparty risk premiumCα1,α2,…,αm is normal then the distribution of the variable zemp2 (defined below) should be chi-square χ2 with degrees of freedom equal to the number of the firms in the network m. For a given network, we conduct the first K-S test as follows.The hypotheses are(3.9)H0:thecopulaCα1,α2,…,αmisGaussian,H1:thecopulaCα1,α2,…,αmisnotGaussian.
We follow the methodology discussed in Malevergne and Sornette [10].(1)  Find the empirical distribution F̂αm of αm and calculate Ûm as in the parameter estimation step:(3.10)Û1=(F̂α1(α1)),⋮Ûn=(F̂αn(αm)),
where
(3.11)F̂1(αω,1)=1Ω∑i,k=1Ω1{αi,1<αk,1},⋮F̂n(αω,m)=1Ω∑i,k=1Ω1{αi,m<αk,m}.(2) Calculate the Gaussian variables yω1,yω2,…,yωm via the following transformation:(3.12)yω1=Φ-1(Ûω1)=Φ-1(F̂α1(αω1)),⋮yωm=Φ-1(Ûωn)=Φ-1(F̂αn(αωm)).(3) Determine the covariance matrix Remp of the Gaussian variables y1,y2,…,ym(4) Calculate the variable zemp.2(3.13)zemp.2(ω)=ytRemp-1y.(5) Find the empirical distribution of zemp.2, namely, F̂zemp2(6) Calculate the Kolmogorov distance:(3.14)D=maxzemp2|F̂zemp2(zemp2)-Fχ2(zemp2)|.(7) Verify the result by using the kstest2 (Kolmogorov test to compare the distribution of two samples) command in Matlab.


#### 3.2.2. Second K-S Test: Empirical versust

As previously mentioned, the second K-S test will allow us to compare the empirical distribution of the copula to the theoretical Student'st distribution. We will then compare the results of the second K-S test to those of the first K-S test. We proceed with the second K-S test as follows:(3.15)H0:thecopulaCα1,α2,…,αmisStudent′scopula,H1:thecopulaCα1,α2,…,αmisnotStudent′scopula.Unlike in the Gaussian copula case, there is no analytical expression for the density of the theoretical Kolmogorov-Smirnov statisticzth.2 in the t-copula case. Thus we generate by simulation (over a very large number of scenarios Ω¯) the distribution F̂zth.2. The overall algorithm is as follows.(1) Generate the “theoretical” data (xω1,xω2,…,xωm) by Conditional Monte Carlo (CMC) simulation following the method adapted from that of Aas et al. [11] and find the empirical distribution F̂xm of the theoretical data. Please refer to the appendix for the Conditional Monte Carlo (CMC) simulation algorithm.(2) Calculate Ûm:
(3.16)Û1th=(F̂x1(x1)),⋮Ûmth=(F̂xm(xm)).(3) Calculate the variables yω1th,yω2th,…,yωnth(3.17)yth(ω,1)=Φ-1(Û1th(ω)),⋮yth(ω,m)=Φ-1(Ûmth(ω)).(4) Determine the covariance matrix Rth of the variables y1th,y2th,…,ymth.(5) Calculate the variable zth2:(3.18)zth.2(ω)=(yth)tRth-1yth.(6) Find the theoretical distribution of zth.2, namely, F̂zth2:(7) Calculate the Kolmogorov distance
(3.19)D=maxzemp2|F̂zemp2(zemp2)-Fzth2(zemp2)|.(8) Verify the result by using the kstest2 command in Matlab.


## 3.1. Parameter Estimation

Assuming that the dependence of the firms' counterparty risk premia is the Student's t-copula, 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)θ1 of the marginals,(2)
infer the parameter(s)θ2 of the Student's t-copula.Our analyses focus on the second step: estimating the parameter(s) of the copulaθ2. We bypass the first step because we take the empirical distribution of the marginals as given. The Student's t-copula has two parameters, namely, the correlation matrix R 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 ν:(3.1)ν̂=argmaxν∑ω=1ΩlncR,ν(F1(αω,1),F2(αω,2),…,Fn(αω,m)),
where
(3.2)cR,ν(u1,u2,…,um)=|R|-1/2Γ((ν+2)/2)Γ(ν/2)(Γ(ν/2)Γ((ν+1)/2))2(1+(1/ν)ςtR-1ς)-(ν+2)/2∏j=1m(1+ςj2/ν)-(ν+1)/2,ςj=tν-1(uj).Letm 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:
(3.3)α11,α21,…,αΩ1,α12,α22,…,αΩ2,⋮α1m,α2m,…,αΩm.(2)
Calculate, for each scenarioω,
(3.4)Ûω1=F̂1(αω,1)=1Ω∑i,k=1Ω1{αi,1<αk,1},Ûω2=F̂2(αω,2)=1Ω∑i,k=1Ω1{αi,2<αk,2},⋮Ûωm=F̂2(αω,m)=1Ω∑i,k=1Ω1{αi,m<αk,m}.(3)
Transform the variableÛ into the variable ζ for each scenario ω:
(3.5)ζω1=tν-1(Ûω1),ζω2=tν-1(Ûω2),⋮ζωm=tν-1(Ûωm).(4)
Calculate the covariance matrix ofζ:
(3.6)Remp=Cov(ζω1,ζω2,…,ζωm).(5)
Maximize over the degrees of freedomν the log likelihood of the Student's t-copula density:
(3.7)maxν∑ω=1Ωlnc(F1(αω,1),F2(αω,2),…,F2(αω,n))=maxln(|Remp|-1/2Γ((ν+2)/2)Γ(ν/2)(Γ(ν/2)Γ((ν+1)/2))2(1+(1/ν)ςitRemp-1ςi)-(ν+2)/2∏j=1n(1+ςi,j2/ν)-(ν+1)/2),
where
(3.8)ςω,1=tν-1(Ûω1),ςω,2=tν-1(Ûω2),⋮ςω,n=tν-1(Ûωn).


## 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 goodness-of-fit test we use to compare the empirical distribution with the hypothesized cumulative distribution is called the Kolmogorov-Smirnov test. We perform two K-S tests: (1) empirical versus normal and (2) empirical versus Student's t.


### 3.2.1. First K-S Test: Empirical versus Normal

The first K-S test rests on the fact that if the copula of the counterparty risk premiumCα1,α2,…,αm is normal then the distribution of the variable zemp2 (defined below) should be chi-square χ2 with degrees of freedom equal to the number of the firms in the network m. For a given network, we conduct the first K-S test as follows.The hypotheses are(3.9)H0:thecopulaCα1,α2,…,αmisGaussian,H1:thecopulaCα1,α2,…,αmisnotGaussian.
We follow the methodology discussed in Malevergne and Sornette [10].(1)  Find the empirical distribution F̂αm of αm and calculate Ûm as in the parameter estimation step:(3.10)Û1=(F̂α1(α1)),⋮Ûn=(F̂αn(αm)),
where
(3.11)F̂1(αω,1)=1Ω∑i,k=1Ω1{αi,1<αk,1},⋮F̂n(αω,m)=1Ω∑i,k=1Ω1{αi,m<αk,m}.(2) Calculate the Gaussian variables yω1,yω2,…,yωm via the following transformation:(3.12)yω1=Φ-1(Ûω1)=Φ-1(F̂α1(αω1)),⋮yωm=Φ-1(Ûωn)=Φ-1(F̂αn(αωm)).(3) Determine the covariance matrix Remp of the Gaussian variables y1,y2,…,ym(4) Calculate the variable zemp.2(3.13)zemp.2(ω)=ytRemp-1y.(5) Find the empirical distribution of zemp.2, namely, F̂zemp2(6) Calculate the Kolmogorov distance:(3.14)D=maxzemp2|F̂zemp2(zemp2)-Fχ2(zemp2)|.(7) Verify the result by using the kstest2 (Kolmogorov test to compare the distribution of two samples) command in Matlab.


### 3.2.2. Second K-S Test: Empirical versust

As previously mentioned, the second K-S test will allow us to compare the empirical distribution of the copula to the theoretical Student'st distribution. We will then compare the results of the second K-S test to those of the first K-S test. We proceed with the second K-S test as follows:(3.15)H0:thecopulaCα1,α2,…,αmisStudent′scopula,H1:thecopulaCα1,α2,…,αmisnotStudent′scopula.Unlike in the Gaussian copula case, there is no analytical expression for the density of the theoretical Kolmogorov-Smirnov statisticzth.2 in the t-copula case. Thus we generate by simulation (over a very large number of scenarios Ω¯) the distribution F̂zth.2. The overall algorithm is as follows.(1) Generate the “theoretical” data (xω1,xω2,…,xωm) by Conditional Monte Carlo (CMC) simulation following the method adapted from that of Aas et al. [11] and find the empirical distribution F̂xm of the theoretical data. Please refer to the appendix for the Conditional Monte Carlo (CMC) simulation algorithm.(2) Calculate Ûm:
(3.16)Û1th=(F̂x1(x1)),⋮Ûmth=(F̂xm(xm)).(3) Calculate the variables yω1th,yω2th,…,yωnth(3.17)yth(ω,1)=Φ-1(Û1th(ω)),⋮yth(ω,m)=Φ-1(Ûmth(ω)).(4) Determine the covariance matrix Rth of the variables y1th,y2th,…,ymth.(5) Calculate the variable zth2:(3.18)zth.2(ω)=(yth)tRth-1yth.(6) Find the theoretical distribution of zth.2, namely, F̂zth2:(7) Calculate the Kolmogorov distance
(3.19)D=maxzemp2|F̂zemp2(zemp2)-Fzth2(zemp2)|.(8) Verify the result by using the kstest2 command in Matlab.


## 3.2.1. First K-S Test: Empirical versus Normal

The first K-S test rests on the fact that if the copula of the counterparty risk premiumCα1,α2,…,αm is normal then the distribution of the variable zemp2 (defined below) should be chi-square χ2 with degrees of freedom equal to the number of the firms in the network m. For a given network, we conduct the first K-S test as follows.The hypotheses are(3.9)H0:thecopulaCα1,α2,…,αmisGaussian,H1:thecopulaCα1,α2,…,αmisnotGaussian.
We follow the methodology discussed in Malevergne and Sornette [10].(1)  Find the empirical distribution F̂αm of αm and calculate Ûm as in the parameter estimation step:(3.10)Û1=(F̂α1(α1)),⋮Ûn=(F̂αn(αm)),
where
(3.11)F̂1(αω,1)=1Ω∑i,k=1Ω1{αi,1<αk,1},⋮F̂n(αω,m)=1Ω∑i,k=1Ω1{αi,m<αk,m}.(2) Calculate the Gaussian variables yω1,yω2,…,yωm via the following transformation:(3.12)yω1=Φ-1(Ûω1)=Φ-1(F̂α1(αω1)),⋮yωm=Φ-1(Ûωn)=Φ-1(F̂αn(αωm)).(3) Determine the covariance matrix Remp of the Gaussian variables y1,y2,…,ym(4) Calculate the variable zemp.2(3.13)zemp.2(ω)=ytRemp-1y.(5) Find the empirical distribution of zemp.2, namely, F̂zemp2(6) Calculate the Kolmogorov distance:(3.14)D=maxzemp2|F̂zemp2(zemp2)-Fχ2(zemp2)|.(7) Verify the result by using the kstest2 (Kolmogorov test to compare the distribution of two samples) command in Matlab.


## 3.2.2. Second K-S Test: Empirical versust

As previously mentioned, the second K-S test will allow us to compare the empirical distribution of the copula to the theoretical Student'st distribution. We will then compare the results of the second K-S test to those of the first K-S test. We proceed with the second K-S test as follows:(3.15)H0:thecopulaCα1,α2,…,αmisStudent′scopula,H1:thecopulaCα1,α2,…,αmisnotStudent′scopula.Unlike in the Gaussian copula case, there is no analytical expression for the density of the theoretical Kolmogorov-Smirnov statisticzth.2 in the t-copula case. Thus we generate by simulation (over a very large number of scenarios Ω¯) the distribution F̂zth.2. The overall algorithm is as follows.(1) Generate the “theoretical” data (xω1,xω2,…,xωm) by Conditional Monte Carlo (CMC) simulation following the method adapted from that of Aas et al. [11] and find the empirical distribution F̂xm of the theoretical data. Please refer to the appendix for the Conditional Monte Carlo (CMC) simulation algorithm.(2) Calculate Ûm:
(3.16)Û1th=(F̂x1(x1)),⋮Ûmth=(F̂xm(xm)).(3) Calculate the variables yω1th,yω2th,…,yωnth(3.17)yth(ω,1)=Φ-1(Û1th(ω)),⋮yth(ω,m)=Φ-1(Ûmth(ω)).(4) Determine the covariance matrix Rth of the variables y1th,y2th,…,ymth.(5) Calculate the variable zth2:(3.18)zth.2(ω)=(yth)tRth-1yth.(6) Find the theoretical distribution of zth.2, namely, F̂zth2:(7) Calculate the Kolmogorov distance
(3.19)D=maxzemp2|F̂zemp2(zemp2)-Fzth2(zemp2)|.(8) 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 t-copula results in a better fit than the Gaussian copula.We show the log-likelihood defined on the right-hand-side of (3.1) as a function of the number of degrees of freedom of the Student's t 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 log-likelihood function is uniformly decreasing. In other words, the degrees of freedom that maximize the log-likelihood are finite numbers in all the cases.We show hereafter a subset of our simulation results. In Figures5 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.Figure 5
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=1.It is worth mentioning again that the log-likelihood function is uniformly decreasing across all samples. As mentioned above, the Gaussian distribution is a Student'st-distribution with an infinite number of degrees of freedom. Thus, from our parameter estimation step, we conclude that any Student's t-distribution results in a better fit than the Gaussian distribution. This  fact will be confirmed when we perform hypothesis testing for goodness-of-fit check.


#### 4.1.1. “Full” or Triangular Networks

See Figures5–16.Figure 6
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=5.Figure 7
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=20.Figure 8
Inference on the triangular network with 3 firms. Marginal distribution is exponential.Figure 9
Inference on the triangular network with 4 firms. Marginal distribution is exponential.Figure 10
Inference on the triangular network with 2, 3, and 4 firms. Marginal distribution is standard uniform.Figure 11
Inference on the triangular network with 2 firms. Marginal distribution is uniform with support (5,10).Figure 12
Inference on the triangular network with 2 firms. Distribution is gamma.Figure 13
Inference on the triangular network with 3 firms. Marginal distribution is gamma.Figure 14
Inference on the triangular network with 4 firms. Marginal distribution is gamma.Figure 15
Inference on the triangular network. Marginal distribution is quadratic.Figure 16
Inference on the triangular network with 4 firms. Marginal distribution is cubic.Figure 17
Inference on the series network. Marginal distribution is exponential withE[νi]=1.Figure 18
Inference on the star network. Marginal distribution is exponential withE[νi]=1.


#### 4.1.2. Other Networks

See Figures17 and 18.


### 4.2. Hypothesis Testing

#### 4.2.1. First K-S Test: Empirical versus Normal

We analyze the Kolmogorov statisticD which represents the maximum distance between the empirical distribution F̂zemp2 and the χ2-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 statisticD 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 D. Note that our Kolmogorov statistic D 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 (χ2, 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 K-S Test: Empirical versust

Note that the copula is likely to be nonGaussian, is it a Student's t-copula? Which copula, Gaussian or Student'st, is more likely? The second K-S test will allow us to compare the empirical distribution of the copula to the theoretical Student's t 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'st 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 t.


#### 4.2.3. P-Value

For networks of 2–5 firms, theP-values indicate that we cannot reject the null hypothesis in neither the first nor the second Kolmogorov-Smirnov tests. However, the P-values resulting from the first K-S test (which tests whether the empirical distribution comes from the normal distribution) are lower than the P-values resulting from the second K-S test (which tests whether the empirical distribution comes from the Student's t distribution) in all cases. For networks of more than 10 firms, the P-values indicate that we can reject the null hypothesis in the first K-S test but not the second K-S 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 t distribution.


#### 4.2.4. Test Power

It is well known that the one-sided Kolmogorov-Smirnov test is superior in many aspects to the traditional goodness-of-fit test: see, for instance, Massey [12]. The latter paper provides a figure, from which the power of that test is given. For two-sided Kolmogorov-Smirnov 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.


## 4.1. Parameter Estimation

The principal result that we obtain is that, for counterparty risk premia, any Student's t-copula results in a better fit than the Gaussian copula.We show the log-likelihood defined on the right-hand-side of (3.1) as a function of the number of degrees of freedom of the Student's t 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 log-likelihood function is uniformly decreasing. In other words, the degrees of freedom that maximize the log-likelihood are finite numbers in all the cases.We show hereafter a subset of our simulation results. In Figures5 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.Figure 5
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=1.It is worth mentioning again that the log-likelihood function is uniformly decreasing across all samples. As mentioned above, the Gaussian distribution is a Student'st-distribution with an infinite number of degrees of freedom. Thus, from our parameter estimation step, we conclude that any Student's t-distribution results in a better fit than the Gaussian distribution. This  fact will be confirmed when we perform hypothesis testing for goodness-of-fit check.


### 4.1.1. “Full” or Triangular Networks

See Figures5–16.Figure 6
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=5.Figure 7
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=20.Figure 8
Inference on the triangular network with 3 firms. Marginal distribution is exponential.Figure 9
Inference on the triangular network with 4 firms. Marginal distribution is exponential.Figure 10
Inference on the triangular network with 2, 3, and 4 firms. Marginal distribution is standard uniform.Figure 11
Inference on the triangular network with 2 firms. Marginal distribution is uniform with support (5,10).Figure 12
Inference on the triangular network with 2 firms. Distribution is gamma.Figure 13
Inference on the triangular network with 3 firms. Marginal distribution is gamma.Figure 14
Inference on the triangular network with 4 firms. Marginal distribution is gamma.Figure 15
Inference on the triangular network. Marginal distribution is quadratic.Figure 16
Inference on the triangular network with 4 firms. Marginal distribution is cubic.Figure 17
Inference on the series network. Marginal distribution is exponential withE[νi]=1.Figure 18
Inference on the star network. Marginal distribution is exponential withE[νi]=1.


### 4.1.2. Other Networks

See Figures17 and 18.


## 4.1.1. “Full” or Triangular Networks

See Figures5–16.Figure 6
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=5.Figure 7
Inference on the triangular network with 2 firms. Marginal distribution is exponential withE[νi]=20.Figure 8
Inference on the triangular network with 3 firms. Marginal distribution is exponential.Figure 9
Inference on the triangular network with 4 firms. Marginal distribution is exponential.Figure 10
Inference on the triangular network with 2, 3, and 4 firms. Marginal distribution is standard uniform.Figure 11
Inference on the triangular network with 2 firms. Marginal distribution is uniform with support (5,10).Figure 12
Inference on the triangular network with 2 firms. Distribution is gamma.Figure 13
Inference on the triangular network with 3 firms. Marginal distribution is gamma.Figure 14
Inference on the triangular network with 4 firms. Marginal distribution is gamma.Figure 15
Inference on the triangular network. Marginal distribution is quadratic.Figure 16
Inference on the triangular network with 4 firms. Marginal distribution is cubic.Figure 17
Inference on the series network. Marginal distribution is exponential withE[νi]=1.Figure 18
Inference on the star network. Marginal distribution is exponential withE[νi]=1.


## 4.1.2. Other Networks

See Figures17 and 18.


## 4.2. Hypothesis Testing

### 4.2.1. First K-S Test: Empirical versus Normal

We analyze the Kolmogorov statisticD which represents the maximum distance between the empirical distribution F̂zemp2 and the χ2-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 statisticD 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 D. Note that our Kolmogorov statistic D 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 (χ2, 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 K-S Test: Empirical versust

Note that the copula is likely to be nonGaussian, is it a Student's t-copula? Which copula, Gaussian or Student'st, is more likely? The second K-S test will allow us to compare the empirical distribution of the copula to the theoretical Student's t 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'st 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 t.


### 4.2.3. P-Value

For networks of 2–5 firms, theP-values indicate that we cannot reject the null hypothesis in neither the first nor the second Kolmogorov-Smirnov tests. However, the P-values resulting from the first K-S test (which tests whether the empirical distribution comes from the normal distribution) are lower than the P-values resulting from the second K-S test (which tests whether the empirical distribution comes from the Student's t distribution) in all cases. For networks of more than 10 firms, the P-values indicate that we can reject the null hypothesis in the first K-S test but not the second K-S 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 t distribution.


### 4.2.4. Test Power

It is well known that the one-sided Kolmogorov-Smirnov test is superior in many aspects to the traditional goodness-of-fit test: see, for instance, Massey [12]. The latter paper provides a figure, from which the power of that test is given. For two-sided Kolmogorov-Smirnov 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.


## 4.2.1. First K-S Test: Empirical versus Normal

We analyze the Kolmogorov statisticD which represents the maximum distance between the empirical distribution F̂zemp2 and the χ2-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 statisticD 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 D. Note that our Kolmogorov statistic D 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 (χ2, 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 K-S Test: Empirical versust

Note that the copula is likely to be nonGaussian, is it a Student's t-copula? Which copula, Gaussian or Student'st, is more likely? The second K-S test will allow us to compare the empirical distribution of the copula to the theoretical Student's t 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'st 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 t.


## 4.2.3. P-Value

For networks of 2–5 firms, theP-values indicate that we cannot reject the null hypothesis in neither the first nor the second Kolmogorov-Smirnov tests. However, the P-values resulting from the first K-S test (which tests whether the empirical distribution comes from the normal distribution) are lower than the P-values resulting from the second K-S test (which tests whether the empirical distribution comes from the Student's t distribution) in all cases. For networks of more than 10 firms, the P-values indicate that we can reject the null hypothesis in the first K-S test but not the second K-S 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 t distribution.


## 4.2.4. Test Power

It is well known that the one-sided Kolmogorov-Smirnov test is superior in many aspects to the traditional goodness-of-fit test: see, for instance, Massey [12]. The latter paper provides a figure, from which the power of that test is given. For two-sided Kolmogorov-Smirnov 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 t-copula 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 Kolmogorov-Smirnov hypothesis testing—confirms our conclusion.Table 1(a)Full (Triangular) NetworkP-valueNo. of firmsNo. of revenue scenariosFirst K-S TestSecond K-S Test250.99581.00002100.99211.00002150.88201.00002200.61071.00002300.77621.00002500.81101.00002700.75621.00002900.53421.0000350.67791.00003100.80001.00003150.87381.00003200.73181.00003300.61071.00003500.60191.00003700.57651.00003900.14001.0000550.95411.00005100.44311.00005150.07181.00005200.52721.00005300.13861.00005500.21841.00005700.17531.00005900.03051.00001050.00211.000010100.00821.0000101501.8654e-0091.0000102001.4725e-0051.0000103009.6681e-0091.0000105005.6889e-0101.0000107002.3889e-0201.0000109007.6473e-0241.000050502.1647e-0231.0000501006.6643e-0191.0000501503.3610e-0371.0000502002.1543e-0271.0000503007.2496e-0671.0000505001.2755e-0691.0000507001.1159e-1071.0000509006.2242e-1311.0000100502.0685e-0171.00001001001.5506e-0451.00001001506.3211e-0591.00001002002.3167e-0551.00001003006.1185e-0741.00001005003.0720e-0901.00001007009.3211e-1281.00001009002.6643e-1831.0000200504.0089e-0151.00002001001.5506e-0451.00002001501.6517e-0281.00002002001.4175e-0891.00002003008.0437e-1191.00002005001.0424e-0971.00002007001.6714e-1251.00002009006.3297e-2031.0000(b)Series NetworkP-valueNo. of firmsNo. of revenue scenariosFirst K-S TestSecond K-S Test250.84091.00002100.96101.00002150.27371.00002200.95961.00002300.99571.00002500.95721.00002700.93421.00002900.84181.0000350.99581.00003100.09941.00003150.70551.00003200.52721.00003300.09341.00003500.60191.00003700.83841.00003900.77021.0000550.84091.00005100.89381.00005150.80001.00005200.37671.00005300.63971.00005500.85671.00005700.87481.00005900.01731.00001050.67791.000010100.79421.000010150.52721.000010200.63971.000010300.37621.000010500.22081.000010700.16721.000010900.03491.000050502.1647e-0231.0000501009.1220e-0091.0000501502.6562e-0071.0000502007.6668e-0101.0000503004.7542e-0051.000050500.05551.0000507008.9159e-0041.000050900.00251.0000100502.9719e-0091.00001001001.5506e-0451.00001001505.8193e-0071.00001002005.4170e-0091.00001003009.6737e-0091.00001005001.3138e-0091.00001007008.7499e-0101.00001009009.0278e-0071.0000200501.0799e-0081.00002001003.6964e-0121.00002001501.3332e-0241.00002002001.4175e-0891.00002003001.4008e-0341.00002005001.2495e-0211.00002007003.3716e-0141.00002009003.6571e-0221.0000(c)Star NetworkP-valueNo. of firmsNo. of revenue scenariosFirst K-S TestSecond K-S Test12502.9719e-0091.0000121009.1220e-0091.0000121501.9937e-0131.0000122002.4000e-0111.0000123008.6264e-0231.0000125001.1180e-0341.0000127001.3228e-0461.0000129001.2521e-0421.0000



---

*Source: 546547-2010-05-19.xml*