DIGITAL OPTION PRICING BASED ON COPULAS WITH STOCHASTIC SIMULATION

The Pure and Applied Mathematics.
2015.
Aug,
22(3):
299-313

- Received : July 14, 2015
- Accepted : August 14, 2015
- Published : August 31, 2015

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

In this paper, we show the effectiveness of copulas by comparing the correlation of market data of year 2010 with those of years 2006-2009 and investigate copula functions as pricing methods of digital and rainbow options through real market data. We propose an accurate method of pricing rainbow options by using the correlation coefficients obtained from the copula functions depending on strike prices between assetes instead of simple traditional correlation coefficients.
Definition 2.1
(
[6]
). A two-dimensional subcopula is a function
C'
with the following properties:
Note that for every (
u
,
v
) in
DomC'
, 0 ≤
C'
(
u
,
v
) ≤ 1 so that
Ran
C'
is also a subset of
I
(i.e. [0, 1]).
Definition 2.2
(
[6]
). A two-dimensional copula is a function
C
whose domain
I
^{2}
with the following properties:
Lemma 2.3
(
[6]
).
Let H be a joint distribution with margins F and G, then there is a unique subcopula
C
'
such that
Lemma 2.4
(
[6]
).
Let
C
'
be a subcopula. Then exists a copula C s.t
Theorem 2.5
(
[6]
). (Sklar’s Theorem)
Let H be a joint distribution function with margins F and G then there exists a copula C such that
If
F
and
G
are continuous, then
C
is unique. That is,
C
is uniquely determined on
Ran F
×
Ran G
. Conversly, if
C
is a copula and
F
and
G
are distribution functions, then
H
is a joint distribution function.
Corollary 2.6
(
[6]
).
Let
Then
x
is given by the
survival function
[6]
The
joint survival function
for the joint distribution function
H
is given by
Suppose that the copula of
X
and
Y
is
C
. Then
with defining a function
(
survival copula
) from
I
^{2}
to
I
by
Note that
where
F
^{ −1}
denotes the inverse cumulative distribution function and
u
is an uniform variable defined over (0, 1).
Since
can be written as below:
From the fact that
F_{X}
(
x
) =
P
[
X
≤
x
] and
, we have
and
Hence the upper tail dependence can be rewritten as:
Lower tail dependence is symmetrically defined:
Therefore we obtain the lower tail dependence as below:
Definition 2.7
(
[6]
). Let
φ
be a continuous, strictly decreasing function from
I
to [0, ∞) such that
φ
(1) = 0. We define that
φ
^{[−1]}
: the pseudo-inverse of
φ
with Dom
φ
^{[−1]}
= [0, ∞) and Ran
φ
^{[−1]}
=
I
where
Note that
φ
^{[−1]}
is continuous and nonincreasing on [0, ∞], and strictly decreasing on [0,
φ
(0)]. Furthermore,
φ
^{[−1]}
(
φ
(
u
)) =
u
on
I
, and
Finally, if
φ
(0) = ∞, then
φ
^{[−1]}
=
φ
^{−1}
.
Lemma 2.8
(
[6]
).
Let φ be a continuous, strictly decreasing function from I to
[0, ∞]
such that
φ
(1) = 0,
and let
φ
^{[−1]}
be the pseudo-inverse of φ
.
Let C be the function from I
^{2}
to I given by
Then C satisfies the boundary conditions for a copula
.
i.e., C
(
u
, 1) =
u
,
C
(1,
v
) =
v
and
C
(
u
, 0) = 0 =
C
(0,
v
).
2.5.2.
Important Examples of Archimedean copulas Let φ be a continuous, strictly decreasing function from I to
[0, ∞]
such that φ
(1) = 0,
and let φ
^{[−1]}
be the pseudoinverse of φ
[6]
.
Let
C
be the function from
I
^{2}
to
I
given by
Examples of bivariate Archimedean copulas are as the followings
[5]
:
A
and
B
has same exercise date
T
and strike prices
K
_{A}
and
K
_{B}
respectively, then the price of this digital put option[
DP
] in a complete market is
where
Q
(
K
_{A}
,
K
_{B}
) is the joint risk-neutral probability that both A and B market indexes are below the corresponding strike prices. In order to compare the price of our bivariate product with that of the univariate ones, it would be great if we could write the price as
Consider a bivariate digital call option[
DC
]. Differently from the digital put option, it pays one unit of currency if both A and B indexes are above the strike levels
K
_{A}
and
K
_{B}
. The relevant probability in this case is
Like
DP
, we can represent this as
where
is a survival copula. Using survival copulas, it can be rewritten as belows
For applying copula functions, we can use Archimedean copulas to obtain
where
φ
is the generator function of Archimedean copula functions. Finally we may price our bivariate claim using
Algorithm 3.1
(
[4]
). :
For each Archimedean copula we need the follwongs with Kendall’s
τ
:
In Gumbel copula, we have a generator function
φ
(
t
)=(− ln
t
)
^{θ}
:
Therefore, we can obtain
u
and
v
as below:
In Clayton copula, with a generator function
φ
(
t
) =
t
^{−θ}
− 1:
Similarly, we have
3.2. Digital Option Pricing
Consider two univariate digital put options with the maturity
T
and strike prices
K
_{1}
and
K
_{2}
respectively. If we denote by
X
_{T}
and
Y
_{T}
the prices of the underlying assets at maturity then we have
The price of the bivariate digital put option can be evaluated by the following formula.
According to Sklar’s theorem, we have
Assume that asset value random processes X and Y have lognormal distributions under a constant interest rate
r
and constant volatilities
σ_{X}
and
σ_{Y}
respectively as in the Black-Scholes model
[1]
, then we have the indivisual prices of Put Options,
and Call Options,
where
and
with
X
_{0}
and
Y
_{0}
, are initial values of underlying assets. Hence the value of digital put option can be written as
and the value of digital call option should be
Consequently, we show the price of digital put option using Clayton copula among Archimedean copula functions as the following.
Similarly, we have for the call option value as below,
where
θ
is a Clayton copula parameter and can be estimated using Kendall’s
τ
[2]
.
X
and
Y
at maturity
T
over strike price
K
. Rainbow call option price under the risk-neutral measure is shown as below
[9]
.
where
with
ρ
which is the correlation coefficient between two assets and
with two asset volatilities
σ_{X}
and
σ_{Y}
respectively.
It can be shown that the rainbow call option value turns out to be equal to the value of an European call option if two assets are definitely same i.e
X
≡
Y
.
Proposition 3.1.
The rainbow call option formula becomes the formula for European call option values if two underlying assets are definitely same
.
Proof
. First, with
X
=
Y
, we have
since
σ_{X}
=
σ_{Y}
. This shows
as
ρ
→ 1. Now from the first part of three terms of the rainbow formula,
where
as
ρ
→ 1 since
σ
goes to 0 as
ρ
→ 1.
From the second part, similarly
Hence the sum of first two terms becomes
Now we invetigate the last part of the formula.
This completes the proof. ☐
The correlation between two assets is important to evaluate rainbow options especially depending on strike prices and it can be estimated from the copula functions established earlier on the relationship between assets. We propose an accurate method of pricing rainbow options by using the correlation coefficients obtained from the copula functions depending on strike prices between asstes instead of simple traditional correlation coefficients.
Distributions of Samsung Electronics & Rftech
For the sake of choosing apt copula functions, the calculated values of upper and lower tail dependences which determine the sensitivity of heavy tail are given by the formulae in Section 2.4 as 0.115 and 0.991 respectively for Samsung-Rftech and 0.182 and 0.992 respectively for Samsung-Displaytech. Therefore, we can use rather Clayton copula for adequate relationship because of the lower tail dependence value which is not close to zero. The result after performing simulations with Gumbel and Clayton copulas for Samsung-Rftech is shown in
Figure 2
.
Gumbel(Left) & Clayton(Right) for Samsung-Rftech
The correlation coefficients between two companies found by copula functions can be applied to expect the behavior of financial stability of firms. Dynamic fluctuation in curve of correlation coefficients between firms represents the active economic behavior while damping to less fluctuatedly static curve in correlation coefficinets leads us to expect the financial problem of one of companies compared. This is another advantage in applications of copula functions.
In
Figure 3
, we see the obvious result explained as the above. First two correlation coefficient curves each row include the 3-4 year coefficient curve for the first figure and the last year curve of 2009 or 2010 for the second. The third figure each row shows the dynamics of scaled stock prices of companies and it is worthy of notice that the company Sanyang was defaulted in 2009.
Samsung-Displaytech(1st row) & Samsung-Sanyang(2nd row)
Copula parameters between companies
As well as parameters for correlation between firms except almost independent relationship LG-CJ used to copula models we investigate, it is not difficult to find the volatility of each company for option pricing using real log return data. For example, Samsung Electronics’ volatility is
σ
_{ss}
= 0.0215614 with starting stock price
S
_{0}
= 659000 and Rftech has
σ
_{rf}
= 0.0328137 as a volatility with
S
_{0}
= 7080. Now, the digital put option value for Samsung-Rftech is calculated using copula functions like Gumbel and Clayton, and
Figure 4
shows option values depending on strike prices of Samsung and Rftech under Clayton copula function.
Digital Put Option Price of Samsung and Rftech
Conversely, the price of rainbow options can be estimated from the correlation coefficients depending on strike prices of these two firms which was calculated using proper copula functions. This is another accurate method to evaluate prices of rainbow options not using simple traditional correlation coefficients.

