# Improved Maximum Likelihood Estimation of Heston Model and Pricing Efficiency Test: Hong Kong Hang Seng Index Option

Authors: Huan Wang; Bin Song; Dongmei Guo
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

This paper selects improved maximum likelihood method to conduct parameter estimation of Heston model, and results show that the share option pricing performance of Hang Seng Index is better and pricing error of at-the-money options is the smallest. By comparing parameter estimation of samples in different intervals, it has been found that parameter estimated results of two-year market data are obviously inferior to estimated effect of one-year data.


---

## Body

## 1. Introduction

This paper adopts Maximum Likelihood Estimation (MLE) to estimate Heston model parameters and test pricing efficiency using Hang Seng Index Option data. Empirical results will be given for reference for promotion of China’s stock index options, provide the scientific evidence for transaction pricing and hedging of issuers and investors, and lay a foundation for stable and healthy development of stock index options market.Cox and Ross [1] elicited CEV (Constant Elasticity of Variance) options pricing model on the condition that variance of underlying asset price yield rate was proportional to share price. Then Hull and White [2] assuming variance of asset price was subject to Geometric Brownian motion showed that option price was equal to mean value of Black-Scholes option prices of average variance of yield rates within duration of options. And Scott [3] and E. M. Stein and J. C. Stein [4] developed O-U (Ornstein-Uhlenbeck) process with mean reversion to describe variance of assets return. Then Heston [5] disproved assumption that risk premium of volatility and correlation coefficients were equal to zero based on O-U (Ornstein-Uhlenbeck) process. And Bates [6] brought random jumps into random volatility model and obtained analytic solution of European style option price through Fourier reverse transformation of characteristic equation complex value. Then the results were promoted by Bakshi et al. [7] by adding random interest rate. In another study, Carr and Wu [8] described asset price process using Levy Process with specified tail index and skewness parameter. However, many studies focus on Heston Model [9, 10]; for example, Ahlip and Rutkowski employed Heston model to describe exchange rate changes and CIR interest rate model to describe interest rate changes and finally they obtained analytic solution; Leung adopted Heston model to solve analytic pricing problem of back-looking options of floating exercise price. Different from the assumption that volatility is subject to some function or process in the above models, implied volatility is volatility value deducted reversely by substituting market price of options into Black-Scholes options pricing equation. For example, Dupire [11] put forward local volatility model based on implied volatility and reconstruct volatility of underlying asset in virtue of quoted price of Black-Scholes options and regularization method and then obtaining analytic formula of volatility.Heston model has been widely applied to industries and academic fields; the reason is that the characteristic equation method is employed to obtain analytic solution of options price while improving the assumption of constant volatility, so it is a major breakthrough in researches on random volatility models. Since Heston model is a nonlinear and non-Gaussian state-dependent space defined and its likelihood function form is quite complicated, thus it makes it very difficult to estimate parameters of Heston model. But there have been some studies on parameter estimation of Heston model. For example, Basawa and Prakasa Rao [12] elicited a function of drift coefficients and diffusion coefficients about unknown parameters by nonlinear Kalman filtering and conducted a MLE on model parameters based on continuous observed values of unknown parameters. Jacquier et al. [13] conducted parameter estimation through Markov Chain Monte Carlo. Christoffersen and Jacobs [14] adopted nonlinear least squares method to conduct parameter estimation of random volatility models. Ruiz [15] first proposed Pseudo-Maximum Likelihood (QML) to estimate parameters of random volatility models for the first time. But Ait-Sahalia [16] developed an Approximate Maximum Likelihood (AML) function and tried to solve incompatible problem of estimated results caused by directly applying AML to disperse observation. And Aït-Sahalia and Kimmel [17] gathered price data of S&P 500 Index Options market and used AML function to conduct parameter estimation of several random volatility models and obtained good estimation results. This paper, based on research of Ait-Sahalia, adopts MLE to conduct parameter estimation and empirical analysis.


## 2. Heston Model

In this section, we will briefly review Heston model.Assume that Heston model satisfies the following conditions:(
1
) Stock price follows stochastic process related to current volatility. (
2
) Short selling is allowed in stock market without transaction cost. (
3
) There is no stock dividend within the term of validity of options.In Heston model, underlying asset price follows stochastic differential equation:(1)
d
S
=
μ
S
d
t
+
Y
S
d
W
S
,where S is underlying asset price, Y is instantaneous variance of underlying asset price, and d
W
S is a Brownian movement. The variance Y derived by(2)
d
Y
=
k
θ
-
Y
d
t
+
σ
Y
d
W
Y
,where κ is mean reversion parameter, θ is long-term volatility level, and σ is volatility of Y. W
S has correlation ρ with W
Y:(3)
d
W
S
d
W
Y
=
ρ
d
t
.By using Ito Lemma and no arbitrage argument, Heston arrived at the following partial differential equation:(4)
∂
C
∂
t
+
S
2
Y
2
∂
2
C
∂
S
2
+
γ
S
∂
C
∂
S
-
γ
C
+
k
θ
-
V
-
λ
Y
∂
C
∂
Y
+
σ
2
Y
2
∂
2
C
∂
Y
2
+
ρ
σ
S
Y
∂
2
C
∂
S
∂
Y
=
0
.λ represents the price of volatility risk; Heston made an argument about analytic solution of the model that it had a form which is similar to Black-Scholes formula:(5)
C
S
t
,
Y
t
,
t
,
T
=
S
t
P
1
-
K
e
-
r
T
-
t
P
2
.P
1 and P
2 are defined by the following Fourier Transform:(6)
P
j
x
,
Y
t
,
T
,
K
=
1
2
+
1
π
∫
0
∞
R
e
e
-
i
ϕ
ln
⁡
K
f
j
x
,
Y
t
,
T
,
ϕ
i
ϕ
d
ϕ
,
j
=
1,2
,
x
=
ln
⁡
S
t
,where(7)
f
j
x
,
Y
t
,
T
,
ϕ
=
exp
⁡
C
T
-
t
,
ϕ
+
D
T
-
t
,
ϕ
Y
t
+
i
ϕ
x
,
C
T
-
t
,
ϕ
=
γ
ϕ
i
γ
+
a
σ
2
b
j
-
ρ
σ
ϕ
i
+
d
γ
-
2
ln
⁡
1
-
g
e
d
r
1
-
g
,
D
T
-
t
,
ϕ
=
b
j
-
ρ
σ
ϕ
i
+
d
σ
2
1
-
e
d
r
1
-
g
e
d
r
,
g
=
b
j
-
ρ
σ
ϕ
i
+
d
b
j
-
ρ
σ
ϕ
i
-
d
,
d
=
b
j
-
ρ
σ
ϕ
i
2
-
σ
2
2
μ
j
ϕ
i
-
ϕ
2
,
j
=
1,2
,
μ
1
=
1
2
,
μ
2
=
1
2
,
a
=
k
θ
,
b
1
=
k
+
λ
-
ρ
σ
,
b
2
=
k
+
λ
.


