Cva Derivatives
Cva Derivatives
Cva Derivatives
1
{>
}
+()(
1
)1
{
1
<
=1
.
Figure 4.1: Structure of a CDS contract (Herbertsson, 2012).
15
Figure 4.2: Different scenarios with no default before time T, and default before time T (Herbertsson, 2012).
By dividing the default leg with the premium leg, we get the T year CDS spread () as
() =
1
{T}
()(1 )
(
1
{>
}
+()(
1
)1
{
1
<
=1
. (4.1)
The CDS spread is expressed in basis points (bp) per annum and the expectation is under the risk neutral
measure, is the recovery in the case of default and () is the discount factor,
() = ()
0
where
is the short term risk free interest rate and a deterministic function
of time,
= ().
However, the equation above can be simplified when we assume that:
1. The interest rate is a deterministic function of time, () = ()
0
, as well as
independent of the default time . (() is the risk free interest rate as a function of time .)
2. The credit loss = 1 is constant where is the recovery rate;
Then the Equation (4.1) can be rewritten as
16
() =
(1 )
()
1 (
) +
(
1
)()
=1
(4.2)
where () = [ ] is the distribution function of the default time for the obligor and
= [()].
Also when we assume that the interest rate
= (). Then
and
= [()] = ()
0
. Moreover, let be constant, and let
() =
()
, then we have
1
{}
()(1 ) = 1
{}
()(1 )
() = (1 ) ()
()
0
.
Furthermore
(
)
1
4
1
{>
}
= (
)
1
4
1
{>
}
= (
)
1
4
[ >
] = (
)
1
4
(1 (
))
and
()(
1
)1
{
1
<<
}
= 1
{
1
<<
}
(
1
)()
()
0
= ()(
1
)
()
is the accrued premium, which is a final payment done by A to B at the default time. The size of this
premium is related to the time interval between the last payment and the default. Hence we can write
the Equation (4.2) as
() =
(1 ) ()
()
0
(
)
1
4
1 (
) + ()(
1
)
()
4
=1
. (4.3)
The CDS-spread formula in Equation (4.3) can be simplified if we make two more assumptions:
1. Ignore the accrued premium term in the premium leg.
2. Assume that the loss is paid at time
=
4
, at the end of each quarter instead of at time when
default occurs in the interval
1
4
,
4
, that is, when
1
4
<
4
.
Then Equation (4.2) can be rewritten into the following simplified expression (Herbertsson, 2012):
() =
(1 ) (
)((
) (
1
))
4
=1
(
)
4
=1
(1 (
))
1
4
(4.4)
where () = [ ].
17
Hence, from Equation (4.4) we see that in order to calculate the CDS spread (), we need a model for
the default time , more specific, we need and explicit expression for the default distribution
() = [ ]. Thus, in the next section, we will therefore discuss one such model for , namely a so
called intensity based model.
4.2 Intensity based model
In this section we will study the so called intensity based model, which is needed when calculating the
CDS spread.
To start with, assume that we have a probability measure , and
)
>0
to be a d-dimensional stochastic process i.e.
=
,1
,
,2
,
,3
,
,
where is an integer and
,1
,
,2
,
,3
,
,
typically models different
kind of economic or financial factors. Hence, in the function :
() = (
)
>0
. Then one can define the random variable as
(Herbertsson, 2012):
= inf 0: (
)
1
0
. (4.4)
Hence, is the first time the increasing process (
0
reaches the random level
1
, see in the
Figure 4.4 (Herbertsson, 2012)
Figure 4.3: The construction of via
.
From Equation (4.4), we can further derive (Herbertsson, 2012)
18
[
> ] = exp (
0
(4.5)
and
() = (
)exp (
0
(4.6)
where
)
>0
. The intuitive meaning of this is as follow.
Consider a single obligor with default time , and we assume as discussed before that
is a stochastic
process and
[ [, +)|
> (4.7)
see also in Figure 4.4. Thus, the probability of having a default in the small time period [, +)
conditional on the information
. To be more specific,
is
denoted as the default intensity of , with respect to the information
(Herbertsson, 2012)
Figure 4.4: Given the information, will arrive in [, +) with probability
.
Finally, one can prove that the construction of in Equation (4.4) leads to the relationship in Equation
(4.7), that is
t
is the intensity of the random variable .
The Intensity () can for example be considered in three different cases, namely:
Intensity () can be a deterministic constant;
Intensity (
0
=
that is, is exponentially distributed with parameter .
When (
0
= ()
0
(4.8)
and
() = () ()
0
. (4.9)
Another important case of deterministic default intensity () is a so called piecewise constant default
intensity. Let
1
,
2
, ,
1
0 <
1
2
1
<
1
<
(4.10)
for some positive constants
1
,
2
, ,
0
(4.11)
20
Figure 4.5: The intensity () as a piecewise constant function.
which together with Equation (4.10) yields that
[ > ] =
0 <
1
1
1
(
1
)
2
1
<
1
=1
1
<
where we define
0
as
0
= 0.
An example of stochastic default intensity
is a Vasicek-process, given by
= (
) +
where
is a Brownian motion under risk-neutral measure P. Then we get that the survival probability
up to time is given by (Bjrk, 2009)
[ > ] = (() ()
0
)
for
() =
1
() = (() )
2
2
2
2
4
()
2
.
Moreover, since
() =
[>]
() is given by
21
() =
2
2
2
+
0
(
1)
2
2
2
()()
0
.
4.3 Valuation of an Interest Rate Swap
In this section we will explain the features of an interest rate swap, first in a general way followed by a
more theoretical manner.
An interest rate swap is an agreement between two parties, A and B, to exchange a fixed leg interest
rate cash flow stream for a floating equivalent under a given period (see Figure 4.6). One party pays the
fixed stream while the other pays the floating. The floating rate is based on an interest rate, for example
the LIBOR. Depending on movements in the underlying interest rate the value of the swap contract
changes with time. A common reason to enter into an interest rate swap contract is to manage risks
related to interest rates (Asgharian, et al., 2007).
An interest rate swap contract specifies the following properties:
Swap rate (the annual fixed interest rate cash flow stream)
The interest rate on which the floating rate is taken from
Maturity
Payment frequency
Notional amount
A simple illustration of the payments streams is displayed in Figure 4.6.
Figure 4.6: Structure of an interest rate swap.
In the following two subsections we are discussing the valuation of an interest rate swap. The interest
rate swap can be viewed as a portfolio of generalized Forward Rate Agreements (FRA). Hence to get a
proper understanding of the value of the Interest rate swap we start by introducing the FRA.
4.3.1 Forward Rate Agreement
In this section we will closely follow the setup and notation from (Brigo, et al., 2006).
22
Forward rates are interest rates that can be locked in today for an investment in a future time period.
The forward rate can be defined through a prototypical Forward Rate Agreement (FRA). An FRA is an
over-the-counter contract between two parties that determines the rate of interest to be paid or
received on an obligation beginning at a future start date. It is characterized by three time instants:
t - the time at which the rate is considered
1
- the expiry date
2
- the time of maturity
where t T
1
T
2
.
The holder of the FRA receives an interest-rate payment at time T
2
for the period between T
1
and T
2
. At
the maturity T
2
, a payment based on the fixed rate K
FRA
is exchanged against a floating payment based
on the spot rate L(T
1
, T
2
). To put it in a simple way, the contract allows one to lock in the interest rate
between time T
1
and time T
2
at a value of K
FRA
, for a contract with simply compounded rates. This
means that the expected cash flows must be discounted from T
2
to T
1
. At time T
2
one receives
(T
1
, T
2
)K
FRA
N units of cash and simultaneously pays the amount (T
1
, T
2
)L(T
1
, T
2
)N. Here N is the
contracts nominal value and (T
1
, T
2
) denotes the year fraction for the contract period [T
1
, T
2
]. Thus,
the value of the FRA, at time T
2
can, for the seller of the FRA (fixed rate receiver), be expressed as (Brigo,
et al., 2006)
N (T
1
, T
2
)(K
FRA
L(T
1
, T
2
)). (4.12)
Further, L(T
1
, T
2
) can also be written as
L(T
1
, T
2
) =
1 P(T
1
, T
2
)
(T
1
, T
2
)P(T
1
, T
2
)
(4.13)
and this enables us to rewrite Equation (4.12) as
N (T
1
, T
2
) K
FRA
1 P(T
1
, T
2
)
(T
1
, T
2
)P(T
1
, T
2
)
= N(T
1
, T
2
)K
FRA
1
P(T
1
, T
2
)
+1. (4.14)
To find the value of the FRA at time t, the cash flow exchanged in Equation (4.14) must be discounted
back to time t, that is, we want to compute the quantity
N P(t, T
2
) (T
1
, T
2
)K
FRA
1
P(T
1
, T
2
)
+1. (4.15)
To do this we first note that according to classical, no arbitrage interest rate theory, the implied forward
rate between time t and T
2
can be derived from two consecutive zero coupon bonds due to the equality
(Filipovic, 2009)
P(t, T
2
) = P(t, T
1
)P(T
1
, T
2
). (4.16)
23
This gives that P(t, T
1
) =
P(t,T
2
)
(T
1
,T
2
)
. Hence, by using Equation (4.16) we then get that
N P(t, T
2
) (T
1
, T
2
)K
FRA
1
P(T
1
, T
2
)
+1
= N[P(t, T
2
)(T
1
, T
2
)K
FRA
P(t, T
1
) +P(t, T
2
)].
Thus, we have that the value of the above FRA contract at time t is given by
FRA(t, T
1
, T
2
, (T
1
, T
2
), N, K
FRA
) = N[P(t, T
2
)(T
1
, T
2
)K
FRA
P(t, T
1
) +P(t, T
2
)]. (4.17)
Only one value of K
FRA
gives the FRA a value of 0 at time t. By solving for this value of K
FRA
, we get that
the appropriate FRA rate to use in the contract is the simply compounded forward interest rate
prevailing at time t for the expiry T
1
> t at maturity T
2
> T
1
, which is defined as
F
s
(t; T
1
, T
2
) =
P(t, T
1
) P(t, T
2
)
(t, T
2
)P(t, T
2
)
=
1
(T
1
, T
2
)
P(t, T
1
)
P(t, T
2
)
1. (4.18)
Rewriting the value of Equation (4.17) in terms of the simply compounded forward interest rate in
Equation (4.18) gives
FRA(t, T
1
, T
2
, (T
1
, T
2
), N, K
FRA
) = N P(t, T
2
)(T
1
, T
2
)K
FRA
F
s
(t; T
1
, T
2
). (4.19)
4.3.2 Interest rate swap
In this subsection we introduce the interest rate swap which is a generalization of the FRA. A
prototypical payer interest rate swap exchanges cash flows between two indexed legs, starting from a
future time. At every time point T
i
, within one, by the swap contract specified period T
+1
, T
, the
fixed leg pays N
i
K
IRS
, where K
IRS
is a fixed interest rate, N is the nominal value and
i
is the year
fraction between T
i1
and T
i
, (
i
= T
i
T
i1
). The floating leg pays N
i
L(T
i1
, T
i
) corresponding to
the interest rate L(T
i1
, T
i
) resetting at the preceding instant T
i1
for the maturity given by T
i
. For
simplicity, in this case, we are considering that the fixed-rate payments and floating-rate payments
occur at the same dates and with the same year fractions. Hence cash flows only take place at the date
of the coupons T
+1
, T
+2
, T
+3
T
. The individual who pays the fixed leg and receives the floating, B
in Figure 4.6, is the payer while the individual on the opposite side is termed the receiver, A in Figure 4.6.
The discounted payoff at time t < T
i=+1
N
i
(L(T
i1
, T
i
) K
IRS
)
and the discounted payoff at time t < T
i=+1
N
i
K
IRS
L(T
i1
, T
i
).
24
Seeing this last contract, from As side, as a portfolio of FRAs, every individual FRA can be valued using
the Formulas (4.17) and (4.19). This implies that the value of the interest rate swap
reciever
(t), is given
by (see also in Brigo, et al., 2006)
reciever
(t) = FRA(t, T
i1
, T
i
,
i
, N, K)
i=+1
reciever
(t) = N
i
P(t, T
i
)K
IRS
F
s
(t; T
i1
, T
i
)
i=+1
.
So using Equation (4.18) in the above expression implies that
reciever
(t) = N
i
K
IRS
P(t, T
i
)
i
P(t, T
i
)
(T
i1
, T
i
)
P(t, T
i1
)
P(t, T
i
)
1
i=+1
which can be simplified into
reciever
(t) = N
i
K
IRS
P(t, T
i
) P(t, T
i1
) P(t, T
i
)
i=+1
.
The sum above can be separated into two sums
N
i
K
IRS
P(t, T
i
)
i=+1
+ N P(t, T
i
) P(t, T
i1
)
i=+1
.
The second sum of the two, can be simplified
N P(t, T
i
) P(t, T
i1
)
i=+1
= N Pt, T
N P(t, T
).
This is because of when adding up the terms from = +1 to = all terms in the sum cancel out
apart from N Pt, T
andN P(t, T
reciever
(t) = N P(t, T
) +N Pt, T
+N
i
K
IRS
P(t, T
i
)
i=+1
. (4.20)
If we want to look at the value of the swap from the side of the payer instead of the receiver, then the
value is simply obtained by changing the sign of the cash flows
reciever
(t) =
payer
(t).
We can write the total value of the swap at t T
payer
(t) = NP(t, T
) Pt, T
K P(t, T
i
)
i=+1
. (4.21)
To clarify the interpretation of Formula (4.21), its features will now be further discussed. The interest
rate swap can be decomposed into two legs, a floating and a fixed. These two legs can be seen as two
fundamental prototypical contracts. The floating leg, N P(t, T
+K P(t, T
i
)
i=+1
in Formula (4.21), can be
seen as a coupon bearing bond. This gives that the Interest rate swap could be seen as an agreement for
exchanging the floating rate note (floating leg) for the coupon bearing bond (fixed leg).
A coupon bearing bond is a contract that guarantees a payment of a deterministic amount of cash at
future times T
+1
, T
+2
, T
+3
T
K
IRS
+N for i = .Here K
IRS
is the fixed interest rate and N is the nominal bond value.
Discounting the cash flows back to present time t from the payment times T
i
, the value of the coupon
bearing bond is (Brigo, et al., 2006)
NPt, T
+K P(t, T
i
)
i=+1
.
where N K P(t, T
i
)
i=+1
are the future discounted cash flows from the coupon payments, and
N Pt, T
+K
P(t, T
i
)
i=+1
is the second part of Formula (4.21).
The floating leg in the interest rate swap, N P(t, T
, T
+1
, T
+2
T
1.
.
Moreover, the note pays a cash flow at T
consisting of the reimbursement of the notional value. The
value of a floating rate note is obtained by changing the sign of the formula for the
reciever
(t) in
Equation (4.20) with a fixed leg of zero and then adding it to the present value of the cash flow paid at
time T
) N Pt, T
0 +N Pt, T
= N P(t, T
).
The reason for this convenient final formula is that the entire floating rate note can be replicated
through a self-financing portfolio. This shows that a floating rate note is always equal to N units of cash
at its reset dates, in our case every quarter of a year. In other words, the floating rate is always equal to
its notional amount when t=T
i
. One could express this as a floating rate note always trades at par
(Bjrk, 2009).
26
The forward swap rate K
,
(t) at time t for the sets of times T and year fraction is the rate in the fixed
leg of the interest rate swap in Formula (4.21) that makes the swap a fair contract at time t (Brigo, et al.,
2006). The formula is displayed below
K
,
(t) =
P(t, T
) Pt, T
i=+1
P(t, T
i
)
. (4.22)
In our calculations we assume that the interest rate swap contract is written at time = T
, this reduces
Equation (4.22) to
K
,
(t) =
1 Pt, T
i=+1
P(t, T
i
)
.
4.4 The CIR Model
In this thesis we will use a model for the short term interest rate given by the CIR model (Cox, et al.,
1985) on the term structure of interest rate.
The CIR model is developed to provide a general framework for determining the term structure of
interest rates. It can be used for the pricing of derivative products and risk free securities. The model has
its origin in Vasiceks model from 1977 but with the change that it introduces a square root term in the
diffusion coefficient on the instantaneous short rate dynamics. The CIR model have been highly
recognized and grown to be a benchmark within interest rate theory. Reasons for this is its analytical
compliance and, different from Vasiceks model, that it assumes the short term instantaneous interest
rates to always be positive, hence making it more handy to use. However the assumption of an always
positive short term instantaneous interest rates is not always true in the financial climate of today. The
CIR model specifies an instantaneous interest rate
= (
) +
where is a strictly positive parameter illustrating the speed of adjustment to the long term mean .
The term
is the standard deviation factor of the process and is often called the volatility. The
condition 2 >
2
puts a positive restriction on
, hence
(, ) denote the discounted value of the bilateral financial contract at time t from the investors
point of view. Note that
(, ) is the sum of all future discounted cash flows to the investor in the
contract over the period t to T. So due to the counterparty risk,
(, ) as (, ). Hence (, ) =
(, ) if there
would be no risk of default.
The Net Present Value (NPV) of the risk-free cash flow is defined as
(, ) = [(, )|
].
Here
is denoted as all the information available at time t. With the above definitions in mind now we
can define
(, ) as
(, ) = 1
{>}
(, ) +1
{}
(, ) +(, ) (, )
+
(, )
+
.
In the formula above we can see, if no default happens, (, ) =
(, )|
] is always non-negative, and the cash flows associated with the scenario where default is
an option are always smaller compared to the corresponding of the default free version (Brigo, 2008).
28
[
(, )|
] = [(, )|
] (1 )1
{}
(, )()
+
|
.
The CVA can be decomposed into three terms, the loss given default, the expected exposure and the
probability of default. Closer focus will be put on each of these terms in Section 5.4, 5.5 and 5.6.
5.2 CVA under Basel III
In this thesis we are calculating the CVA using the advanced formula for CVA which will be introduced in
the coming subsections. Hence, what we are calculating is one CVA value which later can be used as an
input in a VaR-model to calculate the CVA-capital charge defined in Basel III. Nevertheless, when
calculating the CVA-capital charge which is a new feature to Basel III the CVA calculations must be based
on the formulas we use and introduce below in the Section 5.2.3. Hence as an important addition we
will in this section explain the CVA-capital charge which our results can be used to calculate. Also we will
introduce the terms expected exposure, probability of default and loss given default which is included in
our CVA calculations.
5.2.1 CVA Capital Charge
Due to worsening in the counterpartys credit quality, the CVA capital charge was added to the CCR
capital charge to measure the risk of Mark-to-Market (MtM) losses. There exist two different methods
to calculate the CVA capital charge, the standard and the advanced, which one to use depends on what
method a bank is approved for in calculating capital charge for counterparty default risk and certain
interest rate risk. Banks that use the advanced method to calculate CVA capital charge need to have
IMM approval for counterparty credit risk as well as the approval to use the market risk internal models
approach for the specific interest-rate risk of bonds (Basel Committee on Banking Supervision, 2011).
5.2.2 Standard Model
For the banks that do not have the IMM approval for counterparty credit risk, they must calculate its
CVA capital charge using the following formula stated in the Section 104 of Basel III (Basel Committee on
Banking Supervision, 2011)
= 2.33
where
2
= 0.5
=1
2
+0.75
=1
where
29
is the one-year risk horizon (in units of a year), h=1
is the exposure at default of counterparty i who does not granted the approval for
IMM, and the non-IMM banks the exposure should be discounted by using the factor
[1 (0.05
)]/(0.05
).
is the full notional of one or more index CDS used for hedge the CVA risk, and should be
discounted by using the factor [1 (0.05
)]/(0.05
).
is the weight applicable to index hedges. The bank must map indices according to one of
the weights
(the
and
needs to be
add together when there are several positions).
is the maturity of the index hedge ind, which is the notional weight average maturity in
case of several hedge position.
5.2.3 Advanced method
In order to calculate a VaR on CVA, banks uses their own specific interest rate risk VaR model for bonds
by modeling changes in the CDS-spreads of counterparties. The VaR is calculated on the aggregated
CVAs of the banks OTC derivatives. The model will not measure the sensitivity of CVA to changes in
other market factors because of the restriction for the model to changes in the counterparties credit
spread. The CVA capital charge is composed by the sum of a stressed and a non-stressed VaR
component, and due to its property, the calibration of the expected exposure is normally done using the
credit spread calibration of the worst one-year-period contained in the three-year-period.
In Basel III it reads that, the CVA for CVA capital charge calculations must for every counterparty be
based on the formula below, regardless of which accounting valuation method a bank uses to determine
CVA. This is also the formula we use for our CVA calculations. The following formula is taken from
Section 98 of Basel III (Basel Committee on Banking Supervision, 2011)
= (
1
+
=1
0;
where,
is the loss given default of the counterparty, based on the spread of a market
instrument of the counterparty, which must be assessed instead of an internal estimate.
30
, and
0
= 1.
5.3 Exposure
Counterparty Credit Exposure (CCE), or simply exposure, is the amount a company could lose if its
counterparty defaults (Cesari, 2012). The Mark-to-Market (MtM) value of the OTC contract, which states
the current market value of an asset, is the value we need to use to calculate the counterparty credit
exposure. The MtM value can be either positive or negative; hence there is an asymmetry of potential
losses. For example, suppose the MtM value is negative at default, then the company will owe its
counterparty; while if it is the opposite case when the MtM value is positive, the counterparty will be
unable to make future commitments and lose the amount of the MtM value. With the property that the
company loses if the MtM is positive and gain nothing when the MtM is negative, we can define the
expression of the exposure as:
= max(0,
)
where
is the MtM value of the contract at time t and represents the different scenarios. In Figure
5.1, we display 10 simulations of exposure for interest rate swaps with the length of 10 years.
Figure 5.1: Ten simulations of Exposure for interest rate swaps.
31
5.3.1 Quantitative measure of Exposure
There are several ways to quantify exposure and according to the definition by the Basel Committee on
Banking Supervision, three of them are widely used:
Expected Exposure (EE): it is the expectation of loss given default based on a zero recovery rate.
In other word, it is the average exposure after considering different scenarios. Thus, it can be
calculated using the formula
=
1
max0,
=1
where is the number of simulated interest rate paths, and
)].
Potential Future Exposure (PFE): the definition for PFE is very much similar to Value at Risk (VaR),
which is the worst possible exposure with a certain confidence interval (usually 99%). What is
different between the two is that the PFE has a longer time horizon, usually years; while the
VaRs time horizon is accounted in days
Expected Positive Exposure (EPE): is slightly different to the previous two concepts since it is
defined to be the time average of the
(0,
=1
where t is the time now and M is the maturity. What needs to be mention is that the
. (5.1)
This term comes from the approximation (see also in Herbertsson 2012)
()
()
.
Thus
[ > ] =
()
and this implies that
[
1
<
] = [
] [
1
] = [ >
1
] [ >
]
=
()
()
where
1
<
. The max in Equation (5.1) guarantees that the value of PD is non-negative, since in
this formula () may be a function of
.
The term
is the CDS-spread. In our thesis we have created three different CDS-spreads by simulating
three different risk scenarios; High; medium and Low. This structure is inspired by Brigos article in
2008.The complete table with the CDS-spreads will be shown in Section 7.
33
6. Calculating Expected Exposure
In this section we explain the steps taken to calculate our CVA values for an interest rate swap in a CIR-
framework.
6.1 The interest rate path
To model CVA for an interest rate swap in a CIR-framework we need to model the expected exposure.
To get a value for the expected exposure in the swap contract we start by simulating interest rate paths
following a CIR-process.
For the simulation we take support in an algorithm described in (Broadie, 2006). If the interest rate
follows a CIR-process including the parameters , , and
, the value of
given
, where < , is
given by multiplying a scaling factor with a non-central chi-squared distributed random variable.
2
() (6.1)
In Equation (6.1), denotes the non-centrality chi-squared parameter, d is the degrees of freedom and
is a scaling factor. The parameters are dependent on the distance between s and t, , , and the
initial value of the short term rate
0
.
=
4
(t)
2
(1
(t)
)
=
4
=
2
1
(t)
4
.
The CIR-model is estimated with the restriction that 2 >
2
, and if this holds, then it means that d will
always be bigger than one.
When d >1 a non-central chi-squared random variable can be written as a sum of a non-central chi-
squared random variable with one degree of freedom and a normal chi-squared random variable with d-
1 degrees of freedom (Broadie, 2006). So for this case, with an estimation that is restricted subject to
d >1, the interest rate simulation for every period is modeled as follow:
1
2
() +
1
2
where
1
2
() = +
2
and Z is a standard normal distributed random variable.
34
Figure 6.1: CIR-process path simulation with, = 0.1, = 0.03, = 0.1,
= 0.02.
Figure 6.1 show 10 simulated interest rate paths over a period of ten years when the interest rate
follows a CIR process with the parameters = 0.1, = 0.03, = 0.1,
0
= 0.02.
6.2 Bond Price
Subject to assumptions of continuous trading and absence of transaction costs, it is possible to model
the arbitrage-free price of a default-free zero coupon bond (, ) in the following way (Bjrk, 2009)
(, ) = (, )
(,)r
(, ) =
()
(
()
1) +
(, ) =
()
1
(
()
1) +
where
=
2
+2
2
, =
+
2
, =
2
2
.
35
Figure 6.2: Simulation of Bond Prices.
In Figure 6.2 we show the bond price (, ) as a function of time t, given the realization of the interest
rate
as in Figure 6.1.
6.3 The Swap Contract
The value of the swap contract
payer
(t) is calculated at close time points up until maturity. It is
calculated with the formula below
payer
(t) = NP(t, T
) Pt, T
K P(t, T
i
)
i=+1
where the calculation is based on the bond prices P(t, T) from the previous subsection. In this
subsection we will let = (), i.e. the index will change as t runs from 0 to T
()1
<
()
.
To ease the interpretation of the formula, the
payer
(t) can be decomposed into its individual terms.
The factor N is the notional amount, hence the scaling factor for the interest rate swap contract. The
derivation of the term N P(t, T
) =
N, when = () and
()1
<
()
. This will be explained in more detail below.
36
The term NPt, T
K P(t, T
i
)
i=+1
is all the future discounted cash flows at time t from the
fixed leg payment stream. Here N K P(t, T
i
)
i=+1
are the future discounted cash flows from the
coupon payments, this sum decreases in size as the contract approaches maturity. The term N Pt, T
is the discounted reimbursement of the notional value in the fixed leg, this goes to N as the contract
approaches maturity. Here K is the swap rate that makes the contract fair at time t = 0, or in other
words the swap rate that makes the present value of the floating and fixed leg equal at time t = 0.
Finally, is the year fraction between T
i1
and T
i
(Brigo, et al., 2006) (Bjrk, 2009).
The formula above is used to calculate the value of the interest rate swap at time t from t = 0 to
t = T
with close steps. The swap rate K is calculated once, at t = 0. The formula above calculates the
value of an interest rate swap at time t for t T
i=+1
= 0, N P(t, T
) = N and N Pt, T
=
1
max0,
=1
For this we use the method explained in Section 5. It is displayed in Figure 6.4.
38
Figure 6.4: the Expected Exposure.
Now we have calculated all the terms needed for our CVA value.
39
7. Results
Our CVA is calculated using the following formula with a notional amount of 1 and a value of the other
terms as specified later in this section.
= (
1
+
=1
0;
For the CDS we have created CDS spreads by simulating 3 different scenarios; High; medium and Low.
This setup is inspired by Brigos article in 2008 (Brigo, 2008).
Table 7.1: Different CDS spreads representing different risk scenarios.
The discount rate
.
Instead of using a flat interest rate for the discount term
follows the simulated Interest rate. Figure 7.3 illustrate a comparison between the
results.
42
Figure 7.3: Comparison of CVA values as a function of using different discount rate scenarios, where the upper ones
discount rate is dependent on the CIR simulation.
43
Figure 7.4: Comparison of CVA as a function of different discount rates.
From Figure 7.4, it can be concluded that the two different discount scenarios seem to give relatively
similar CVA values.
8. Discussion and conclusion
By looking at the model risk for CVA on an interest rate swap in the CIR-framework and by investigating
how sensitive the CVA-value is with respect to the underlying parameters, we can conclude that the
level of the underlying parameters has a big impact on the final CVA value.
The estimated value of the EE will change with respect to , , ,
0
, driving the simulated interest rate
path. The parameters that affect the CVA value the most are and . That the CVA calculation is the
most sensitive to changes in and is fairly intuitive considering the meaning of the included
parameters. The long term mean of the process will have a great impact on the level of the CIR-
process, in our case the interest rate. Also considering that is the volatility of the process, an increase
in this parameter will lead to an increased uncertainty about the future interest rate, hence, a bigger risk
will be associated with a high . The level of
0
will measured over a longer period not have as big of an
impact on the interest rate. Further, if is high it means that the process quickly will return to its mean,
however the speed of the return to the long term mean does clearly not have a dramatic impact on the
CVA.
Even if the CVA value is small and a small change of its value doesnt seem too intimidating, it could
have big effects. Since we have done our estimations with a notional of one, the real CVA must be scaled
with the proper amount of the notional. It is also discovered by our research that the different CDS
levels have a great impact on the final level of the CVA.
44
Comparing our work to previous research from Brigo (2008) and Douglas (2012), we can conclude that
the results from our CVA calculations give similar CVA values as these papers. With similar input
variables our output matches the results of previous efforts to calculate CVA for an interest rate swap.
They as us, emphasize the importance of being aware of the volatility in the CVA and account for it in a
sound manner. The scientific contribution of our work is the finding that parameters , , ,
0
have an
impact on the CVA for an interest rate swap in a CIR framework.
Lastly, since there is significant insecurity in the modeling of every term in the CVA model, it is very
difficult to say with confidence that the estimated CVA value calculated in our model is accurate.
Nevertheless if the CVA is not taken seriously, there will be a great chance that what happened in the
financial crisis of 2008-2009 might happen again. During that crisis losses could be traced to volatility in
CVA. So now when the importance of accurately measuring the CVA has been recognized, we believe
that the research within this field will grow further. To conclude, it is evident that the final CVA value is
sensitive to changes in the underlying parameters , , ,
0
as well as to variations in the other terms
included in the CVA model, which indicates an evident model risk.
45
Bibliography
Asgharian, Hossein and Norden, Lars. 2007. Rntebrande Instrument, Studentlitteratur 1st edition. s.l. :
Studentlitteratur AB, 2007.
Bank for International. 2009. History of the Basel Committee and its Membership. Bank for
International Settlements. [Online] August 2009. http://www.bis.org/bcbs/history.pdf.
Bank for International Settlements. 2011. Capital treatment for bilateral counterparty credit risk
finalised by the Basel Committee. Bank for International Settlements. [Online] June 1, 2011.
http://www.bis.org/press/p110601.htm.
Bank For International Settlements. 2005. The Application of Basel II to Trading Activiteis and the
Treatment of Double Default Effects. Bank for International Settlements. [Online] Bank For International
Settlements, April 2005. http://www.bis.org/publ/bcbs111.htm.
Basel Committee on Banking Supervision. 2011. Basel III: A global regulatory framework for more
resilient banks and banking systems. Bank for International Settlements. [Online] June 2011.
http://www.bis.org/publ/bcbs189.pdf.
Bjrk, Tomas. 2009. Arbitrage Theory in Continuous Time, 3rd edition. New York : Oxford University
Press, 2009.
Bodie, Zvi, Kane, Alex and Marcus, Alan. 2012. Essentials of Investments, 9th edition. s.l. : Richard D.
Irwin, Inc., 2012.
Brigo, Damiano and Chourdakis, Kyriakos. 2009. Counterparty Risk For Credit Default Swaps: Impact of
Spread Volatility And Default Correlation. International Journal of Theoretical and Applied Finance. 2009,
Vol. 12, 07, pp. 1007-1026.
Brigo, Damiano and Mercurio, Fabio. 2006. Interest Rate Models - Theory and Practice, 2nd edition.
Berlin Heidelberg : Springer-Verlag, 2006.
Brigo, Damiano. 2008. COUNTERPARTY RISK FOR CREDIT DEFAULT SWAPS. International Journal of
Theoretical and Applied Finance. 2, 2008, Vol. 12, 07, pp. 1007-1026.
Broadie, Mark. 2006. Exact Simulation of Stochastic Volatility and Other Affine Jump Diffusion Processes.
Operations Research. 2, 2006, Vol. 54, 2, pp. 217-231.
Calabresi, Massimo. 2009. Banks Balk at Selling Toxic Assets. TIME Magazine. [Online] TIME Magazine,
April 16, 2009. http://www.time.com/time/business/article/0,8599,1891614,00.html.
Caruana, Jaime. 2011. Regulatory reform: remaining challenges. Bank for International Settlements.
[Online] Bank for International Settlements, July 7, 2011. http://www.bis.org/speeches/sp110711.htm.
46
Cesari, Giovanni. 2012. Modelling, Pricing, and Hedging Counterparty Credit Exposure, 1st edition. New
York : Springer-Verlag, 2012.
Cox, C. John, Ingersoll, E. Jonathan and Ross, A. Stephen. 1985. A Theory of the Term Structure of
Interest Rates. Econometrica. 1, 1985, Vol. 53, 2, pp. 385-408.
Douglas, Rohan. 2012. Compareing Alternate Methods for Calculating CVA Capital Charges under Basel
III. s.l. : Global Association of Risk Professionals, 2012.
Engelmann, Bernd. 2011. The Basel II Risk Parameters, 2nd edition. New York : Springer-Verlag, 2011.
Filipovic, Damir. 2009. Term - Structure Models, A Graduate Course. Berlin Heidelberg : Springer-Verlag,
2009.
Fusar, Peter C. 2008. The Professional Risk Managers' Guide, 7th edition. s.l. : McGraw-Hill, 2008.
Herbertsson, Alexander. 2012. Credit Risk Modelling Lecture Notes. Gothenburg : University of
Gothenburg, 2012.
Hull, John. 2012. Options, Futures, and Other Derivatives, 8th edition. s.l. : Prentice Hall, Inc., 2012.
International Accounting Standards Board. 2009. International Accounting Standards. IFRS. [Online]
January 2009. http://www.ifrs.org/Meetings/Documents/IASBJan09/Derec0901b02Cobs.pdf.
Kozul, Natasha. 2011. Mastering Investment Banking Securities. s.l. : Financial Times/ Prentice Hall, 2011.
Moles, Peter and Terry, Nicholas. 2012. The Handbook of International Financial Terms. s.l. : Oxford
University Press, 2012.
Saunders, Anthony. 2012. Financial Markets and Institutions, 5th edition. s.l. : McGraw-Hill Publishing
Company, 2012.
Schnbucher, Philipp J. 2003. Credit Derivatives Pricing Models: Models, Pricing and Implementation.
s.l. : Wiley,Chichester, 2003.
Schuermann, Til. 2004. What Do We Know About Loss Given Default? s.l. : Wharton Financial Institutions
Center, 2004.
Standard & Poor's. 2009. Understanding Standard & Poor's Rating Definitions. s.l. : Standard & Poor's,
2009.
47
Appendix
48