1. INTRODUCTION

Copula functions represent a methodology which has recently become the most significant new tool to handle in a flexible way the co-movement between markets, risk factors and other relevant variables studied in finance. While the tool is borrowed from the theory of statistics, it has been gathering more and more popularity both among academics and practitioners in the field of finance principally because of the intimate increase of volatility and erratic behavior of financial markets
[2]
. These new developments have caused standard tools of financial mathematics. The need to reach effective diversification has led to new investment products, bound to exploit the credit risk features of the assets. It is particularly for the evaluation of these new products, such as securitized assets and basket credit derivatives that the need to account for co-movement among non-normally distributed variables has become an unavoidable task.
Copula functions have been applied to the solution of these problems in mathematical finance. In fact, the use of copula functions enables the task of specifying the marginal distribution to be decoupled from the dependence structure of variables. For instance, the fact that the forward-futures spread is more significant for long maturity is shown by evaluating the spread based on copula models
[8]
. For applying copula functions, we have to settle three problems which are 1)what is more proper in Gaussian or Archimedean copula functions and 2)how to determine marginal distributions about real data, finally 3)how to select kinds of Archimedean copula functions
[4]
. In this paper, they are used in numerical examples that illustrate the use of the copula functions in the digital and rainbow option pricing.
This paper is organized as follows. Section 2 provides an overview of copula functions. In Section 3, copula algorithm method as well as digital option pricing and rainbow option pricing is established and shown for numerical experiments. In Section 4, stochastic simulations are performed and investigated through copula models. The conclusion of this work is in Section 5.
2. PRELIMINARIES