## 3. The Estimation Method

In the paper, we choose maximum likelihood method to estimate Heston model parameters. First, we determine a likelihood function for the state variables, representing conditional probability density of state variables. Based on the closed form approximation method of Ait-Sahalia [16], we calculate transferring density likelihood function of state variables and obtained likelihood function of its conditional probability density according to Markov property of state variable, making MLE method be applied to parameter estimation of models in a simpler way.


### 3.1. Likelihood Function of Multivariate Diffusion Model

We present the equation in full generality:(8)
d
X
t
=
μ
X
t
;
θ
d
t
+
σ
X
t
;
θ
d
W
t
,where X
t is an m-vector of state variables and μ
(
X
t
,
θ
) is an m-dimensional function, σ
(
X
t
,
θ
) is m
×
m matrix-valued function, W
t is an m-dimensional canonical Brownian motion, and θ is model parameter. Equation (8) can express multivariate diffusion models, including relevancy between variables in X
t, like Heston model with correlation coefficient ρ between two Brownian movements. The first element in X
t usually represents underlying asset price S
t, X
t can be expressed as X
t
=
S
t
;
Y
t
′, and Y
t is (m
-
1)-dimensional vector. On the condition that X
t
=
x
0, conditional probability density of X
t
+
Δ
=
x will be called transfer equation, and it can be expressed as p
X
(
x
∣
x
0
,
Δ
). As we suppose that discrete data is of Markov characteristics, according to Bayes rules, simple form of log-likelihood function can be easily obtained:(9)
l
n
θ
,
Δ
=
∑
i
=
1
n
l
X
X
i
Δ
∣
X
i
-
1
Δ
,
Δ
,where l
X
=
ln
⁡
p
X. Since analytic formula of p
X and ln
⁡
p
X cannot be obtained, then Ait-Sahalia [16] offered the analytic form of log-likelihood function of multivariate diffusion with the form of a power series: (10)
l
X
J
Δ
,
x
∣
x
0
;
θ
=
-
m
2
ln
⁡
2
π
Δ
-
D
v
x
,
θ
+
C
X
-
1
x
∣
x
0
;
θ
Δ
+
∑
k
=
0
J
C
X
k
x
∣
x
0
;
θ
Δ
k
k
!
,where D
v
(
x
;
θ
)
=
1
/
2
ln
⁡
(
det
⁡
[
v
(
x
;
θ
)
]
) and v
(
x
)
=
σ
x
σ
′
(
x
). This equation can be calculated to arbitrary order J. C
X
(
k
) is corresponding to Taylor expansion nearby x
0 at order j
k. Ait-Sahalia [16] indicates j
k
=
2
(
J
-
k
). C
X
(
k
) calculated in this way will be expressed as C
X
(
j
k
,
k
). Then (10) can be expressed as(11)
l
~
X
J
Δ
,
x
∣
x
0
;
θ
=
-
m
2
ln
⁡
2
π
Δ
-
D
v
x
,
θ
+
C
X
j
-
1
,
-
1
x
∣
x
0
;
θ
Δ
+
∑
k
=
0
J
C
X
j
k
,
k
x
∣
x
0
;
θ
Δ
k
k
!
.According to the method of Ait-Sahalia, we can use Kolmogorov equation to solveC
X
(
j
k
,
k
).


### 3.2. Likelihood Function Calculation of Heston Model

Write Heston model in a multivariate model framework, and reexpress it as the following matrix form; consider it first under measureQ; then(12)
d
X
t
=
d
S
t
Y
t
=
r
-
d
S
t
κ
′
γ
′
-
Y
t
d
t
+
1
-
ρ
2
Y
t
S
t
0
ρ
Y
t
S
t
σ
Y
t
β
d
W
1
Q
t
W
2
Q
t
.Then conduct logarithmetics of stock price to obtain the following form:(13)
d
S
t
Y
t
=
r
-
d
-
1
2
Y
t
κ
′
γ
′
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
Q
t
W
2
Q
t
.Convert it under measureP; (12) will be transformed into(14)
d
S
t
Y
t
=
r
-
d
-
1
2
Y
t
+
b
Y
t
κ
γ
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
P
t
W
2
P
t
,
a
=
r
-
d
,
b
=
λ
1
1
-
ρ
2
+
λ
2
ρ
-
1
2
,
κ
=
κ
′
-
λ
2
σ
,
γ
=
κ
+
λ
2
σ
κ
γ
′
.Stock dividend is not considered; then Heston model underP measure can be simplified as(15)
d
S
t
Y
t
=
r
+
λ
1
-
ρ
2
-
1
2
Y
t
κ
γ
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
P
t
W
2
P
t
.According to this model and the method about likelihood function analytic formula proposed by Ait-Sahalia [16], we can calculate coefficient C
X
(
j
k
,
k
) in likelihood function of conditional probability density of state variables in Heston model. In Heston model, likelihood function of conditional probability density of state variable X
t is(16)
l
X
J
Δ
,
x
∣
x
0
;
θ
=
-
m
2
ln
⁡
2
π
Δ
-
D
v
x
,
θ
+
C
X
-
1
x
∣
x
0
;
θ
Δ
+
∑
k
=
0
J
C
X
k
x
∣
x
0
;
θ
Δ
k
k
!
.Using this likelihood expression formula, we input the observed value[
S
t
,
Y
t
] and adopt Matlab to compile MLE function to estimate the parameters.


### 3.3. Input Variable in the Likelihood Function

