PRICING FORWARD-FUTURES SPREAD BASED ON COPULAS WITH STOCHASTIC SIMULATION
PRICING FORWARD-FUTURES SPREAD BASED ON COPULAS WITH STOCHASTIC SIMULATION
The Pure and Applied Mathematics. 2014. Feb, 21(1): 77-93
• Received : January 14, 2014
• Accepted : February 17, 2014
• Published : February 28, 2014 PDF e-PUB PubReader PPT Export by style
Article
Author
Metrics
Cited by
TagCloud
YUQI, PU
DEPARTMENT OF MATHEMATICS, SUNGKYUNKWAN UNIVERSITY, SUWON, GYEONGGI 440-746, KOREAEmail address:magie609@gmail.com
SEKI, KIM
DEPARTMENT OF MATHEMATICS, SUNGKYUNKWAN UNIVERSITY, SUWON, GYEONGGI 440-746, KOREAEmail address:skim@skku.edu

Abstract
This paper focuses on computational contractual distinctions as an explanation for the spread between a forward contract and a similar futures contract which is derived and investigated. We evaluate this spread by constructing a time series model, which was established based on copula functions, and also show that the forward-futures spread is more significant for long maturity.
Keywords
1. INTRODUCTION
Forward and futures contracts are two seemingly similar financial instruments. Both are obligations to buy or sell a certain asset at an agreed price at a specified date in the future. Yet, the two instruments differ  . Whereas the forward is only settled at maturity, the futures contract is settled continuously through its lifespan. In other words, with the forward contract, cash is exchanged only at maturity, whereas with the futures contract, it is exchanged daily between contract initiation and maturity. This difference in cash flows has an impact on the value of the contracts.
It is our goal to show under which conditions the values of a forward and a futures contract are equal, and, if they differ, but how much. The difference in value between a forward contract and a similar futures contract is usually called forward-futures spread . Cox, Ingersoll, and Ross  analyze the relationship between forward and futures prices in a framework with stochastic interest rate. Grinblatt and Jegadeesh  compare the futures price with the forward price implied by the term structure of spot London Interbank Offered rates. Murawski  develop intuitive expressions for the forward-futures spread taking into account the probability of counterparty default. In this work, the methodology which is suited in pricing this spread is provided by generating stochastic processes with a given structure prescribed in terms of copula. The analysis of copula functions is started from the approach suggested by Nelsen  . Roughly speaking, the copula is a function that joins the marginal distributions of random variables in a multivariate distribution describing the joint behavior. Or statistically speaking, the copula connects a multivariate distribution to its marginals in such a way that it captures the entire dependence structure in the multivariate distribution. The importance of copulas in modeling the distribution of a multivariate random variable is justified by the Sklar’s theorem: any multivariate distribution can be expressed as its copula function evaluated at its marginal distribution function; and any copula function when evaluated at any marginal distributions is a multivariate distribution. Hence the information in the joint distribution is decomposed into those in the marginal distributions and that in the copula function. Consequently copulas allow one to model the marginal distributions and the dependence structure of a multivariate random variable separately. Moreover, the copula measure of dependence is invariant to any increasing transformation of individual series.
In the copula approach to time series modeling, the finite dimensional distributions of the time series are generated by copulas. Chen and Fan  and Pumi and Lopes  exhibit the simulation of time series based on some prescribed copula. By coupling different marginal distributions with different copula functions, copula-based time series models are able to model the dependence structure and the marginal behavior of a time series separately, allowing for a wide variety of marginal behaviors (such as skewness and fat tailed-ness) and dependence properties (such as asymmetric dependence and positive tail dependence). This separate modeling of the temporal dependence and the marginal behavior is particularly important when it is known that the dependence structure and the marginal properties of a time series are affected by different exogenous variables, which can easily modeled via the parametric copula approach.
This paper is organized as follows. Section 2 provides an overview of copula functions, including the measure of concordance. In Section 3, the forward contract and the similar futures contract as well as the difference between them are investigated. Section 4 establishes the copula-based time series model and Section 5 gives theevaluation of forward-futures spread based on this model. The conclusion of this paper is in Section 6.
2. COPULA FUNCTIONS
In this section we give the precise definition and important fundamental properties as well as some examples that are necessary for this work. We then present the measure of dependence, Kendall’s tau and also provide the relationships between the tau and copulas  .
- 2.1. Copula and Conditional Distribution
Definition 2.1. A two-dimensional copula is a function C : [0, 1] × [0, 1] → [0, 1] with the following property:
1. C(u, 0) = 0 = C(0, v) for every u, v ∈ [0, 1]. 2. C(u, 1) = u and C(1, v) = v for every u, v ∈ [0, 1]. 3. For every u1, u2, v1, v2 ∈ [0, 1] with u1u2 and v1v2:
C(u2, v2) − C(u2, v1) − C(u1, v2) + C(u1, v1) ≥ 0.
As such, it can present the joint distribution function of two standard uniform random variables U 1 , U 2 :
C ( u 1 , u 2 ) = Pr ( U 1 u 1 , U 2 u 2 )
We can use this feature in order to re-write via copulas the joint distribution function of two (even non-uniform) random variables. The most interesting fact about copulas in this sense is the Sklar’s theorem  .
Theorem 2.1 (Sklar’s Theorem). Let H ( x, y ) be a joint distribution function with continuous marginals F 1 and F 2 . Then there exists a copula C with
H ( x, y ) = C ( F 1 ( x ), F 2 ( y ))
for every PPT Slide
Lager Image
. If F 1 and F 2 are continuous, then C is unique. Otherwise , C is uniquely determined on Range F 1 × Range F 2 . On the other hand, if C is a copula and F 1 and F 2 are distribution functions, then the function H ( x , y ) = C ( F 1 ( x ), F 2 ( y )) is a joint distribution function with margins F 1 and F 2 .
The theorem is very important, because it provides a way to analyze the dependence structure of bivariate distributions without studying marginal distributions. Copulas are thus bivariate uniform distributions, which describe the dependence structure of random variables. Indeed, for many financial applications, the problem is not to use a given bivariate distribution but consists in finding a convenient distribution to describe some stylized facts, for example the relationships between different asset returns. Beside as the literature mentioned, strictly increasing transformation of the underlying random variables result in the transformed variables having the same copula, which is very important for practical applications.
Theorem 2.2. Consider a copula C ( u , v ), then
max( u + v − 1, 0) ≤ C ( u , v ) ≤ min ( u , v ).
As a consequence of Sklar’ theorem, if X and Y are random variables with a joint distribution function H and margins F and G , respectively, then for all x , y R , we have
max( F ( x ) + G ( y ) − 1, 0) ≤ H ( x , y ) ≤ min( F ( x ), G ( y ))
The above bounds are called the Fréchet-Hoeffding bounds for joint distribution functions H with margins F and G .
As we have known, dependency is an important concept inferring outcomes from a random variable based on the knowledge of a related factor. For an illustration, consider two uniform random variables U and V with a known copula C and U is observed. The goal is to deduce the conditional distribution which can then be used for predicting or estimating V . Assuming sufficient regularity, we obtain for the conditional cumulative distribution function PPT Slide
Lager Image
Hence, the conditional cumulative distribution function is just the partial derivative function of copula and can be derived directly from the copula itself. The next theorem will give an explanation to the existence of the partial derivative of copulas.
Theorem 2.3. Let C be a copula. For any v ∈ [0, 1], the partial derivative ∂C/∂u exists for almost all u, and for such v and u , PPT Slide
Lager Image
Similarly, for any u ∈ [0, 1], the partial derivative ∂C/∂v exists for almost all v, and for such v and u , PPT Slide
Lager Image
Furthermore, the functions u ∂C ( u , v )/ ∂u and are v ∂C ( u , v )/ ∂v defined and nondecreasing almost everywhere on [0, 1].
2.2. Some important Copulas It is time to present a few examples of copulas, especially for a particular class of copulas, the so-called Archimedean, which is particularly easy to handle and will be used in this paper.
2.2.1. Perfect dependence and independence The following three specific copulas are worth mentioning: the product copula, the minimum copula and the maximum copula . First, the product copula PPT Slide
Lager Image
related to the case of no dependence. As a direct consequence of Sklar’s theorem we obtain that random variables are independent if and only if their copula is the product copula. As already mentioned, the Frechet-Hoeffding bounds are also related to copulas, where the maximum copula is the Frechet-Hoeffding upper bound which is given by PPT Slide
Lager Image
Maximum copula refers to the case of perfect positive dependence. The other extreme is given by the minimum copula and as already mentioned, which is the Frechet-Hoeffding lower bound PPT Slide
Lager Image
This case refers to perfect negative dependence.
2.2.2. Archimedean copulas Archimedean copulas may be constructed using a generating function 𝜙 : [0, 1] → [0, ∞] continuous, strictly decreasing, convex and such that 𝜙(1) = 0. Given such a function 𝜙 a copula may be generated computing PPT Slide
Lager Image
where 𝜙 [-1] is the pseudo-inverse of 𝜙, defined as PPT Slide
Lager Image
Among Archimedean copulas, we are going to consider only the one-parameter ones, which are constructed using a generator 𝜙 θ ( t ), indexed by the parameter θ . Now three kinds of Archimedean copulas will be introduced as follows.
Gumbel Copula The bivariate Gumbel copula or Gumbel-Hougaard copula is given by the function PPT Slide
Lager Image
The parameter θ may take all values in the interval [1, ∞). For θ = 1 we discover the product copula, while for θ → ∞ the Gumbel copula tends to the maximum copula so that the Gumbel copula interpolates between independence and perfect positive dependence.
Clayton Copula A second example is the Clayton copula, given by the function PPT Slide
Lager Image
where θ ∈ [−1, 0)∪(0, ∞). For the limits as θ → 0, we obtain the product copula, while for θ → ∞, the Clayton copula arrives as the maximum copula. For θ = −1, we obtain the minimum copula. Thus, like the Gumbel copula, the Clayton copula interpolates between certain dependence structures, namely perfect negative dependence, independence and perfect positive dependence.
Frank Copula Another example is the Frank copula, given by PPT Slide
Lager Image
where θ ∈ (−∞, 0)∪(0, ∞). For θ → 0 Frank copula reduces to the product copula, while it can reach the maximum copula and the minimum copula for θ → ∞ and θ → −∞,respectively. Therefore, same with the Clayton copula, the Frank copula also interpolates between certain dependence structures, namely perfect negative dependence, independence and perfect positive dependence.
2.3. Measure of Concordance While a copula function fully describes the dependence structure between two random variables, we are often interested in simple summary statistics that indicate the extent and the direction of co-movement. The most basics of these measures are concordance measures, which assume the value of −1 for perfect negative dependence (corresponding to the lower Fréchet bound or the minimum copula), 0 for independence (corresponding to product copula) and 1 for perfect positive dependence (corresponding to the upper Fréchet bound or the maximum copula). Roughly speaking, two random variables are said to be concordant if large values of one come with large values of the other and small values of one come with small values of the other. Two samples ( xi , yi ) and ( xj , yj ) of a continuous random vector ( X , Y ) are concordant if ( xi xj )( yi yj ) ≥ 0. This means that either xi xj and yi yj or xi xj and yi yj . Concordance measures are completely determined by the underlying copula and independent of the marginal distributions unlike Pearson’s Correlation.
The two most common concordance measures are Kendall’s tau and Spearman’s rho. For two pairs of i.i.d. random vectors ( X 1 , Y 1 ) and ( X 2 , Y 2 ), Kendall’s tau has an interpretation as the difference between the probability of a joint outcome with the same sign Pr [(( X 1 X 2 )( Y 1 Y 2 ) > 0)] and the probability of a joint outcome with different sign Pr [(( X 1 X 2 )( Y 1 Y 2 ) < 0)], in other words the difference between the probability of the variables either rising or falling together (concordance) minus the probability of them moving in different directions (discordance). Similarly, Spearman’s rho measures the correlation between rank ordered data. The standard Pearson correlation coefficient is not a concordance measure and in particular, in contrast to Kendall’s tau and Spearman’s rho, it will be changed but non-monotonic transformations of the data since it is not independent of the marginal distributions.
Both measures have a lot of properties have a lot of properties in common. Being obviously measures with values in [−1, 1], they take the value 0 for independent variables (while there might also be non-independent random variables with zero rank correlation) and return 1 (−1) for the perfect positive dependence (perfect negative dependence) case. Moreover, they can directly be derived from the unique copula C that describes the dependence between two random variables X and Y . PPT Slide
Lager Image PPT Slide
Lager Image
In the following, the expression of τ as a function of the parameter θ of the copula are given in the following table, PPT Slide
Lager Image
where PPT Slide
Lager Image
dt and n is an integer.
3. FORWARD AND FUTURES CONTRACTS
Forward contract and futures contract are fundamentally similar. Both contracts are agreements between two parties to trade a specific good or asset at a future date at a strike price determined today. These agreements facilitate the transfer of risk between contracting parties, permitting the creation of hedging strategies  .
In spite of their conceptual similarity, the implementation of these two contracts is quite different. Participants in futures markets trade highly standardized contracts by open outcry on organized exchanges. Clearing associations for the exchanges control default risk by assuming the opposite side of each transaction, guaranteeing payment. To reduce traders’ incentive to default, the exchanges mark positions to market by settling contract gains and losses daily.
In contrast to the futures market, the market for forward contracts is mostly over the counter market. In the absence of a clearinghouse, trader reputation helps control default risk. Unlike futures contracts, no marking-to-market occurs: the parties in the contract make no payments until maturity of the contract. The forward market for financial interest rate instruments departs even further from the futures market format: a trader can use the cash market for interest rate instruments to establish a forward position. For example, a trader can secure a forward position for a three-month T-bill beginning three months in the future by buying a T-bill with a six-month maturity and selling a T-bill with a three-month maturity.
Some notions for forward contract and futures contract are given in the following.
• T: maturity date of the forward and futures contract,
• S(t): price at time t of the good or asset on which the contracts are written,
• R(t): the interest rate at time t,
• B(t,T): the price at time t of a zero-coupon bond paying 1 at time T,
• ForS(t,T): forward price at time t with maturity date T,
• FutS(t,T): futures price at time t with maturity date T.
3.1. Forward Contract As already mentioned, forward contract is an agreement to buy or sell an asset at a certain future time for a certain price K . It is traded in the over-the counter market, usually between two financial institutions or between a financial institution and one of its clients. One of the parties to a forward contract assumes a long position and agrees to buy the underlying asset on a certain specified future date for a certain specified price. The other party assumes a short position and agrees to sell the asset on the same date for the same price.
Now, we are going to give the price formula of forward contract. First we define the discount process PPT Slide
Lager Image
, then the price at time t of a zero-coupon bond paying 1 at time T can be written as PPT Slide
Lager Image
where PPT Slide
Lager Image
is the risk-neutral expectation.
From the fact that the discounted asset price is a martingale under the risk-neutral measure PPT Slide
Lager Image
, PPT Slide
Lager Image
This value should be zero and we have the following theorem.
Theorem 3.1. Assume that zero-coupon bonds of all maturities can be traded. Then the price of forward contract can be written as PPT Slide
Lager Image
3.2. Futures Contract Like a forward contract, a future contract is an agreement between two parties to buy or sell an asset at a certain time in the future for a certain price. Unlike forward contracts, futures contracts are normally traded on an exchange. To make trading possible, the exchange specifies certain standardized features of the contract. As the two parties to the contract do not necessarily know each other, the exchange also provides a mechanism that gives the two parties a guarantee that the contract will be honored. Hence the futures price of an asset should be an adapted stochastic process FutS ( t , T ) with two properties.
(1) The futures price agrees with the asset price on the delivery date i.e.
FutS ( T , T ) = S ( T ).
(2) The value of holding the futures contract over a period of time and receiving the cash flows associated with this position is zero: PPT Slide
Lager Image
The unique process satisfying these two properties can be defined as follows.
Definition 3.1. The futures price of an asset whose value at time T is S ( t ) is given by the formula PPT Slide
Lager Image
3.3. Forward-Futures Spread If the interest is a constant r , then B ( t , T ) = e r(Tt) and PPT Slide
Lager Image PPT Slide
Lager Image
In this case, the forward and futures prices agree.
When interest rates vary unpredictably (as they do in the real world), forward and futures prices are in theory no longer the same. We compare ForS ( t , T ) and FutS ( t , T ) with a random interest rate. In this case, from the fact that PPT Slide
Lager Image
we have PPT Slide
Lager Image
, and the forward-futures spread can be written as PPT Slide
Lager Image
where PPT Slide
Lager Image
denotes the covariance of D ( T ) and S ( T ) under the riskneutral measure. If the interest rate is nonrandom, this covariance is zero and the futures price agrees with the forward price.
Theorem 3.2. Under the risk-neutral measure with random interest rate, the price of forward-futures spread can be written as PPT Slide
Lager Image
4. TIME SERIES MODEL BASED ON COPULA
Let { Xt } be a stationary Markov process of order one. Then statistical properties are completely determined by the joint distribution of X t−1 , Xt , and H ( X t−1 , Xt ). By Sklar’s theorem, one can express H ( X t−1 , Xt ) in terms of the marginal distribution of X and the copula function of X t−1 and Xt . As a result, the statistical properties of a stationary Markov process { Xt } are completely determined by its marginal distribution and the copula of X t−1 and Xt . This suggests the copula approach as an alternative approach to modeling a stationary Markov process: instead of specifying the joint distribution of X t−1 and Xt , one specifies the marginal distribution of Xt and the copula function of X t−1 and Xt . The advantage of copula approach is that one has the freedom to choose the marginal distribution and the copula function separately; the former characterizes the marginal behavior such as the fat-tailed of the time series { Xt }, while the latter characterizes the temporal dependence property such as nonlinear, asymmetric dependence, of the time series.
Let PPT Slide
Lager Image
be a sample from a stationary first order Markov process and let F be the true invariant distribution, Suppose that F is absolutely continuous. Let Cθ be the true parametric copula associated to the pair ( X t−1 , Xt ) up to an unknown parameter θ . In this setting , if one writes Ut = F ( Xt ) and consider PPT Slide
Lager Image
, the transformed process is still a stationary Markov process, the joint distribution of ( U t−1 , Ut ) will be Cθ ( u t−1 , ut ) which can be written as according to Sklar’s theorem: PPT Slide
Lager Image
As already mentioned, the conditional distribution of Xt also can be written as PPT Slide
Lager Image
After getting the conditional distribution of Xt , one can generate the time-series by using inverse function PPT Slide
Lager Image PPT Slide
Lager Image
where ξt is the random number which follows U (0, 1)
Now it is very easy to simulate a time series from a copula-based Markov model. Let PPT Slide
Lager Image
be the conditional distribution of Ut given U t−1 = u t−1 . To generate a series PPT Slide
Lager Image
from a copula-based time series model, one may proceed as follows.
1. Choose the starting value u0 ∈ [0, 1]. 2. Generate n i.i.d. U(0, 1) random variables (r1, r2,···, rn). 3. For each i = 1, 2, ··· , n, set 4. For each i = 1, 2, ··· , n, set xi = F−1(ui). 5. The desired sample is (x1, x2, ··· , xn).
The above algorithm can be used to generate a wide variety of dependence structures.
In this part, we evaluate marginal densities of the stock and stochastic interest rate to estimate the forward-futures spread between the forward price and the futures price with stochastic simulations.
5.1. Computation of all Statistics and Marginal Densities The data on which our investigation is based are the historical stock price (SP500) and interest rate (LIBOR) ranging from 2 Jan. 1986 until 29 Dec. 2006, totally 5298 observations. The following table presents the sample moments of these time series, and figure 1 displays these two time series. For convenience, we denote the stock price as S ( t ) and the interest rate as R ( t ). PPT Slide
Lager Image PPT Slide
Lager Image
Scatter plots for marginal distribution
Form the historical data we obtain the marginal density function of each time series reported in figure 2 and figure 3 , respectively, and the following table explains the marginals for them. PPT Slide
Lager Image PPT Slide
Lager Image
Marginal density for stock price PPT Slide
Lager Image
Marginal density for interest rate
We have an intuitive impression of the dependence structure for each time series by graphical statistical analysis previously. And also we can have a further understand of such a dependence structure through calculating the concordance measure, Kendall’s tau. However, recalling that these statistics are particularly case of copulas’ application, the copula notion extends the concordance notion to the continuous framework. Moreover, even if the two concordance measures mentioned above allow to measure the association degree between two random variables, these statistics provide less information than a copula. In the following table, we report the estimated parameters θ for the families of Archimedean copulas described previously, as well as the measure concordance, Kendall’s tau. PPT Slide
Lager Image
5.2. Computation of Prices of Spread In what follows we will use Gumbel copula for simulating the time series. Here again, we write the form of Gumbel copula, and discover its partial derivative function. PPT Slide
Lager Image PPT Slide
Lager Image
We simulate 2000 times for S ( t ) and R ( t ) based on the above algorithm, then we get 2000 paths for each process. We take average after calculating the covariance of S ( t ) and D ( t ) and the bond price B (0, t ). PPT Slide
Lager Image PPT Slide
Lager Image
where Di ( t ), Si ( t ) and Bi (0, T ) is the value of D ( T ), S ( T ) and B (0, T ) for the i th path.
We vary T between 0.25 and 10, i.e. between three months and ten years, the values of the forward-futures spread are shown in the following table and figure 4 . PPT Slide
Lager Image PPT Slide
Lager Image