The analysis of copula functions started from the approach suggested by Nelsen. Cherubibi, Luciano, and Vecchiato applied copula functions to finance fields. In this section, we introduce the fundamental concept of copula and some knowledge needed.
- 2.1. Copula

- 1.Dom C'=S1×S2whereS1andS2are subsets ofIcontaining 0 and 1.
- 2.C'is grounded and 2-increasing.
- 3. For everyuinS1and everyvinS2,

PPT Slide

Lager Image

- 1. For everyu,vinI

PPT Slide

Lager Image

- 2. For everyu1,u2,v1,v2inIsuch thatu1≤u2andv1≤v2

PPT Slide

Lager Image

- 2.2. Sklar’s Theorem

- 1.DomC' =RanF×Ran G
- 2.H(x, y) =C'(F(x),G(y))for all x, y.

PPT Slide

Lager Image

PPT Slide

Lager Image

- H, F, G, andC'be joint distribution functions,
- F(−1)andG(−1)be quasi-inverses of F and G respectively.

PPT Slide

Lager Image

- 2.3. Survival Copulas

In many applications, the random variables of interest represent the lifetimes of individuals or objects in some population. The probability of an individual living or surviving beyond time
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.4. Upper and Lower tail dependence

Tail dependence can be measured and relates the amount of dependence in the upper right quadrant tail or lower left one of a bivariate distribution
[6]
.
Upper tail dependence is defined as:
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.5. Copula Functions

- 2.5.1.Archimedean copula

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.6. Digital Option

In finance, a digital[binary] option is a type of the option in which the payoff is either some fixed amount of some asset or nothing at all
[3]
. To give an intuitive grasp of the use of copula functions in finance, consider a very simple product, a bivariate digital option. This option pays one unit of currency if two stocks or indexes are above or below a pair of strike price levels. Options like these are very often used in structured finance, particularly index-linked products
[2]
.
For example, assume the option of companies
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

3. METHOD ALGORITHM

- 3.1. Copula Algorithm

- 1.Simulate two independent U(0, 1)random variates s and q.
- 2.Set,whereKcopulais the distribution function C(u,v).
- 3.Set u=φ−1(sφ(t))and v=φ−1((1 −s)φ(t))

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 3.3. Rainbow Option Pricing

A rainbow option is a derivative whose value is depending on two or more underlying securities or events. Rainbow options are usually calls or puts on the best or worst of underlying asstes, or options which pay the best or worst of assets. The number of assets underlying the option is called the number of colors of the rainbow. We consider a rainbow call option corrected formula based on the minimum value of two underlying assets
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

4. NUMERICAL RESULTS

- 4.1. Information of Stock Price