In the following, we divide two steps to talk about the input variable in the likelihood function. First, reduce spatial dimension of input variable; second, determine the replacement of momentary volatility.(1
) Reduce Spatial Dimension of Input Variable. We examine each of those steps to reduce spatial dimension of input variables. Firstly, we limit options duration T
i
-
t to a fixed value to reduce a dimension. For example, the options are selected which will be due three months later among listed options on each trading day; secondly, we consider selecting a fixed option value which can be expressed as S
/
K. For a call option, in-the-money options refer to options with S
/
K being greater than 1, at-the-money options refer to options with S
/
K being 1, and out-of-the-money options refer to options with S
/
K being less than 1. Fixed exercise prices K are easy way to reducing a dimension. However, set an exercise price of at-the-money option around share price or stock index in each trading day, and with the fluctuation of stock price S, the distribution of exercise price K is not machine-parsable, so it cannot be guaranteed that option contract of exercise price K exists in each trading day within sampling period when data is being collected. Based on this listing method, it is more feasible for us to set an option price of rational value state; thirdly, assume the option price is homogeneous. Payoff of European call option is a linear function about stock price and exercise price. The price of one European call option can be expressed as [
E
(
S
T
)
-
K
]
+
·
e
-
∫
r
(
t
)
d
t, so it can be considered that(17)
C
T
,
α
S
T
,
α
K
=
α
S
T
-
α
K
+
=
α
S
T
-
K
+
=
α
C
T
,
S
T
,
K
.Generally, the price of one option is not of homogeneity except for some additional limitations for the model. For example,(18)
σ
1
X
t
σ
1
′
X
t
=
φ
11
Y
t
S
t
2
,
σ
1
X
t
σ
i
′
X
t
=
φ
1
i
Y
t
S
t
=
φ
i
1
Y
t
S
t
,
i
>
1
σ
i
X
t
σ
j
′
X
t
=
φ
i
j
Y
t
,
i
>
1
,
j
>
1
μ
i
Q
X
t
=
ψ
i
Y
t
,
i
>
1
.So far, we have successfully reduced three dimensions of input variables.(2
) Replacement of Instantaneous Volatility
Y
t. The instantaneous volatility Y
t of stock price is an unobservable variable. Generally there are two methods to solve this problem.One method is to modify input variables. Convert input variables of likelihood function from observable variables into unobservable variables.AsY
t is unobservable, we convert it into an observable variable, such as C
t. We have obtained expression formula of conditional probability density log-likelihood function (16) adopted in parameter estimation. This expression formula requires input variable to be S
t
;
Y
t
′. According to MLE, a specific observed value is input through all forms of parameter combination θ. However, during the process of searching the optimal θ, with continuous variation of θ, X
t value will change, which indicates that, in a strict sense, the parameter estimation method which inputs [
S
t
;
Y
t
] to conduct maximum searching of likelihood function is unscientific.Hence, we can try to convert input variable from unobservable state variableX
t
=
S
t
;
Y
t
′ into observable asset price vector G
t
=
S
t
;
C
t
′, thus converting input variable S
t
;
Y
t
′ into S
t
;
C
t
′. In fact, what we obtain from the preceding part of this paper is likelihood function about X
t
=
S
t
;
Y
t
′. A series of volatility Y
t can be used to calculate a series of option price for one set of parameters. By comparing it with C
t, select the most appropriate volatility Y
t of this set of parameters, so there will be corresponding X
t
=
S
t
;
Y
t
′ for this set of parameters.In actual operation of parameter estimation, it is necessary to make certain adjustment of likelihood function aboutX
t. Then X
t needs to be multiplied by a Jacobian term, and it will be a partial differential equation of state variable X
t with respect to S
t and C
t. Put these differential expressions in a matrix, and Jacobian term is the determinant of this matrix, which is in the following form:
(19)
J
t
=
det
⁡
∂
S
t
∂
S
t
∂
S
t
∂
Y
t
1
⋯
∂
S
t
∂
Y
t
N
∂
C
t
1
∂
S
t
∂
C
t
1
∂
Y
t
1
⋯
∂
C
t
1
∂
Y
t
N
⋮
⋮
⋱
⋮
∂
C
t
N
∂
S
t
∂
C
t
N
∂
Y
t
1
⋯
∂
C
t
N
∂
Y
t
N
=
det
⁡
1
0
⋯
0
∂
C
t
1
∂
S
t
∂
C
t
1
∂
Y
t
1
⋯
∂
C
t
1
∂
Y
t
N
⋮
⋮
⋱
⋮
∂
C
t
N
∂
S
t
∂
C
t
N
∂
Y
t
1
⋯
∂
C
t
N
∂
Y
t
N
=
det
⁡
∂
C
t
1
∂
Y
t
1
⋯
∂
C
t
1
∂
Y
t
N
⋮
⋱
⋮
∂
C
t
N
∂
Y
t
1
⋯
∂
C
t
N
∂
Y
t
N
.The following is the process of converting likelihood function of state variableX
t into that of state variable G
t through Jacobian term. Suppose that transition density function of X
t is p
X
(
Δ
,
x
∣
x
0
;
θ
) and that of G
t is p
G
(
Δ
,
g
∣
g
0
;
θ
); write observable asset price vector G
t as function of state variable X
t:(20)
G
t
+
Δ
=
f
X
t
+
Δ
;
θ
.Then its inverse function to expresses that state is the function of observed asset prices:(21)
X
t
+
Δ
=
f
-
1
G
t
+
Δ
;
θ
.Given thatG
t
=
g
0, conditional probability density of G
t
+
Δ
=
g is(22)
p
G
Δ
,
g
∣
g
0
;
θ
=
d
e
t
∂
f
f
-
1
g
;
θ
∂
x
-
1
p
X
Δ
,
f
-
1
g
;
θ
∣
f
-
1
g
0
;
θ
;
θ
=
J
t
Δ
,
g
∣
g
0
;
θ
-
1
p
X
Δ
,
f
-
1
g
;
θ
∣
f
-
1
g
0
;
θ
;
θ
,where the multiplier J
t
(
Δ
,
g
∣
g
0
;
θ
) will be calculated by referring to (19).Observable asset prices have the property of Markov, and they can be obtained using Bayes Rule. Log-likelihood function for discrete data of vectorg
t at dates t
0
,
t
1
,
…
,
t
n can be expressed as(23)
l
n
θ
≡
n
-
1
∑
i
=
1
n
l
G
t
i
-
t
i
-
1
,
g
t
i
∣
g
t
i
-
1
;
θ
,where(24)
l
G
Δ
,
g
∣
g
0
;
θ
=
ln
⁡
p
G
Δ
,
g
∣
g
0
;
θ
=
ln
⁡
J
t
Δ
,
g
∣
g
0
;
θ
+
l
X
Δ
,
f
-
1
g
,
θ
∣
f
-
1
g
0
,
θ
;
θ
.Then log-likelihood function about observable asset prices can be obtained.Next we introduce the alternative method to replace momentary volatility by Black-Scholes implied volatility.Input variable of this method is[
S
t
,
,
σ
i
m
v
o
l
t
]. The evolution of stock price does not conform to the assumptions of Black-Scholes model in this paper; because of operability of empirical inspection, we will use implied volatility of short-term at-the-money options calculated by Black-Scholes formula as a substitute variable of momentary volatility of stock price. This method has greatly simplified complexity of parameter estimation, especially calculating Jacobian term according to when input variables are converted.Here, the reasons for selecting implied volatility of short-term at-the-money options to replace momentary volatility are as follows: firstly, based on preceding researches, as options are close to due, implied volatility of short-term at-the-money options is converged to momentary volatility of stock price; secondly, short-term at-the-money options are the most active option kind with the best mobility, so price of this kind of options will be less influenced by factors such as market microstructure; thirdly, compared with other options, prices of at-the-money options are the most sensitive to change of volatility, and a slight change of volatility will result in a major change of at-the-money options price. Correspondingly, if the observation error on options price is small, then the error of volatility will be less, which can be ignored.There is one shortcoming of this method. It does not take some fluctuating factors underQ measure. While in practical process of searching the best θ, with continuous change of θ, X
t also changes not as invariable as input observable value [
S
t
,
,
σ
i
m
v
o
l
t
] is.


## 3.1. Likelihood Function of Multivariate Diffusion Model

We present the equation in full generality:(8)
d
X
t
=
μ
X
t
;
θ
d
t
+
σ
X
t
;
θ
d
W
t
,where X
t is an m-vector of state variables and μ
(
X
t
,
θ
) is an m-dimensional function, σ
(
X
t
,
θ
) is m
×
m matrix-valued function, W
t is an m-dimensional canonical Brownian motion, and θ is model parameter. Equation (8) can express multivariate diffusion models, including relevancy between variables in X
t, like Heston model with correlation coefficient ρ between two Brownian movements. The first element in X
t usually represents underlying asset price S
t, X
t can be expressed as X
t
=
S
t
;
Y
t
′, and Y
t is (m
-
1)-dimensional vector. On the condition that X
t
=
x
0, conditional probability density of X
t
+
Δ
=
x will be called transfer equation, and it can be expressed as p
X
(
x
∣
x
0
,
Δ
). As we suppose that discrete data is of Markov characteristics, according to Bayes rules, simple form of log-likelihood function can be easily obtained:(9)
l
n
θ
,
Δ
=
∑
i
=
1
n
l
X
X
i
Δ
∣
X
i
-
1
Δ
,
Δ
,where l
X
=
ln
⁡
p
X. Since analytic formula of p
X and ln
⁡
p
X cannot be obtained, then Ait-Sahalia [16] offered the analytic form of log-likelihood function of multivariate diffusion with the form of a power series: (10)
l
X
J
Δ
,
x
∣
x
0
;
θ
=
-
m
2
ln
⁡
2
π
Δ
-
D
v
x
,
θ
+
C
X
-
1
x
∣
x
0
;
θ
Δ
+
∑
k
=
0
J
C
X
k
x
∣
x
0
;
θ
Δ
k
k
!
,where D
v
(
x
;
θ
)
=
1
/
2
ln
⁡
(
det
⁡
[
v
(
x
;
θ
)
]
) and v
(
x
)
=
σ
x
σ
′
(
x
). This equation can be calculated to arbitrary order J. C
X
(
k
) is corresponding to Taylor expansion nearby x
0 at order j
k. Ait-Sahalia [16] indicates j
k
=
2
(
J
-
k
). C
X
(
k
) calculated in this way will be expressed as C
X
(
j
k
,
k
). Then (10) can be expressed as(11)
l
~
X
J
Δ
,
x
∣
x
0
;
θ
=
-
m
2
ln
⁡
2
π
Δ
-
D
v
x
,
θ
+
C
X
j
-
1
,
-
1
x
∣
x
0
;
θ
Δ
+
∑
k
=
0
J
C
X
j
k
,
k
x
∣
x
0
;
θ
Δ
k
k
!
.According to the method of Ait-Sahalia, we can use Kolmogorov equation to solveC
X
(
j
k
,
k
).


## 3.2. Likelihood Function Calculation of Heston Model

Write Heston model in a multivariate model framework, and reexpress it as the following matrix form; consider it first under measureQ; then(12)
d
X
t
=
d
S
t
Y
t
=
r
-
d
S
t
κ
′
γ
′
-
Y
t
d
t
+
1
-
ρ
2
Y
t
S
t
0
ρ
Y
t
S
t
σ
Y
t
β
d
W
1
Q
t
W
2
Q
t
.Then conduct logarithmetics of stock price to obtain the following form:(13)
d
S
t
Y
t
=
r
-
d
-
1
2
Y
t
κ
′
γ
′
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
Q
t
W
2
Q
t
.Convert it under measureP; (12) will be transformed into(14)
d
S
t
Y
t
=
r
-
d
-
1
2
Y
t
+
b
Y
t
κ
γ
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
P
t
W
2
P
t
,
a
=
r
-
d
,
b
=
λ
1
1
-
ρ
2
+
λ
2
ρ
-
1
2
,
κ
=
κ
′
-
λ
2
σ
,
γ
=
κ
+
λ
2
σ
κ
γ
′
.Stock dividend is not considered; then Heston model underP measure can be simplified as(15)
d
S
t
Y
t
=
r
+
λ
1
-
ρ
2
-
1
2
Y
t
κ
γ
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
P
t
W
2
P
t
.According to this model and the method about likelihood function analytic formula proposed by Ait-Sahalia [16], we can calculate coefficient C
X
(
j
k
,
k
) in likelihood function of conditional probability density of state variables in Heston model. In Heston model, likelihood function of conditional probability density of state variable X
t is(16)
l
X
J
Δ
,
x
∣
x
0
;
θ
=
-
m
2
ln
⁡
2
π
Δ
-
D
v
x
,
θ
+
C
X
-
1
x
∣
x
0
;
θ
Δ
+
∑
k
=
0
J
C
X
k
x
∣
x
0
;
θ
Δ
k
k
!
.Using this likelihood expression formula, we input the observed value[
S
t
,
Y
t
] and adopt Matlab to compile MLE function to estimate the parameters.


## 3.3. Input Variable in the Likelihood Function