From 2005, cellular phone industry is not only for simple communication but also seeking for a wide range of development such as game, videotelephony, up-to-date smartphone. These days, one of top companies in cellular phone industry is Samsung Electronics. Therefore, it can be selected with relevant companies which are cooperative firms of Samsung Electronics. These are Abico, Displaytech, Sanyang, Rftech, Mtechvision. For the purpose of comparing non-subsidiary company or from different industry, LG Electronics, SK Hynix and CJ Corporation are investigated.
We consider the relationship between them and analyze how effect volatilities of stock prices make to pricing of financial derivatives with the relationship like correlation coefficients between them. We choose total 5 years(from 2006 to 2010) including around 1200 discrete data each company and compare data from the first 4 years with those of the last year. The volatility of stock prices during 2006-2009 analyzed and estimated correlations between them are applied to investigate the financial stability condition of companies in 2010.
- 4.2. Distribution and Correlation

We use the estimated distributions of selected companies from the real log-return data to find the parameters of relationship between them. The distributions of Samsung Electronics & Rftech obtained are approximately normal as shown in
Figure 1
.
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 4.3. Option Pricing

From real log-return data, we obtain estimated parameters for copula functions to establish copula models as
Table 1
.
Copula parameters between companies

PPT Slide

Lager Image

PPT Slide

Lager Image

5. Conclusion

Upper and lower tail dependences which determine the sensitivity of heavy tail give the right choice for proper copula functions to show the adequate relationship between companies. We can select Clayton, Gumbel, or Frank copula for better representation of correlation between firms. We show the effectiveness of copulas by comparing the correlation of market data of year 2010 with those of years 2006-2009 and investigate copula functions as pricing methods of digital options through real market data. The correlation coefficients between two companies found by copula functions can be applied to expect the behavior of financial stability of firms. Dynamic fluctuation in curve of correlation coefficients between firms represents the active economic behavior while damping to less fluctuatedly static curve in correlation coefficinets leads us to expect the financial problem of one of companies compared. This is an important advantage in applications of copula functions. Stochastic simulations for digital options are performed under the condition with adequately chosen copula functions and estimated parameters from real market data. Moreover, the accurate method of pricing rainbow options is proposed using the correlation coefficients obtained from the copula functions depending on strike prices between assetes instead of simple traditional correlation coefficients. This application can be extended to evaluate credit risk features and credit derivatives.
Acknowledgements

The corresponding author was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education, Science and Technology(2010-0010694).

Black F.
,
Scholes M.
1973
The Pricing of Options and Corporate Liabilities
Journal of Political Economics
81
637 -
654
** DOI : 10.1086/260062**

Cherubini U.
,
Luciano E.
,
Vecchiato W.
2004
Copula Methods in Finance
John Wiley & Sons

Hull John C.
2009
Options, Futures, and Other Derivatives
Seventh edition
Prentice Hall

Melchiori Mario R.
2003
Which Archimedean Copula is the right one?
The YieldCurve.com e-Journal

Naimy Viviane Y.
2012
Gaussian Copula vs Loans Loss Assessment: A Simplified and EasyTo-Use Model
Journal of Business Case Studies
8
533 -
542

Nelsen Roger B.
2006
An Introduction to Copulas
Second Edition
Springer Verlag
New York

Philipp Arbenz
2013
Bayesian Copulae Distributions, with Application to Operational Risk Management-Some Comments
Methodology and Computing in Applied Probability
15
(1)
105 -
108
** DOI : 10.1007/s11009-011-9224-0**

Pu Yuqi
,
Kim Seki
2014
Pricing Forward-Future Spread based on Copulas with Stochastic Simulation
The Pure and Applied Mathematics
21
77 -
93
** DOI : 10.7468/jksmeb.2014.21.1.77**

Stulz Rene M.
Option on the Minimum or Maximum of Two Risky Assets
Journal of Financial Economics
10
161 -
185
** DOI : 10.1016/0304-405X(82)90011-3**

Citing 'DIGITAL OPTION PRICING BASED ON COPULAS WITH STOCHASTIC SIMULATION
'

@article{ SHGHCX_2015_v22n3_299}
,title={DIGITAL OPTION PRICING BASED ON COPULAS WITH STOCHASTIC SIMULATION}
,volume={3}
, url={http://dx.doi.org/10.7468/jksmeb.2015.22.3.299}, DOI={10.7468/jksmeb.2015.22.3.299}
, number= {3}
, journal={The Pure and Applied Mathematics}
, publisher={Korean Society of Mathematical Education}
, author={KIM, M. S.
and
KIM, SEKI}
, year={2015}
, month={Aug}