In the following, we divide two steps to talk about the input variable in the likelihood function. First, reduce spatial dimension of input variable; second, determine the replacement of momentary volatility.(1
) Reduce Spatial Dimension of Input Variable. We examine each of those steps to reduce spatial dimension of input variables. Firstly, we limit options duration T
i
-
t to a fixed value to reduce a dimension. For example, the options are selected which will be due three months later among listed options on each trading day; secondly, we consider selecting a fixed option value which can be expressed as S
/
K. For a call option, in-the-money options refer to options with S
/
K being greater than 1, at-the-money options refer to options with S
/
K being 1, and out-of-the-money options refer to options with S
/
K being less than 1. Fixed exercise prices K are easy way to reducing a dimension. However, set an exercise price of at-the-money option around share price or stock index in each trading day, and with the fluctuation of stock price S, the distribution of exercise price K is not machine-parsable, so it cannot be guaranteed that option contract of exercise price K exists in each trading day within sampling period when data is being collected. Based on this listing method, it is more feasible for us to set an option price of rational value state; thirdly, assume the option price is homogeneous. Payoff of European call option is a linear function about stock price and exercise price. The price of one European call option can be expressed as [
E
(
S
T
)
-
K
]
+
·
e
-
∫
r
(
t
)
d
t, so it can be considered that(17)
C
T
,
α
S
T
,
α
K
=
α
S
T
-
α
K
+
=
α
S
T
-
K
+
=
α
C
T
,
S
T
,
K
.Generally, the price of one option is not of homogeneity except for some additional limitations for the model. For example,(18)
σ
1
X
t
σ
1
′
X
t
=
φ
11
Y
t
S
t
2
,
σ
1
X
t
σ
i
′
X
t
=
φ
1
i
Y
t
S
t
=
φ
i
1
Y
t
S
t
,
i
>
1
σ
i
X
t
σ
j
′
X
t
=
φ
i
j
Y
t
,
i
>
1
,
j
>
1
μ
i
Q
X
t
=
ψ
i
Y
t
,
i
>
1
.So far, we have successfully reduced three dimensions of input variables.(2
) Replacement of Instantaneous Volatility
Y
t. The instantaneous volatility Y
t of stock price is an unobservable variable. Generally there are two methods to solve this problem.One method is to modify input variables. Convert input variables of likelihood function from observable variables into unobservable variables.AsY
t is unobservable, we convert it into an observable variable, such as C
t. We have obtained expression formula of conditional probability density log-likelihood function (16) adopted in parameter estimation. This expression formula requires input variable to be S
t
;
Y
t
′. According to MLE, a specific observed value is input through all forms of parameter combination θ. However, during the process of searching the optimal θ, with continuous variation of θ, X
t value will change, which indicates that, in a strict sense, the parameter estimation method which inputs [
S
t
;
Y
t
] to conduct maximum searching of likelihood function is unscientific.Hence, we can try to convert input variable from unobservable state variableX
t
=
S
t
;
Y
t
′ into observable asset price vector G
t
=
S
t
;
C
t
′, thus converting input variable S
t
;
Y
t
′ into S
t
;
C
t
′. In fact, what we obtain from the preceding part of this paper is likelihood function about X
t
=
S
t
;
Y
t
′. A series of volatility Y
t can be used to calculate a series of option price for one set of parameters. By comparing it with C
t, select the most appropriate volatility Y
t of this set of parameters, so there will be corresponding X
t
=
S
t
;
Y
t
′ for this set of parameters.In actual operation of parameter estimation, it is necessary to make certain adjustment of likelihood function aboutX
t. Then X
t needs to be multiplied by a Jacobian term, and it will be a partial differential equation of state variable X
t with respect to S
t and C
t. Put these differential expressions in a matrix, and Jacobian term is the determinant of this matrix, which is in the following form:
(19)
J
t
=
det
⁡
∂
S
t
∂
S
t
∂
S
t
∂
Y
t
1
⋯
∂
S
t
∂
Y
t
N
∂
C
t
1
∂
S
t
∂
C
t
1
∂
Y
t
1
⋯
∂
C
t
1
∂
Y
t
N
⋮
⋮
⋱
⋮
∂
C
t
N
∂
S
t
∂
C
t
N
∂
Y
t
1
⋯
∂
C
t
N
∂
Y
t
N
=
det
⁡
1
0
⋯
0
∂
C
t
1
∂
S
t
∂
C
t
1
∂
Y
t
1
⋯
∂
C
t
1
∂
Y
t
N
⋮
⋮
⋱
⋮
∂
C
t
N
∂
S
t
∂
C
t
N
∂
Y
t
1
⋯
∂
C
t
N
∂
Y
t
N
=
det
⁡
∂
C
t
1
∂
Y
t
1
⋯
∂
C
t
1
∂
Y
t
N
⋮
⋱
⋮
∂
C
t
N
∂
Y
t
1
⋯
∂
C
t
N
∂
Y
t
N
.The following is the process of converting likelihood function of state variableX
t into that of state variable G
t through Jacobian term. Suppose that transition density function of X
t is p
X
(
Δ
,
x
∣
x
0
;
θ
) and that of G
t is p
G
(
Δ
,
g
∣
g
0
;
θ
); write observable asset price vector G
t as function of state variable X
t:(20)
G
t
+
Δ
=
f
X
t
+
Δ
;
θ
.Then its inverse function to expresses that state is the function of observed asset prices:(21)
X
t
+
Δ
=
f
-
1
G
t
+
Δ
;
θ
.Given thatG
t
=
g
0, conditional probability density of G
t
+
Δ
=
g is(22)
p
G
Δ
,
g
∣
g
0
;
θ
=
d
e
t
∂
f
f
-
1
g
;
θ
∂
x
-
1
p
X
Δ
,
f
-
1
g
;
θ
∣
f
-
1
g
0
;
θ
;
θ
=
J
t
Δ
,
g
∣
g
0
;
θ
-
1
p
X
Δ
,
f
-
1
g
;
θ
∣
f
-
1
g
0
;
θ
;
θ
,where the multiplier J
t
(
Δ
,
g
∣
g
0
;
θ
) will be calculated by referring to (19).Observable asset prices have the property of Markov, and they can be obtained using Bayes Rule. Log-likelihood function for discrete data of vectorg
t at dates t
0
,
t
1
,
…
,
t
n can be expressed as(23)
l
n
θ
≡
n
-
1
∑
i
=
1
n
l
G
t
i
-
t
i
-
1
,
g
t
i
∣
g
t
i
-
1
;
θ
,where(24)
l
G
Δ
,
g
∣
g
0
;
θ
=
ln
⁡
p
G
Δ
,
g
∣
g
0
;
θ
=
ln
⁡
J
t
Δ
,
g
∣
g
0
;
θ
+
l
X
Δ
,
f
-
1
g
,
θ
∣
f
-
1
g
0
,
θ
;
θ
.Then log-likelihood function about observable asset prices can be obtained.Next we introduce the alternative method to replace momentary volatility by Black-Scholes implied volatility.Input variable of this method is[
S
t
,
,
σ
i
m
v
o
l
t
]. The evolution of stock price does not conform to the assumptions of Black-Scholes model in this paper; because of operability of empirical inspection, we will use implied volatility of short-term at-the-money options calculated by Black-Scholes formula as a substitute variable of momentary volatility of stock price. This method has greatly simplified complexity of parameter estimation, especially calculating Jacobian term according to when input variables are converted.Here, the reasons for selecting implied volatility of short-term at-the-money options to replace momentary volatility are as follows: firstly, based on preceding researches, as options are close to due, implied volatility of short-term at-the-money options is converged to momentary volatility of stock price; secondly, short-term at-the-money options are the most active option kind with the best mobility, so price of this kind of options will be less influenced by factors such as market microstructure; thirdly, compared with other options, prices of at-the-money options are the most sensitive to change of volatility, and a slight change of volatility will result in a major change of at-the-money options price. Correspondingly, if the observation error on options price is small, then the error of volatility will be less, which can be ignored.There is one shortcoming of this method. It does not take some fluctuating factors underQ measure. While in practical process of searching the best θ, with continuous change of θ, X
t also changes not as invariable as input observable value [
S
t
,
,
σ
i
m
v
o
l
t
] is.


## 4. Empirical Analysis: Based on Hong Kong Hang Seng Index Options

On one hand, empirical analysis can examine applicability of Heston model to Hong Kong market, estimate relevant parameters, and provide investors of Hang Seng Index Options with a reference for price forecasting of future transactions. On the other hand, Hong Kong market is closed related with domestic market of China, so pricing method of Hong Kong market may have good reference significance for promoting Hushen 300 Index Options in China.


### 4.1. Data Selection

This section selects day-degree data of Hang Seng Index Options price in Hong Kong Exchange. The data set consists of continuous 500 transaction days spanning from November 24th, 2011, to November 29th, 2013, and complete 247 transaction days spanning from December 3rd, 2012, to November 29th, 2013.In order to reduce dimensions to match number of parameters and input variables, we need to fix due dates of options selected in the sample. However, among option contract product sequences in Hong Kong Exchange, options that can be transacted each day include the following:(1)
Hang Seng Index Options expire at the end of each of the recent three month.
(2)
Hang Seng Index Options expire on the last transaction day at the end of March, June, September, and December within the year.Hence, it cannot be guaranteed that durations of options of each transaction day in the sample are identical. We can only defer the expiration date of sample options of each month by one month and guarantee that duration of selected option in the sample will be with limits. Specifically, we select data with its expiration date being 60–90 days. For example, for January 1st, we select the option contract that will expire at the end of March, so its duration will be about 90 days, but, for March 20th, we select the option contract which will become due at the end of May, so its duration will be about 70 days. For empirical research, we consider expiration dates of the two as identical.For convenience, we only select call option contract among these data. Since the trading volume of in-the-money options of Hong Kong market is small, and some only have single-digit trading volume, thus we select at-the-money options and some out-of-the-money options to conduct an empirical analysis with big trading volumes. Based on variation situation of everyday Hang Seng Index, we respectively find out exercise priceK of corresponding contracts approximate to at-the-money option S
t
/
K
≈
1, out-of-the-money option S
t
/
K
≈
0.95, and S
t
/
K
≈
0.92 as well as closing price C
t of these option contracts.


### 4.2. Analysis of Empirical Results

(
1) Parameter Estimation Results and Analysis. Table 1 lists estimation results selecting at-the-money option price data and out-of-the-money option price data within the sample.Table 1
Heston model estimation values and standard deviation of at-the-money call option and two out-of-the-money options based on three-month duration.
Heston model parameters
(
S
t
/
K
)
=
1Parameter estimation values
(
S
t
/
K
)
=
1Standard deviation
(
S
t
/
K
)
=
0.95Parameter estimation value
(
S
t
/
K
)
=
0.95Standard deviation
(
S
t
/
K
)
=
0.92Parameter estimation value
(
S
t
/
K
)
=
0.92Standard deviation
κ
0.4264
0.2784
0.4373
0.2086
0.4459
0.2113
γ
0.0946
0.0181
0.0915
0.0163
0.0909
0.0168
ρ
−0.6284
0.0867
−0.7011
0.0715
−0.6812
0.0790
λ
3.4543
3.0332
3.4280
2.9312
3.3788
2.7003
σ
0.4822
0.0433
0.3565
0.0333
0.3421
0.0315It can be seen from Table1 that, among three samples, κ obtained through the estimation basically stabilizes within 0.42–0.45. Parameter κ is mean reversion speed of volatility square which represents the speed of volatility and its square reversing to mean value when they deviate from the long-term mean value. Aït-Sahalia and Kimmel [17] conducted an empirical analysis using S&P 500 and its stock index option data and obtained that estimated value of reversion rate κ is 5.07, which means that, compared with American market, the mean reversion rate of Hong Kong market is slower, which indicates that when volatility in Hong Kong stock market deviates from mean value upwards or downwards, it will reverse to the mean level at a slow speed. Long-term mean value of volatility square is around 0.09, which indicates that underlying asset price volatility will move around 30% within a long period, and correlation coefficient ρ of two Brownian movements will be around −0.6 to −0.7, and λ stabilizes around 3.4, but variable range estimated value of volatility σ of stock index volatility is slightly larger, from 0.34 to 0.48.In the model,ρ is the risk source of stock index risk source and stock index volatility square, that is, the correlation coefficient between two Brownian movements d
W
1 and d
W
2; the above empirical results show that ρ is a negative value and remains within −0.6 to −0.7, which indicates that HSI stock index yield rate has a negative correlation with its volatility. This result has overthrown complete negative correlation or complete noncorrelation between underlying asset price and its volatility in other random volatility models to exactly explained “leverage effect” observed in the market.On the whole, standard deviation of estimated variable value of at-the-money option sample(
S
t
/
K
)
=
1 is higher than the two out-of-the-money option samples (
S
t
/
K
)
=
0.95 and (
S
t
/
K
)
=
0.92. We try to expand data length within the sample and select 500-day at-the-money option price data within the sample to conduct parameter estimation, and the results are as shown in Table 2.Table 2
Parameter estimation results of 500-day at-the-money option price data.
Heston model parameters
κ
γ
ρ
λ
σ
Estimated values of Heston model parameters
1.1409
12.8701
−0.9957
−175.9290
0.0000It can be seen that, among 500-day parameter estimation results, interest rater has negative values, σ is 0, and reversed mean γ of volatility mean reversion equation reaches as high as 12.8701, which severely does not conform to current situation. This phenomenon can be explained from the angle of parameter estimation method. MLE method adopts closing price of at-the-money call option contracts of the current day as data of this day, so the bigger the selected sample is, the bigger the time span will be. If 500-day data in the sample is selected, sample span should be more than two years. As options market is not fixed, so when MLE method is used to conduct parameter estimation of options pricing research, larger sample size does not mean more ideal results. Instead, a proper sample size should be selected to guarantee adequacy of sample data and avoid estimation deviation caused by changes of multiperiod price characteristics.(
2) Analysis of Options Pricing Error. Here, according to parameter estimation results of data in option sample under different value states, we calculate options pricing errors 15 days beyond the sample, and Figures 1–3, respectively, express options pricing errors under different value states.Figure 1
Options pricing error of at-the-money option(
S
t
/
K
)
=
1 price.Figure 2
Options pricing error of out-of-the-money option(
S
t
/
K
)
=
0.95 price.Figure 3
Options pricing error of in-the-money call option(
S
t
/
K
)
=
0.92 price.From Figures1–3, we can see that it has a higher accuracy to estimate model parameters and deduce options price using at-the-money options, compared with the two out-of-the-money options, and pricing errors of 15 option contracts outside the sample are within 10%, and there are 7 ones with their errors within 3% and 11 ones with errors within 6%, which indicates that the model can be better priced at-the-money call options with expiration date being about 3 months. Table 3 reflects the average pricing error, where(25)
B
i
a
s
i
=
C
i
m
o
d
e
l
t
-
C
i
m
a
r
k
e
t
t
,
P
¯
a
b
s
o
l
u
t
e
=
a
v
e
r
a
g
e
∑
t
=
1
T
∑
i
=
1
n
t
B
i
a
s
i
t
C
i
m
a
r
k
e
t
t
,
P
¯
a
c
t
u
a
l
=
a
v
e
r
a
g
e
∑
t
=
1
T
∑
i
=
1
n
t
B
i
a
s
i
t
C
i
m
a
r
k
e
t
t
.Table 3
Options pricing error.
(
S
t
/
K
)
=
1
(
S
t
/
K
)
=
0.95
(
S
t
/
K
)
=
0.92
Average absolute error (P
¯
a
b
s
o
l
u
t
e)
3.66%
9.23%
10.07%
Average relative error (P
¯
a
c
t
u
a
l)
3.30%
−1.18%
3.35%Other origins of parameter estimation error and model pricing error include simplification of real measure. In order to simplify the model and reduce number of parameters to be estimated, for the convenience of parameter estimation, the model adopted in the empirical analysis is the model form in risk neutral measure (Q measure); real measure (P measure) is not taken into consideration when parameters are being set, so price of risk factors in P measure has not been estimated out; for example, Heston model in (17) is(26)
d
S
t
Y
t
=
r
-
d
-
1
2
Y
t
+
λ
1
-
ρ
2
Y
t
κ
γ
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
P
t
W
2
P
t
.(
3) Parameter Sensitivity Analysis of Options Pricing. In order to inspect reliance of options pricing effect on accuracy of parameter estimation, we, respectively, analyze sensitivity of each parameter of the model for options price estimation value.From Figures4–8, it is easy to see the parameters monotonicity of options price. Parameters κ and γ have positive correlation with estimations of options price, and overestimation of parameters will result in overestimation of options price; parameters λ and σ have negative correlation with estimations of options price, and overestimation of parameters will result in underestimation of options price. In particular, influence of correlation coefficient ρ on estimation of options price differs in different ranges of ρ; when ρ is less than 0.4, overestimation of ρ will result in underestimation of options price. When ρ is greater than 0.4, the result is opposite.Figure 4
Sensitivity ofκ.Figure 5
Sensitivity ofγ.Figure 6
Sensitivity ofρ.Figure 7
Sensitivity ofλ.Figure 8
Sensitivity ofσ.Calculate point elasticity of at-the-money options of parameters estimation, representing sensitivity of options price to each parameter shown in Table4.Table 4
Analysis of sensitivity of options price to parameters.
Parameters
Parameters estimation value
Sensitivity of options price to parameters
κ
0.4264
0.0264
γ
0.0946
0.0398
ρ
−0.6284
0.0259
λ
3.4543
0.1089
σ
0.4822
0.0551Point elasticities of at-the-money options of parameters estimation remain between 0.02 and 0.11, which indicate that sensitivities of options price to parameters are small and that the influences of errors on options price are not great.


## 4.1. Data Selection

This section selects day-degree data of Hang Seng Index Options price in Hong Kong Exchange. The data set consists of continuous 500 transaction days spanning from November 24th, 2011, to November 29th, 2013, and complete 247 transaction days spanning from December 3rd, 2012, to November 29th, 2013.In order to reduce dimensions to match number of parameters and input variables, we need to fix due dates of options selected in the sample. However, among option contract product sequences in Hong Kong Exchange, options that can be transacted each day include the following:(1)
Hang Seng Index Options expire at the end of each of the recent three month.
(2)
Hang Seng Index Options expire on the last transaction day at the end of March, June, September, and December within the year.Hence, it cannot be guaranteed that durations of options of each transaction day in the sample are identical. We can only defer the expiration date of sample options of each month by one month and guarantee that duration of selected option in the sample will be with limits. Specifically, we select data with its expiration date being 60–90 days. For example, for January 1st, we select the option contract that will expire at the end of March, so its duration will be about 90 days, but, for March 20th, we select the option contract which will become due at the end of May, so its duration will be about 70 days. For empirical research, we consider expiration dates of the two as identical.For convenience, we only select call option contract among these data. Since the trading volume of in-the-money options of Hong Kong market is small, and some only have single-digit trading volume, thus we select at-the-money options and some out-of-the-money options to conduct an empirical analysis with big trading volumes. Based on variation situation of everyday Hang Seng Index, we respectively find out exercise priceK of corresponding contracts approximate to at-the-money option S
t
/
K
≈
1, out-of-the-money option S
t
/
K
≈
0.95, and S
t
/
K
≈
0.92 as well as closing price C
t of these option contracts.


## 4.2. Analysis of Empirical Results

(
1) Parameter Estimation Results and Analysis. Table 1 lists estimation results selecting at-the-money option price data and out-of-the-money option price data within the sample.Table 1
Heston model estimation values and standard deviation of at-the-money call option and two out-of-the-money options based on three-month duration.
Heston model parameters
(
S
t
/
K
)
=
1Parameter estimation values
(
S
t
/
K
)
=
1Standard deviation
(
S
t
/
K
)
=
0.95Parameter estimation value
(
S
t
/
K
)
=
0.95Standard deviation
(
S
t
/
K
)
=
0.92Parameter estimation value
(
S
t
/
K
)
=
0.92Standard deviation
κ
0.4264
0.2784
0.4373
0.2086
0.4459
0.2113
γ
0.0946
0.0181
0.0915
0.0163
0.0909
0.0168
ρ
−0.6284
0.0867
−0.7011
0.0715
−0.6812
0.0790
λ
3.4543
3.0332
3.4280
2.9312
3.3788
2.7003
σ
0.4822
0.0433
0.3565
0.0333
0.3421
0.0315It can be seen from Table1 that, among three samples, κ obtained through the estimation basically stabilizes within 0.42–0.45. Parameter κ is mean reversion speed of volatility square which represents the speed of volatility and its square reversing to mean value when they deviate from the long-term mean value. Aït-Sahalia and Kimmel [17] conducted an empirical analysis using S&P 500 and its stock index option data and obtained that estimated value of reversion rate κ is 5.07, which means that, compared with American market, the mean reversion rate of Hong Kong market is slower, which indicates that when volatility in Hong Kong stock market deviates from mean value upwards or downwards, it will reverse to the mean level at a slow speed. Long-term mean value of volatility square is around 0.09, which indicates that underlying asset price volatility will move around 30% within a long period, and correlation coefficient ρ of two Brownian movements will be around −0.6 to −0.7, and λ stabilizes around 3.4, but variable range estimated value of volatility σ of stock index volatility is slightly larger, from 0.34 to 0.48.In the model,ρ is the risk source of stock index risk source and stock index volatility square, that is, the correlation coefficient between two Brownian movements d
W
1 and d
W
2; the above empirical results show that ρ is a negative value and remains within −0.6 to −0.7, which indicates that HSI stock index yield rate has a negative correlation with its volatility. This result has overthrown complete negative correlation or complete noncorrelation between underlying asset price and its volatility in other random volatility models to exactly explained “leverage effect” observed in the market.On the whole, standard deviation of estimated variable value of at-the-money option sample(
S
t
/
K
)
=
1 is higher than the two out-of-the-money option samples (
S
t
/
K
)
=
0.95 and (
S
t
/
K
)
=
0.92. We try to expand data length within the sample and select 500-day at-the-money option price data within the sample to conduct parameter estimation, and the results are as shown in Table 2.Table 2
Parameter estimation results of 500-day at-the-money option price data.
Heston model parameters
κ
γ
ρ
λ
σ
Estimated values of Heston model parameters
1.1409
12.8701
−0.9957
−175.9290
0.0000It can be seen that, among 500-day parameter estimation results, interest rater has negative values, σ is 0, and reversed mean γ of volatility mean reversion equation reaches as high as 12.8701, which severely does not conform to current situation. This phenomenon can be explained from the angle of parameter estimation method. MLE method adopts closing price of at-the-money call option contracts of the current day as data of this day, so the bigger the selected sample is, the bigger the time span will be. If 500-day data in the sample is selected, sample span should be more than two years. As options market is not fixed, so when MLE method is used to conduct parameter estimation of options pricing research, larger sample size does not mean more ideal results. Instead, a proper sample size should be selected to guarantee adequacy of sample data and avoid estimation deviation caused by changes of multiperiod price characteristics.(
2) Analysis of Options Pricing Error. Here, according to parameter estimation results of data in option sample under different value states, we calculate options pricing errors 15 days beyond the sample, and Figures 1–3, respectively, express options pricing errors under different value states.Figure 1
Options pricing error of at-the-money option(
S
t
/
K
)
=
1 price.Figure 2
Options pricing error of out-of-the-money option(
S
t
/
K
)
=
0.95 price.Figure 3
Options pricing error of in-the-money call option(
S
t
/
K
)
=
0.92 price.From Figures1–3, we can see that it has a higher accuracy to estimate model parameters and deduce options price using at-the-money options, compared with the two out-of-the-money options, and pricing errors of 15 option contracts outside the sample are within 10%, and there are 7 ones with their errors within 3% and 11 ones with errors within 6%, which indicates that the model can be better priced at-the-money call options with expiration date being about 3 months. Table 3 reflects the average pricing error, where(25)
B
i
a
s
i
=
C
i
m
o
d
e
l
t
-
C
i
m
a
r
k
e
t
t
,
P
¯
a
b
s
o
l
u
t
e
=
a
v
e
r
a
g
e
∑
t
=
1
T
∑
i
=
1
n
t
B
i
a
s
i
t
C
i
m
a
r
k
e
t
t
,
P
¯
a
c
t
u
a
l
=
a
v
e
r
a
g
e
∑
t
=
1
T
∑
i
=
1
n
t
B
i
a
s
i
t
C
i
m
a
r
k
e
t
t
.Table 3
Options pricing error.
(
S
t
/
K
)
=
1
(
S
t
/
K
)
=
0.95
(
S
t
/
K
)
=
0.92
Average absolute error (P
¯
a
b
s
o
l
u
t
e)
3.66%
9.23%
10.07%
Average relative error (P
¯
a
c
t
u
a
l)
3.30%
−1.18%
3.35%Other origins of parameter estimation error and model pricing error include simplification of real measure. In order to simplify the model and reduce number of parameters to be estimated, for the convenience of parameter estimation, the model adopted in the empirical analysis is the model form in risk neutral measure (Q measure); real measure (P measure) is not taken into consideration when parameters are being set, so price of risk factors in P measure has not been estimated out; for example, Heston model in (17) is(26)
d
S
t
Y
t
=
r
-
d
-
1
2
Y
t
+
λ
1
-
ρ
2
Y
t
κ
γ
-
Y
t
d
t
+
1
-
ρ
2
Y
t
0
ρ
Y
t
σ
Y
t
β
d
W
1
P
t
W
2
P
t
.(
3) Parameter Sensitivity Analysis of Options Pricing. In order to inspect reliance of options pricing effect on accuracy of parameter estimation, we, respectively, analyze sensitivity of each parameter of the model for options price estimation value.From Figures4–8, it is easy to see the parameters monotonicity of options price. Parameters κ and γ have positive correlation with estimations of options price, and overestimation of parameters will result in overestimation of options price; parameters λ and σ have negative correlation with estimations of options price, and overestimation of parameters will result in underestimation of options price. In particular, influence of correlation coefficient ρ on estimation of options price differs in different ranges of ρ; when ρ is less than 0.4, overestimation of ρ will result in underestimation of options price. When ρ is greater than 0.4, the result is opposite.Figure 4
Sensitivity ofκ.Figure 5
Sensitivity ofγ.Figure 6
Sensitivity ofρ.Figure 7
Sensitivity ofλ.Figure 8
Sensitivity ofσ.Calculate point elasticity of at-the-money options of parameters estimation, representing sensitivity of options price to each parameter shown in Table4.Table 4
Analysis of sensitivity of options price to parameters.
Parameters
Parameters estimation value
Sensitivity of options price to parameters
κ
0.4264
0.0264
γ
0.0946
0.0398
ρ
−0.6284
0.0259
λ
3.4543
0.1089
σ
0.4822
0.0551Point elasticities of at-the-money options of parameters estimation remain between 0.02 and 0.11, which indicate that sensitivities of options price to parameters are small and that the influences of errors on options price are not great.


## 5. Conclusion

This paper applies MLE to Heston model and conducts an empirical analysis with market price data of Hong Kong Hang Seng Index Options, and results show that Heston model has a better pricing effect on Hang Seng Index Options. Compared with other options, this model has presented the least pricing error of at-the-money options. Also, we find out that market characteristics have changed in Hong Kong market within two-year period. The parameters estimation effect obviously decreases using two-year market data, compared with estimation effect of about one-year data. This means sample length of MLE should not be too long. And parameters sensitivity analysis of option price indicates that point elasticities of at-the-money options remain between 0.02 and 0.11, which shows that options price is not that sensitive to parameters and reflects better applicability of MLE.



---

*Source: 7549537-2016-08-07.xml*